| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 1,825,970 | 2 bp→CT | 100% | coding (742‑743/1422 nt) | PP_1624 → | group II intron‑encoding maturase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 1,825,970 | 0 | T | C | 88.9% | 18.8 / ‑3.6 | 9 | S248M (TCG→CTG) | PP_1624 | group II intron‑encoding maturase |
| Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (6/2); total (7/2) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
| * | NC_002947 | 1,825,971 | 0 | C | T | 88.9% | 19.2 / ‑3.7 | 9 | S248M (TCG→CTG) | PP_1624 | group II intron‑encoding maturase |
| Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (6/2); total (7/2) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CGGGATCATGTCGGGCGGTGTCGTCAGCCCACGGCAGGAGGGGACGCCGCAAGGCGGCCCGCTCTCGCCGTTGCTGTCGAATATCCTGCTCGACGAACTCGACCGCGAGCTGGAGCGGCGGGGCCATCGCTTCGTACGTTATGCCGATGATGCGAACATCTATGTACGCAGTCCCCGGGCCGGCGAACGGGTG > NC_002947/1825906‑1826098 || cGGGATCATGTCGGGCGGTGTCGTCAGCCCACGGCAGGAGGGGACGCCGCAAGGCGGCCCGCTTTCGCCGTtgctgtcgaatatcctgctcgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttatctttgac > 2:6761340‑M1/1‑71 (MQ=255) aaggcggaggggacgCCGCAAGGCGGCCC‑CTCCTGCCGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtcgtga < 2:1899671‑M1/145‑112 (MQ=255) gggACGCCGCAAGGCGGCCC‑CTCCTGCCGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtcgtgattgacccga > 2:4873168‑M1/1‑30 (MQ=255) ggACGCCGCAAGGCGGCCC‑CTCCTGCCGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtcgtgattgacccgat > 1:3163552‑M1/1‑29 (MQ=255) ggACGCCGCAAGGCGGCCC‑CTCCTGCCGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtcgtgattgacccgat > 1:4232015‑M1/1‑29 (MQ=255) ggACGCCGCAAGGCGGCCC‑CTCCTGCCGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtcgtgattgacccgat > 1:811486‑M1/1‑29 (MQ=255) gACGCCGCAAGGCGGCCC‑CTCCTGCCGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtcgtgattgacccgatc > 2:3477570‑M1/1‑28 (MQ=255) gACGCCGCAAGGCGGCCC‑CTCCTGCCGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctt > 1:1238954‑M1/1‑28 (MQ=255) gACGCCGCAAGGCGGCCC‑CTCCTGCCGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctt < 2:1238954‑M1/101‑74 (MQ=255) || CGGGATCATGTCGGGCGGTGTCGTCAGCCCACGGCAGGAGGGGACGCCGCAAGGCGGCCCGCTCTCGCCGTTGCTGTCGAATATCCTGCTCGACGAACTCGACCGCGAGCTGGAGCGGCGGGGCCATCGCTTCGTACGTTATGCCGATGATGCGAACATCTATGTACGCAGTCCCCGGGCCGGCGAACGGGTG > NC_002947/1825906‑1826098 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |