Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,983,700 | 0 | T | G | 47.4% | 9.3 / 12.3 | 19 | T397P (ACC→CCC) | ECB_03713 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (0/10); new base G (9/0); total (9/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.08e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTGC > REL606/3983625‑3983775 | gAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGt < 1:990365/76‑1 (MQ=255) cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCTCCCCCaaa > 2:866175/1‑76 (MQ=255) cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCCCCCCCaaa > 1:1913472/1‑76 (MQ=255) ggggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCGCCCGCaaaa > 2:950984/1‑75 (MQ=255) ggggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCCCCCCCaaaa > 1:590553/1‑75 (MQ=255) gggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAAtt < 1:5177/76‑1 (MQ=255) ggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAAttt < 2:849893/76‑1 (MQ=255) gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTg < 2:525963/76‑1 (MQ=255) gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCCCCCCCAAAtttt > 1:143549/1‑75 (MQ=255) gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCCCCCCCAAAATTg > 1:1805598/1‑76 (MQ=255) tACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGGCCCCCCCAAAttttaa > 1:1107950/1‑73 (MQ=255) tGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTTCCCCCAAATTT‑TCCATCGCTTgaga > 2:1497495/1‑76 (MQ=255) tGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGgaga < 1:285760/76‑1 (MQ=255) gATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAAATTGACGATCGCTGGGGAt > 1:1578069/1‑76 (MQ=255) ttgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGa < 2:1243470/76‑1 (MQ=255) gtgcATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGcc < 2:479993/75‑1 (MQ=255) tACGGTTGATTTTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGaa < 1:210123/76‑1 (MQ=255) tACGGCTGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAATTTCGAGATCGCCGGAGGTGAGGCCGACGCGaa > 1:76828/1‑76 (MQ=255) gCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcca < 1:364873/76‑1 (MQ=255) tGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCCCAAATTTGACGATCGCTGGGGATGACGCCACCGCGGAGccacc > 1:166002/1‑76 (MQ=255) tGTTATTGGTGGCATCGGCA‑GGGGGGCCTCCCACAAAATTTACGAACCCCCGAGAAGACGCCATCGCGACGccacc > 2:1560272/1‑76 (MQ=255) ttGGTGGCATCGGCA‑GGGGGGCCTCCCACAAAATTGACGATCGCTGGAGATGACGCCATCGCCAAGccacccaccc > 2:1899861/1‑74 (MQ=255) ggtggCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCt < 1:232425/76‑1 (MQ=255) cGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTc < 1:690445/76‑1 (MQ=255) tGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTgc > 2:1344732/1‑76 (MQ=255) tGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTgc > 1:1475809/1‑76 (MQ=255) | GAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCA‑GGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTGC > REL606/3983625‑3983775 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |