Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,089,197 A→T 6.2% intergenic (‑69/+14) yeeZ ← / ← yoeB putative epimerase/toxin of the YoeB‑YefM toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,089,1970AT 6.2% 331.9 / 9.6 96intergenic (‑69/+14)yeeZ/yoeBputative epimerase/toxin of the YoeB‑YefM toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base A (46/44);  new base T (2/4);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 104 total aligned reads displayed.

TTGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAC  >  NC_000913/2089051‑2089344
                                                                                                                                                  |                                                                                                                                                    
ttGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTat‑g                                                                                                                                                   >  4:175981/1‑149 (MQ=255)
ttGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTat‑g                                                                                                                                                   >  1:508098/1‑149 (MQ=255)
  gCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTat‑gaa                                                                                                                                                 <  1:595143/149‑1 (MQ=255)
  gCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTat‑gaa                                                                                                                                                 <  2:626206/149‑1 (MQ=255)
      gCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAt‑gaatgaa                                                                                                                                             <  2:132789/149‑1 (MQ=255)
      gCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCATa‑gaatgaa                                                                                                                                             <  1:164399/149‑1 (MQ=255)
       ccTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑gaatgaat                                                                                                                                            <  2:144263/149‑1 (MQ=255)
       ccTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑gaatgaat                                                                                                                                            <  4:404078/149‑1 (MQ=255)
        cTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCATA‑gaatgaatt                                                                                                                                           <  3:463025/149‑2 (MQ=255)
         tCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCt                                                                                                                                          >  4:434748/1‑149 (MQ=255)
         tCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCt                                                                                                                                          <  1:228401/149‑1 (MQ=255)
            cgcTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGtt                                                                                                                                       <  2:543833/149‑1 (MQ=255)
            cgcTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTC‑ATAgaatgaattgtca                                                                                                                                       <  3:113196/149‑6 (MQ=255)
             gcTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTc                                                                                                                                      <  2:508098/149‑1 (MQ=255)
                   aGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCaataat                                                                                                                                <  2:609069/149‑1 (MQ=255)
                      gACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGTGACCTTTTTCATCCGTTACCTCCTGATACACCATAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATGTGTTCAATAATGAt                                                                                                                             >  3:633712/1‑149 (MQ=255)
                      gACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGAt                                                                                                                             >  2:327999/1‑149 (MQ=255)
                       aCATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATa                                                                                                                            >  1:426391/1‑149 (MQ=255)
                         aTCGCCAGCGGCATGCCT‑ACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACg                                                                                                                         >  2:63507/1‑149 (MQ=255)
                          tCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACg                                                                                                                         >  1:454029/1‑149 (MQ=255)
                           cGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGa                                                                                                                        <  3:295447/149‑1 (MQ=255)
                             ccAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTat‑gaat                                                                                                                                                <  2:38441/123‑1 (MQ=255)
                             ccAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTat‑gaat                                                                                                                                                >  1:38441/1‑123 (MQ=255)
                               aGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCACCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACAtg                                                                                                                    <  2:454029/149‑1 (MQ=255)
                                   gCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCATA‑gaatgaatt                                                                                                                                           <  1:747376/122‑2 (MQ=255)
                                   gCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCATA‑gaatgaatt                                                                                                                                           >  2:747376/1‑121 (MQ=255)
                                     aTGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCaa                                                                                                              >  1:364158/1‑149 (MQ=255)
                                           aaCCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGa                                                                                                                              <  3:659091/127‑1 (MQ=255)
                                           aaCCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGa                                                                                                                              >  4:659091/1‑127 (MQ=255)
                                           aaCCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCa                                                                                                        >  3:473983/1‑149 (MQ=255)
                                            accaccCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAg                                                                                                       >  3:753629/1‑149 (MQ=255)
                                            accaccCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAg                                                                                                       <  4:290052/149‑1 (MQ=255)
                                               accCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCaataat                                                                                                                                <  4:409981/121‑1 (MQ=255)
                                                 ccTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAAt                                                                                                  <  1:680961/149‑1 (MQ=255)
                                                     aCCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAAtcgtc                                                                                              <  2:85081/149‑1 (MQ=255)
                                                     aCCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAAtcgtc                                                                                              <  4:131195/149‑1 (MQ=255)
                                                      cccGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAAgcgtct                                                                                             <  4:516808/149‑2 (MQ=255)
                                                        cGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGt                                                                                           >  1:319711/1‑149 (MQ=255)
                                                         gACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTa                                                                                          >  2:35852/1‑149 (MQ=255)
                                                                 cGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGta                                                                                  <  4:472235/149‑1 (MQ=255)
                                                                  gACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGtat                                                                                 >  2:418167/1‑149 (MQ=255)
                                                                   aCCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTAAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGtata                                                                                >  3:654818/1‑149 (MQ=255)
                                                                     cTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATAcc                                                                              >  4:676598/1‑149 (MQ=255)
                                                                      tttttCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCa                                                                             >  2:434801/1‑149 (MQ=255)
                                                                          tCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACg                                                                         >  1:508814/1‑149 (MQ=255)
                                                                           cATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACgg                                                                        >  3:605124/1‑149 (MQ=255)
                                                                             tCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTg                                                                      >  2:202468/1‑149 (MQ=255)
                                                                               cGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGct                                                                    >  3:700743/1‑149 (MQ=255)
                                                                                  tATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGctcct                                                                 <  1:63507/149‑1 (MQ=255)
                                                                                      tcCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGt                                                             >  1:333797/1‑149 (MQ=255)
                                                                                      tcCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGt                                                             >  3:218814/1‑149 (MQ=255)
                                                                                       ccTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTa                                                            <  1:149038/149‑1 (MQ=255)
                                                                                        cTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTaa                                                           <  3:536281/149‑1 (MQ=255)
                                                                                        cTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTaa                                                           >  1:738730/1‑149 (MQ=255)
                                                                                              cacCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTc                                                     >  1:459087/1‑149 (MQ=255)
                                                                                               acCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCg                                                    <  4:146521/149‑1 (MQ=255)
                                                                                                ccTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCgg                                                   >  1:189838/1‑149 (MQ=255)
                                                                                                 cTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCggg                                                  >  2:103577/1‑149 (MQ=255)
                                                                                                  ttAGATCTATAAGGCTACGCTAGCGTATCATAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGa                                                 >  3:232430/1‑149 (MQ=255)
                                                                                                  ttAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCgg                                                   >  2:645726/1‑147 (MQ=255)
                                                                                                  ttAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCgg                                                   <  1:645726/147‑1 (MQ=255)
                                                                                                   tAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGAc                                                <  1:410250/149‑1 (MQ=255)
                                                                                                      aTCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAg                                             >  1:207456/1‑149 (MQ=255)
                                                                                                         tataAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGaagccaga                                          >  4:585303/1‑142 (MQ=255)
                                                                                                          ataAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAAc                                         >  2:764454/1‑149 (MQ=255)
                                                                                                           taAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAAcc                                        >  1:771634/1‑149 (MQ=255)
                                                                                                            aaGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCATTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGACACCt                                       <  2:587158/149‑1 (MQ=255)
                                                                                                            aaGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCt                                       <  1:133450/149‑1 (MQ=255)
                                                                                                                 tACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACaa                                  <  4:103151/149‑1 (MQ=255)
                                                                                                                  aCGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACaaa                                 <  1:327999/149‑1 (MQ=255)
                                                                                                                    gCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAAtt                               >  2:531471/1‑149 (MQ=255)
                                                                                                                        gCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGt                           >  1:496254/1‑149 (MQ=255)
                                                                                                                          gTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCt                                       >  2:520164/1‑135 (MQ=255)
                                                                                                                          gTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCt                                       <  1:520164/135‑1 (MQ=255)
                                                                                                                            aTCAAAACTGACAATTCATTCTAT‑GAATTAATCTGTTAAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGATTACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATGATGTTTCa                       >  4:149847/1‑149 (MQ=255)
                                                                                                                            aTCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCa                       <  1:495932/149‑1 (MQ=255)
                                                                                                                             tCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTc                                                     <  1:479853/118‑1 (MQ=255)
                                                                                                                             tCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTc                                                     >  2:479853/1‑118 (MQ=255)
                                                                                                                             tCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAg                      <  4:605124/149‑1 (MQ=255)
                                                                                                                             tCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAg                      <  2:279458/149‑1 (MQ=255)
                                                                                                                              cAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAgg                     <  2:319711/149‑1 (MQ=255)
                                                                                                                              cAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAgg                     <  2:480078/149‑1 (MQ=255)
                                                                                                                              cAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAgg                     <  1:661096/149‑1 (MQ=255)
                                                                                                                               aaaaCTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAggg                    >  3:589956/1‑149 (MQ=255)
                                                                                                                                aaaCTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTAt                             <  1:702685/139‑1 (MQ=255)
                                                                                                                                aaaCTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTAt                             >  2:702685/1‑139 (MQ=255)
                                                                                                                                aaaCTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAgggg                   >  2:424607/1‑149 (MQ=255)
                                                                                                                                 aaCTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGt                  >  1:548706/1‑149 (MQ=255)
                                                                                                                                  aCTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGtt                 <  2:189838/149‑1 (MQ=255)
                                                                                                                                   cTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTc                >  2:127642/1‑149 (MQ=255)
                                                                                                                                    tGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCt               <  2:207456/149‑1 (MQ=255)
                                                                                                                                      aCAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTgg             <  2:281607/149‑1 (MQ=255)
                                                                                                                                       cAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGc            <  3:140814/149‑1 (MQ=255)
                                                                                                                                       agattcattcATA‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACcgc                                                                                     <  1:519513/74‑1 (MQ=38)
                                                                                                                                         attcattcATA‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGc            >  2:590741/1‑147 (MQ=255)
                                                                                                                                                 tAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAc  >  3:164631/1‑149 (MQ=255)
                                                                                                                                                 tAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAc  <  2:738730/149‑1 (MQ=255)
                                                                                                                                                 tAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAc  >  1:178879/1‑149 (MQ=255)
                                                                                                                                                  at‑gaatgaATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGc                                                                     >  2:409000/1‑81 (MQ=255)
                                                                                                                                                  at‑gaatgaATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGc                                                                     <  1:409000/81‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTGCCAGCCTCGCGCTGAAAGTGACATCGCCAGCGGCATGCCTAACCACCCTAACCCGACAATTGCGACCTTTTTCATCCGTTATCTCCTGATACACCTTAGATCTATAAGGCTACGCTAGCGTATCAAAACTGACAATTCATTCTAT‑GAATGAATCTGTTCAATAATGATAACGACATGCTGCAATGAGCAGTGAATCGTCGGTAACCGCGTATACCAGACGGTGCTCCTCTGTAATGCGTCGGGACCAGAAACCTGACAAATTATGTTTCAGGGGTTCTGGCTTCCCCTTAC  >  NC_000913/2089051‑2089344

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑