Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,978 A→C 17.1% intergenic (‑188/+120) 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,9780AC17.1% 594.9 / 102.1 216intergenic (‑188/+120)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (70/100);  new base C (29/2);  total (114/102)
Fisher's exact test for biased strand distribution p-value = 2.02e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.07e-02

Only 100 of 223 total aligned reads displayed.

GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705879‑3706077
                                                                                                   |                                                                                                         
gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                           <  1:2352659/100‑1 (MQ=255)
     aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                      <  2:2931050/100‑1 (MQ=255)
       tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                    >  2:2540644/1‑100 (MQ=255)
           gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                >  1:314620/1‑100 (MQ=255)
           gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                >  1:540358/1‑100 (MQ=255)
           gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGccgcatcc                                                                                                <  2:3152318/100‑7 (MQ=37)
            ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGccgagatcg                                                                                               >  1:2998049/1‑93 (MQ=37)
              ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                             <  1:498298/100‑1 (MQ=255)
                   gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                       >  1:1003086/1‑100 (MQ=255)
                      cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                                                                                    <  1:2152923/100‑1 (MQ=255)
                       tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                                   <  1:2684622/100‑1 (MQ=255)
                          aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                                >  2:435474/1‑100 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                                <  1:66700/99‑1 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                               <  1:1750611/100‑1 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                               >  2:554841/1‑100 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                               >  1:1188652/1‑100 (MQ=255)
                             cAGAGGCACAGCGCAACCGGGTATACAGATTGTTGACAAAGTGCGCTTAGTTTATGCCGGATGCAGCGTTAACGCCTTATCCGGCCTACAAAATAGTgaa                                                                             >  2:2696329/1‑98 (MQ=18)
                              agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaa                                                                            <  2:1755901/100‑1 (MQ=255)
                                  gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                        <  2:1984590/100‑1 (MQ=255)
                                   cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                       >  1:1203667/1‑100 (MQ=255)
                                   cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                       >  1:2312687/1‑100 (MQ=255)
                                   cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                       >  2:801970/1‑100 (MQ=255)
                                    acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                      <  1:2227/100‑1 (MQ=255)
                                     caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                     <  1:1490724/100‑1 (MQ=255)
                                     caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                     >  2:322727/1‑100 (MQ=255)
                                      aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                                    <  1:2260519/100‑1 (MQ=255)
                                       gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                   <  2:1717408/100‑1 (MQ=255)
                                           cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                               >  2:713589/1‑100 (MQ=255)
                                             ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCa                                                             >  1:592661/1‑100 (MQ=255)
                                              cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                           <  1:815204/100‑1 (MQ=255)
                                               gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                              <  1:2540662/97‑1 (MQ=255)
                                               gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGa                                                          <  1:2726085/100‑1 (MQ=255)
                                                 gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTACAAAATCGTGCAAATTCAATACATTGCAGaa                                                         >  2:1000886/1‑99 (MQ=255)
                                                    taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtttt                                                   >  1:96292/1‑100 (MQ=255)
                                                     aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTc                                                  <  2:2553549/100‑1 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                               <  1:1412898/100‑1 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑g                                               <  2:2911768/100‑1 (MQ=255)
                                                         ggcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAg                                            >  1:2269458/8‑100 (MQ=18)
                                                         aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAg                                            <  1:3241177/100‑1 (MQ=255)
                                                          ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGG                                           >  1:3036586‑M1/1‑100 (MQ=255)
                                                            gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCC                                         <  2:3057609‑M1/100‑1 (MQ=255)
                                                              tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTG                                       >  1:2947864‑M1/1‑100 (MQ=255)
                                                               gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGA                                      >  2:1112880‑M1/1‑100 (MQ=255)
                                                               gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGA                                      <  1:754944‑M1/100‑1 (MQ=255)
                                                                 cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATa                                    <  2:496698/100‑1 (MQ=255)
                                                                   aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCGTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAg                                  <  2:2192780/100‑1 (MQ=255)
                                                                   aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAg                                  <  1:566298/100‑1 (MQ=255)
                                                                       gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGa                              >  1:1636400/1‑100 (MQ=255)
                                                                        cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAg                             <  2:314618/100‑1 (MQ=255)
                                                                           gttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTCCGT‑‑GTAGGCCTGATAAGCGTAgcgc                          >  1:1381591/10‑100 (MQ=18)
                                                                            ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                         <  1:2319709‑M2/100‑2 (MQ=255)
                                                                            atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                        >  1:459037/9‑100 (MQ=18)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCatc                       <  1:1243170‑M2/100‑4 (MQ=255)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCatc                       <  1:1823918‑M2/100‑4 (MQ=255)
                                                                                ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCatcag                     <  1:2027936‑M2/100‑6 (MQ=255)
                                                                                gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                     >  1:3210370/5‑100 (MQ=17)
                                                                                gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAg                     >  1:877645/5‑100 (MQ=18)
                                                                                 taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                    >  2:141498/4‑100 (MQ=17)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                   >  2:2064445/4‑100 (MQ=21)
                                                                                 tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCatcagg                    <  2:56686‑M2/100‑7 (MQ=255)
                                                                                  ntgCCGGAAGCGGCGTTAACGCCTTATCCGGCCTACAAAATAGTGCAAATACAATACATTGCAGAATCTCC‑‑GTAGGCCGGATAAGCGGAGCGCatcaggc                   <  2:1296808‑M2/99‑8 (MQ=255)
                                                                                  aggccGGATAAGGCGTTAACGACGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCatcaggc                   <  2:1491592‑M2/98‑8 (MQ=255)
                                                                                  aTGCCGGATGCGGCGTTAACGCCTTATCCGGTCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATgcgctccgcttatcaggc                   <  2:1618080‑M1/100‑19 (MQ=255)
                                                                                  aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATgcgctccgcttatcaggc                   <  1:1992686‑M1/100‑19 (MQ=255)
                                                                                  aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATgcgctccgcttatcaggc                   <  1:1118506‑M1/100‑19 (MQ=255)
                                                                                  aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATgcgctccgcttatcaggc                   <  1:1186655‑M1/100‑19 (MQ=255)
                                                                                   tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                  >  1:1677161/1‑100 (MQ=255)
                                                                                   tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                  >  1:346974/1‑100 (MQ=255)
                                                                                   tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                  >  1:1574772/1‑100 (MQ=255)
                                                                                   tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                  >  1:2514116/1‑100 (MQ=255)
                                                                                   tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                  >  2:334523/1‑100 (MQ=255)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                 >  1:1546718/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                 >  1:1544029/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:628195/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:1954102/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:2496663/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                  >  2:478083/2‑100 (MQ=21)
                                                                                    gtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                 >  2:940961/10‑99 (MQ=21)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                 >  1:1073161/1‑100 (MQ=255)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                 >  1:2426375/1‑99 (MQ=18)
                                                                                     ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTGC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                <  1:73055/100‑1 (MQ=255)
                                                                                      cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATgcgctccgcttatcaggcctac               <  1:2065379‑M1/100‑23 (MQ=255)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt              >  2:125044/7‑100 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt              >  2:2154284/7‑100 (MQ=14)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt              >  2:2741940/7‑100 (MQ=17)
                                                                                        gataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcg            >  2:1654229/6‑97 (MQ=21)
                                                                                        gataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt             <  1:478130/95‑1 (MQ=21)
                                                                                        gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGACgcgctccgcttatcaggcctacga             <  1:41579‑M1/100‑25 (MQ=255)
                                                                                         ataaggcgTTTACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttta            >  2:852679/5‑99 (MQ=18)
                                                                                         aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg            <  2:1203661/100‑1 (MQ=255)
                                                                                         aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg            <  1:2862721/100‑1 (MQ=255)
                                                                                          tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc           <  2:592659/100‑1 (MQ=255)
                                                                                           gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc          <  1:1936026/100‑1 (MQ=255)
                                                                                            cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct         <  2:3131500/100‑1 (MQ=255)
                                                                                                gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc     >  2:1570183/1‑100 (MQ=255)
                                                                                                gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc     <  1:1725137/100‑1 (MQ=255)
                                                                                                gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATGCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc     <  1:2282249/100‑1 (MQ=255)
                                                                                                 ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc    <  2:1351307/100‑1 (MQ=255)
                                                                                                  tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa   >  2:1455542/1‑100 (MQ=255)
                                                                                                   aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt  >  2:1878284/1‑100 (MQ=255)
                                                                                                   |                                                                                                         
GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAT  >  NC_000913/3705879‑3706077

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: