Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,931,972 | G→A | 18.2% | R202W (CGG→TGG) | rluA ← | bifunctional tRNA pseudouridine(32) synthase TrmQ/ribosomal large subunit pseudouridine synthase RluA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,931,972 | 0 | G | A | 18.2% | 46.6 / 4.7 | 22 | R202W (CGG→TGG) | rluA | bifunctional tRNA pseudouridine(32) synthase TrmQ/ribosomal large subunit pseudouridine synthase RluA |
Reads supporting (aligned to +/- strand): ref base G (13/5); new base A (2/2); total (15/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.65e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01 |
CCCGCTACCACAGGTATTGCACAGCCCCTGAGCGCAGTGCGGCCCTTGCAGGAGCGGCTTTAGCCGCGAAGAGGCCGGTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCATGCAGGCACAGCCGCTCATGGGCCGCCAGGGCTTGCGGGTTGGCGTACAGGCGATCACCCAGCAGCGGGTGCCCGATCGACAGCATGTGCACG > NC_002947/1931846‑1932105 | cccGCTACCACAGGTATTCCACAGCCCCTGAGCGCCGGGCGGCCCTTCCCGGCGCGGCTTTCGCCGCGAAGAGGCCGGTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGa < 1:346255/136‑1 (MQ=255) aGCCCCTGAGCGCAGTGCGGCCCTTGCAGGAGCGGCTTTAGCCGCGAAGAGGCCGGTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGc > 1:103603/1‑136 (MQ=255) cTGAGCGCAGTGCGGCCCTTGCAGGAGCGGCTTTAGCCGCGAAGAGGCCGGTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGa > 1:133159/1‑136 (MQ=255) cTGAGCGCAGTGCGGCCCTTGCAGGAGCGGCTTTAGCCGCGAAGAGGCCGGTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGa > 1:257333/1‑136 (MQ=255) cgcAGTGCGGCCCTTGCAGGAGCGGCTTTAGCCGCGAAGAGGCCGGTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCAt > 2:343888/1‑136 (MQ=255) cgcAGTGCGGCCCTTGCAGGAGCGGCTTTAGCCGCGAAGAGGCCGGTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCAt > 2:101285/1‑136 (MQ=255) cGGCCCTTGCAGGAGCGGCTTTAGCCGCGAAGAGGCCGGTACAGGCCCTCCAAAATCAGAACGCCGCAGGGCACACGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGACGCAAGCAGGca > 1:314067/1‑136 (MQ=255) cTTTAGCCGCGAAGAGGCCGGTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTaaa > 2:99054/1‑94 (MQ=255) cTTTAGCCGCGAAGAGGCCGGTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTaaa < 1:99054/94‑1 (MQ=255) cTTTAGCCGCGAAGAGGCCGGTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCATGCAGGCACAGCCGCTCATGGgccgc > 1:136186/1‑136 (MQ=255) aaGAGGCCGGTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCATGCAGGCACAGCCGCTCATGGGCCGCCAGGGCTTGCg < 2:480455/136‑1 (MQ=255) agagGCCGGTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCACTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCATGCAGGCACAGCCGCTCATGGGCCGCCAGGGCTTGCgg > 1:377345/1‑136 (MQ=255) agagGCCGGTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCACTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCATGCAGGCACAGCCGCTCATGGGCCGCCAGGGCTTGCgg > 1:374982/1‑136 (MQ=255) gagGCCGGTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCATGTGGGCACAGCCGCACATGGGCCGCCAGGCCTTGCggg > 2:313876/1‑136 (MQ=255) ggTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCATGCAGGCACAGCCGCTCATGGGCCGCCAGGGCTTGCg > 1:500677/1‑128 (MQ=255) ggTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCATGCAGGCACAGCCGCTCATGGGCCGCCAGGGCTTGCg < 2:500677/128‑1 (MQ=255) ggCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCACTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCATGCAGGCACAGCCGCTCATGGGCCGCCAGGGCTTGCGGGTTGGCGTACAgg < 2:374982/136‑1 (MQ=255) ggCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCACTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCATGCAGGCACAGCCGCTCATGGGCCGCCAGGGCTTGCGGGTTGGCGTACAgg < 2:377345/136‑1 (MQ=255) aaaaTCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCATGCAGGCACAGCCGCTCATGGGCCGCCAGGGCTTGCGGGTTGGCGTACAGGCGATCAcc > 1:321668/1‑136 (MQ=255) aaaaTCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCATGCAGGCACAGCCGCTCATGGGCCGCCAGGGCTTGCGGGTTGGCGTACAGGCGATCAcc > 1:240493/1‑136 (MQ=255) ggCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCATGCAGGCACAGCCGCTCATGGGCCGCCAGGGCTTGCGGGTTGGCGTACAGGCGATCACCCAGCAGCGGGTGCCCGAt < 1:378735/136‑1 (MQ=255) ccGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCATGCAGGCACAGCCGCTCATGGGCCGCCAGGGCTTGCGGGTTGGCGTACAGGCGATCACCCAGCAGCGGGTGCCCGATCGACAGCATGTGCAcg > 1:502880/1‑136 (MQ=255) | CCCGCTACCACAGGTATTGCACAGCCCCTGAGCGCAGTGCGGCCCTTGCAGGAGCGGCTTTAGCCGCGAAGAGGCCGGTACAGGCCCTGCAAAATCAGAACGGCGCAGGGCACTCGAACTTCAACCGCTGCCCGGACACCGGGTGGGTAAAGCTCAGCATCGAGGCATGCAGGCACAGCCGCTCATGGGCCGCCAGGGCTTGCGGGTTGGCGTACAGGCGATCACCCAGCAGCGGGTGCCCGATCGACAGCATGTGCACG > NC_002947/1931846‑1932105 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |