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.

GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑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)
                   gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                     >  1:3143893/1‑100 (MQ=255)
                    cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                                    <  2:3783812/100‑1 (MQ=255)
                     gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                   >  2:2650893/1‑100 (MQ=255)
                     gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                   >  1:759743/1‑100 (MQ=255)
                       tGCACGCAGAGGCACAGCGCCACCGCGTATACCGATTATTGACAACGTGCGCATTGTATATACAGGATCAGTCATTAACGCCTTATCCGGACTAC‑aaact                                                                                 >  1:2932332/1‑98 (MQ=255)
                        gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                >  1:1343124/1‑100 (MQ=255)
                        gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                >  2:4289383/1‑100 (MQ=255)
                         cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                               >  1:152041/1‑100 (MQ=255)
                           cgatctggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                             <  1:556371/94‑1 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                             >  2:3311832/1‑100 (MQ=255)
                            gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                            <  2:1161768/100‑1 (MQ=255)
                            gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                            <  2:1562290/100‑1 (MQ=255)
                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                           <  2:4410172/100‑1 (MQ=255)
                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                           <  1:3809608/100‑1 (MQ=255)
                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                           <  2:1373877/100‑1 (MQ=255)
                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCCCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                           <  1:3082814/100‑1 (MQ=255)
                                 ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTgagatcg                                                                       >  2:556373/1‑94 (MQ=255)
                                  gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                      <  1:715987/100‑1 (MQ=255)
                                   cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                     >  1:3290612/1‑100 (MQ=255)
                                     caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                   <  2:1412288/100‑1 (MQ=255)
                                     caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                   >  2:854330/1‑100 (MQ=255)
                                         gcCACCGGTTATACAGATTGTTGACAAAGTGCCCTTTGTTTATGCCGGATGTGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                 <  2:3355360/98‑1 (MQ=255)
                                         gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                               >  2:2606558/1‑100 (MQ=255)
                                            aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                            <  1:2591134/100‑1 (MQ=255)
                                                ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑aa                                                    >  2:2171039/1‑100 (MQ=255)
                                                  tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCCGATTCGGCGTTAACGACTTATCCGGCCTAC‑AAAATCGTGCAGATTCAATACATTGCA‑G‑‑‑AAt‑t                                                 >  1:1484584/1‑100 (MQ=18)
                                                    taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAt‑ttt                                               >  1:1074708/1‑100 (MQ=255)
                                                         ggcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAg                                        >  2:2017649/8‑100 (MQ=17)
                                                          ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAgg                                       <  1:247495/100‑1 (MQ=255)
                                                           tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGc                                      >  2:4266395/1‑100 (MQ=255)
                                                                    aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGc                             <  1:1976622/100‑1 (MQ=255)
                                                                           tttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑G‑‑‑GAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAgcgc                      >  2:209628/10‑100 (MQ=18)
                                                                           tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcgc                      <  2:2377869/100‑1 (MQ=255)
                                                                            ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                     <  2:1735579/100‑1 (MQ=255)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                   <  1:1690462/100‑1 (MQ=255)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                   >  2:2231739/1‑100 (MQ=255)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                   >  1:141880/1‑100 (MQ=255)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                   <  2:2062359/100‑1 (MQ=255)
                                                                               tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                  <  2:488108/100‑1 (MQ=255)
                                                                               tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                  >  1:2461817/1‑100 (MQ=255)
                                                                               tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑G‑‑‑AATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                 >  2:695297/6‑100 (MQ=18)
                                                                               tgtaggCCGGATAAGGCGTTCACGCCGCATCAGGCCTAC‑AAAGACCTGCAAATTCAATATATTGCA‑G‑‑‑GAT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                  >  2:2981499/6‑100 (MQ=17)
                                                                                 taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATG‑‑‑GGC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                >  2:2704805/4‑100 (MQ=17)
                                                                                 taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑G‑‑‑GAA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAgg                >  1:2266878/4‑100 (MQ=18)
                                                                                 tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                >  1:3266858/1‑100 (MQ=255)
                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCA‑G‑‑‑AAA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc               >  1:3205087/3‑100 (MQ=18)
                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGCA‑ATTCAAC‑TT‑‑‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                >  2:4188952/3‑100 (MQ=21)
                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAAGCGTGCAAATTCAATACATTGCA‑T‑‑‑GGG‑CCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc               >  1:4033404/3‑100 (MQ=17)
                                                                                  aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc               >  2:380383/1‑100 (MQ=255)
                                                                                   ggctggataaggcgTTCACGCCGCATCCAGCCTAC‑AAAATTGTGCAAAATCAATGGATTGCA‑C‑‑‑AGC‑TAGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa              >  2:2329378/11‑100 (MQ=25)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATG‑‑‑GGC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa              >  2:3933407/2‑100 (MQ=17)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑A‑‑‑TCT‑TCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCa              >  1:1574412/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAAATTCAATATATTGCA‑G‑‑‑AGA‑TTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa              >  1:4098772/2‑100 (MQ=11)
                                                                                    gtcggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCA‑G‑‑‑GAA‑CGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa             >  2:3706644/10‑99 (MQ=21)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt             >  1:3624096/1‑100 (MQ=255)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt             >  1:944810/1‑100 (MQ=255)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑T‑‑‑GAG‑CCAC‑‑GTAGTCCTGATAAGCGTAGCGCATCAGGCag             >  1:1008699/1‑99 (MQ=17)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑C‑‑‑CAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa             >  2:168247/1‑99 (MQ=21)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑G‑‑‑AGA‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa             >  2:654224/1‑99 (MQ=17)
                                                                                    gCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa              >  1:3517417/1‑99 (MQ=17)
                                                                                     ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑G‑‑‑AATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt           >  2:2534255/9‑100 (MQ=18)
                                                                                     ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑G‑‑‑AAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt            >  1:1376719/9‑98 (MQ=18)
                                                                                     ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑G‑‑‑AAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt            >  1:1177795/9‑98 (MQ=18)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑C‑‑‑CAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat            >  2:2852309/9‑98 (MQ=21)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑G‑‑‑GAG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt           >  2:1286656/9‑100 (MQ=18)
                                                                                     ccGGATAAGGCGGTCACTCCGCATCCGACCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑G‑‑‑GAG‑GTAT‑‑GTATGCCTGATAAGCGTAGCGCATCAGGCAtt            >  1:4302543/1‑99 (MQ=18)
                                                                                      cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑G‑‑‑AGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt           >  2:4026745/8‑100 (MQ=17)
                                                                                      cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑‑‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaatt           >  2:2769333/8‑97 (MQ=18)
                                                                                      cggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑G‑‑‑GAG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttta          >  2:831849/8‑99 (MQ=21)
                                                                                      cggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt           >  2:183644/8‑100 (MQ=18)
                                                                                       ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑G‑‑‑AATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc         >  2:3501717/7‑98 (MQ=21)
                                                                                       ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑G‑‑‑AATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc         >  2:4089109/7‑98 (MQ=21)
                                                                                       ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑G‑‑‑AATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc         >  1:3421876/7‑98 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑G‑‑‑GAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt          >  1:3852956/7‑100 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑G‑‑‑GAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt          >  1:2287956/7‑100 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑G‑‑‑GAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt          >  1:3330045/7‑100 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑A‑‑‑GAT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgct          >  2:320232/7‑96 (MQ=17)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCACTTAAATTCAATATATTGCA‑C‑‑‑CAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattc          >  1:2987483/7‑99 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑G‑‑‑GAG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat         >  2:3028480/7‑98 (MQ=21)
                                                                                       ggataaggcgGTTACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCA‑A‑‑‑TCT‑CCGT‑‑GTAGGCATGATAAGCGTAGCGAATCAGGccat            >  1:4190648/7‑95 (MQ=17)
                                                                                         aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg        <  2:2461829/100‑1 (MQ=255)
                                                                                             ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttttttttt    >  1:1803177/1‑95 (MQ=255)
                                                                                             ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt    <  2:3556852/100‑1 (MQ=255)
                                                                                             ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct     <  1:3746961/99‑1 (MQ=255)
                                                                                               cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg  >  1:1301512/1‑100 (MQ=255)
                                                                                               cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg  <  1:991962/100‑1 (MQ=255)
                                                                                                   |                                                                                                       
GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑G‑‑‑AAT‑TTTC‑‑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: ‑