Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,978 A→C 30.4% 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,9780AC30.4% 529.1 / 257.0 323intergenic (‑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 (114/89);  new base C (76/2);  total (231/92)
Fisher's exact test for biased strand distribution p-value = 2.39e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03

Only 100 of 325 total aligned reads displayed.

GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTG  >  NC_000913/3705879‑3706073
                                                                                                   |                                                                                                       
gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                         <  1:3546962/100‑1 (MQ=255)
  aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                       <  2:939979/100‑1 (MQ=255)
    tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                     <  1:1865137/100‑1 (MQ=255)
     aTTGCTGCCGGATAGCGCTGCACGCAGAGGCACAGCGCAACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                    >  1:2722073/1‑100 (MQ=255)
      ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                   <  2:2694274/100‑1 (MQ=255)
        gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                 >  1:2176118/1‑100 (MQ=255)
         ctgcCGGATGGCGCTGCACGCGGAGGCACAGCGCCACGGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                >  1:2363313/1‑100 (MQ=255)
         ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                >  2:1943922/1‑100 (MQ=255)
            ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                             >  2:1244450/1‑100 (MQ=255)
              ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                           >  2:2179141/1‑100 (MQ=255)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                          <  1:1244446/100‑1 (MQ=255)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                          <  2:4139052/100‑1 (MQ=255)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                          <  1:2237018/100‑2 (MQ=37)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTCTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                          <  2:3113597/100‑1 (MQ=255)
                   gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                     >  1:3143893/1‑100 (MQ=255)
                    cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                    <  2:3783812/100‑1 (MQ=255)
                     gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                   >  2:2650893/1‑100 (MQ=255)
                     gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                   >  1:759743/1‑100 (MQ=255)
                       tGCACGCAGAGGCACAGCGCCACCGCGTATACCGATTATTGACAACGTGCGCATTGTATATACAGGATCAGTCATTAACGCCTTATCCGGACTACaaact                                                                                 >  1:2932332/1‑98 (MQ=255)
                        gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                                >  1:1343124/1‑100 (MQ=255)
                        gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                                >  2:4289383/1‑100 (MQ=255)
                         cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                               >  1:152041/1‑100 (MQ=255)
                           cgatctggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                             <  1:556371/94‑1 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                             >  2:3311832/1‑100 (MQ=255)
                            gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                            <  2:1161768/100‑1 (MQ=255)
                            gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                            <  2:1562290/100‑1 (MQ=255)
                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                           <  2:4410172/100‑1 (MQ=255)
                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                           <  1:3809608/100‑1 (MQ=255)
                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                           <  2:1373877/100‑1 (MQ=255)
                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCCCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                           <  1:3082814/100‑1 (MQ=255)
                                 ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTgagatcg                                                                       >  2:556373/1‑94 (MQ=255)
                                  gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                      <  1:715987/100‑1 (MQ=255)
                                   cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                     >  1:3290612/1‑100 (MQ=255)
                                     caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                   <  2:1412288/100‑1 (MQ=255)
                                     caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                   >  2:854330/1‑100 (MQ=255)
                                         gcCACCGGTTATACAGATTGTTGACAAAGTGCCCTTTGTTTATGCCGGATGTGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                 <  2:3355360/98‑1 (MQ=255)
                                         gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                               >  2:2606558/1‑100 (MQ=255)
                                            aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                            <  1:2591134/100‑1 (MQ=255)
                                                ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑aa                                                    >  2:2171039/1‑100 (MQ=255)
                                                  tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCCGATTCGGCGTTAACGACTTATCCGGCCTACAAAATCGTGCAGATTCAATACATTGCAG‑‑‑AAtt                                                 >  1:1484584/1‑100 (MQ=18)
                                                    taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AAtttt                                               >  1:1074708/1‑100 (MQ=255)
                                                         ggcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAg                                        >  2:2017649/8‑100 (MQ=17)
                                                          ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAgg                                       <  1:247495/100‑1 (MQ=255)
                                                           tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGc                                      >  2:4266395/1‑100 (MQ=255)
                                                                    aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGc                             <  1:1976622/100‑1 (MQ=255)
                                                                           tttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAG‑‑‑GACCAT‑‑GTAGGCCTGATAAGCGTAgcgc                      >  2:209628/10‑100 (MQ=18)
                                                                           tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAgcgc                      <  2:2377869/100‑1 (MQ=255)
                                                                            ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                     <  2:1735579/100‑1 (MQ=255)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                   <  1:1690462/100‑1 (MQ=255)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                   >  2:2231739/1‑100 (MQ=255)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                   >  1:141880/1‑100 (MQ=255)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                   <  2:2062359/100‑1 (MQ=255)
                                                                               tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                  <  2:488108/100‑1 (MQ=255)
                                                                               tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                  >  1:2461817/1‑100 (MQ=255)
                                                                               tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCG‑‑‑AATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                 >  2:695297/6‑100 (MQ=18)
                                                                               tgtaggCCGGATAAGGCGTTCACGCCGCATCAGGCCTACAAAGACCTGCAAATTCAATATATTGCAG‑‑‑GATACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                  >  2:2981499/6‑100 (MQ=17)
                                                                                 taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATG‑‑‑GGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                >  2:2704805/4‑100 (MQ=17)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAG‑‑‑GAAACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAgg                >  1:2266878/4‑100 (MQ=18)
                                                                                 tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                >  1:3266858/1‑100 (MQ=255)
                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAG‑‑‑AAATCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc               >  1:3205087/3‑100 (MQ=18)
                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                >  2:4188952/3‑100 (MQ=21)
                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAAGCGTGCAAATTCAATACATTGCAT‑‑‑GGGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc               >  1:4033404/3‑100 (MQ=17)
                                                                                  aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc               >  2:380383/1‑100 (MQ=255)
                                                                                   ggctggataaggcgTTCACGCCGCATCCAGCCTACAAAATTGTGCAAAATCAATGGATTGCAC‑‑‑AGCTAGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa              >  2:2329378/11‑100 (MQ=25)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATG‑‑‑GGCCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa              >  2:3933407/2‑100 (MQ=17)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCAA‑‑‑TCTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCa              >  1:1574412/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAG‑‑‑AGATTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa              >  1:4098772/2‑100 (MQ=11)
                                                                                    gtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAG‑‑‑GAACGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa             >  2:3706644/10‑99 (MQ=21)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt             >  1:3624096/1‑100 (MQ=255)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt             >  1:944810/1‑100 (MQ=255)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCAT‑‑‑GAGCCAC‑‑GTAGTCCTGATAAGCGTAGCGCATCAGGCag             >  1:1008699/1‑99 (MQ=17)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCAC‑‑‑CATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa             >  2:168247/1‑99 (MQ=21)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAG‑‑‑AGATTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa             >  2:654224/1‑99 (MQ=17)
                                                                                    gCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa              >  1:3517417/1‑99 (MQ=17)
                                                                                     ccggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCG‑‑‑AATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt           >  2:2534255/9‑100 (MQ=18)
                                                                                     ccggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAG‑‑‑AAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt            >  1:1376719/9‑98 (MQ=18)
                                                                                     ccggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAG‑‑‑AAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt            >  1:1177795/9‑98 (MQ=18)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCAC‑‑‑CATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat            >  2:2852309/9‑98 (MQ=21)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAG‑‑‑GAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt           >  2:1286656/9‑100 (MQ=18)
                                                                                     ccGGATAAGGCGGTCACTCCGCATCCGACCTACAAATCATGCAAATCCAAAATATTGCAG‑‑‑GAGGTAT‑‑GTATGCCTGATAAGCGTAGCGCATCAGGCAtt            >  1:4302543/1‑99 (MQ=18)
                                                                                      cggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAG‑‑‑AGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt           >  2:4026745/8‑100 (MQ=17)
                                                                                      cggataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑AATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaatt           >  2:2769333/8‑97 (MQ=18)
                                                                                      cggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAG‑‑‑GAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttta          >  2:831849/8‑99 (MQ=21)
                                                                                      cggataaggcgTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt           >  2:183644/8‑100 (MQ=18)
                                                                                       ggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCG‑‑‑AATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc         >  2:3501717/7‑98 (MQ=21)
                                                                                       ggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCG‑‑‑AATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc         >  2:4089109/7‑98 (MQ=21)
                                                                                       ggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCG‑‑‑AATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc         >  1:3421876/7‑98 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAG‑‑‑GACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt          >  1:3852956/7‑100 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAG‑‑‑GACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt          >  1:2287956/7‑100 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAG‑‑‑GACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt          >  1:3330045/7‑100 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAA‑‑‑GATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct          >  2:320232/7‑96 (MQ=17)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCACTTAAATTCAATATATTGCAC‑‑‑CATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattc          >  1:2987483/7‑99 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAG‑‑‑GAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat         >  2:3028480/7‑98 (MQ=21)
                                                                                       ggataaggcgGTTACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAA‑‑‑TCTCCGT‑‑GTAGGCATGATAAGCGTAGCGAATCAGGccat            >  1:4190648/7‑95 (MQ=17)
                                                                                         aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg        <  2:2461829/100‑1 (MQ=255)
                                                                                             ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttttttttt    >  1:1803177/1‑95 (MQ=255)
                                                                                             ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt    <  2:3556852/100‑1 (MQ=255)
                                                                                             ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct     <  1:3746961/99‑1 (MQ=255)
                                                                                               cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg  >  1:1301512/1‑100 (MQ=255)
                                                                                               cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg  <  1:991962/100‑1 (MQ=255)
                                                                                                   |                                                                                                       
GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAG‑‑‑AATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTG  >  NC_000913/3705879‑3706073

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: