| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 540,089 | T→C | 14.3% | C822R (TGC→CGC) ‡ | 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 | 540,089 | 0 | T | C | 14.3% | 62.3 / 4.1 | 28 | C822R (TGC→CGC) ‡ | rpoB | DNA‑directed RNA polymerase subunit beta |
| Reads supporting (aligned to +/- strand): ref base T (11/13); new base C (2/2); total (13/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01 | |||||||||||
GTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGAAATCACTGCGGACATCCCGAACGTGGGTGAAGCCGCACTGAACAAG > NC_002947/539943‑540232 | gTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACCTGCGCATCGCGTTCATGGCGTGGAACGGGTTCAACTTCGAAGACTCCATCTGcc < 2:4737/150‑1 (MQ=255) gCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTc < 1:335516/150‑1 (MQ=255) ggTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGAAGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGa < 1:26939/150‑1 (MQ=255) aaGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGtgtt > 2:337235/1‑150 (MQ=255) gCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGAt > 1:311853/1‑145 (MQ=255) gCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGAt < 2:311853/145‑1 (MQ=255) aCGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATcc < 2:356650/141‑1 (MQ=255) aCGGTCCGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATcc > 1:356650/1‑141 (MQ=255) aCGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAgg < 1:28468/150‑1 (MQ=255) tccGTCCACCGACCTCGGTGAACTGGCGCTTGGTCAGAACATGCGCGTAGCGTTCATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACt > 1:94495/2‑150 (MQ=255) cccGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACt > 1:59276/1‑150 (MQ=255) gACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCtgtgtg < 1:264069/150‑1 (MQ=255) gggTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGcgcg < 1:337235/150‑1 (MQ=255) gTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGcc > 2:284/1‑150 (MQ=255) tCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGccc > 1:280410/1‑150 (MQ=255) cAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCa < 2:280410/150‑1 (MQ=255) gTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACTCCAAGCTCGGccc < 2:262834/132‑1 (MQ=255) gTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGccc > 1:262834/1‑132 (MQ=255) tGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAgg > 2:95183/1‑132 (MQ=255) tGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAgg < 1:95183/132‑1 (MQ=255) gTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGaa > 2:255341/1‑130 (MQ=255) gTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGaa < 1:255341/130‑1 (MQ=255) acAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAg < 2:94495/149‑1 (MQ=255) cAACTTCGAAGACTCCATCTGCCTGTCCGATCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGAAATCACTGCGGACATCCCGAACGTGGGTGAAGc < 1:284/150‑1 (MQ=255) cTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGAAATCACTGCGGACATCCCGAACGTGGGTGAAGCCGc > 2:237982/1‑150 (MQ=255) tCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCGCGGCACAGAGGAAATCACTGCGGACATCCCGAACGTGGGTGAAGCCGCACTGAACAAg > 1:140366/1‑150 (MQ=255) aTCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGAAATCACTGCg < 2:334668/111‑1 (MQ=255) aTCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGAAATCACTGCg > 1:334668/1‑111 (MQ=255) | GTCCGCTGGTGAGCAAGGGTGACAAGGTTCAGCGTGGTGACATCATGGCCGACGGCCCGTCCACCGACATGGGTGAGCTGGCTCTGGGTCAGAACATGCGCATCGCGTTCATGGCGTGGAACGGCTTCAACTTCGAAGACTCCATCTGCCTGTCCGAGCGTGTTGTTCAGGAAGATCGCTTCACCACCATCCACATTCAGGAACTGACCTGTGTGGCGCGTGACACCAAGCTCGGCCCAGAGGAAATCACTGCGGACATCCCGAACGTGGGTGAAGCCGCACTGAACAAG > NC_002947/539943‑540232 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |