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‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116
                                                                                                      |                                                                                                         
gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATAc                                                                                                           >  1:2019660/1‑100 (MQ=255)
   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACAtt                                                                                                        >  1:950880/1‑100 (MQ=255)
              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTc                                                                                         <  2:3074507/100‑1 (MQ=255)
               cAGATTGTTGACAAANTGCGNTTTGTTNATGCCGGNTGCGGCGTTAACGNCTTNTCCGGCCT‑‑ACAAAATCGTGCAAATNCAANACATTGCAGAATTTTc                                                                                         >  2:2221156/1‑99 (MQ=17)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                                                      <  2:647631/100‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                                                      >  1:1918879/1‑100 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                                                      <  1:1469346/100‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                                                      <  1:3068839/100‑1 (MQ=255)
                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAgg                                                                                  <  1:3513273/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGcc                                                                                <  2:740877/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGcc                                                                                <  1:3311824/100‑1 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGa                                                                             >  2:1263608/1‑100 (MQ=255)
                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATa                                                                           >  1:93419/1‑100 (MQ=255)
                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGc                                                                        <  1:3213095/100‑1 (MQ=255)
                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGAAGAATTTTC‑‑GTGGGACTGATatg                                                                         >  2:3579623/1‑97 (MQ=17)
                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCg                                                                       <  2:2801349/100‑1 (MQ=255)
                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCg                                                                       >  2:1485102/1‑100 (MQ=255)
                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCg                                                                       >  1:3020061/1‑100 (MQ=255)
                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCg                                                                       >  1:115850/1‑100 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTTCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                  <  1:2451639/100‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                  >  1:446404/1‑100 (MQ=255)
                                   agttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAAAACTGTGCAAATTCAATATATTGCAGGAAACAC‑‑GTAGGCCTGATAAGCGAAg                                                                    >  2:754061/11‑100 (MQ=2)
                                    tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcgc                                                                 >  2:328768/1‑100 (MQ=255)
                                         gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                           >  1:2192925/5‑100 (MQ=18)
                                         gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAACAGTGCAAATTCAATATATTGCAATCTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAg                                                            >  2:2763245/5‑100 (MQ=18)
                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                            >  1:2535929/5‑100 (MQ=18)
                                          taggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCTAATTCAATATATTGCAGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                           >  1:3458623/4‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTTACGCCGCATCTGGCCT‑‑ACAAATTCTTCCAAATTCAATATATTGCAGAAATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  2:1822416/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCagatcgga                                                         >  2:1140318/3‑93 (MQ=18)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:1847401/3‑100 (MQ=17)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAACAGTGCAAATTCAATATATTGCAATCTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc                                                          >  2:488448/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACCTTGCTAATTCAATATATTGCAGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                                                          >  1:549715/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAACATGCAAATTCAATAGATTGCAGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:2900952/3‑100 (MQ=17)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  1:3553624/1‑100 (MQ=255)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                        >  1:2154372/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:2813960/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:2065962/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:2437847/2‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  2:2463804/2‑99 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:13546/2‑100 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                        >  1:1986198/1‑99 (MQ=17)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAGTTCGTGCAAATTCAATAAATTGCAATATGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcgat                                                       >  2:582636/1‑97 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCATGAAAATTCAATACATTGCAAGATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCac                                                        >  1:736990/1‑99 (MQ=17)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATCTCCGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                        >  2:2319068/1‑99 (MQ=18)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAATCATGCAAATCCAAAATATTGCAGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                       >  1:1666181/1‑100 (MQ=21)
                                             gCCGGATAAGGCGTTAACGCCGCATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACACTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                        >  1:2869837/1‑100 (MQ=18)
                                              ncGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTTGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                       >  1:135915/2‑100 (MQ=255)
                                              ccggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAACAGTGCAAATTCAATATATTGCAATCTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaag                                                       >  1:1770171/9‑98 (MQ=18)
                                              ccggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAGATCGTGCAAATTCAATATATTGCATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                       >  2:2357918/9‑98 (MQ=17)
                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                       >  1:1554762/1‑100 (MQ=255)
                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                       >  1:2078441/1‑100 (MQ=255)
                                               cggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg                                                     >  1:1821163/8‑99 (MQ=21)
                                               cggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAACAGTGCAAATTCAATATATTGCAATCTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaagc                                                      >  2:1405910/8‑97 (MQ=18)
                                               cggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAGATCGTGCAAATTCAATATATTGCATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCACATCAGGCagt                                                       >  1:2959286/8‑97 (MQ=17)
                                               cGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                                                      <  2:2670477/100‑1 (MQ=255)
                                                ggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                    >  1:1051619/7‑98 (MQ=21)
                                                ggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGCCAtttgc                                                    >  2:694816/7‑98 (MQ=18)
                                                ggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:1127222/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:1005414/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:1740579/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:2485018/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:1726801/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                     >  1:2577577/7‑96 (MQ=25)
                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAATCATGAAAATTCAATACATTGCAAGATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                                                     >  1:1778313/7‑96 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAACATGCAAATTCAATAGATTGCAGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                     >  1:825580/7‑100 (MQ=17)
                                                 gataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAGATCGTGCAAATTCAATATATTGCATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                    >  1:2710035/6‑100 (MQ=17)
                                                 gataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAATCATGAAAATTCAATACATTGCAAGATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgctc                                                    >  1:3251252/6‑95 (MQ=17)
                                                   taaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAGCACGCAAACTCAATATATTGCAGAGATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttagc                                                  >  1:40017/4‑97 (MQ=17)
                                                    gcggcgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                    >  2:3405976/1‑97 (MQ=255)
                                                     cggcgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                <  2:1171781/100‑1 (MQ=255)
                                                        cgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                             <  2:3140650/100‑1 (MQ=255)
                                                           tAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                          >  1:620544/1‑100 (MQ=255)
                                                            aaCGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGACTTTTC‑‑GTAGGCCTGATAAGAGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                          >  2:1497602/1‑99 (MQ=255)
                                                             aCGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                        <  2:1957047/100‑1 (MQ=255)
                                                             aCGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                        >  1:3377901/1‑100 (MQ=255)
                                                              cGCCTTATCCGGCCT‑‑ACAAGATCGTGAAAATTCAATACATTGCAGCATTTTC‑‑GTAGGCCGGATCAGCGGAACGCATTAGGCATTtttgcctttggta                                          >  2:3022143/1‑94 (MQ=17)
                                                              cGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                        >  2:3341122/1‑99 (MQ=255)
                                                                  ttATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                   <  1:1135717/100‑1 (MQ=255)
                                                                   tATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt                                  >  1:2451083/1‑100 (MQ=255)
                                                                   tATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt                                  <  2:1085753/100‑1 (MQ=255)
                                                                       cGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              >  1:2373034/1‑100 (MQ=255)
                                                                          ccT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                           <  1:2288731/100‑1 (MQ=255)
                                                                          ccT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                           >  1:2146171/1‑100 (MQ=255)
                                                                            t‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                         <  1:1321418/100‑1 (MQ=255)
                                                                               aCAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                        >  1:1383554/1‑100 (MQ=255)
                                                                                cAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                       <  2:1660389/100‑1 (MQ=255)
                                                                                cAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                       <  1:3123319/100‑1 (MQ=255)
                                                                                         tGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               >  1:3596236/1‑100 (MQ=255)
                                                                                         tGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  2:1287548/100‑1 (MQ=255)
                                                                                         tGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  1:1203237/100‑1 (MQ=255)
                                                                                          gCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              >  2:2921043/1‑100 (MQ=255)
                                                                                           cAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  2:3552922/98‑1 (MQ=255)
                                                                                           cAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa             <  1:101148/100‑1 (MQ=255)
                                                                                            aaaTTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt            <  2:1578034/100‑1 (MQ=255)
                                                                                              aTTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag          <  2:3051039/100‑1 (MQ=255)
                                                                                               ttCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga         >  2:2745605/1‑100 (MQ=255)
                                                                                                 cAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag        <  1:1261290/99‑1 (MQ=255)
                                                                                                 cAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga       <  1:1542700/100‑1 (MQ=255)
                                                                                                  aaTACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc      <  1:1927093/100‑1 (MQ=255)
                                                                                                   aTACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACg     >  1:1086813/1‑100 (MQ=255)
                                                                                                      |                                                                                                         
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑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: