Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 11,649 A→G 13.8% S183G (AGC→GGC)  dnaN → DNA polymerase III subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294711,6490AG13.8% 64.7 / 3.5 29S183G (AGC→GGC) dnaNDNA polymerase III subunit beta
Reads supporting (aligned to +/- strand):  ref base A (13/12);  new base G (2/2);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

TGCGCCGCCTGATCGAGCGAACCAGCTTCGCCATGGCCCAGCAGGATGTGCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGTGCCGCGTAAAGGTATCCTGGAGTTGGCGCGCCTGCTGACCGATCCTGAAGGTATGGTCAGCATCGTGCTGGGCCAGCATCACATCCGCGCCACTACTGGTGAGTT  >  NC_002947/11506‑11794
                                                                                                                                               |                                                                                                                                                 
tGCGCCGCCTGATCGAGCGAACCAGCTTCGCCATGGCCCAGCAGGATGTGCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCAATGAgcgcgc                                                                                                                                             <  1:201332/150‑1 (MQ=255)
    ccgccTGATCGAGCGAACCAGCTTCGCCATGGCCCAGCAGGATGTGCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGAt                                                                                                                                         >  2:258116/1‑150 (MQ=255)
      gccTGATCGAGCGAACCAGCTTCGCCATGGCCCAGCAGGATGTGCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATCGgcgcgctcatcg                                                                                                                                       >  2:444853/1‑144 (MQ=255)
             cGAGCGAACCAGCTTCGCCATGGCCCAGCAGGATGTGCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATCGGCGCGCTCATCGAGCaaag                                                                                                                                <  1:444853/150‑4 (MQ=255)
                 cGAACCAGCTTCGCCATGGCCCAGCAGGATGTGCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGAt                                                                                                                            <  1:281827/150‑1 (MQ=255)
                  gAACCAGCTTCGCCATGGCCCTGCAGGATGTTCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGATc                                                                                                                           <  1:311821/150‑1 (MQ=255)
                  gAACCAGCTTCGCCATGGCCCAGCAGGATGTGCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGATc                                                                                                                           <  2:332448/150‑1 (MQ=255)
                      cAGCTTCGCCATGGCCCAGCAGGATGTGCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGATCGaca                                                                                                                       >  1:534745/1‑150 (MQ=255)
                      cAGCTTCGCCATGGCCCAGCAGGATGTGCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGATCGaca                                                                                                                       >  2:96876/1‑150 (MQ=255)
                                     ccAGCAGGATGTGCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGGGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGTGcc                                                                                                        >  2:40663/1‑150 (MQ=255)
                                     ccAGCAGGATGTGCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGTGcc                                                                                                        >  2:529306/1‑150 (MQ=255)
                                        gcagGATGTGCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAgc                                                                                                                                                 <  1:336665/106‑1 (MQ=255)
                                        gcagGATGTGCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAgc                                                                                                                                                 >  2:336665/1‑106 (MQ=255)
                                             ttgttCGCGATTCACTTAACGGCATGCTGCTTGCGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCCCCGCCTGGCGCTTTGCTCGATGAGCGCGCAGATTGAGCAGGAGGATCGACACCAGGTCATCGTGCCGCGTAAAg                                                                                                <  1:538759/149‑1 (MQ=255)
                                              tgtgCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGTGCCGCGTAAAgg                                                                                               >  1:258527/1‑150 (MQ=255)
                                                           cTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGTGCCGCGTAAAGGTATCCTGGAGTTg                                                                                  >  1:528324/1‑150 (MQ=255)
                                                             cAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGTGCCACGTAAAGGTATCCTGGAGTTGgc                                                                                >  2:184519/1‑150 (MQ=255)
                                                                        tgctTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGcc                                                                                                                                            >  2:408326/1‑79 (MQ=255)
                                                                        tgctTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGcc                                                                                                                                            <  1:408326/79‑1 (MQ=255)
                                                                                   tCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGTGCCACGTAAAGGTATCCTGGAGTTGGCGCGCCTGCTGACCGATCCTGaa                                                          <  1:184519/150‑1 (MQ=255)
                                                                                                           tCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGTGCCGCGTAAAGGTATCCTGGAGTTGGCGCGCCTGCTGACCGATCCTGAAGGTATGGTCAGCGTCGTGCTGGGc                                  >  1:476620/1‑150 (MQ=255)
                                                                                                                 gACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGTGCCGCGTAAAGGTATCCTGGAGTTGGCGCGCCTGCTGACCGATCCTGAAGGTATGGTCAGCATCGTGCTGGGCCAGCAt                            >  2:379883/1‑150 (MQ=255)
                                                                                                                     gTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGt                                                                                                           <  1:90211/67‑1 (MQ=255)
                                                                                                                     gTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGt                                                                                                           >  2:90211/1‑67 (MQ=255)
                                                                                                                                        cTCGATGGGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGTGCCGCGTAAAGGTATCCTGGAGTTGGCGCGCCTGCTGACCGATCCTGAAGGTATGGTCAGCATCGTGCTGGGCCAGCATCACATCCGCGCCACTACTGGTGa     <  1:40663/150‑1 (MQ=255)
                                                                                                                                        cTCGATGAGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGTGCCGCGTAAAGGTATCCTGGAGTTGGCGCGCCTGCTGACCGATCCTGAAGGTATGGTCAGCATCGTGCTGGGCCAGCATCACATCCGCGCCACTACTGGTGa     <  1:191356/150‑1 (MQ=255)
                                                                                                                                        cTCGATGAGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGTGCCGCGTAAAGGTATCCTGGAGTTGGCGCGCCTGCTGACCGATCCTGAAGGTATGGTCAGCATCGTGCTGGGCCAGCATCACATCCGCGCCACTACTGGTGa     >  1:220844/1‑150 (MQ=255)
                                                                                                                                         tCGATGAGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGTGCCGCGTAAAGGTATCCTGGAGTTGGCGCGCCTGCTGACCGATCCTGAAGGTATGGTCAGCATCGTGCTGGGCCAGCATCACATCCGCGCCACTACTGGTGAg    <  1:258116/150‑1 (MQ=255)
                                                                                                                                           gATGAGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGTGCCGCGTAAAGGTATCCTGGAGTTGGCGCGCCTGCTGACCGATCCTGAAGGTATGGTCAGCATCGTGCTGGGCCAGCATCACATCCGCGCCACTACTGGCGAGtt  <  2:534745/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
TGCGCCGCCTGATCGAGCGAACCAGCTTCGCCATGGCCCAGCAGGATGTGCGCTATTACCTCAACGGCATGCTGCTTGAGGTTTCCCGCAACACGTTGCGCGCTGTATCCACCGACGGTCACCGCCTGGCGCTTTGCTCGATGAGCGCGCCGATCGAGCAGGAGGATCGACACCAGGTCATCGTGCCGCGTAAAGGTATCCTGGAGTTGGCGCGCCTGCTGACCGATCCTGAAGGTATGGTCAGCATCGTGCTGGGCCAGCATCACATCCGCGCCACTACTGGTGAGTT  >  NC_002947/11506‑11794

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