Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 21.7% intergenic (‑194/+114) 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,9840TG21.7% 556.2 / 152.8 220intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (73/99);  new base G (44/4);  total (117/103)
Fisher's exact test for biased strand distribution p-value = 2.37e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01

Only 100 of 222 total aligned reads displayed.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGT  >  NC_000913/3705886‑3706081
                                                                                                  |                                                                                                    
tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                     >  1:771545/1‑100 (MQ=255)
    gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                 >  1:314620/1‑100 (MQ=255)
    accGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                 >  1:333656/2‑100 (MQ=255)
       ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                              <  1:498298/100‑1 (MQ=255)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                        >  1:1003086/1‑100 (MQ=255)
               cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaaa                                                                                     <  1:818695/100‑1 (MQ=255)
                tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                    <  1:2684622/100‑1 (MQ=255)
                tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                    >  1:2521787/1‑100 (MQ=255)
                 gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                   >  2:1917645/1‑100 (MQ=255)
                   aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                                 >  2:3156420/1‑100 (MQ=255)
                   aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                                 >  2:435474/1‑100 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                >  1:1188652/1‑100 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                <  1:1750611/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                >  2:554841/1‑100 (MQ=255)
                      cAGAGGCACAGCGCTACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                              <  2:1323579/100‑1 (MQ=255)
                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                              <  2:2039549/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)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                        >  1:1203667/1‑100 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                        >  1:2312687/1‑100 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                        >  2:801970/1‑100 (MQ=255)
                             acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCaa                                                                       <  1:2227/100‑1 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                      <  1:942234/100‑1 (MQ=255)
                               aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCaa                                                                       <  1:1888486/98‑1 (MQ=255)
                               aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                                     <  1:2260519/100‑1 (MQ=255)
                               aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                                     <  1:549344/100‑1 (MQ=255)
                                 cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACa                                                                   >  2:1110791/1‑100 (MQ=255)
                                    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                >  2:1658610/1‑100 (MQ=255)
                                    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                >  2:713589/1‑100 (MQ=255)
                                       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGAGGCGTTAACGCCTTATCCGGCCTAC‑AAAATAGTACAAGTTCAATATATTGCA‑g                                                            >  2:2549675/1‑100 (MQ=17)
                                        gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga                                                           <  1:2726085/100‑1 (MQ=255)
                                           tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAt‑t                                                       >  2:1187417/1‑100 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCg                                                   <  2:2911768/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCg                                                   <  1:1412898/100‑1 (MQ=255)
                                                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTa                                                 >  1:996289/1‑100 (MQ=255)
                                                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTa                                                 >  1:71821/1‑100 (MQ=255)
                                                    tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGC                                              <  1:2988547‑M1/100‑1 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCC                                             <  1:2453805‑M1/100‑1 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCC                                             >  1:2298712‑M1/1‑100 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCC                                             <  1:42093‑M1/100‑1 (MQ=255)
                                                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGA                                          <  1:754944‑M1/100‑1 (MQ=255)
                                                           aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCtaataa                                       >  2:2980698/1‑100 (MQ=255)
                                                            aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCGTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAg                                      <  2:2192780/100‑1 (MQ=255)
                                                            aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAg                                      <  1:566298/100‑1 (MQ=255)
                                                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCgg                                   >  2:1209692/1‑100 (MQ=255)
                                                                gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGa                                  >  1:1636400/1‑100 (MQ=255)
                                                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAg                                 <  2:314618/100‑1 (MQ=255)
                                                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAg                                 <  1:2768602/100‑1 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAgcg                               >  2:2536623/1‑100 (MQ=255)
                                                                    gttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCA‑ATCT‑CCGTGTAGGCCTGATAAGCGTAgcgc                              >  1:1381591/10‑100 (MQ=18)
                                                                     atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCAt                            >  1:459037/9‑100 (MQ=18)
                                                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCatc                           <  1:1243170‑M2/100‑4 (MQ=255)
                                                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCatc                           <  1:1823918‑M2/100‑4 (MQ=255)
                                                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCatc                           >  2:1255170‑M2/1‑97 (MQ=255)
                                                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCatcag                         <  1:2027936‑M2/100‑6 (MQ=255)
                                                                          taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAgg                        >  2:141498/4‑100 (MQ=17)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGc                       >  2:2064445/4‑100 (MQ=21)
                                                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCatcagg                        <  2:56686‑M2/100‑7 (MQ=255)
                                                                           aggccGGATAAGGCGTTAACGACGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCatcaggc                       <  2:1491592‑M2/98‑8 (MQ=255)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGA‑TCATGTAGGCCTGATAAGCGTAGCGCATCAGGc                       >  1:2929041/3‑100 (MQ=14)
                                                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATgcgctccgcttatcaggc                       <  1:1118506‑M1/100‑19 (MQ=255)
                                                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATgcgctccgcttatcaggc                       <  1:1186655‑M1/100‑19 (MQ=255)
                                                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATgcgctccgcttatcaggc                       <  1:1992686‑M1/100‑19 (MQ=255)
                                                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                      >  1:1677161/1‑100 (MQ=255)
                                                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                      >  1:1574772/1‑100 (MQ=255)
                                                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                      >  1:651762/1‑100 (MQ=255)
                                                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                      >  2:334523/1‑100 (MQ=255)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                     >  1:1546718/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                     >  1:1544029/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                      >  2:628195/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                      >  2:478083/2‑100 (MQ=21)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GGAA‑ACACGTAGGCCTGATAAGCGAAGCGCATCAGGCa                      >  1:2505895/2‑100 (MQ=18)
                                                                             gtcggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCA‑GGAA‑CGATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                     >  2:940961/10‑99 (MQ=21)
                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                     >  1:1073161/1‑100 (MQ=255)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                     >  1:2377400/1‑99 (MQ=25)
                                                                               cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaatc                   >  2:1718712/8‑97 (MQ=21)
                                                                               cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATgcgctccgcttatcaggcctac                   <  1:2065379‑M1/100‑23 (MQ=255)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                  >  2:3245867/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGA‑TCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                  >  2:2154284/7‑100 (MQ=14)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGA‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                  >  2:2741940/7‑100 (MQ=17)
                                                                                 gataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                 <  1:478130/95‑1 (MQ=21)
                                                                                  ataaggcgTTTACGCCTTATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGGGC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttta                >  2:852679/5‑99 (MQ=18)
                                                                                   tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc               <  2:592659/100‑1 (MQ=255)
                                                                                    gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc              <  2:306408/100‑1 (MQ=255)
                                                                                    gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc              >  2:1652162/1‑100 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct             >  1:1487863/1‑100 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct             <  2:3131500/100‑1 (MQ=255)
                                                                                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt            >  1:2129788/1‑100 (MQ=255)
                                                                                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt            >  2:513034/1‑100 (MQ=255)
                                                                                       gcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATgcgctccgcttatcaggcctacgaaaattc           <  1:3235040‑M1/100‑31 (MQ=255)
                                                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc         <  2:1807535/100‑1 (MQ=255)
                                                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc         <  1:1725137/100‑1 (MQ=255)
                                                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc         >  2:1570183/1‑100 (MQ=255)
                                                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa       >  2:1455542/1‑100 (MQ=255)
                                                                                            aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑TAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa       >  1:317301/1‑99 (MQ=255)
                                                                                            aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt      <  1:388807/100‑1 (MQ=255)
                                                                                              cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa    >  2:244358/1‑100 (MQ=255)
                                                                                               gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg   <  2:1595774/100‑1 (MQ=255)
                                                                                                ccTTATCCGGCCTAC‑AAAATCGATCAAATACAATAAATTGCA‑GAAT‑ATTCGTAGGCCTGATAAGCGGAGCGCATCAAGCATTTTTGCCTTTGCTAtta    >  2:2597756/1‑96 (MQ=18)
                                                                                                  |                                                                                                    
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGT  >  NC_000913/3705886‑3706081

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