Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 538,830 | T→C | 12.9% | L402P (CTG→CCG) | rpoB → | DNA‑directed RNA polymerase subunit beta |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 538,830 | 0 | T | C | 12.9% | 75.5 / 3.2 | 31 | L402P (CTG→CCG) | rpoB | DNA‑directed RNA polymerase subunit beta |
Reads supporting (aligned to +/- strand): ref base T (15/12); new base C (2/2); total (17/14) | |||||||||||
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.41e-01 |
ACACCACCAGCAACCAACTGGAAGCGCTGGTCGAGATCTACCGCATGATGCGTCCAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACATCGTCGAGGTTCTGAAGACCCTGGTCGATATCCGTAACGGCAAAGGCATCGTCGACGACATC > NC_002947/538701‑538975 | acacCACCAGCAACCAACTGGAAGCGCTGGTCGAGATCTACCGCATGATGCGTCCAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGa > 2:129033/1‑150 (MQ=255) ccaccaTCAACCAACTGGAAGCGCTGGTCGAGATCTACCGCATGATGCGTCCAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGt < 1:503923/150‑1 (MQ=255) ccaccaGCAACCAACTGGAAGCGCTGGTCGAGATCTACCGCATGATGCGTCCAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGt > 2:379165/1‑150 (MQ=255) caacaGCAACCAACTGGAAGCGCTTGTCGAGATCTACCGCATGATGCGTCCAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGtt > 2:647683/4‑150 (MQ=255) caGCAACCAACTGGAAGCGCTGGTCGAGATCTACCGCATGATGCGTCCAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCaa < 2:199313/150‑1 (MQ=255) caGCAACCAACTGGAAGCGCTGGTCGAGATCTACCGCATGATGCGTCCAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCaa > 1:74278/1‑150 (MQ=255) cTGGAAGCGCTGGTCGAGATCTACCGCATGATGCGTCCAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATc > 2:643014/1‑150 (MQ=255) tGGAAGCGCTGGTCGAGATCTACCGCATGATGCGTCCAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCg > 1:198988/1‑150 (MQ=255) ggTCGAGATCTACCGCATGATGCGTCCAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCga < 1:379165/150‑1 (MQ=255) aTCTACCGCATGATGCGTCCAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACATCCTGTTCTTCAGTGCCGAGCGTTACTAGCTGTCCGCCGGAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAg < 2:90076/150‑1 (MQ=255) ccGCATGATGCGTCCAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCCGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTc > 1:760756/1‑150 (MQ=255) gCATGATGCGTCCAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCgg > 1:550511/1‑150 (MQ=255) tCCAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGTCGTGCTGAg > 2:307420/1‑150 (MQ=255) cAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTa < 1:331399/120‑1 (MQ=255) cAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTa > 2:331399/1‑120 (MQ=255) aaccaaGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACAtcgt > 1:477959/1‑150 (MQ=255) caaGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACATCGTCGa < 1:69913/150‑1 (MQ=255) cGCAGCCGAGACCCTGTGCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGAGCTTCTGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACATCGTCGAGGTTCt < 1:647683/150‑1 (MQ=255) aGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCCGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACATCGTCGAGGTTCTGaa < 2:760756/150‑1 (MQ=255) aGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCCGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACATCGTCGAGGTTCTGaa > 2:448170/1‑150 (MQ=255) ccTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACATCGTCGAGGTTCTGAAGACCCTGGt > 1:458356/1‑150 (MQ=255) ccTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACATCGTCGAGGTTCTGAAGACCCTGGt < 2:458356/150‑1 (MQ=255) cTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACATCGTCGAGGTTCTGAAGACCCTGGTc > 1:634918/1‑150 (MQ=255) tCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCCGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACATCGTCGAGGTTCTGAAGACCCTGGTCGata < 1:448170/150‑1 (MQ=255) gTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCgg > 2:131455/1‑89 (MQ=255) gTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCgg < 1:131455/89‑1 (MQ=255) gCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACATCGTCGAGGTTCTGAAGACCCTGGTCGATATCCGTAACGGCAAAGGCAtc < 1:643014/150‑1 (MQ=255) gTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACATCGTCGAGGTTCTGAAGACCCTGGTCGATATcc < 1:363647/126‑1 (MQ=255) gTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACATCGTCGAGGTTCTGAAGACACTGGTCGATATcc > 2:363647/1‑126 (MQ=255) aCGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACATCGTCGAGGTTCTGAAGACCCTGGTCGATATCCGTAACGGCAAAGGCATCGTCGACGACa < 2:634918/150‑1 (MQ=255) gATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACATCGTCGAGGTTCTGAAGACCCTGGTCGATATCCGTAACGGCAAAGGCATCGTCGACGACATc > 1:485672/1‑150 (MQ=255) | ACACCACCAGCAACCAACTGGAAGCGCTGGTCGAGATCTACCGCATGATGCGTCCAGGCGAGCCGCCAACCAAGGACGCAGCCGAGACCCTGTTCAACAACCTGTTCTTCAGTGCCGAGCGTTACGATCTGTCCGCCGTAGGCCGCATGAAGTTCAACCGTCGTATCGGTCGTACCGAGATCGAGGGTTCGGGCGTGCTGAGCAAGGAAGACATCGTCGAGGTTCTGAAGACCCTGGTCGATATCCGTAACGGCAAAGGCATCGTCGACGACATC > NC_002947/538701‑538975 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |