Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,681,655 G→A 21.0% L44L (CTG→CTA recJ → ssDNA‑specific exonuclease RecJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,681,6550GA21.0% 35.3 / 4.0 19L44L (CTG→CTArecJssDNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  ref base G (7/8);  new base A (2/2);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGT  >  NC_002947/1681535‑1681780
                                                                                                                        |                                                                                                                             
agcTCGCCCGCTGCCGCAGACCCTGCCATTTCTGGTTAACCTGACACCCTTGCTGACCCGCCTTTACGCCGGGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACcagcag                                                                                                                <  2:33982/134‑1 (MQ=255)
        cGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGCGCTGGACAAGAGCCTGGCCCGCCTACTGCCCTACCAGCAGCTCAAggg                                                                                                        >  1:214140/1‑136 (MQ=255)
        cGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTACTGCCCTACCAGCAGCTCAAggg                                                                                                        >  1:212641/1‑136 (MQ=255)
              cgccgACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCCCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGa                                                                                                  >  1:286569/1‑136 (MQ=255)
                ccgACCCTGCCATTTCTGGGTAACCTGCCCCCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAgg                                                                                                <  2:209648/136‑1 (MQ=255)
                                                           gCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcgccagcg                                                     >  2:150740/1‑136 (MQ=255)
                                                                  cgccgcGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTc                                              <  2:186595/136‑1 (MQ=255)
                                                                          gcggcgTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCgg                                      >  1:104628/1‑136 (MQ=255)
                                                                                           gCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGGCGAAGCCATCGCCCCGCGCCAGCGGATCCTCATCGTCGGCGACGTCGCTGCCgccc                     >  1:39654/1‑133 (MQ=255)
                                                                                             cGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATcc                                                >  1:201446/1‑107 (MQ=255)
                                                                                             cGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATcc                                                <  2:201446/107‑1 (MQ=255)
                                                                                             cGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGt                   >  1:172812/1‑136 (MQ=255)
                                                                                             cGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGt                   >  1:240737/1‑136 (MQ=255)
                                                                                             cGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGt                   <  2:250888/136‑1 (MQ=255)
                                                                                                     aCAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGAATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGc           <  1:109323/136‑1 (MQ=255)
                                                                                                        agagCCTGGCACGCCTGCGGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCTAAGCCATCCACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAg        <  1:123986/136‑1 (MQ=255)
                                                                                                             cTGGCACGCCTACTGCCCTCCCAGCAGCTCACGGGCATCGCGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGCCCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCg   <  2:214140/136‑1 (MQ=255)
                                                                                                             cTGGCACGCCTACTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCg   <  2:212641/136‑1 (MQ=255)
                                                                                                              gggCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGt  <  1:150740/135‑1 (MQ=255)
                                                                                                                        |                                                                                                                             
ACCTCGCCCGCTGCCGCCGACCCTGCCATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGT  >  NC_002947/1681535‑1681780

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: