Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,983,700 | 0 | T | G | 31.6% | 63.9 / 19.2 | 38 | T397P (ACC→CCC) | ECB_03713 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (3/23); new base G (12/0); total (15/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.68e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTGC > REL606/3983626‑3983775 | aaGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTg < 1:1345308/76‑1 (MQ=255) gtgtAAATGCGGGGTTGCCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGcc < 1:724560/76‑1 (MQ=255) tAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTcc < 2:431073/76‑1 (MQ=255) tAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTcc > 1:407183/1‑76 (MQ=255) aaaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGcccccc > 2:548073/1‑76 (MQ=255) aaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCca < 2:1015661/76‑1 (MQ=255) aaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTcccc > 2:1509914/1‑75 (MQ=255) tGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCca‑ca < 1:1328323/76‑1 (MQ=255) cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCC‑Ca‑aa > 2:197872/1‑76 (MQ=255) cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCGCCCA‑Ca‑aa > 1:803552/1‑76 (MQ=255) tACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCCGGGGGGCCCCCCC‑CACAAttttc > 1:1196745/1‑74 (MQ=255) tACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGAc < 1:647894/76‑1 (MQ=255) tACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGAc < 1:1217667/76‑1 (MQ=255) accaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACCA‑AAtttacc > 1:312592/1‑73 (MQ=255) ccaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCCCCCA‑CC‑AATTTGACGa > 2:1987221/1‑76 (MQ=255) caccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGAt < 2:416309/76‑1 (MQ=255) cGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCCCCCC‑CA‑AATTTTACGATccc > 1:1106462/1‑74 (MQ=255) gAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCt < 2:1648975/76‑1 (MQ=255) gTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGgag < 2:1198281/76‑1 (MQ=255) aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATg < 1:1177084/76‑1 (MQ=255) aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATg < 2:355051/76‑1 (MQ=255) ttgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGa < 2:1152770/76‑1 (MQ=255) ttgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCA‑CA‑AATTTTACGATCGCCGGGGATGa > 1:1277509/1‑76 (MQ=255) tgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGAc < 2:927801/76‑1 (MQ=255) tgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCC‑CA‑AAATTTTAGAGCGCTGGGGAAGAc > 1:201034/1‑76 (MQ=255) tgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGAc > 2:292176/1‑76 (MQ=255) gcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACg < 1:906933/76‑1 (MQ=255) ttgcTTCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACGcc < 1:1428903/76‑1 (MQ=255) ttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACGcc < 2:663335/76‑1 (MQ=255) cATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACGCCATc < 1:930441/76‑1 (MQ=255) tCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCC‑CA‑AATTTTACCATCGCTGGAGAAGACCCCAtacc > 1:1910680/1‑73 (MQ=255) ccTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCCCCCC‑CA‑AAATTGACGATCGCTGGGGATGACGCCACcgcg > 1:1367162/1‑76 (MQ=255) cTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccac > 1:1466011/1‑76 (MQ=255) tGTTATTGGTGGCATCGGCAGGGGGGCCTCCCA‑CA‑AATTTGATGATCGCTGGGGATGACGCCATCGCGAAGccccc > 1:127551/1‑76 (MQ=255) gTTATTGGTGGCATCGGCAGGGGGGCCTCCCC‑CA‑AATTTGACGAACGCTGGGGATGACGCCATCGCGAAGccaccc > 2:1386180/1‑76 (MQ=255) tATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccacccac < 1:1257862/76‑1 (MQ=255) ttGGTGGCATCGGCAGGGGGGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCacac > 2:1423811/1‑76 (MQ=255) ggtggCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCt < 1:1914650/76‑1 (MQ=255) aTCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGtt < 1:548073/76‑1 (MQ=255) cGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTc < 2:1087927/76‑1 (MQ=255) gCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAg > 1:1879300/1‑76 (MQ=255) aGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGCTTCAGTa > 1:333142/1‑76 (MQ=255) gaggTGCCTTCCA‑CA‑AAGTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACACTCACCTGCGGTTTCAGTaa < 2:264702/74‑1 (MQ=255) ggTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAtt < 1:462989/76‑1 (MQ=255) tGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTgc > 1:436138/1‑76 (MQ=255) | AAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCA‑CA‑AATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTGC > REL606/3983626‑3983775 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |