Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,824,771 C→A 5.6% intergenic (+139/+35) PP_5102 → / ← trmB hypothetical protein/tRNA (guanosine(46)‑N7)‑methyltransferase TrmB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,824,7710CA 5.6% 344.4 / 2.4 113intergenic (+139/+35)PP_5102/trmBhypothetical protein/tRNA (guanosine(46)‑N7)‑methyltransferase TrmB
Reads supporting (aligned to +/- strand):  ref base C (51/54);  new base A (4/3);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

GCCGCGGCTGAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGTGCCGTCTGCC  >  NC_002947/5824622‑5824920
                                                                                                                                                      |                                                                                                                                                     
gccgcGGCTGAGCGAATGACCCCTTGTAGGAGCGGATGCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTc                                                                                                                                                       >  1:1352566/1‑150 (MQ=255)
 ccgcGGCTGAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCa                                                                                                                                                      >  1:2292983/1‑150 (MQ=255)
       cTGAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCgg                                                                                                                                                <  1:2282270/150‑1 (MQ=255)
       cTGAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCgg                                                                                                                                                >  1:1200323/1‑150 (MQ=255)
         gaccGAATGACCCCTCGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAATGAAACAACCGCCCTGAAAAGGGCGGtt                                                                                                                                              <  1:266117/147‑1 (MQ=18)
         gAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGtt                                                                                                                                              >  1:218411/1‑150 (MQ=255)
         gAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGtt                                                                                                                                              <  1:714728/150‑1 (MQ=255)
         gAGCGAAGGCCCCCTCGTGGGAGCGGAGTCAACGGCGAAAGGGCCGGCCCGCCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAACCCCGCTCTTACAAAGAAAAAACCGCCCTTTTCAGGGCGGtt                                                                                                                                              <  1:1147241/150‑1 (MQ=21)
           gCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGtttt                                                                                                                                            >  1:2206482/1‑150 (MQ=255)
           gCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGtttt                                                                                                                                            >  1:2046494/1‑150 (MQ=255)
            cGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGTGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTa                                                                                                                                           <  1:1920599/150‑1 (MQ=255)
            cGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTa                                                                                                                                           <  1:188437/150‑1 (MQ=255)
             gatccaCCCCTTGTAGGAGCGGAGCAACCCGCACAAGGTGCGGCCCTGCTCGCCAATCTCCACCCCCCCCTCACCGTGGTCGCCGCTACCATATTCGCGGTTAAACCCGCTCCTACAAAGAAAAAACCGCCCTGAAAAGGGCGGtttttt                                                                                                                                          <  1:1873767/145‑3 (MQ=255)
              aaTGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTAcc                                                                                                                                         >  1:1189919/1‑150 (MQ=255)
                  aCCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGaaa                                                                                                                                     >  2:575181/1‑150 (MQ=255)
                     ccTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGaaaa                                                                                                                                    <  2:1648956/148‑1 (MQ=255)
                     ccTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGaaaa                                                                                                                                    >  1:1648956/1‑148 (MQ=255)
                     ccTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAgg                                                                                                                                  <  1:981486/150‑1 (MQ=255)
                     ccTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAgg                                                                                                                                  <  2:2206482/150‑1 (MQ=255)
                       ttGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGCCCAACCCGCTCATGGCTCTTGTAAACCCGCACCAACAAAAAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGcc                                                                                                                                <  2:1834308/150‑1 (MQ=16)
                          tAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTa                                                                                                                             <  2:1783398/150‑1 (MQ=255)
                            ggAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATc                                                                                                                           >  1:353289/1‑150 (MQ=255)
                               gCGGATTCATCCGCCAAAGGGCCGTCGCTGCGGGCCTATTTCCACCCGCTCTTCACCGTGGGCGCCGCTTCGGTATTCGCGGGTAAACCCGCTTCTACAAAGAAAAAACCGCTCTTTTCATGGCGGTTTTACCGAAAAGGCCGTATCAGc                                                                                                                        <  2:181817/150‑1 (MQ=21)
                               gCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGc                                                                                                                        <  2:1973425/150‑1 (MQ=255)
                               gCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGc                                                                                                                        >  1:678021/1‑150 (MQ=255)
                                   aTTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCt                                                                                                                    >  2:1598610/1‑150 (MQ=255)
                                        cccGCAAAAGGGCCGGCCATCCTGCCCAATCTCCGCCCGCGCTTCACCGTGGGAGCCGCTACCGTATTCGCGGGAAACACCGCGCCAACAAAGAAAAAACCGCCCTTTTCAGCGCGGTTTTACCGCTAAGGCCGTATCAGCCGCTCAATc                                                                                                               <  2:371149/149‑1 (MQ=25)
                                        aCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTGAAAAGGGCGGttt                                                                                                                                             <  1:1954824/120‑1 (MQ=18)
                                        aCCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGGATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCGTTTAGGGGCGGGTTTACCGAAAAGGCCGTATCAGCCGCtgcatc                                                                                                               >  1:251082/1‑145 (MQ=25)
                                        aCCGCGAAAGCGCCGCCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCCCTCCTACAAACAAAAACCCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATc                                                                                                               <  1:1802924/150‑1 (MQ=25)
                                         ccGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATcc                                                                                                              <  2:353289/150‑1 (MQ=255)
                                           gcgAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTAt                                                                                                                            >  1:1225694/1‑134 (MQ=255)
                                           gcgAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTAt                                                                                                                            <  2:1225694/134‑1 (MQ=255)
                                              aaaGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCtt                                                                                                         >  1:1235055/1‑150 (MQ=255)
                                                 gggCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTtctc                                                                                                      <  2:1719060/150‑1 (MQ=255)
                                                  ggccggccCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCg                                                                                                     <  1:575181/150‑1 (MQ=255)
                                                    ccggccCTGCTGGCCCATCTCCACCCGCCCTTATGCGGGGACTTCTTTATCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGaa                                                                                                   <  1:1319740/150‑1 (MQ=37)
                                                       gccCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACtt                                                                                                >  1:933047/1‑150 (MQ=255)
                                                            gctgGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGt                                                                                           <  1:1510268/150‑1 (MQ=255)
                                                            gctgGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGt                                                                                           <  2:2175739/150‑1 (MQ=255)
                                                              tgGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCGCGAACTTCAGGTcc                                                                                         >  2:605977/1‑150 (MQ=255)
                                                               ggCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCGCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCCATCCACCTTCTCGAACTTCAAGTccc                                                                                        >  1:1543178/1‑150 (MQ=255)
                                                                  cAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCcaca                                                                                     >  1:374398/1‑150 (MQ=255)
                                                                  cAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCcaca                                                                                     >  1:676328/1‑150 (MQ=255)
                                                                  cAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCcaca                                                                                     >  2:1175332/1‑150 (MQ=255)
                                                                  cAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTGTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCcaca                                                                                     >  1:187258/1‑150 (MQ=255)
                                                                     tctcCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACAccc                                                                                  >  1:2012630/1‑150 (MQ=255)
                                                                        ccACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGt                                                                               >  2:2177132/1‑150 (MQ=255)
                                                                         cACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTg                                                                              <  2:933047/150‑1 (MQ=255)
                                                                          aCCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGc                                                                             <  2:2292983/150‑1 (MQ=255)
                                                                            ccGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTCTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGccc                                                                           <  1:1759151/150‑1 (MQ=255)
                                                                             cGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGcccc                                                                          >  1:801321/1‑150 (MQ=255)
                                                                              gCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCaa                                                                                                                 <  1:1169988/110‑1 (MQ=255)
                                                                              gCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCCCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCaa                                                                                                                 >  2:1169988/1‑110 (MQ=255)
                                                                              cctccccACCGCGGGCGCCCCCACCGTATTGGCGGGTACACCCGCTCCTACAAAGAAAAAACCGCCCTGTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCGCCATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCa                                                                         <  2:1450228/144‑1 (MQ=255)
                                                                                ccTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGc                                                                       <  1:1545298/150‑1 (MQ=255)
                                                                                       cGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGc                                                                                                                                 <  2:1964983/85‑1 (MQ=255)
                                                                                       cGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGc                                                                                                                                 >  1:1964983/1‑85 (MQ=255)
                                                                                          gggCGCCGCTACCGTATTCGCGGGTAAACCCGCTCGTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTCCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACcgc                                                             <  2:676328/150‑1 (MQ=255)
                                                                                                cgcTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGtt                                                       <  1:2258701/150‑1 (MQ=255)
                                                                                                   tACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGa                                                    >  1:1364695/1‑150 (MQ=255)
                                                                                                    aCCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCgcggc                                                          >  2:587904/1‑143 (MQ=255)
                                                                                                    aCCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCgcggc                                                          <  1:587904/143‑1 (MQ=255)
                                                                                                     ccGTATTCCCGTGTAAACGCGTTCCTACAACGAAAAACCCGCCCTTTTCAGGGCGGTTTTACCTAAATGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCATGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAAc                                                  <  2:343886/150‑1 (MQ=255)
                                                                                                         aTTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCAc                                                                                                            <  1:125037/88‑1 (MQ=255)
                                                                                                         aTTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCAc                                                                                                            >  2:125037/1‑88 (MQ=255)
                                                                                                         aTTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTgg                                              <  2:801321/150‑1 (MQ=255)
                                                                                                              gcgGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAGCTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGAt                                          >  1:616010/1‑150 (MQ=255)
                                                                                                              gcgGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGAt                                          >  2:635746/1‑150 (MQ=255)
                                                                                                              gcgGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGAt                                          >  1:1523273/1‑150 (MQ=255)
                                                                                                                 ggTAACCCCGCTCTTACAAAGAAAACATCTCCCTTTTCAGGGCGGCTTTACCGAAAAGGCCGTATCAGCCGCGCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCTAACTTGGTGATCgg                                       <  2:280535/150‑1 (MQ=255)
                                                                                                                 ggTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCTTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCgg                                       <  1:917295/150‑1 (MQ=255)
                                                                                                                       cccGCTCCTACAAAGAAAAACACGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTg                                               <  1:712191/136‑1 (MQ=255)
                                                                                                                       cccGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTg                                               >  2:712191/1‑136 (MQ=255)
                                                                                                                       cccGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGctc                                 >  1:1639396/1‑150 (MQ=255)
                                                                                                                        ccGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGcc                                                                            >  1:314783/1‑106 (MQ=255)
                                                                                                                        ccGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGcc                                                                            <  2:314783/106‑1 (MQ=255)
                                                                                                                        ccGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGctct                                >  1:1306439/1‑150 (MQ=255)
                                                                                                                        ccGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGctct                                >  2:891070/1‑150 (MQ=255)
                                                                                                                         cGCTCCTACAAAGAAACAACCGCCCTTTTTAGGTCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCAAGCAATTCCCCGCGGCGTTCGAACTTGGTGATCGGGCGCTCtt                               <  2:573429/150‑1 (MQ=255)
                                                                                                                               tACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGgc                         >  1:2022534/1‑150 (MQ=255)
                                                                                                                                  aaaGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGcgcg                      >  1:2130751/1‑150 (MQ=255)
                                                                                                                                  aaaGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGcgcg                      >  2:1036766/1‑150 (MQ=255)
                                                                                                                                   aaGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCgg                     <  2:1523273/150‑1 (MQ=255)
                                                                                                                                      gtaaaaCCGCCCTGAAAAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGcac                  >  1:261447/3‑150 (MQ=255)
                                                                                                                                        aaaaCCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACTGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACAt                >  2:1219361/1‑150 (MQ=255)
                                                                                                                                         aaaCCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATa               >  2:300087/1‑150 (MQ=255)
                                                                                                                                          aaCCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAg              <  1:891070/150‑1 (MQ=255)
                                                                                                                                            ccGCCCTTTTCAGGGGCGTCTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACGTCAAGTCCCACACCCCGTGCCCCAGCCGCTCGCCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGt            <  2:1543178/150‑1 (MQ=255)
                                                                                                                                             cGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGTCCCAGCCGCTCAcc                                                               <  2:885908/98‑1 (MQ=255)
                                                                                                                                             cGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGTCCCAGCCGCTCAcc                                                               >  1:885908/1‑98 (MQ=255)
                                                                                                                                             cGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCgg                                       >  1:933942/1‑122 (MQ=255)
                                                                                                                                             cGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCgg                                       <  2:933942/122‑1 (MQ=255)
                                                                                                                                               cccTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGACCTTGGTGATCGGGCGCTCGTCGGGGCGCGGCACATAGGTGcc         >  1:1682756/1‑150 (MQ=255)
                                                                                                                                                 cTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGc                                                                       <  1:1785971/86‑1 (MQ=255)
                                                                                                                                                 cTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGc                                                                       >  2:1785971/1‑86 (MQ=255)
                                                                                                                                                  ttttCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACATTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAAc                                                  <  1:167853/106‑1 (MQ=255)
                                                                                                                                                    ttCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGTGCCGTCTg    >  2:1702388/1‑150 (MQ=255)
                                                                                                                                                      cAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGTGCCGTCTGcc  <  2:2076479/150‑1 (MQ=255)
                                                                                                                                                      cAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGTGCCGTCTGcc  <  2:1189919/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCCGCGGCTGAGCGAATGACCCCTTGTAGGAGCGGATTCAACCGCGAAAGGGCCGGCCCTGCTGGCCAATCTCCACCCGCCCTTCACCGTGGGCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACAAAGAAAAAACCGCCCTTTTCAGGGCGGTTTTACCGAAAAGGCCGTATCAGCCGCTCAATCCACCTTCTCGAACTTCAAGTCCCACACCCCGTGCCCCAGCCGCTCACCGCGGCGTTCGAACTTGGTGATCGGGCGCTCTTCGGGGCGCGGCACATAGGTGCCGTCTGCC  >  NC_002947/5824622‑5824920

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