Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,983,700 | 0 | T | G | 20.8% | 125.1 / 11.8 | 48 | T397P (ACC→CCC) | ECB_03713 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (10/28); new base G (10/0); total (20/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.82e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAA > REL606/3983627‑3983771 | aGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTg > 1:2082758/1‑75 (MQ=255) tgtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGcc < 1:1889313/75‑1 (MQ=255) tAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGGA‑GGGGGGcccc > 2:265513/1‑73 (MQ=255) tAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGcccc > 1:1208748/1‑73 (MQ=255) aaaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTcc > 1:2093727/1‑75 (MQ=255) aaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTccc < 1:1307981/75‑1 (MQ=255) aaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTccc < 1:1448934/75‑1 (MQ=255) aaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCTccc > 2:70023/1‑75 (MQ=255) aTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTcccc > 2:1849206/1‑74 (MQ=255) tGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCcac < 2:1998824/75‑1 (MQ=255) tGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCcac < 2:1060995/75‑1 (MQ=255) gggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCCCCCCCaaaa > 1:2125979/1‑74 (MQ=255) ggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCCCCAAAt > 2:685762/1‑75 (MQ=255) ttACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCTCCCCCAAATTTg > 2:1454795/1‑75 (MQ=255) ttACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCCCCCCCAAAtttt > 1:151914/1‑74 (MQ=255) caccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGa < 1:438067/75‑1 (MQ=255) accGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGAt < 2:1850200/75‑1 (MQ=255) ccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATc < 2:854432/75‑1 (MQ=255) ccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCTCCCCCACAATTGACGaac > 2:460540/1‑73 (MQ=255) cGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCg < 2:1415437/75‑1 (MQ=255) gggTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTg < 1:1546567/75‑1 (MQ=255) gggTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTg < 2:2088403/75‑1 (MQ=255) gTGATTGCTTGCATCCTACGGCTGTTGTTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGga < 1:632545/75‑1 (MQ=255) tGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCCCCCACAAAATTGACGAACGCCGgag > 2:340670/1‑75 (MQ=255) aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGGCTCCCCCAAATTTGACCAACGCCGgtgag > 1:1799046/1‑74 (MQ=255) aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGAt < 1:298741/75‑1 (MQ=255) aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGAt < 1:241222/75‑1 (MQ=255) aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGAt < 1:513151/75‑1 (MQ=255) ttgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATg < 1:1984350/75‑1 (MQ=255) ttgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGAGCTCCCCCAAATTTTACGAACCCCGGGGATg > 1:1318947/1‑75 (MQ=255) tgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGa < 2:1874297/75‑1 (MQ=255) tgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCTCCCCCAAATTTGACGATCGCTGGggagga > 2:1683943/1‑75 (MQ=255) gcttgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGAc < 1:665238/75‑1 (MQ=255) ttgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGc < 2:812161/75‑1 (MQ=255) ttgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCCAATTTGACGATCGCTGGGGAGGACGc > 1:1838207/1‑75 (MQ=255) ttgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAATTTGAAGATCGCTGGAGAAGACGc > 2:1623939/1‑75 (MQ=255) ttgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAAATTTAGGCTCGGTGGGGAT‑ACGcc > 1:2110128/1‑75 (MQ=255) tgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGcc < 2:149254/75‑1 (MQ=255) gcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCCCCCCCAAAATTTACGATCGCTGGGGATGACGccc > 2:344027/1‑74 (MQ=255) aTCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCACAAAATTTACGAACGCTGGAGATGAAGccccc > 2:1009484/1‑72 (MQ=255) cTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCCCCCCCAAAATTTACGATCGCTGGAGATGACGCCACcgcg > 2:1551188/1‑75 (MQ=255) tACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAATTTTACGATCGCTGGGGAGGAAGCCCTCCCGa > 1:2017615/1‑75 (MQ=255) cGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCCCAAAATTGACGATCGCTTGGGATGGAGCCAACGCCAAg > 1:445308/1‑75 (MQ=255) ggCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGc < 2:232111/75‑1 (MQ=255) ggCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGc < 1:917019/75‑1 (MQ=255) gCTGTTATTGGTGGCATCGGCA‑GGGGTGCCCCCCACAAATTTGACGATCGCTGGAGATGACGCCCTCGCGAAGcc > 1:959881/1‑75 (MQ=255) cTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGGTGAAGCCAACGCGAAggcc > 2:711959/1‑72 (MQ=255) gtgaTTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccacc < 1:32089/72‑1 (MQ=255) gTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccacc < 1:657552/75‑1 (MQ=255) ttATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAATTTTACGATCGCTGGGGAAGAAGCGCTCCCGGAGccaccc > 1:860486/1‑75 (MQ=255) ttATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccaccc > 1:235423/1‑75 (MQ=255) ttATTGGTGGCATCGGCA‑GGGGGGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccaccc > 1:288480/1‑75 (MQ=255) tATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccaccca > 2:184351/1‑75 (MQ=255) tATTGGTGGCATCGGCA‑GGGGGGCCTCCCACAAATTTGACGATCGCTGGGGATGACGCCCTCGCGAAGccacccc > 1:50175/1‑74 (MQ=255) ttGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAAATTGACGATCGCTGGGGATGACGCCATCGCGGAGCCACCccca > 2:1431142/1‑75 (MQ=255) ggtggCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACAcc < 2:733290/75‑1 (MQ=255) gtggCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCt < 1:1055376/75‑1 (MQ=255) cATCGGCA‑GGGGGGCCTCCCACAAATTTGACGATCCCTGGAGATGACGCCATCGCGCAGCCACCCACCCCTCCgg > 1:1651555/1‑75 (MQ=255) cGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGttt < 1:789578/75‑1 (MQ=255) cGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGttt < 1:344027/75‑1 (MQ=255) ggCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTc > 1:1952901/1‑75 (MQ=255) gCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCa > 2:1202291/1‑75 (MQ=255) cA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGTCGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAg < 1:954671/75‑1 (MQ=255) ggggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa > 2:1772218/1‑75 (MQ=255) ggggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa > 2:381342/1‑75 (MQ=255) ggggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa > 1:496826/1‑75 (MQ=255) gggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTaa > 1:1828066/1‑75 (MQ=255) | AGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAA > REL606/3983627‑3983771 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |