Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 3,784,658 | 0 | A | C | 55.6% | 8.8 / 19.3 | 18 | V567G (GTG→GGG) | dmdC_4 | 3‑methylmercaptopropionyl‑CoA dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (0/8); new base C (10/0); total (10/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.29e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGAT‑GGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGA > ACB143/3784529‑3784754 | ggACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGAT‑GGCCGCCCCCGGcc > 8:296344/1‑136 (MQ=255) aCCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGGGGGCGCGTGGAT‑GGCCGCcccccggccg > 2:321051/1‑130 (MQ=255) aCCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGGGGGCGCCTCGGA‑GGCCCCCCCCAGcccg > 8:135879/1‑134 (MQ=255) aCCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGGCCCGGCGCT‑GGCCGCCCCCAGCCgg > 2:59932/1‑136 (MQ=255) aTGCGCTCCCGCCCCTCGGCACGCCAGCCGCCGCCGGGGTTCATGACCCCCTCGCACGTGGGTTTCACACCGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCTAG‑GGCCGCCACCTGCCGGTCGAGGTc < 8:343411/136‑1 (MQ=255) aTGCGCTACCGCCCCTCGGCACGCCCGCCGCCGCCGGTGTTCATGACGACCTCGCACGGGGGTTTCAGAACGATTCGCAAATGAATGGCACCGGCGTGGCGGCCGCGGCGAT‑GGCCGCCACCTGCCGGTCGAGGTc < 8:284942/136‑1 (MQ=255) aCCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGGGGCT‑GGCCGCCCCCTGCCCGGCGAgggcgggccg > 5:358815/1‑133 (MQ=255) cACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGCTTCGGCAATGAATGGCAACGGCGTGGCGGCCCCGGCGGG‑GGGCGCCCCCGGGCGGTGGCGGGCGGGCCGCCGGGAGcccggg > 1:221453/1‑132 (MQ=255) ccgcagccgGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGAT‑GGCCGCCCCCTGCCGGTCGAGGTCGGGCAGCAGGTATCCGAGGCAGAAccc > 2:286727/1‑136 (MQ=255) ggTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCCGCCGCGTCGCT‑GGCCGCCCCCGGCCGGTCGAGGGCGGGCCGCAGGGAGCCGCGGAACCACCCCAcccctt > 2:39773/1‑131 (MQ=255) gTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGCT‑GGCCGCCCCCAGCCCGTCGAGGTCGGGCCGCACGGATCCGCGGAAGAACCCAACGCATTCCGGCGTTTACCGCCGCTgcgc > 3:367994/1‑136 (MQ=255) ggTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCCGCCGT‑GGCCCCCCCCAGCCCGTCGCGGGGGGGGCGCAGGGATCCGAGGGAGAACCCAACCCGTTACAGCGTTTACCCGCGCTGcgcccc > 6:346935/1‑134 (MQ=255) ttCAGAACGATTCGGCAATGAATGGCAACGGCGTGGGGGCCGCGGCGAG‑GGGCGCCCCCACCCGGGCGGGGGCGGGGCGCCGGGGGGCGAGGAAGAACCCCCACCCTTGCTGCTTTTTACCGCGCGGCGCcccccc > 2:136349/1‑136 (MQ=255) gAACGATTCGGCAATGAATGCCAACGGCGGGGCGGCCGCGGCGAG‑GGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACt < 4:68253/136‑1 (MQ=255) caaCGATGCGGCAATGAAGGCCACCGGCGGGGCGGCCGCGTCGAT‑GGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCACTTATATGCGCCGCCCCCACt < 1:51365/135‑1 (MQ=255) aCGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGAT‑GGCCGCCACCTGCCGGTCGATGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTgg < 1:200427/136‑1 (MQ=255) cgtttCCGCAATGTATGGCCACGGGGGGGGGGGCGGGGCGTTGGGCCCCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTgg < 1:136349/133‑1 (MQ=255) ttCGGCAATTAATGGCAACGGGGTGGCGGCCGCGGCGAT‑GGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAAc < 1:286727/136‑1 (MQ=255) atgaatgGCAACGGCGTGGCGGCCGCGTCGAT‑GGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGa < 6:358815/136‑1 (MQ=255) | GGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGCACGTGGGTTTCAGAACGATTCGGCAATGAATGGCAACGGCGTGGCGGCCGCGTCGAT‑GGCCGCCACCTGCCGGTCGAGGTCGGGCAGCAGGTAGCCGAGGAAGAACCCAACGCTTTCCTGCTTTTCCCGGCGCTGCGCCGCCCCCACTGGCAACGCCTGGA > ACB143/3784529‑3784754 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |