Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 206,867 A→C E581A (GAG→GCG)  dnaE → DNA polymerase III alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913206,8670AC94.5% 525.7 / 5.7 166E581A (GAG→GCG) dnaEDNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (3/6);  new base C (89/67);  total (93/73)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

GCCGGACTGGTGAAGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGAGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGGCGGTATTCCAGCTTGAATC  >  NC_000913/206770‑206963
                                                                                                 |                                                                                                
 ccGGGCTGGTGAAGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCgcc                                                                                               >  2:3382048/1‑100 (MQ=255)
 ccGGACTGGTGAAGTTCGACTTCCTTGGTTTGCGTACGCTGACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGACGAATAGCGCgcc                                                                                               >  2:3200412/1‑100 (MQ=255)
 ccGGACTGGTGAAGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCgcc                                                                                               >  1:2356987/1‑100 (MQ=255)
  cGGACTGGTGAAGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCgccg                                                                                              <  2:1048597/100‑1 (MQ=255)
     aCTGGTGAAGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCgccgccn                                                                                           <  1:1495512/100‑2 (MQ=255)
       tGGTGAAGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCt                                                                                         <  2:1298895/100‑1 (MQ=255)
       tGGTGAAGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCt                                                                                         <  1:3504387/100‑1 (MQ=255)
        ggTGAAGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTg                                                                                        >  1:2828975/1‑100 (MQ=255)
        ggTGAAGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTg                                                                                        >  2:2253855/1‑100 (MQ=255)
        ggTGAAGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTg                                                                                        >  2:3774321/1‑100 (MQ=255)
        ggTGAAGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTg                                                                                        >  2:844810/1‑100 (MQ=255)
         gTGAAGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTgg                                                                                       >  2:1176220/1‑100 (MQ=255)
         gTGAAGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTgg                                                                                       >  1:2382174/1‑100 (MQ=255)
          tGAAGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGa                                                                                      <  1:3451120/100‑1 (MQ=255)
           gAAGTTCGACTTCCCTGGTTTGCGTACGCGCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGat                                                                                     <  1:3636711/100‑1 (MQ=255)
            aaGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGata                                                                                    <  1:973280/100‑1 (MQ=255)
             aGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGatat                                                                                   >  2:2372971/1‑100 (MQ=255)
              gTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATc                                                                                  >  1:2753651/1‑100 (MQ=255)
              gTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATc                                                                                  >  1:1203664/1‑100 (MQ=255)
                 cGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCt                                                                               >  1:2613175/1‑100 (MQ=255)
                 cGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGCCGAAGAATGGCGCGCCGCCGCTGGATATCGCt                                                                               >  2:2263370/1‑100 (MQ=255)
                   aCTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGc                                                                             >  1:1811271/1‑100 (MQ=255)
                      tCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGAt                                                                          >  2:1372803/1‑100 (MQ=255)
                      tCCGTGGATTGCGTACGCTCACCATAAGCAACTGGGCGCTGGATATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGAt                                                                          >  2:1806012/1‑100 (MQ=255)
                      nccTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCCACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATGGCTGCGAt                                                                          >  1:3252282/2‑100 (MQ=255)
                       ccTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATc                                                                         >  1:1697843/1‑100 (MQ=255)
                        cTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATcc                                                                        >  1:4781/1‑100 (MQ=255)
                        cTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGAGCCGCCGCTGGATATCGCTGCGATcc                                                                        <  1:380356/100‑1 (MQ=255)
                         ttGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATccc                                                                       <  1:2618818/100‑1 (MQ=255)
                         ttGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATccc                                                                       <  2:3579060/100‑1 (MQ=255)
                          tGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGAGCCGCCGCTGGATATCGCTGCGATCCCg                                                                      <  2:1241697/100‑1 (MQ=255)
                            gTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCt                                                                    <  2:2355665/100‑1 (MQ=255)
                            gTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCt                                                                    >  2:3288803/1‑100 (MQ=255)
                            gTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCACTGGATATCGCTGCGATCCCGCt                                                                    >  1:2878479/1‑100 (MQ=255)
                             tttGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGAGCCGCCGCTGGATATCGCTGCGATCCCGCTg                                                                   <  1:507046/100‑1 (MQ=255)
                             tttGCGGACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTg                                                                   <  2:3054302/100‑1 (MQ=255)
                                gCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGgat                                                                >  1:91075/1‑100 (MQ=255)
                                  gTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGgatga                                                              <  2:3237114/100‑1 (MQ=255)
                                  gTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGgatga                                                              >  1:1606735/1‑100 (MQ=255)
                                     cGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATaa                                                           <  2:1131295/100‑1 (MQ=255)
                                     cGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATaa                                                           >  2:3437555/1‑100 (MQ=255)
                                     cGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATaa                                                           >  2:912287/1‑100 (MQ=255)
                                       cTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGCTATCGCTGCGATCCCGCTGGATGATAAGa                                                         >  1:771743/1‑100 (MQ=255)
                                       cTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGa                                                         <  1:3950216/100‑1 (MQ=255)
                                         cACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTTGATATCGCTGCGATCCCGCTGGATGATAAGaaa                                                       <  2:580082/100‑1 (MQ=255)
                                         cACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGaaa                                                       <  2:3952207/100‑1 (MQ=255)
                                          aCCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGaaaa                                                      <  1:341833/100‑1 (MQ=255)
                                          aCCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGaaaa                                                      >  1:1918742/1‑100 (MQ=255)
                                             atcatcAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCt                                                   <  1:2990730/100‑1 (MQ=255)
                                               catcaACTGGGCGCTGGAGATGATCAACTAGCGGCGGGCGAAGAAGGGCGCTCCGCCGCTGGATACCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTc                                                 <  2:3528616/100‑1 (MQ=255)
                                               catcaACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTc                                                 <  1:3319698/100‑1 (MQ=255)
                                               catcaACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGAGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTc                                                 <  2:1696986/100‑1 (MQ=255)
                                               cagccaCTGGGCGCTGGAGATGATCAACAACCGGCGGGCGAAGAATGGCGCGCCGCTGCTGGATCTCCCTGCGATCCCGCGGGATGATATGAAAAGCTTc                                                 >  2:653046/6‑100 (MQ=255)
                                                atcaACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCg                                                >  1:3940780/1‑100 (MQ=255)
                                                atcaACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCg                                                <  1:1648320/100‑1 (MQ=255)
                                                 tcaACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGa                                               >  2:1853282/1‑100 (MQ=255)
                                                  caACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGAc                                              <  1:1780216/100‑1 (MQ=255)
                                                  caACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGAc                                              <  1:1861258/100‑1 (MQ=255)
                                                  caACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGAc                                              >  1:2989805/1‑100 (MQ=255)
                                                   aaCTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACa                                             >  2:453428/1‑100 (MQ=255)
                                                   aaCTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACa                                             <  1:2430442/100‑1 (MQ=255)
                                                    aCTGGGCGGTGGAGATGATCAACAAGCGGCGGGCGAAGAAGGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACAt                                            <  2:1530938/100‑1 (MQ=255)
                                                    aCTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACAt                                            >  2:1475086/1‑100 (MQ=255)
                                                    aCTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACAt                                            <  1:3140301/100‑1 (MQ=255)
                                                       gggCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACAtgct                                         <  1:1769971/100‑1 (MQ=255)
                                                       gggCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACAtgct                                         >  2:282596/1‑100 (MQ=255)
                                                              ggAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGc                                  >  2:2165382/1‑100 (MQ=255)
                                                              ggAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGc                                  >  1:2634202/1‑100 (MQ=255)
                                                              ggAGATGATCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGc                                  >  1:2157582/1‑100 (MQ=255)
                                                                 gatgatCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTGCGACATGCTGCAACGCACg                               >  2:2442044/1‑100 (MQ=255)
                                                                  gatctCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCgg                              <  1:1281751/96‑1 (MQ=255)
                                                                   ngatCAACAAGCGGCGGGCGAAGAATGGCGAGCCGCCGCTGGATATCGCTGCGATACCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGa                             >  1:1666275/2‑100 (MQ=255)
                                                                    gatCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGan                            <  1:4000437/100‑2 (MQ=255)
                                                                    gatCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGaa                            <  2:261794/100‑1 (MQ=255)
                                                                    gatCAACAAGCGGCGGGCGAAGAATGGCGAGCCGCCGCTTGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGaa                            >  2:1188312/1‑100 (MQ=255)
                                                                    gatCAACAAGCGACGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATACCGCTGGATGATAAGAAAAGCTTCGACATGATGAAACGCTCGGaa                            >  2:771321/1‑100 (MQ=255)
                                                                     atCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGaaa                           <  1:1952669/100‑1 (MQ=255)
                                                                      tCAACAAGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAAc                          >  1:3279611/1‑100 (MQ=255)
                                                                        aacaaGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATCTCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCa                        >  2:3072616/1‑100 (MQ=255)
                                                                        aacaaGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCa                        >  1:1468351/1‑100 (MQ=255)
                                                                         acaaGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCAc                       <  2:2124355/100‑1 (MQ=255)
                                                                          caaGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCAcg                      >  1:3577118/1‑100 (MQ=255)
                                                                            aGCGGCGGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCAcggc                    <  2:2482781/100‑1 (MQ=255)
                                                                             gcggcgGGCGAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCAcggcg                   >  1:2252556/1‑100 (MQ=255)
                                                                                 cgggcgAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACTCTCGGAAACCACGGCGGTAt               <  1:3187200/100‑1 (MQ=255)
                                                                                  gggcgAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGGCGGTAtt              <  2:3858731/100‑1 (MQ=255)
                                                                                    gcgAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGGCGGTATTcc            <  2:2678243/100‑1 (MQ=255)
                                                                                    gcgAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGGCGGTATTcc            >  1:1375495/1‑100 (MQ=255)
                                                                                    gcgAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGCCGGTATTcc            >  2:2044149/1‑100 (MQ=255)
                                                                                     cgAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGGCGGTATTccc           >  2:547135/1‑99 (MQ=255)
                                                                                     cgAAGAATGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGGCGGTATTCCa           >  2:74965/1‑100 (MQ=255)
                                                                                        agaaTGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGGCGGTATTCCAGCt        <  2:1287630/100‑1 (MQ=255)
                                                                                         gaaTGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGGCGGTATTCCAGCtt       >  2:1274055/1‑100 (MQ=255)
                                                                                         gaaTGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGGCGGTATTCCAGCtt       <  1:2699747/100‑1 (MQ=255)
                                                                                         gaaTGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGGCGGTATTCCAGCtt       >  2:3331114/1‑100 (MQ=255)
                                                                                          aagggCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACTGCGGTATTCCAGCTTg      <  2:15934/97‑1 (MQ=255)
                                                                                           aTGGCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGGCGGTATTCCAGCTTGa     >  2:2801482/1‑100 (MQ=255)
                                                                                             ggCGCGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGGCGGTATTCCAGCTTGAAt   <  1:875955/100‑1 (MQ=255)
                                                                                              gcgcgcCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGGCGGTATTCCAGCTTGAATc  <  2:2577493/100‑1 (MQ=255)
                                                                                              gcgcgcCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGGCGGTATTCCAGCTTGAATc  >  1:3756551/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GCCGGACTGGTGAAGTTCGACTTCCTTGGTTTGCGTACGCTCACCATCATCAACTGGGCGCTGGAGATGATCAACAAGCGGCGGGCGAAGAATGGCGAGCCGCCGCTGGATATCGCTGCGATCCCGCTGGATGATAAGAAAAGCTTCGACATGCTGCAACGCTCGGAAACCACGGCGGTATTCCAGCTTGAATC  >  NC_000913/206770‑206963

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: