Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,378,320 G→A 22.2% G405G (GGC→GGT dnaK ← chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,378,3200GA22.2% 31.9 / 4.9 18G405G (GGC→GGTdnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base G (7/7);  new base A (2/2);  total (9/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 = 9.31e-01

GCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGAC  >  NC_002947/5378196‑5378442
                                                                                                                            |                                                                                                                          
gCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCAACACGCCGCCCATGGtttg                                                                                                                 >  2:268485/1‑135 (MQ=255)
    aCGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATa                                                                                                             <  1:325427/136‑1 (MQ=255)
     cGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATAc                                                                                                            >  2:353942/1‑136 (MQ=255)
                      gTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTg                                                                                           <  1:353942/136‑1 (MQ=255)
                         aaTGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACg                                                                                        >  1:11370/1‑136 (MQ=255)
                                       actGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAGGCGGGTGAACACgccaccga                                                                                                                        <  2:299464/89‑3 (MQ=255)
                                       actGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAGGCGGGTGAACACgccaccga                                                                                                                        >  1:299464/2‑88 (MQ=255)
                                          ggTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTc                                                                       <  2:254107/136‑1 (MQ=255)
                                                 gtcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCACGGTCAGCGGGTTGACGTCCAGCAGCACGACGTCTTTCACg                                                                >  2:203521/1‑136 (MQ=255)
                                                      cggTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAGGCGGGTGAACACGCCACCGAGCGTTTCAATGCCGAGCGAAAGCGGGGTCACGTCGAGCAGCAGAACGTCCTTGACGTc                                                              >  1:112988/3‑133 (MQ=255)
                                                      cggTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAGGCGGGTGAACACGCCACCGAGCGTTTCAATGCCGAGCGAAAGCGGGGTCACGTCGAGCAGCAGAACGTCCTTGACGTc                                                              <  2:112988/131‑1 (MQ=255)
                                                      cggTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAGGCGGGTGAACACGCCACCGAGCGTTTCAATGCCGAGCGAAAGCGGGGTCACGTCGAGCAGCAGAACGTCCTTGACGTc                                                              <  2:110558/131‑1 (MQ=255)
                                                      cggTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAGGCGGGTGAACACGCCACCGAGCGTTTCAATGCCGAGCGAAAGCGGGGTCACGTCGAGCAGCAGAACGTCCTTGACGTc                                                              >  1:110558/3‑133 (MQ=255)
                                                      cAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCAcc                                                           <  2:271845/136‑1 (MQ=255)
                                                        gTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCgg                                                         <  2:127606/136‑1 (MQ=255)
                                                        gTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCgg                                                         <  2:335675/136‑1 (MQ=255)
                                                                      cGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCAccc                                           >  1:400764/1‑136 (MQ=255)
                                                                                 tCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGTTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgcagc                                <  1:203521/136‑1 (MQ=255)
                                                                                                  ttctCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCtt               >  1:319177/1‑136 (MQ=255)
                                                                                                               cgGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGAc  >  2:510256/1‑136 (MQ=255)
                                                                                                                            |                                                                                                                          
GCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGAC  >  NC_002947/5378196‑5378442

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