Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,983,700 | 0 | T | G | 36.2% | 61.9 / 29.3 | 47 | T397P (ACC→CCC) | ECB_03713 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (5/25); new base G (17/0); total (22/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.61e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATT > REL606/3983625‑3983773 | gAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGt < 2:1388578/76‑1 (MQ=255) aaGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTg < 2:2049675/76‑1 (MQ=255) tgtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCt < 2:1841831/76‑1 (MQ=255) tgtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCt < 1:725839/76‑1 (MQ=255) gtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGcccc > 2:988727/1‑74 (MQ=255) tAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑cc > 1:1080541/1‑76 (MQ=255) aaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCca < 2:2005985/76‑1 (MQ=255) tGCGGGGTTACCACCGAGGGTGATTGCTTTCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCcaca < 2:380159/76‑1 (MQ=255) tGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCccca > 1:1901980/1‑76 (MQ=255) cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCC‑CCCACcaa > 2:866122/1‑76 (MQ=255) ggggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAAt < 1:1399006/76‑1 (MQ=255) ggggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAAt < 2:845101/76‑1 (MQ=255) gTTACCGCCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTg < 1:936439/76‑1 (MQ=255) gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTg < 1:487681/76‑1 (MQ=255) gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCCAAAttt > 2:818248/1‑76 (MQ=255) ttACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGa < 1:939069/76‑1 (MQ=255) ttACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCCCAAAATTGa > 1:508611/1‑76 (MQ=255) tgccaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGAc < 1:734059/74‑1 (MQ=255) tACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCC‑CCCCCAAATTTTAc > 2:2093046/1‑76 (MQ=255) ccaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGa < 1:1429080/76‑1 (MQ=255) ccaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCCCCAAATTGACGa > 2:1227361/1‑76 (MQ=255) gAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCCCAAAATTTGCGATAGCt > 1:1891331/1‑76 (MQ=255) aGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCCCAAAATTGACGATAGCTg > 2:720936/1‑76 (MQ=255) gggTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTgg < 2:2101361/76‑1 (MQ=255) gggTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCACAAATTTTACGATCGCTgg > 1:1510530/1‑76 (MQ=255) tGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCCCCAAATTTACGATAGCCGggga > 1:796538/1‑76 (MQ=255) gATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCCCAAATTTGACGATCGCTGgggaa > 1:423283/1‑75 (MQ=255) gATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCCCAAATTTTACGATCGCTGGAGAt > 1:186118/1‑76 (MQ=255) ttgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGa < 2:527265/76‑1 (MQ=255) tgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCCGGGGAAGaa > 2:264522/1‑75 (MQ=255) tgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCCCAAAATTTACGATCGCCGGGGGTGAc > 1:1508003/1‑76 (MQ=255) cttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCCCAAAATTTACGATCGCTGGGGATGACGc > 1:1061090/1‑76 (MQ=255) tttcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCCCAAAATTTACGATCGCTGGGGGAGACGcc > 1:1262925/4‑76 (MQ=255) ttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCCCCAAATTTACGAACGCTGGGGATGACGcc > 1:1432781/1‑76 (MQ=255) gcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCACAAAATTTACGATCGCTGGGGATGACGcccc > 2:1186431/1‑74 (MQ=255) gcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCC‑CCCACCAATTTTACGAACCCTGGGGAAGACGcccc > 2:676414/1‑74 (MQ=255) aTCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATcg < 2:1177854/76‑1 (MQ=255) aTCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATcg < 2:1972977/76‑1 (MQ=255) aTCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCC‑CCCCCAAAATTTACGATCGCTGGGGAGGACGCCCTcg > 1:1269452/1‑76 (MQ=255) tCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCC‑CCCACAAAATTTACCATCGCCGGGGAAGACGCCATcgc > 2:1667487/1‑76 (MQ=255) cTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCACAAAATTTACGATCGCTGGGGATGACGCCATCGCGa > 1:848432/1‑76 (MQ=255) tACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGaa < 1:562883/76‑1 (MQ=255) tACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGaa < 1:40962/76‑1 (MQ=255) tGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccacc < 1:1363051/76‑1 (MQ=255) gTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccaccc < 2:1522220/76‑1 (MQ=255) gTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccaccc > 1:175592/1‑76 (MQ=255) ttATTGGTGGCATCGGCAGGGGGGCCT‑CCCACAAAATTTACGATCGCTGGAGATGACGCCCTCGCGAAGccacccc > 1:1645277/1‑75 (MQ=255) tggttgCATCGGCAGGGGGGCCT‑CCCACAAATTTTTCGATCGCTGGAGATGAAGCCATCGCCACGCCACCCCCAcc > 2:2197031/1‑76 (MQ=255) tggtggCATCGGCAGGGGGGCCT‑CCCACAAATTTTACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACAcc > 2:186750/1‑76 (MQ=255) ggtggCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCt < 1:2132580/76‑1 (MQ=255) ggtggCATCGGCAGGGGGGCCT‑CCCACAAATTTGACGATCGCTGGGGATGACGCCATCGCGAAGACACCCACAccc > 2:1042360/1‑75 (MQ=255) ggtggCATCGGCAGGGGGGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCTACAccc > 2:1360603/1‑75 (MQ=255) gtggCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTg < 1:1393252/76‑1 (MQ=255) ggCATCGGCAGGGGTGCCT‑CCCACAAAATTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCcgcg > 1:1242701/1‑76 (MQ=255) tCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGttt < 2:796538/76‑1 (MQ=255) ggCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCa < 2:35607/76‑1 (MQ=255) cAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGt > 1:679939/1‑76 (MQ=255) aGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa > 1:576623/1‑76 (MQ=255) ggggTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTaa > 1:2026808/1‑76 (MQ=255) ggTGCCT‑CCCACAAATCTGACGATCGCTGAAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAtt < 2:365717/76‑1 (MQ=255) | GAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATT > REL606/3983625‑3983773 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |