Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 20.2% intergenic (‑181/+127) 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,9710CA20.2% 521.9 / 109.1 187intergenic (‑181/+127)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (61/88);  new base A (34/4);  total (95/92)
Fisher's exact test for biased strand distribution p-value = 3.75e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

Only 100 of 207 total aligned reads displayed.

AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705872‑3706070
                                                                                                   |                                                                                                       
 gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                        >  2:2227/1‑100 (MQ=255)
 gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                        >  2:1936014/1‑100 (MQ=255)
 gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                        >  2:1722095/1‑100 (MQ=255)
  ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                       <  2:2590641/100‑1 (MQ=255)
    aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATAGAGATTGTTGTCAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                     >  2:2282233/1‑100 (MQ=255)
     aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                     <  1:1673237/99‑1 (MQ=255)
     aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                    >  2:3241149/1‑100 (MQ=255)
     aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                    <  1:244426/100‑1 (MQ=255)
       gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                  <  1:2265874/100‑1 (MQ=255)
       gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                  <  1:2352659/100‑1 (MQ=255)
        aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                 <  2:1306513/100‑1 (MQ=255)
         aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                <  2:361079/100‑1 (MQ=255)
           tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGAGTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                              <  1:588548/100‑1 (MQ=255)
            aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                             <  1:2317900/100‑1 (MQ=255)
              tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                           >  2:2540644/1‑100 (MQ=255)
                ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                         >  1:222900/1‑100 (MQ=255)
                  gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                       >  1:540358/1‑100 (MQ=255)
                  gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                       >  1:314620/1‑100 (MQ=255)
                   ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGccgagatcg                                                                                      >  1:2998049/1‑93 (MQ=37)
                     ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                    <  1:1655425/100‑1 (MQ=255)
                     ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                    <  1:498298/100‑1 (MQ=255)
                     ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                    <  2:314862/100‑1 (MQ=37)
                          gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                              >  1:1003086/1‑100 (MQ=255)
                          gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                              <  2:2910889/100‑1 (MQ=255)
                              tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                          >  1:2521787/1‑100 (MQ=255)
                                 aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                       >  2:3156420/1‑100 (MQ=255)
                                 aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                       >  2:435474/1‑100 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                      <  1:2084365/100‑1 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                      <  1:1750611/100‑1 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                      >  1:1188652/1‑100 (MQ=255)
                                   gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                     <  1:1396191/100‑1 (MQ=255)
                                   gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                     <  2:1846287/100‑1 (MQ=255)
                                    cAGAGGCACAGCGCAACCGGGTATACAGATTGTTGACAAAGTGCGCTTAGTTTATGCCGGATGCAGCGTTAACGCCTTATCCGGCCTAC‑AAAATAGTgaa                                                                    >  2:2696329/1‑98 (MQ=18)
                                      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                                                            <  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)
                                              gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                          <  2:756339/100‑1 (MQ=255)
                                                  cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                      >  2:713589/1‑100 (MQ=255)
                                                  cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                      >  2:1658610/1‑100 (MQ=255)
                                                    ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                    <  1:1722107/100‑1 (MQ=255)
                                                     cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAg                                                   <  1:815204/100‑1 (MQ=255)
                                                        gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTAC‑AAAATCGTGCAAATTCAATACATTGCAGaa                                                 >  2:1000886/1‑99 (MQ=255)
                                                           taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAt‑ttt                                            >  1:96292/1‑100 (MQ=255)
                                                            aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTc                                           <  2:2553549/100‑1 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑g                                        <  1:554843/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)
                                                                   gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCC                                  <  2:3057609‑M1/100‑1 (MQ=255)
                                                                   gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCC                                  <  1:42093‑M1/100‑1 (MQ=255)
                                                                     tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTG                                >  1:2947864‑M1/1‑100 (MQ=255)
                                                                     tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTG                                >  2:679678‑M1/1‑100 (MQ=255)
                                                                      gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGA                               >  2:1112880‑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)
                                                                         aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCtaataa                            >  2:2980698/1‑100 (MQ=255)
                                                                                   atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                 >  1:459037/9‑100 (MQ=18)
                                                                                    tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCAt                 >  1:2764559/8‑100 (MQ=18)
                                                                                     gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCatc                >  2:1255170‑M2/1‑97 (MQ=255)
                                                                                     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:2027936‑M2/100‑6 (MQ=255)
                                                                                       ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCatcag              >  1:1209698‑M2/1‑95 (MQ=255)
                                                                                       gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCAGAGA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  1:3210370/5‑100 (MQ=17)
                                                                                       gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAg              >  1:877645/5‑100 (MQ=18)
                                                                                        taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAGAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:141498/4‑100 (MQ=17)
                                                                                        tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCatcagg             <  1:1609638‑M2/100‑7 (MQ=255)
                                                                                         aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGAAA‑CCTT‑‑GTAGGCCTGATAAGCGTACCGCATCAGGc            >  2:2132410/3‑100 (MQ=18)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc            >  2:1433463/3‑100 (MQ=25)
                                                                                         aTGCCGGATGCGGCGTTAACGCCTTATCCGGTCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATgcgctccgcttatcaggc            <  2:1618080‑M1/100‑19 (MQ=255)
                                                                                         aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATgcgctccgcttatcaggc            <  1:1118506‑M1/100‑19 (MQ=255)
                                                                                         aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATgcgctccgcttatcaggc            <  1:739240‑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:346974/1‑100 (MQ=255)
                                                                                          ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt          >  1:758786/2‑100 (MQ=18)
                                                                                          ggCCGGATAAGGCGTTTACGACGCATCCGACCTAC‑AAAATCGTGCTAATTCAATATATTGCAGAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:1523100/2‑100 (MQ=17)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:2496663/2‑100 (MQ=18)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTGGGCT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:478083/2‑100 (MQ=21)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCAGGAA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa           >  1:2719297/2‑100 (MQ=18)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCAGGAA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa           >  1:2505895/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)
                                                                                           gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  2:2926657/1‑99 (MQ=18)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa          >  1:2377400/1‑99 (MQ=25)
                                                                                            ccGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt         >  2:2161661/1‑98 (MQ=18)
                                                                                            ccGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCATGAG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt         >  2:62520/1‑98 (MQ=17)
                                                                                                aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  2:1203661/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)
                                                                                                   cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  <  2:3131500/100‑1 (MQ=255)
                                                                                                   |                                                                                                       
AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705872‑3706070

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