Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,852,459 A→T 17.4% intergenic (‑130/+167) recR ← / ← ybaB recombination protein RecR/nucleoid‑associated broad specificity regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,852,4590AT17.4% 47.1 / 4.8 23intergenic (‑130/+167)recR/ybaBrecombination protein RecR/nucleoid‑associated broad specificity regulator
Reads supporting (aligned to +/- strand):  ref base A (12/7);  new base T (2/2);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

GCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCG  >  NC_002947/4852324‑4852604
                                                                                                                                       |                                                                                                                                                 
gCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAagc                                                                                                                                                >  8:72156/1‑139 (MQ=255)
gCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAagc                                                                                                                                                <  7:72156/139‑1 (MQ=255)
 cTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCTTCGCAGccctttcgc                                                                                                                                     <  8:213041/149‑7 (MQ=39)
                 caCAATTCAATAAGCGCCGGTTAAACCCCCGGGGGGGCGGTGGCGTCAAATCCCCTCCCCAAAAGCCTGCAAGGTCTGGTTTTTTGTGGCCGCGGCCTTGTGTCGCTAATGGGCTGCTAAGCAGCCCCGGCACTCTGTGCTGCGAAGCt                                                                                                                     <  8:317965/148‑1 (MQ=255)
                         aaTAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCTTCGCAGccctttc                                                                                                                                       <  8:103630/123‑5 (MQ=38)
                         aaTAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCTTCGCAGccctttc                                                                                                                                       >  7:103630/1‑119 (MQ=38)
                             agagGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCgggg                                                                                                         >  3:362675/1‑149 (MQ=255)
                                  cgttttATACCTCGGGGGGGGGGTGGGTTAAAATCACCTCCCCAAACGCCTTCATGGTTTGTTTTTTTGTGGGAGCGGCTTTGTGGCGCTAATGGGCTGCTATGCTGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCt                                                                                                    <  4:362675/146‑1 (MQ=255)
                                               gggggggcgggggcCTAAAACCCCCCCCCCCAACGCCGGCCCGGTCTGCTTTTTTGGGGCCGCGCCCTGGGGGCGCCAACGGGCGGCGACGCCCCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTc                                                                                       <  3:68671/147‑1 (MQ=255)
                                                                                         ggTCTGCTTTTTTGTGGAAGCGGCCTTGTGTCGCGAAAGGGCTGCGAagca                                                                                                                                               >  5:366316/1‑51 (MQ=255)
                                                                                         ggTCTGCTTTTTTGTGGAAGCGGCCTTGTGTCGCGAAAGGGCTGCGAagca                                                                                                                                               <  6:366316/51‑1 (MQ=255)
                                                                                                                ccTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAAc                      >  5:197440/1‑149 (MQ=255)
                                                                                                                         gcgAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCCCtgt             >  4:276538/1‑149 (MQ=255)
                                                                                                                         gcgAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgt             >  1:177636/1‑149 (MQ=255)
                                                                                                                           gccggggCTGCTTCGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttg           >  8:97855/5‑149 (MQ=255)
                                                                                                                           gAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttg           >  5:366425/1‑149 (MQ=255)
                                                                                                                             aaGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttgt          >  4:340973/1‑148 (MQ=255)
                                                                                                                                 gCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGTCACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcg     >  1:128621/1‑149 (MQ=255)
                                                                                                                                   tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGa                                                       <  7:273612/97‑1 (MQ=255)
                                                                                                                                   tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGa                                                       >  8:273612/1‑97 (MQ=255)
                                                                                                                                   tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACt               <  6:294172/137‑1 (MQ=255)
                                                                                                                                   tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACt               >  5:294172/1‑137 (MQ=255)
                                                                                                                                    gCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcggcg  >  6:7384/1‑149 (MQ=255)
                                                                                                                                       |                                                                                                                                                 
GCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCG  >  NC_002947/4852324‑4852604

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