Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 303,104 G→T 36.4% R423L (CGG→CTG)  envZ → osmolarity sensor protein EnvZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947303,1040GT36.4% 11.0 / 5.9 11R423L (CGG→CTG) envZosmolarity sensor protein EnvZ
Reads supporting (aligned to +/- strand):  ref base G (2/5);  new base T (2/2);  total (4/7)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

GGCACGGGGCGGCAAGGGCACCGGGCTGGGGCTGGCGATCGTCAAGCGGATTGCTGCGCAACATGGCGGGAACGTGGAGTTGCGTAACCGGTCTGGCGGCGGGATCGAGGCGCGGGTCAGGTTGCCGCTGGGGCTGCTGCTGCCGCGTAATGCTGTGTGAGGCTCACTGGCCTCTTCACGGCTAAAGCCGCACCTACAGGGATTGCGCAGTACCTGTGGC  >  NC_002947/302991‑303210
                                                                                                                 |                                                                                                          
ggCACGGGGCGGCAAGGGCACCGGGCTGGGGCTGGCGATCGTCAAGCGGATTGCTGCGCAACATGGCGGGAACGTGGAGTTGCGTAACCGGTCTGGCGGCGGGATCGAGGCGCGGGTCAGGTTGCCGCTGGGgctg                                                                                      <  1:86486/136‑1 (MQ=255)
   aCGGGGCGGCAAGGGCACCGGGCTGGGGCTGGCGATCGTCAAGCGGATTGCTGCGCAACATGGCGGGAACGTGGAGTTGCGTAACCGGTCTGGCGGCGGGATCGAGGCGCGGGTCAGGTTGCCGCTGGGgctgctg                                                                                   <  2:191769/136‑1 (MQ=255)
           gcAAGGGCACCGGGCTGGGGCTGGCGATCGTCAAGCGGATTGCTGCGCAACATGGCGGGAACGTGGAGTTGCGTAACCGGTCTGGCGGCGGGATCGAGGCGCGGGTCAGGTTGCCGCTGGGGCTGCTGCTGCcgcg                                                                           <  1:375124/136‑1 (MQ=255)
                  cACCGGGCTGGGGCTGGCGATCGTCAAGCGGATTGCTGCGCAACATGGCGGGAACGTGGAGTTGCGTAACCGGTCTGGCGGCGGGATCGAGGCGCTGGTCAGGTTGCCGCTGTGGCTGCTGCTGCCGCGTAATGCt                                                                    >  2:2716/1‑136 (MQ=255)
                  cACCGGGCTGGGGCTGGCGATCGTCAAGCGGATTGCTGCGCAACATGGCGGGAACGTGGAGTTGCGTAACCGGTCTGGCGGCGGGATCGAGGCGCTGGTCAGGTTGCCGCTGTGGCTGCTGCTGCCGCGTAATGCt                                                                    >  2:79980/1‑136 (MQ=255)
                      gggctggggctgGCGATCGTCAAGCGGATTGCTGCGCAACATGGCGGGAACGTGGAGTTGCGTAACCGGTCTGGCGGCGGGATCGAGGCGCGGGTCAGGTTGCCGCTGGGGCTGCTGCTGCCGCGTAATGCtgtgt                                                                <  2:69613/136‑1 (MQ=255)
                                     aTCGTCAAGCGGATTGCTGCGCAACATGGCGGGAACGTGGAGTTGCGTAACCGGTCTGGCGGCGGGATCGAGGCGCGGGTCAGGTTGCCGCTGGGGCTGCTGCTGCCGCGTAATGCTGTGTGAGGCTCACTGGcca                                                 >  2:393804/1‑135 (MQ=255)
                                                                    ggAACGTGGAGTTGCGTAACCGGTCTGGCGGCGGGATCGAGGCGCGGGTCAGGTTGCCGCTGGGGCTGCTGCTGCCGCGTAATGCTGTGTGAGGCTCACTGGCCTCTTCACGGCTAAAGCCGCACCTACAGGGAtt                  <  1:234458/136‑1 (MQ=255)
                                                                      aaCGTGGAGTTGCGTAACCGGTCTGGCGGCGGGATCGAGGCGCTGGTCAGGTTGCCGCTGTGGCTGCTGCTGCCGCGTAATGCTGTGTGAGGCTCACTGGCCTCTTCACGGCTAAAGCCGCACCTACAGGGATTgc                <  1:2716/136‑1 (MQ=255)
                                                                      aaCGTGGAGTTGCGTAACCGGTCTGGCGGCGGGATCGAGGCGCTGGTCAGGTTGCCGCTGTGGCTGCTGCTGCCGCGTAATGCTGTGTGAGGCTCACTGGCCTCTTCACGGCTAAAGCCGCACCTACAGGGATTgc                <  1:79980/136‑1 (MQ=255)
                                                                                    tAACCGGTCTGGCGGCGGGATCGAGGCGCGGGTCAGGTTGCCGCTGGGGCTGCTGCTGCCGCGTAATGCTGTGTGAGGCTCACTGGCCTCTTCACGGCTAAAGCCGCACCTACAGGGATTGCGCAGTACCTGTGGc  >  1:172473/1‑136 (MQ=255)
                                                                                                                 |                                                                                                          
GGCACGGGGCGGCAAGGGCACCGGGCTGGGGCTGGCGATCGTCAAGCGGATTGCTGCGCAACATGGCGGGAACGTGGAGTTGCGTAACCGGTCTGGCGGCGGGATCGAGGCGCGGGTCAGGTTGCCGCTGGGGCTGCTGCTGCCGCGTAATGCTGTGTGAGGCTCACTGGCCTCTTCACGGCTAAAGCCGCACCTACAGGGATTGCGCAGTACCTGTGGC  >  NC_002947/302991‑303210

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