| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,647,315 | C→A | 20.8% | L162L (CTC→CTA) | dbpA → | 23S rRNA specific ATP‑dependent RNA helicase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,647,315 | 0 | C | A | 20.8% | 62.0 / 3.0 | 23 | L162L (CTC→CTA) | dbpA | 23S rRNA specific ATP‑dependent RNA helicase |
| Reads supporting (aligned to +/- strand): ref base C (10/8); new base A (2/2); total (12/11) | |||||||||||
| 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.80e-01 | |||||||||||
TCGCTTCGCTGGAGCACGGCGCGCACATCATCGTCGGCACCCCAGGGCGCATCCAGCAGCACCTGGACAAAGGCACCCTGGTCCTGGACGGGTTGAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTCGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCCTGCTGTTCTCGGCCACCTACCCGGCCGGCATCAAGCAACTGGCCGCCGACTTCATGCGCAACCCT > NC_002947/5647185‑5647444 | tcgcttcgctGGAGCACGGCGCGCACATCATCGTCGGCACCCCAGGGCGCATCCAGCAGCACCTGGACAAAGGCACCCTGGTCCTGGACGGGTTGAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTCGACAt < 2:202100/136‑1 (MQ=255)tcgcttcgctGGAGCACGGCGCGCACATCATCGTCGGCACCCCAGGGCGCATCCAGCAGCACCTGGACAAAGGCACCCTGGTCCTGGACGGGTTGAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTCGACAt < 2:199450/136‑1 (MQ=255) ttcgctGGAGCACGGCGCGCACATCATCGTCGGCACCCCAGGGCGCATCCAGCAGCACCTGGACAAAGGCACCCTGGTCCTGGACGGGTTGAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTc < 1:193441/127‑1 (MQ=255) ttcgctGGAGCACGGCGCGCACATCATCGTCGGCACCCCAGGGCGCATCCAGCAGCACCTGGACAAAGGCACCCTGGTCCTGGACGGGTTGAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTc > 2:193441/1‑127 (MQ=255) gcACATCATCGTCGGCACCCCAGGGCGCATCCAGCAGCACCTGGACAAAGGCACCCTGGTCCTGGACGGGTTGAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTCGACATGGGCTTCTTCGACGCCATCGcc > 2:420568/1‑136 (MQ=255) aGGGCGCATCCAGCAGCACCTGGACAAAGGCACCCTGGTCCTGGACGGGTTGAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTCGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCg > 1:78914/1‑136 (MQ=255) cagcagCACCTGGACAAAGGCACCCTGGTCCTGGACGGGTTGAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTCg < 1:308976/79‑1 (MQ=255) cagcagCACCTGGACAAAGGCACCCTGGTCCTGGACGGGTTGAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTCg > 2:308976/1‑79 (MQ=255) cagcagCACCTGGACAAAGGCACCCTGGTCCTGGACGGGTTGAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTCGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCgccgcc < 1:358059/136‑1 (MQ=255) cagcagCACCTGGACAAAGGCACCCTGGTCCTGGACGGGTTGAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTCGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCgccgcc > 2:411083/1‑136 (MQ=255) tGGACAAAGGCACCCTGGTCCTGGACGGGTTGAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTAGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCctgct > 1:514923/1‑136 (MQ=255) tGGACAAAGGCACCCTGGTCCTGGACGGGTTGAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTAGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCctgct > 1:113343/1‑136 (MQ=255) gccTGGACGGGTTGAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTCGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCCTGCTGTTCTCGGCCACCTAccc < 1:40630/135‑1 (MQ=255) gAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTCGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCCTGCTGTTCTCGGCCACCTACCCGGCCGGCATCAAg > 1:70573/1‑136 (MQ=255) ggtggtGGAGGTAGCCGATCGCTTGCTTGAGTTGGGGTTCTTGGTGGCGATCGCGAGGTTGATCGGCTAGACCCGGTCGCGGCGCCAGACCCTGCTGTTCTCGGCCACCTAGCCGGCCGGCATCATGCAACTGgcc < 2:319808/136‑1 (MQ=255) ggTGCTGGACGAAGCCGACCGCATGCTCGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCCTGCTGTTCTCGGCCACCTACCCTGCCGGCATCAAGCAACTGgcc < 2:435167/136‑1 (MQ=255) ggTGCTGGACGAAGCCGACCGCACGCTCGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCCTGCTGTTCTCGGCCACCTACCCTGCCGGCATCAAGCAACTGgcc > 1:435167/1‑136 (MQ=255) gAAGCCGACCGCATGCTCGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCCTGCTGTTCTCGGCCACCTACCCGGCCGGCATCAAGCAACTGGCCGCCGACTTCa > 1:411129/1‑136 (MQ=255) gAAGCCGACCGCATGCTCGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCCTGCTGTTCTCGGCCACCTACCCGGCCGGCATCAAGCAACTGGCCGCCGACTTCa > 1:319808/1‑136 (MQ=255) gAAGCCGACCGCATGCTCGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCCTGCTGTTCTCGGCCACCTACCCGGCCGGCATCAAGCAACTGGCCGCCGACTTCa > 2:36916/1‑136 (MQ=255) aaGCCGACCGCATGCTCGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCCTGCTGTTCTCGGCCACCTACCCGGCCGGCATCAAGCAACTGGCCGCCGACTTCAt < 1:420568/136‑1 (MQ=255) gCATGCTCGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCCTGCTGTTCTCGGCCACCTACCCGGCCGGCATCAAGCAACTGgcc > 1:462466/1‑116 (MQ=255) gCATGCTCGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCCTGCTGTTCTCGGCCACCTACCCGGCCGGCATCAAGCAACTGgcc < 2:462466/116‑1 (MQ=255) cATGCTAGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCCTGCTGTTCTCGGCCACCTACCCGGCCGGCATCAAGCAACTGGCCGCCGACTTCATGCGCAACCCt < 2:113343/136‑1 (MQ=255) cATGCTAGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCCTGCTGTTCTCGGCCACCTACCCGGCCGGCATCAAGCAACTGGCCGCCGACTTCATGCGCAACCCt < 2:514923/136‑1 (MQ=255) | TCGCTTCGCTGGAGCACGGCGCGCACATCATCGTCGGCACCCCAGGGCGCATCCAGCAGCACCTGGACAAAGGCACCCTGGTCCTGGACGGGTTGAACACCCTGGTGCTGGACGAAGCCGACCGCATGCTCGACATGGGCTTCTTCGACGCCATCGCCAGCATCATCGGCAAGACCCCGTCGCGCCGCCAGACCCTGCTGTTCTCGGCCACCTACCCGGCCGGCATCAAGCAACTGGCCGCCGACTTCATGCGCAACCCT > NC_002947/5647185‑5647444 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |