Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 417,321 | 0 | T | G | 29.0% | 49.2 / 10.3 | 31 | T27P (ACA→CCA) | ECB_00378 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (7/15); new base G (9/0); total (16/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.16e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGGCGTCGGGACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTCATACATAATACGGCTA > REL606/417249‑417395 | tGGCGTCGGGACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGtttt < 1:634155/76‑1 (MQ=255) cGTCGGGACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGgtttttc > 2:753438/1‑74 (MQ=255) cGGGACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAgg > 1:1050435/1‑76 (MQ=255) cGGGACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTCCCgg > 2:87066/1‑76 (MQ=255) gACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCAt < 1:960047/76‑1 (MQ=255) aCTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCAtt < 2:570540/76‑1 (MQ=255) aCTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCAtt < 1:879333/76‑1 (MQ=255) ntAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTTCCGGCAttt > 1:269754/2‑75 (MQ=255) tAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTAt < 1:1155481/76‑1 (MQ=255) aGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTTCAGGCATTATg > 1:32500/1‑76 (MQ=255) tttGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTTCCGGCATTATGAcc > 2:240071/1‑76 (MQ=255) gTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCttt < 1:740559/76‑1 (MQ=255) gTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTCCAGGCATTATGACCTCttt > 1:810060/1‑76 (MQ=255) ttCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTg < 1:669687/76‑1 (MQ=255) cgcgCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACt < 1:659018/76‑1 (MQ=255) cgcgCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACt < 1:188902/76‑1 (MQ=255) cgcATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTCCAGGCATTATGACCTCTTTGATTAACtga > 2:842540/1‑76 (MQ=255) aTTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACtgaatg < 2:559687/76‑1 (MQ=255) tttttCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACtgaatga < 1:12676/76‑1 (MQ=255) ttttCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACtgaatgaa > 1:1131110/1‑76 (MQ=255) ttCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAAt > 1:389569/1‑76 (MQ=255) ttCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAAt < 1:1056233/76‑1 (MQ=255) tCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATc > 1:195162/1‑76 (MQ=255) tCAGCGTATTATCGACATCGGGTAGAGTGGCTGGTTTTCCAGGCATTATGACCTCTTTGATTTACTGAATGAAATc > 2:1076578/1‑76 (MQ=255) gTTTTATCGATATCGGGTAGAGTGGCTGGTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAgg > 2:689/1‑76 (MQ=255) gTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAgg > 1:258436/1‑76 (MQ=255) ttttATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAggg < 2:827524/76‑1 (MQ=255) ttATCGACATCGGGTAGAGTGGCTGGTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAg > 2:877373/1‑76 (MQ=255) cGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGa < 2:259673/76‑1 (MQ=255) gACATCGGGTAGAGTGGCTGGTTTTCCAGGCATTATGACCTCTTTTATTAACTGAATGAAATCGTAGGGAGTCGAt > 2:213584/1‑76 (MQ=255) gggTAGAGTGGCTGTTTTTCCAGGCATTATGAGCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTcata < 1:155901/76‑1 (MQ=255) gTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTcataca < 1:134332/76‑1 (MQ=255) tAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTTGATCTcatacat > 1:509742/1‑76 (MQ=255) tAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTcatacat > 1:213975/1‑76 (MQ=255) gTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTCATACATAATACGGCTa > 1:1028062/1‑76 (MQ=255) | TGGCGTCGGGACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTCATACATAATACGGCTA > REL606/417249‑417395 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |