Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,017 C→T 30.8% intergenic (‑227/+81) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,706,0170CT30.8% 350.5 / 165.6 227intergenic (‑227/+81)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (70/74);  new base T (62/2);  total (149/78)
Fisher's exact test for biased strand distribution p-value = 2.79e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 229 total aligned reads displayed.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116
                                                                                                      |                                                                                                         
gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATAc                                                                                                           >  1:2019660/1‑100 (MQ=255)
   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAtt                                                                                                        >  1:950880/1‑100 (MQ=255)
              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTc                                                                                         <  2:3074507/100‑1 (MQ=255)
               cAGATTGTTGACAAANTGCGNTTTGTTNATGCCGGNTGCGGCGTTAACGNCTTNTCCGGCCT‑‑AC‑AAAATCGTGCAAATNCAANACATTGCA‑GA‑AT‑T‑TTc                                                                                         >  2:2221156/1‑99 (MQ=17)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑g                                                                                      <  2:647631/100‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑g                                                                                      >  1:1918879/1‑100 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑g                                                                                      <  1:1469346/100‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑g                                                                                      <  1:3068839/100‑1 (MQ=255)
                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAgg                                                                                  <  1:3513273/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGcc                                                                                <  2:740877/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGcc                                                                                <  1:3311824/100‑1 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGa                                                                             >  2:1263608/1‑100 (MQ=255)
                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATa                                                                           >  1:93419/1‑100 (MQ=255)
                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGc                                                                        <  1:3213095/100‑1 (MQ=255)
                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGAA‑GA‑AT‑T‑TTC‑‑GTGGGACTGATatg                                                                         >  2:3579623/1‑97 (MQ=17)
                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCg                                                                       <  2:2801349/100‑1 (MQ=255)
                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCg                                                                       >  2:1485102/1‑100 (MQ=255)
                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCg                                                                       >  1:3020061/1‑100 (MQ=255)
                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCg                                                                       >  1:115850/1‑100 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTTCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                  <  1:2451639/100‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                  >  1:446404/1‑100 (MQ=255)
                                   agttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAAAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑A‑CAC‑‑GTAGGCCTGATAAGCGAAg                                                                    >  2:754061/11‑100 (MQ=2)
                                    tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAgcgc                                                                 >  2:328768/1‑100 (MQ=255)
                                         gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                           >  1:2192925/5‑100 (MQ=18)
                                         gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑T‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAg                                                            >  2:2763245/5‑100 (MQ=18)
                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                            >  1:2535929/5‑100 (MQ=18)
                                          taggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                           >  1:3458623/4‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTTACGCCGCATCTGGCCT‑‑AC‑AAATTCTTCCAAATTCAATATATTGCA‑GA‑AA‑T‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  2:1822416/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCagatcgga                                                         >  2:1140318/3‑93 (MQ=18)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:1847401/3‑100 (MQ=17)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑T‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc                                                          >  2:488448/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACCTTGCTAATTCAATATATTGCAGGG‑AC‑T‑AT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                                                          >  1:549715/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑T‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:2900952/3‑100 (MQ=17)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  1:3553624/1‑100 (MQ=255)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                        >  1:2154372/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:2813960/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:2065962/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:2437847/2‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑C‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  2:2463804/2‑99 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑C‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:13546/2‑100 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                        >  1:1986198/1‑99 (MQ=17)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑A‑TA‑T‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcgat                                                       >  2:582636/1‑97 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCac                                                        >  1:736990/1‑99 (MQ=17)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑CT‑C‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                        >  2:2319068/1‑99 (MQ=18)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑T‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                       >  1:1666181/1‑100 (MQ=21)
                                             gCCGGATAAGGCGTTAACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACACTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                        >  1:2869837/1‑100 (MQ=18)
                                              ncGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTTGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                       >  1:135915/2‑100 (MQ=255)
                                              ccggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑T‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaag                                                       >  1:1770171/9‑98 (MQ=18)
                                              ccggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑C‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                       >  2:2357918/9‑98 (MQ=17)
                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                       >  1:1554762/1‑100 (MQ=255)
                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                       >  1:2078441/1‑100 (MQ=255)
                                               cggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg                                                     >  1:1821163/8‑99 (MQ=21)
                                               cggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑T‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaagc                                                      >  2:1405910/8‑97 (MQ=18)
                                               cggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑C‑CAC‑‑GTAGGCCTGATAAGCGTAGCACATCAGGCagt                                                       >  1:2959286/8‑97 (MQ=17)
                                               cGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                                                      <  2:2670477/100‑1 (MQ=255)
                                                ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                    >  1:1051619/7‑98 (MQ=21)
                                                ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑AT‑CCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGCCAtttgc                                                    >  2:694816/7‑98 (MQ=18)
                                                ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:1127222/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:1005414/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:1740579/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:2485018/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:1726801/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑A‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                     >  1:2577577/7‑96 (MQ=25)
                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                                                     >  1:1778313/7‑96 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑T‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                     >  1:825580/7‑100 (MQ=17)
                                                 gataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑C‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                    >  1:2710035/6‑100 (MQ=17)
                                                 gataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgctc                                                    >  1:3251252/6‑95 (MQ=17)
                                                   taaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑C‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttagc                                                  >  1:40017/4‑97 (MQ=17)
                                                    gcggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                    >  2:3405976/1‑97 (MQ=255)
                                                     cggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                <  2:1171781/100‑1 (MQ=255)
                                                        cgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                             <  2:3140650/100‑1 (MQ=255)
                                                           tAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                          >  1:620544/1‑100 (MQ=255)
                                                            aaCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑CT‑T‑TTC‑‑GTAGGCCTGATAAGAGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                          >  2:1497602/1‑99 (MQ=255)
                                                             aCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                        <  2:1957047/100‑1 (MQ=255)
                                                             aCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                        >  1:3377901/1‑100 (MQ=255)
                                                              cGCCTTATCCGGCCT‑‑AC‑AAGATCGTGAAAATTCAATACATTGCA‑GC‑AT‑T‑TTC‑‑GTAGGCCGGATCAGCGGAACGCATTAGGCATTtttgcctttggta                                          >  2:3022143/1‑94 (MQ=17)
                                                              cGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                        >  2:3341122/1‑99 (MQ=255)
                                                                  ttATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                   <  1:1135717/100‑1 (MQ=255)
                                                                   tATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt                                  >  1:2451083/1‑100 (MQ=255)
                                                                   tATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt                                  <  2:1085753/100‑1 (MQ=255)
                                                                       cGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              >  1:2373034/1‑100 (MQ=255)
                                                                          ccT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                           <  1:2288731/100‑1 (MQ=255)
                                                                          ccT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                           >  1:2146171/1‑100 (MQ=255)
                                                                            t‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                         <  1:1321418/100‑1 (MQ=255)
                                                                               aC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                        >  1:1383554/1‑100 (MQ=255)
                                                                                c‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                       <  2:1660389/100‑1 (MQ=255)
                                                                                c‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                       <  1:3123319/100‑1 (MQ=255)
                                                                                         tGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               >  1:3596236/1‑100 (MQ=255)
                                                                                         tGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  2:1287548/100‑1 (MQ=255)
                                                                                         tGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  1:1203237/100‑1 (MQ=255)
                                                                                          gCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              >  2:2921043/1‑100 (MQ=255)
                                                                                           cAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  2:3552922/98‑1 (MQ=255)
                                                                                           cAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa             <  1:101148/100‑1 (MQ=255)
                                                                                            aaaTTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt            <  2:1578034/100‑1 (MQ=255)
                                                                                              aTTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag          <  2:3051039/100‑1 (MQ=255)
                                                                                               ttCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga         >  2:2745605/1‑100 (MQ=255)
                                                                                                 cAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag        <  1:1261290/99‑1 (MQ=255)
                                                                                                 cAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga       <  1:1542700/100‑1 (MQ=255)
                                                                                                  aaTACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc      <  1:1927093/100‑1 (MQ=255)
                                                                                                   aTACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACg     >  1:1086813/1‑100 (MQ=255)
                                                                                                      |                                                                                                         
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑