Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,970 G→A 19.8% intergenic (‑180/+128) 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,9700GA19.8% 450.9 / 92.7 187intergenic (‑180/+128)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (62/88);  new base A (33/4);  total (95/92)
Fisher's exact test for biased strand distribution p-value = 8.24e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-02

Only 100 of 206 total aligned reads displayed.

CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069
                                                                                                   |                                                                                                       
 aGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGACGGATgc                                                                                                        <  1:1704478/100‑1 (MQ=255)
  gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                       >  2:2485483/1‑100 (MQ=255)
  gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                       >  2:1936014/1‑100 (MQ=255)
  gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                       >  2:1722095/1‑100 (MQ=255)
    gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                     >  1:2372596/1‑100 (MQ=255)
    gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGACACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                     >  2:1419986/1‑100 (MQ=255)
      aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                    <  1:1673237/99‑1 (MQ=255)
      aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                   >  2:3241149/1‑100 (MQ=255)
        gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                 <  1:2265874/100‑1 (MQ=255)
          aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                               <  2:361079/100‑1 (MQ=255)
             aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                            <  2:2931050/100‑1 (MQ=255)
             aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                            <  1:2317900/100‑1 (MQ=255)
                 ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                        >  1:222900/1‑100 (MQ=255)
                   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                      >  1:314620/1‑100 (MQ=255)
                   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                      >  2:2698935/1‑100 (MQ=255)
                   accGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                      >  1:333656/2‑100 (MQ=255)
                    ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGccgagatcg                                                                                     >  1:2998049/1‑93 (MQ=37)
                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                   <  1:1655425/100‑1 (MQ=255)
                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                   <  2:314862/100‑1 (MQ=37)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                             >  1:1003086/1‑100 (MQ=255)
                                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                      >  2:3156420/1‑100 (MQ=255)
                                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                      >  2:435474/1‑100 (MQ=255)
                                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                      >  2:2660637/1‑100 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                     <  1:2084365/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                     >  1:1188652/1‑100 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                    <  2:1846287/100‑1 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                    <  1:1396191/100‑1 (MQ=255)
                                     cAGAGGCACAGCGCTACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                   <  2:1323579/100‑1 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                   <  2:2251073/100‑1 (MQ=255)
                                       gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                 <  1:2272759/100‑1 (MQ=255)
                                        acgCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCTGCCTAC‑AAAATCGTGCAAAt                                                                <  2:1976780/98‑1 (MQ=255)
                                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                              <  2:1984590/100‑1 (MQ=255)
                                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                              <  2:559423/100‑1 (MQ=255)
                                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                              <  1:1557699/100‑1 (MQ=255)
                                            acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCaa                                                            >  1:1983299/1‑100 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                           <  2:2152544/100‑1 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                           <  1:1490724/100‑1 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                           >  1:1807547/1‑100 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCaa                                                            <  1:1888486/98‑1 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                           <  1:1255176/99‑1 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                          <  1:549344/100‑1 (MQ=255)
                                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                         <  2:756339/100‑1 (MQ=255)
                                                cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACa                                                        >  2:1110791/1‑100 (MQ=255)
                                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                     >  2:713589/1‑100 (MQ=255)
                                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAg                                                  <  1:815204/100‑1 (MQ=255)
                                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGAGGCGTTAACGCCTTATCCGGCCTAC‑AAAATAGTACAAGTTCAATATATTGCAg                                                  >  2:2549675/1‑100 (MQ=17)
                                                        ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGaa                                                >  1:515116/1‑100 (MQ=255)
                                                         gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTAC‑AAAATCGTGCAAATTCAATACATTGCAGaa                                                >  2:1000886/1‑99 (MQ=255)
                                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAt‑t                                             >  2:1187417/1‑100 (MQ=255)
                                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAt‑ttt                                           >  1:96292/1‑100 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑g                                       <  2:1585909/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑g                                       <  2:258807/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑g                                       <  2:2911768/100‑1 (MQ=255)
                                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTa                                     >  1:71821/1‑100 (MQ=255)
                                                                 aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAg                                    <  2:56786/100‑1 (MQ=255)
                                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGG                                   <  1:120160‑M1/100‑1 (MQ=255)
                                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGG                                   >  1:3036586‑M1/1‑100 (MQ=255)
                                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGC                                  >  1:450165‑M1/1‑100 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCC                                 <  1:2453805‑M1/100‑1 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCC                                 <  2:3057609‑M1/100‑1 (MQ=255)
                                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTG                               >  2:679678‑M1/1‑100 (MQ=255)
                                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTG                               >  1:2947864‑M1/1‑100 (MQ=255)
                                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTG                               <  1:1419994‑M1/99‑1 (MQ=255)
                                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATa                            <  2:496698/100‑1 (MQ=255)
                                                                           aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCGTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAg                          <  2:2192780/100‑1 (MQ=255)
                                                                           aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAg                          <  1:566298/100‑1 (MQ=255)
                                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCgg                       >  2:1209692/1‑100 (MQ=255)
                                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                   >  2:2536623/1‑100 (MQ=255)
                                                                                     tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCAt                >  1:2764559/8‑100 (MQ=18)
                                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCatc               <  1:1823918‑M2/100‑4 (MQ=255)
                                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCatca              >  1:2578396‑M2/1‑96 (MQ=255)
                                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCatcag             >  1:1209698‑M2/1‑95 (MQ=255)
                                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCatcag             <  1:2027936‑M2/100‑6 (MQ=255)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCAGAGA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  1:3210370/5‑100 (MQ=17)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAg             >  1:877645/5‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGGAG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:2064445/4‑100 (MQ=21)
                                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCatcagg            <  1:1609638‑M2/100‑7 (MQ=255)
                                                                                          ntgCCGGAAGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATAGTGCAAATACAATACATTGCAGAAT‑CTCC‑‑GTAGGCCGGATAAGCGGAGCGCatcaggc           <  2:1296808‑M2/99‑8 (MQ=255)
                                                                                          aggccGGATAAGGCGTTAACGACGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCatcaggc           <  2:1491592‑M2/98‑8 (MQ=255)
                                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATgcgctccgcttatcaggc           <  1:739240‑M1/100‑19 (MQ=255)
                                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATgcgctccgcttatcaggc           <  1:1118506‑M1/100‑19 (MQ=255)
                                                                                          aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGTTT‑ATTG‑‑GTAGGCCTGATGCGCGGCGCTTATCAGGc           <  1:2180039/100‑1 (MQ=21)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  1:651762/1‑100 (MQ=255)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  1:2514116/1‑100 (MQ=255)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt         >  1:758786/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTGGGCT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:478083/2‑100 (MQ=21)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCAGGAA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa          >  1:2719297/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  2:3142817/2‑100 (MQ=18)
                                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt         >  1:1073161/1‑100 (MQ=255)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt        >  1:2302856/1‑100 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTGGGCT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  1:2426375/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACCAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  1:2341314/1‑99 (MQ=21)
                                                                                             ccGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCATGAG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt        >  2:62520/1‑98 (MQ=17)
                                                                                               ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATgcgctccgcttatcaggcctacg      <  1:1738871‑M1/100‑24 (MQ=255)
                                                                                                gATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGACgcgctccgcttatcaggcctacga     <  1:41579‑M1/100‑25 (MQ=255)
                                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg    <  1:493807/100‑1 (MQ=255)
                                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg    <  1:2862721/100‑1 (MQ=255)
                                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg    >  2:898528/1‑100 (MQ=255)
                                                                                                  tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc   >  1:1544192/1‑100 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc  >  2:1652162/1‑100 (MQ=255)
                                                                                                   |                                                                                                       
CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069

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: ‑