Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 539,720 C→T 26.7% R699C (CGT→TGT)  rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947539,7200CT26.7% 16.7 / 4.9 15R699C (CGT→TGT) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (5/6);  new base T (2/2);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

ACGTCACCCTGATGGACGTTTCGCCGAAGCAGGTTGTTTCCGTTGCAGCGTCGCTGATTCCGTTCCTCGAGCACGACGACGCCAACCGTGCGTTGATGGGTTCGAACATGCAGCGTCAGGCTGTACCGACCCTGCGTGCCGACAAGCCGCTGGTAGGTACCGGCATGGAGCGCAACGTTGCCCGTGACTCCGGTGTCTGCGTGGTTGCTCGCCGCGGTGGTGTGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCC  >  NC_002947/539586‑539851
                                                                                                                                      |                                                                                                                                    
aCGTCACCCTGATGGACGTTTCGCCGAAGCAGGTTGTTTCCGTTGCAGCGTCGCTGATTCCGTTCCTCGAGCACGACGACGCCAACCGTGCGTTGATGGGTTCGAACATGCAGCGTCAGGCTGTACCGACCCTGCGTGCCGACAAGCCGc                                                                                                                       <  2:96450/150‑1 (MQ=255)
                cGTTTCGCCGAAGCAGGTTGTTTCCGTTGCAGTGTCGCGGATCCCGTTCCTCGAGCACGACGACGCCAACCGTGCGTTGATGGGTTCGAACATGCAGCGTCAGGCTGTACCGACCCTGCGTGCCGACAAGCCGCTGGTAGGTACCGGCAt                                                                                                      <  2:40885/150‑1 (MQ=255)
                cGTTTCGCCGAAGCAGGTTGTTTCCGTTGCAGCGTCGCTGATTCCGTTCCTCGAGCACGACGACGCCAACCGTGCGTTGATGGGTTCGAACATGCAGCGTCAGGCTGTACCGACCCTGCGTGCCGACAAGCCGCTGGTAGGTACCGGCAt                                                                                                      <  1:30195/150‑1 (MQ=255)
                                                            cGTTCCTCGAGCACGACGACGCCAACCGTGCGTTGATGGGTTCGAACATGCAGCGTCAGGCTGTACCGACCCTGTGTGCCGACAAGCCGCTGGTGGATGGAAGCAGGTAGAGCAACGTTGCCCGTGACTCCGGTGTCTGCGTGGTTGCTc                                                          >  2:395401/1‑150 (MQ=255)
                                                            cGTTCCTCGAGCACGACGACGCCAACCGTGCGTTGATGGGTTCGAACATGCAGCGTCAGGCTGTACCGACCCTGTGTGCCGACAAGCCGCCGGACAGCAAAAGGGGTGAGCGCAACGTTGCCCGTGACTCCGGTGTCTGCGTGGTTGCTc                                                          >  2:395293/1‑150 (MQ=255)
                                                                       cACGACGACGCCAACCGTGCGTTGATGGGTTCGAACATGCAGCGTCAGGCTGTACCGACCCTGCGTGCCGACAAGCCGCTGGTAGGTACCGGCATGGAGCGCAACGTTGCCCGTGACTCCGGTGTCTGCGTGGTTGCTCGCCGCggtggt                                               <  1:299435/150‑1 (MQ=255)
                                                                            cgacgCCAACCGTGCGTTGATGGGTTCGAACATGCAGCGTCAGGCTGTACCGACCCTGTGTGCCGACAAGCCGCTGGTAGGTACCGGCATGGAGCGCAACGTTGCCCGTGACTCCGGTGCAAACTTGGTTGCTCGCCGCGGTGGTGTGAt                                          <  1:395293/150‑1 (MQ=255)
                                                                            cgacgCCAACCGTGCGTTGATGGGTTCGAACATGCAGCGTCAGGCTGTACCGACCCTGTGTGCCGACAAGCCGCTGGTAGGTACCGGCATGGAGCGCAACGTTGCCCGTGACTCCGGTCGGTACGTTGTTGCTCGCCGCGGTGGTGTGAt                                          <  1:395401/150‑1 (MQ=255)
                                                                                   aaCCGTGCGTTGATGGGTTCGAACATGCAGCGTCAGGCTGTACCGACCCTGCGTGCCGACAAGCCGCTGGTAGGTACCGGCATGGAGCGCAACGTTGCCCGTGACTCCGGTGTCTGCGTGGTTGCTCGCCGCGGTGGTGTGATCGACTCg                                   >  2:409553/1‑150 (MQ=255)
                                                                                                   gTTCGAACATGCAGCGTCAGGCTGTACCGACCCTGCGTGCCGACAAGCCGCTGGTAGGTACCGGCATGGAGCGCAACGTTGCCCGTGACTCCGGTGTCTGCGTGGTTGCTCGCCGCGGTGGTGTGATCGACTCGGTCGACGCCAGCCGTa                   >  1:304566/1‑150 (MQ=255)
                                                                                                   gTTCGAACATGCAGCGTCAGGCTGTACCGACCCGGCGTGCCGACAAGCCGCTGGTAGGTACCGGCATGGAGCGCAACGTTGCCCGTGACTCCGGTGTCTGCGTGGTTGCTCGCCGCGGTGGTGTGATCGACTCGGTCGACGCCAGCCGTa                   >  2:104124/1‑150 (MQ=255)
                                                                                                       gAACATGCAGCGTCAGGCTGTACCGACCCTGCGTGCCGACAAGCCGCTGGTAGGTACCGGCATGGAGCGCAACGTTGCCCGTGACTCCGGTGTCTGCGTGGTTGCTCGCCGCGGTGGTGTGATCGACTCGGTCGACGCCAGCCGTATCgt               >  1:405719/1‑150 (MQ=255)
                                                                                                                    cAGGCTGTACCGACCCTGCGTGCCGACAAGCCGCTGGTAGGTACCGGCATGGAGCGCAACGTTGCCCGTGACTCCGGTGTCTGCGTGGTTGCTCGCCGCGGTGGTGTGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGcc  <  2:304566/150‑1 (MQ=255)
                                                                                                                      ggCTGTACCGACCCTGCGTGCCGACAAGCCGCTGGTAGGTACCGGCATGGAGCGCAACGTTGCCCGTGACTCCGGTGTCTGCGTGGTTGCTCGCCGCGGTGGTGTGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGc   <  1:334779/147‑1 (MQ=255)
                                                                                                                      ggCTGTACCGACCCTGCGTGCCGACAAGCCGCTGGTAGGTACCGGCATGGAGCGCAACGTTGCCCGTGACTCCGGTGTCTGCGTGGTTGCTCGCCGCGGTGGTGTGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGc   >  2:334779/1‑147 (MQ=255)
                                                                                                                                      |                                                                                                                                    
ACGTCACCCTGATGGACGTTTCGCCGAAGCAGGTTGTTTCCGTTGCAGCGTCGCTGATTCCGTTCCTCGAGCACGACGACGCCAACCGTGCGTTGATGGGTTCGAACATGCAGCGTCAGGCTGTACCGACCCTGCGTGCCGACAAGCCGCTGGTAGGTACCGGCATGGAGCGCAACGTTGCCCGTGACTCCGGTGTCTGCGTGGTTGCTCGCCGCGGTGGTGTGATCGACTCGGTCGACGCCAGCCGTATCGTTGTTCGCGTTGCC  >  NC_002947/539586‑539851

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