Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,983,700 | 0 | T | G | 29.6% | 96.8 / 21.4 | 54 | T397P (ACC→CCC) | ECB_03713 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (4/34); new base G (16/0); total (20/34) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.30e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTGC > REL606/3983625‑3983775 | gAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGt < 1:1264894/76‑1 (MQ=255) gtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTc < 2:1817616/76‑1 (MQ=255) gtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTc < 2:714113/76‑1 (MQ=255) tAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCcccc > 2:1955445/1‑76 (MQ=255) aaaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTccc < 1:156477/76‑1 (MQ=255) ntGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTccccc > 1:448962/2‑74 (MQ=255) aTGCGGGGTTACCGCCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCcac < 1:1315620/76‑1 (MQ=255) aTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCcac < 2:1762290/76‑1 (MQ=255) aTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCcac < 2:1561903/76‑1 (MQ=255) gggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAAtt < 2:99221/76‑1 (MQ=255) ggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCCCCAAAtt > 1:1170604/1‑76 (MQ=255) gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTg < 1:1744955/76‑1 (MQ=255) gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCCCCCCCAAAattt > 1:523108/1‑75 (MQ=255) gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAAATTg > 1:84442/1‑76 (MQ=255) tACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGAc < 2:890150/76‑1 (MQ=255) tACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCTCCCCCAAAATTGaa > 2:1333816/1‑75 (MQ=255) accaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCCCCCCCAAAATTGACg > 2:1705314/1‑76 (MQ=255) ccaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGa < 2:1906686/76‑1 (MQ=255) caccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCTCCCCCAAAATTGACCAt > 1:1432355/1‑76 (MQ=255) ccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAAATTTACGCTCg > 2:1995185/1‑76 (MQ=255) cGAGTGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGc < 2:203048/76‑1 (MQ=255) gAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCt < 2:1072763/76‑1 (MQ=255) ggTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGga < 1:395741/76‑1 (MQ=255) tGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAAATTGACGATCGCTGggga > 2:1375281/1‑76 (MQ=255) tGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCCCCCCCAAAATTTACGATCGCTGggga > 1:1878045/1‑76 (MQ=255) gATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCCCCCACAAATTTTACGATCGCTGgagag > 1:1923191/1‑75 (MQ=255) aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATg < 1:1619996/76‑1 (MQ=255) aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCGCCCACAAACTTAAAAAGCGCTGGTGATg > 1:440101/1‑76 (MQ=255) tgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGAc < 2:123225/76‑1 (MQ=255) tgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCCAATTTTGCGATCGCTGGAGGTGAc > 2:575889/1‑76 (MQ=255) tgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAATTTTACGATCGCTGGAGATGAc > 2:2113489/1‑76 (MQ=255) ccngcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAATTTGACGATCGCTGGGGGTGACGc > 2:1499273/4‑76 (MQ=255) gcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCTCCCCCAAAATTTACGATCGCTGGGGGTGACGCCAt > 2:609585/1‑76 (MQ=255) gcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCACAAATTTTACGATCGCTGGGGATGACGCCac > 2:1886920/1‑75 (MQ=255) cATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAAATTTACGATCGCTGGGGGTGACGCCaac > 2:202192/1‑74 (MQ=255) aTCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATcg < 1:1999743/76‑1 (MQ=255) aTCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCACAAATTTGACGATCGCTGGGGATGACGCCCTcg > 2:1878772/1‑76 (MQ=255) tcctccGGTTTTTATTGGTGGCATCGGCG‑GGGGTGCCTCCCACAAA‑TTGACGATCGCTGGAGATGACGCCATcgcg < 1:1031086/76‑1 (MQ=255) cTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGa < 2:66248/76‑1 (MQ=255) tACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCCCCAAATTT‑AAGAACCCCGGAGAAGAAGCCCTCCCGaa > 1:1543162/1‑76 (MQ=255) tACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCTCCCCCAAAATTTACGATCGCCGGGGGTGAAGCCCTCGCGaa > 1:429181/1‑76 (MQ=255) aCGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAg < 2:1257159/76‑1 (MQ=255) aCGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCTCCCACAAAATTTACGATCGCTGGAGATGACGCCATCGCGAAg > 1:1074754/1‑76 (MQ=255) ggCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcc < 2:1253057/76‑1 (MQ=255) cTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAATTTTACGATCGCTGGAGATGACGCCATCGCGAAGccac > 1:2060204/1‑76 (MQ=255) gTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccaccc < 2:1897175/76‑1 (MQ=255) gTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCACAAATTTGACGAACGCTGGGGATGACGCCATCGCGAAGccaccc > 1:1472765/1‑76 (MQ=255) tATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccacccac < 2:1731618/76‑1 (MQ=255) aTTGGTGGCATCGGCA‑GGGGGGCCTCCCACAAATTTGACGATCGCTGGGGGTGGAGCCATCGCGAAGCCCCCcaca > 2:718595/1‑76 (MQ=255) tggtggCATCGGCA‑GGGGGGCCTCCCACAAATTTGACGATCGCTGGGGATGACGCCATCGCGGAGCCACCCACAcc > 2:543406/1‑76 (MQ=255) ggtggCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCt < 1:755461/76‑1 (MQ=255) gtggCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTg < 2:138169/76‑1 (MQ=255) cATCGGCA‑GGGGGGCCTCCCCCAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGt > 1:1858594/1‑76 (MQ=255) aTCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGtt < 1:1283155/76‑1 (MQ=255) aTCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGtt < 1:2010209/76‑1 (MQ=255) aTCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGtt < 1:1375281/76‑1 (MQ=255) aTCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGtt < 1:570962/76‑1 (MQ=255) tCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGttt < 2:440917/76‑1 (MQ=255) ggCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCa > 2:499165/1‑76 (MQ=255) gCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAg < 2:338457/76‑1 (MQ=255) gCA‑GGGGTGCCTCCCACAAATTCGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCCGCGGGTTCAg > 1:593203/1‑76 (MQ=255) a‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa > 1:1941075/1‑76 (MQ=255) a‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa > 1:741892/1‑76 (MQ=255) ggggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTaa < 1:620313/76‑1 (MQ=255) ggggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTaa < 2:124627/76‑1 (MQ=255) gggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAt < 1:1989565/76‑1 (MQ=255) gggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAt < 2:174362/76‑1 (MQ=255) tGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTgc > 1:574092/1‑76 (MQ=255) tGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTgc < 1:1363695/76‑1 (MQ=255) | GAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTGC > REL606/3983625‑3983775 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |