Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,017 C→T 22.9% 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,0170CT22.9% 338.8 / 87.1 166intergenic (‑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 (54/65);  new base T (34/2);  total (99/67)
Fisher's exact test for biased strand distribution p-value = 5.29e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01

Only 100 of 169 total aligned reads displayed.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116
                                                                                                    |                                                                                                       
  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                                                       >  1:354841/1‑100 (MQ=255)
    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                                                     <  1:1547195/100‑1 (MQ=255)
    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                                                     >  2:2631388/1‑100 (MQ=255)
    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                                                     <  1:1694706/100‑1 (MQ=255)
     aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                                                    <  1:1133869/100‑1 (MQ=255)
       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAg                                                                                                  <  2:768951/100‑1 (MQ=255)
         ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAACTTCAATACATTGCAGa‑a                                                                                               >  2:1046636/1‑100 (MQ=255)
             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑At‑ttt                                                                                          >  2:206616/1‑100 (MQ=255)
             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑At‑ttt                                                                                          >  2:241697/1‑100 (MQ=255)
              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTc                                                                                         >  1:1065956/1‑100 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑g                                                                                      <  2:2193861/100‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑g                                                                                      <  1:666549/100‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑g                                                                                      <  1:729138/100‑1 (MQ=255)
                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTa                                                                                    >  2:46913/1‑100 (MQ=255)
                   gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCGGGCCTCC‑AAAGTCGTGCAAGTTCAATACATTGTAGA‑Ag‑ttt                                                                                          >  2:2412171/7‑94 (MQ=14)
                    tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGc                                                                                 <  1:699052/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGcc                                                                                <  2:2429716/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGcc                                                                                <  2:710041/100‑1 (MQ=255)
                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTg                                                                              <  1:300445/100‑1 (MQ=255)
                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTg                                                                              >  1:1220658/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGa                                                                             >  2:2443537/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGa                                                                             >  1:407172/1‑100 (MQ=255)
                         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGAt                                                                            <  1:2648097/100‑1 (MQ=255)
                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATa                                                                           <  2:1750035/100‑1 (MQ=255)
                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGc                                                                        <  1:8254/100‑1 (MQ=255)
                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGc                                                                        <  1:2450090/100‑1 (MQ=255)
                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCg                                                                       <  1:2591667/100‑1 (MQ=255)
                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGGATAT‑AAAATCGTTCAGATTCAATACATTGCAGA‑AT‑TTAC‑‑GTAGTCCTGAt                                                                            <  2:667154/94‑1 (MQ=17)
                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                                                                      <  2:1711901/100‑1 (MQ=255)
                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                                                                      <  1:501753/100‑1 (MQ=255)
                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                                                                    >  2:2447668/1‑100 (MQ=255)
                                   gcttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCAGG‑AG‑CTGC‑‑GTAGGCCTGATAAGCGTAgcg                                                                  >  1:1750756/11‑100 (MQ=18)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATaagcgaagcg                                                                  >  2:308633/1‑100 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGcgcg                                                                  <  1:1602341/100‑1 (MQ=255)
                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                              <  1:1421892/100‑1 (MQ=255)
                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                              >  1:1519874/1‑100 (MQ=255)
                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                            <  1:1105509/100‑1 (MQ=255)
                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                           <  1:1939482/100‑1 (MQ=255)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑A‑AT‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  2:1379932/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                          >  2:101518/3‑100 (MQ=21)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGTTAATTCAATATATTGCAGG‑GA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                                                          >  2:6433/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCAGG‑AT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:2683554/3‑100 (MQ=17)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:2052756/3‑100 (MQ=21)
                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                          >  1:1507151/1‑100 (MQ=255)
                                            ggctggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGCAAAATCAATGGATTGCACA‑GC‑TAGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:840190/11‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:287713/2‑99 (MQ=17)
                                            ggCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAAATTCACTATATTGCAGA‑GA‑TTGC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  2:1641264/2‑100 (MQ=11)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:1350640/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                         >  1:2152905/2‑100 (MQ=21)
                                            agCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                        >  1:225696/2‑100 (MQ=18)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                        >  1:594792/1‑100 (MQ=255)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                       >  2:955987/1‑100 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                        >  1:1011896/1‑99 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑A‑AT‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                        >  1:2395303/1‑98 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAGA‑GT‑TCAC‑‑GTAGGCCTGTTAAGCGTAGCGCATCAGGCag                                                        >  2:537208/1‑99 (MQ=17)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCAAT‑CT‑CCGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                       >  2:87525/9‑98 (MQ=18)
                                               cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatc                                                      >  1:183495/8‑97 (MQ=25)
                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:1768415/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGATCG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:1910347/7‑100 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCAGG‑AG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:2240272/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCCAC‑AAAGTCTTGCAAAATCAATAAATTGCAGG‑AA‑CGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                                                     >  2:2515076/7‑96 (MQ=21)
                                                  ataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAGA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttg                                                   >  2:1107756/5‑100 (MQ=17)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                   <  1:537143/100‑1 (MQ=255)
                                                   taaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCAGAGAT‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttagc                                                  >  2:2112853/4‑97 (MQ=17)
                                                   tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                    >  2:1544002/1‑98 (MQ=255)
                                                   tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                  <  1:1879496/100‑1 (MQ=255)
                                                    gctgcgTTAACGCCTTATCCGGCCAAC‑AAAATCGTGCAAATACAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                          <  2:568403/88‑1 (MQ=255)
                                                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCCGGCATTtttgcctt                                               >  1:493300/1‑100 (MQ=255)
                                                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTACAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                               >  2:1162677/1‑100 (MQ=255)
                                                       gcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCAGA‑GA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAATttagcgttt                                              >  1:2598292/1‑100 (MQ=17)
                                                       gcgTTAACGCCTTATCCCGCCCAC‑ACCATCGTCCAAATCCAATACATCGCACA‑AT‑CTTC‑‑CTCGGCCTCACAACCAGACAGCATCAGGCAATcctgccttt                                              <  2:2294759/100‑1 (MQ=255)
                                                        cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                             >  1:312692/1‑100 (MQ=255)
                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                            <  2:1507955/100‑1 (MQ=255)
                                                          ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                           <  1:2527957/100‑1 (MQ=255)
                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                          >  1:2363959/1‑100 (MQ=255)
                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                          >  2:1454594/1‑100 (MQ=255)
                                                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                        <  2:1199266/100‑1 (MQ=255)
                                                                      ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                               >  1:1107012/1‑100 (MQ=255)
                                                                      ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                               >  1:1167656/1‑100 (MQ=255)
                                                                         gCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                            <  1:2032122/100‑1 (MQ=255)
                                                                          ccTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAACAGGCTTGTATGGc                           >  2:1451755/1‑100 (MQ=255)
                                                                           cTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                          >  1:600509/1‑100 (MQ=255)
                                                                             aC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                        <  1:465553/100‑1 (MQ=255)
                                                                             aC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                        >  2:2497113/1‑100 (MQ=255)
                                                                              c‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                       <  1:2397478/100‑1 (MQ=255)
                                                                                  aaTCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                    <  1:2107805/100‑1 (MQ=255)
                                                                                  aaTCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                    <  1:1619307/100‑1 (MQ=255)
                                                                                   aTCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatt                   <  2:292363/100‑1 (MQ=255)
                                                                                      gTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                <  2:736842/100‑1 (MQ=255)
                                                                                        gCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              >  1:1011631/1‑100 (MQ=255)
                                                                                        gCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              <  2:121432/100‑1 (MQ=255)
                                                                                         cAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa             >  1:2554022/1‑100 (MQ=255)
                                                                                         cAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa             >  2:841699/1‑100 (MQ=255)
                                                                                          aaaTTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt            <  2:554877/100‑1 (MQ=255)
                                                                                           aaTTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa           <  1:332112/100‑1 (MQ=255)
                                                                                            aTTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag          >  2:65050/1‑100 (MQ=255)
                                                                                             ttCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga         <  2:600507/100‑1 (MQ=255)
                                                                                             ttCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga         <  1:2036009/100‑1 (MQ=255)
                                                                                              tCAATACATTCCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAGCAGTCTTGTATGGCTTTTAATTATTCGATagag        <  2:1801475/100‑1 (MQ=255)
                                                                                               cAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga       <  2:2324699/100‑1 (MQ=255)
                                                                                                    |                                                                                                       
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑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: ‑