Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,983,700 | 0 | T | G | 34.2% | 64.8 / 23.1 | 41 | T397P (ACC→CCC) | ECB_03713 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (7/20); new base G (14/0); total (21/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.40e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTG > REL606/3983630‑3983774 | gtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGcccc > 1:1746678/1‑74 (MQ=255) gtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑c > 1:1679771/1‑76 (MQ=255) gtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑c > 1:1294471/1‑76 (MQ=255) aaaTGCGGGGTTACTACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑ccc > 2:118961/1‑76 (MQ=255) aaaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑ccc < 1:158576/76‑1 (MQ=255) aaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCca < 1:1616894/76‑1 (MQ=255) cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACaaa < 1:152113/76‑1 (MQ=255) cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCC‑CCCACcaa > 1:1838118/1‑76 (MQ=255) ggggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCACAAAt > 2:101194/1‑76 (MQ=255) gggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCCCaaaat > 2:86489/1‑74 (MQ=255) ttACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCCCAAAATTGa > 2:939697/1‑76 (MQ=255) ttACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCCCAAAtttta > 2:1408138/1‑74 (MQ=255) tACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCCAAAtttcc > 2:306402/1‑74 (MQ=255) accaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACg < 1:1436645/76‑1 (MQ=255) accGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCCCAAAATTTACGATc > 2:783808/1‑76 (MQ=255) accGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCCCAAAATTTACGATc > 2:1132448/1‑76 (MQ=255) gAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCt < 2:671269/76‑1 (MQ=255) aGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTg < 1:1432926/76‑1 (MQ=255) aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATg < 1:1172467/76‑1 (MQ=255) aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCACAAAATTTACGATCGCTGGGGATg > 1:422500/1‑76 (MQ=255) ttgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGa < 2:724065/76‑1 (MQ=255) ttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCCCCAAATTTTCGATCGCTGGGGATGAAGcc > 1:1144405/1‑76 (MQ=255) tgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCACAAAATTGACGATCGCCGGGGATGACGCCa > 1:109303/1‑76 (MQ=255) tACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCCCAAAATTTACGATCGCCGGGGAAGACGCCATcgcgga > 1:1470000/1‑74 (MQ=255) tACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCACAAAATTTACGATCGCTGGGGATGACGCCCTcgcgga > 1:1043619/1‑74 (MQ=255) tACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCC‑CCCCCAAAATTTTCGGTCGCTGGGGGTGACCCCATtgcgga > 1:2538/1‑74 (MQ=255) ggCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcc < 2:570162/76‑1 (MQ=255) ggCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcc < 2:116294/76‑1 (MQ=255) tGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCCCAAAATTTACGATCGCTGGAGATGAAGCCATCGCGGAGccacc > 2:766235/1‑76 (MQ=255) ttATTGGTGGCATCGGCAGGGGGGCCT‑CCCCCAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccccccc > 2:871422/1‑75 (MQ=255) ttGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCacac < 2:685699/76‑1 (MQ=255) ttGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCacac < 1:1452134/76‑1 (MQ=255) tggtggCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACAcc > 1:517213/1‑76 (MQ=255) ggtggCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACAccc > 1:1293207/1‑75 (MQ=255) tggCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGc < 1:1247134/76‑1 (MQ=255) tggCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGc < 1:1459083/76‑1 (MQ=255) ggCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCg < 1:118961/76‑1 (MQ=255) ggCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCa < 2:1504957/76‑1 (MQ=255) ggCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCa < 1:71590/76‑1 (MQ=255) cAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGt > 2:104083/1‑76 (MQ=255) aGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa > 1:712893/1‑76 (MQ=255) ggggTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTaa < 1:871422/76‑1 (MQ=255) ggTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAtt < 1:306402/76‑1 (MQ=255) ggTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAtt > 1:1321702/1‑76 (MQ=255) ggTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAtt < 1:783808/76‑1 (MQ=255) ggTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCAGCGGTTTCAGTAAtt > 1:1116147/1‑76 (MQ=255) gTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTg > 1:860275/1‑76 (MQ=255) | GTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTG > REL606/3983630‑3983774 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |