Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,357,611 | 0 | A | C | 31.0% | 115.5 / 14.6 | 58 | T184P (ACC→CCC) | ackA | acetate kinase |
Reads supporting (aligned to +/- strand): ref base A (27/13); new base C (0/18); total (27/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.69e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACC > REL606/2357561‑2357660 | cTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCa > 1:859076/1‑51 (MQ=255) cTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCa > 1:2723163/1‑51 (MQ=255) ttACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCAc < 1:417095/51‑1 (MQ=255) ttACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCAc > 1:4411621/1‑51 (MQ=255) ttACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCAc > 1:3848933/1‑51 (MQ=255) ttACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCAc > 1:4479280/1‑51 (MQ=255) aCAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCa < 1:2498569/51‑1 (MQ=255) aCAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCa < 1:6325309/51‑1 (MQ=255) cAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAg > 1:4563222/1‑51 (MQ=255) aaCCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGc > 1:4892416/1‑51 (MQ=255) aCCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGcc < 1:485925/51‑1 (MQ=255) ccTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGTCa < 1:2777700/51‑1 (MQ=255) ccTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCa > 1:4670072/1‑51 (MQ=255) cAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTa > 1:3701400/1‑51 (MQ=255) aaaGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCt < 1:6229299/49‑1 (MQ=255) aaaGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTa < 1:5789612/50‑1 (MQ=255) ggCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATg < 1:4857655/43‑1 (MQ=255) gCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAgg > 1:2616400/1‑51 (MQ=255) cATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGa > 1:1536021/1‑51 (MQ=255) cATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGa > 1:6757381/1‑51 (MQ=255) aTCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGaa > 1:3843547/1‑51 (MQ=255) tCCGTCGTTACGGCGCGCACGGCCCCAGCCACTTCTATGTAACCCAGGAAg < 1:3920140/51‑1 (MQ=255) ccGTCGTTACGGCGCGCACGGCCCCAGCCACTTc < 1:369448/34‑1 (MQ=255) ccGTCGTTACGGCGCGCACGGCCCCAGCCACTTCt < 1:1138249/35‑1 (MQ=255) ccGTCGTTACGGCGCGCACGGCCCCAGCCACTTCTATGTAAcc < 1:6122698/43‑1 (MQ=255) ccGTCGTTACGGCGCGCACGGCCCCAGCCACTTCTATGTAAc < 1:2634366/42‑1 (MQ=255) ccGTCGTTACGGCGCGCACGGCCCCAGCCACTTCTATGTAACCCAGGaa < 1:3958076/49‑1 (MQ=255) cgtcgtTACGGCGCGCACGGCCCCAGCCACTTCTATGTAAc < 1:2023996/41‑1 (MQ=255) cgtcgtTACGGCGCGCACGGCCCCAGCCACTTCTATGTAACCCAGGAAGCg < 1:2514882/51‑1 (MQ=255) gtcgtTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCgg > 1:2684894/1‑51 (MQ=255) tcgtTACGGCGCGCACGGCCCCAGCCACTTCTATGTAACCCAGGaa < 1:2234167/46‑1 (MQ=255) cgtTACGGCGCGCACGGCCCCAGCCACTTCTATGTAACCCAGGAAGCGGCa < 1:6608373/51‑1 (MQ=255) cgtTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCa > 1:117798/1‑51 (MQ=255) gtTACGGCGCGCACGGCACCAGCCACTTCTTTGTAACCCAGGAAGCGGCaa > 1:373992/1‑51 (MQ=255) gtTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCaa > 1:7006820/1‑51 (MQ=255) gtTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCaa > 1:6278914/1‑51 (MQ=255) gtTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCaa > 1:3534799/1‑51 (MQ=255) ttACGGCGCGCACGGCCCCAGCCACTTCTATGTAACCCAGGAAGCgg < 1:2232015/47‑1 (MQ=255) ttACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCaaa > 1:5570777/1‑51 (MQ=255) ggCGCGCACGGCCCCAGCCACTTCTATGTaa < 1:300137/31‑1 (MQ=255) ggCGCGCACGGCCCCAGCCACTTCTATGTAACCCAGGAAGCGGCaa < 1:2535581/46‑1 (MQ=255) ggCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGaa < 1:5550477/39‑1 (MQ=255) ggCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATg < 1:6064503/51‑1 (MQ=255) cgcgcACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCt > 1:4563965/1‑51 (MQ=255) cgcgcACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCt > 1:808016/1‑51 (MQ=255) cgcgcACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCt > 1:714135/1‑51 (MQ=255) cgcgcACGCCACCAGCCACTTCTATGTAACCCAGGAAGCGGCaaa < 1:1049992/45‑1 (MQ=255) gcgcACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTg > 1:3819196/1‑51 (MQ=255) gcACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGaa > 1:1230589/1‑51 (MQ=255) gcACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGaa > 1:5626312/1‑51 (MQ=255) cACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAAc > 1:3605873/1‑51 (MQ=255) cGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACaa > 1:4096675/1‑51 (MQ=255) ggCCCCAGCCACTTCTATGTAACCCAGGaa < 1:3465098/30‑1 (MQ=38) ggCCCCAGCCACTTCTATGTAACCCAGGAAGCGGCaaaa < 1:744339/39‑1 (MQ=255) ggCCCCAGCCACTTCTATGTAACCCAGGAAGCGGCaaa < 1:242515/38‑1 (MQ=255) ggCCCCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACaaa < 1:7185412/51‑1 (MQ=255) ggCCCCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACaa < 1:7293633/50‑1 (MQ=255) cACCAGCCACTTCTATGTAACCCAGGAAGCgg < 1:5075990/32‑1 (MQ=255) cACCAGCCACTTCTATGTAACCCAGGAAGCGGCaaaaa < 1:464630/38‑1 (MQ=255) cACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAAcc > 1:4612126/1‑51 (MQ=255) | CTTACAACCTGTACAAAGAGCACGGCATCCGTCGTTACGGCGCGCACGGCACCAGCCACTTCTATGTAACCCAGGAAGCGGCAAAAATGCTGAACAAACC > REL606/2357561‑2357660 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |