Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 41.4% intergenic (‑195/+113) 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,705,9850TC41.4% 234.0 / 304.8 242intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (70/72);  new base C (89/11);  total (159/83)
Fisher's exact test for biased strand distribution p-value = 3.94e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-02

Only 100 of 246 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084
                                                                                                  |                                                                                                            
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                             >  1:2116179/1‑100 (MQ=255)
  tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGCCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                           >  2:519938/1‑100 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                       <  1:2809485/100‑1 (MQ=37)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                      <  1:1699370/100‑2 (MQ=37)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                     <  2:1683930/100‑1 (MQ=255)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                     <  1:221226/100‑1 (MQ=255)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                     <  2:743460/100‑3 (MQ=37)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑a                                                                                                  <  1:213316/100‑1 (MQ=255)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCAT‑‑a                                                                                                  <  1:2517435/100‑1 (MQ=25)
          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑Ac                                                                                                 >  1:161693/1‑100 (MQ=255)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑Ac                                                                                                 >  2:3460089/1‑99 (MQ=255)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACa                                                                                               <  2:1529982/100‑1 (MQ=255)
             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACaaa                                                                                             >  1:2248387/1‑100 (MQ=255)
                gCACGCAGAGGCACAGCGCCATCGGGTAGACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaaag                                                                                             <  2:2831839/100‑11 (MQ=18)
                gCACGCAGAGGAACAGCGCCACCGGGTATACAGATAGTGGACAAAGTGCGCTTTGTTTATGCCGAATAAGGCGTTAACGCCGCATCCGGcataaacaaa                                                                                              <  2:1717787/99‑10 (MQ=18)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCg                                                                                         <  2:2057814/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTg                                                                                       <  1:1370036/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTg                                                                                       >  2:613041/1‑100 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGc                                                                                      <  1:1974731/100‑1 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGc                                                                                      <  2:3582598/100‑1 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTc                                                                                <  1:3035695/99‑1 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCa                                                                               >  1:3140673/1‑100 (MQ=255)
                            gcaGCGCCATCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCaa                                                                              >  1:3643810/2‑100 (MQ=255)
                            acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCaa                                                                              >  1:1794866/1‑100 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAAt                                                                             <  2:671520/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAAt                                                                             <  2:310388/100‑1 (MQ=255)
                                 gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACAt                                                                         >  1:134892/1‑100 (MQ=255)
                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACAtt                                                                        >  1:950880/1‑100 (MQ=255)
                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAAtt                                                             >  1:2340773/1‑100 (MQ=255)
                                           ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAAttt                                                           >  1:2943577/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                      <  1:1469346/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                      <  2:647631/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                                      <  1:3068839/100‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTa                                                    >  2:3064913/1‑100 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTa                                                    >  2:3407641/1‑100 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGcc                                                <  1:3311824/100‑1 (MQ=255)
                                                    ataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGcc                                                >  1:2819664/5‑100 (MQ=18)
                                                     ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCt                                               <  2:1056185/100‑1 (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                                       >  1:3020061/1‑100 (MQ=255)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCg                                       >  1:115850/1‑100 (MQ=255)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCg                                       <  2:1865346/100‑1 (MQ=255)
                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGGATTTTC‑‑GTAGGCCTGCTAAGCGGa                                     >  1:3262708/1‑100 (MQ=255)
                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGa                                     <  1:805363/100‑1 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAGTCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAg                                    <  2:1587848/99‑1 (MQ=255)
                                                                  agttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAAAACTGTGCAAATTCAATATATTGCAGGAAACAC‑‑GTAGGCCTGATAAGCGAAg                                    >  2:754061/11‑100 (MQ=2)
                                                                   tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAgcgc                                 >  2:328768/1‑100 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                             >  1:1692790/1‑100 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCT‑‑ACAAATTCTTCCAAATTCAATATATTGCAGAAATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                            >  1:1925295/5‑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                          >  1:1833242/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAACAGTGCAAATTCAATATATTGCAATCTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc                          >  2:488448/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACTCCGCATCCGGCCT‑‑ACAAAAGCACGCAAATTCAATATATTGCAGAGATTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:3564957/3‑100 (MQ=11)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:1988866/3‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACCTTGCTAATTCAATATATTGCAGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                          >  1:1032579/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAACATGCAAATTCAATAGATTGCAGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:2900952/3‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                        >  1:2154372/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2437847/2‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCATGAAAATTCAATACATTGCAAGATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:145228/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:2463804/2‑99 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:13546/2‑100 (MQ=18)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAACAGTGCAAATTCAATATATTGCAATCTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaa                        >  1:1772547/1‑99 (MQ=18)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCATGAAAATTCAATACATTGCAAGATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCac                        >  1:736990/1‑99 (MQ=17)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATCTCCGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  2:2319068/1‑99 (MQ=18)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAATCATGCAAATCCAAAATATTGCAGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  2:2222745/1‑100 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAATCATGAAAATTCAATACATTGCAAGATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacg                       >  1:572264/9‑98 (MQ=17)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       >  1:1554762/1‑100 (MQ=255)
                                                                             ccGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat                       >  2:2763196/1‑98 (MQ=21)
                                                                              cggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAGATCGTGCAAATTCAATATATTGCATGAGCCAC‑‑GTATGCCTGATAAGCGTAGAGCATCAGGCa                         >  1:1318405/8‑97 (MQ=17)
                                                                              cggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAGCACGCAAACTCAATATATTGCAGAGATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                      >  1:2044099/8‑97 (MQ=17)
                                                                              cGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                      <  2:2670477/100‑1 (MQ=255)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGCCAtttgc                    >  2:694816/7‑98 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:1127222/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:2485018/7‑100 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:1005414/7‑100 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:1740579/7‑100 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAGATCGTGCAAATTCAATATATTGCATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGAAgttt                     >  1:3546941/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                     >  1:2577577/7‑96 (MQ=25)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAACCTTGCTAATTCAATATATTGCAGGGACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                     >  1:459185/7‑96 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAACATGCAAATTCAATAGATTGCAGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                     >  1:3495876/7‑100 (MQ=17)
                                                                                gataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAGATCGTGCAAATTCAATATATTGCATGAGCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                    >  1:2710035/6‑100 (MQ=17)
                                                                                gataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttg                    >  1:123204/6‑99 (MQ=18)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:1065965/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  2:3140650/100‑1 (MQ=255)
                                                                                          tAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          >  1:620544/1‑100 (MQ=255)
                                                                                          tAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          <  1:3407666/100‑1 (MQ=255)
                                                                                           aaCGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGACTTTTC‑‑GTAGGCCTGATAAGAGGAGCGCATCAGGCATTTTTGCCTTTGCCa          >  2:1497602/1‑99 (MQ=255)
                                                                                           aaCGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATCTTC‑‑GTAGGCCTGATAAGCGGATCGCATCAGGCATTtttgcctttgtcct         >  2:2211785/1‑96 (MQ=255)
                                                                                            aCGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        >  1:3377901/1‑100 (MQ=255)
                                                                                            aCGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  2:1957047/100‑1 (MQ=255)
                                                                                             cGCCTTATCCGGCCT‑‑ACAAGATCGTGAAAATTCAATACATTGCAGCATTTTC‑‑GTAGGCCGGATCAGCGGAACGCATTAGGCATTtttgcctttggta          >  2:3022143/1‑94 (MQ=17)
                                                                                             cGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        >  2:3341122/1‑99 (MQ=255)
                                                                                              gCCTTATCCGGCGT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      <  2:1266239/100‑1 (MQ=255)
                                                                                                 ttATCCGGCTT‑‑ACAAAATCGTGCAAATTCAATAAATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCAATCAGTCt   <  2:1240889/100‑1 (MQ=255)
                                                                                                 ttATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt   <  1:1135717/100‑1 (MQ=255)
                                                                                                  tATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  1:2451083/1‑100 (MQ=255)
                                                                                                  |                                                                                                            
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084

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