Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,983,700 | 0 | T | G | 35.0% | 66.6 / 18.5 | 41 | T397P (ACC→CCC) | ECB_03713 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (6/20); new base G (14/0); total (21/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.34e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTG > REL606/3983625‑3983774 | gAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGt < 2:788070/76‑1 (MQ=255) gAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGAAGGGGt < 1:1150151/76‑1 (MQ=255) aaGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTg < 2:934633/76‑1 (MQ=255) aGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGc < 1:684639/76‑1 (MQ=255) aGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGc < 1:1746855/76‑1 (MQ=255) gtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGccgc > 1:1447410/1‑74 (MQ=255) aaaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑ccc > 2:1335884/1‑76 (MQ=255) gCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCCCaa > 2:1452120/1‑76 (MQ=255) ctgggTTTCCTCCGAGGGTTTTTGCTTGCATCCTTCGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACaaa < 2:760256/74‑1 (MQ=255) cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCCCaaa > 2:748444/1‑76 (MQ=255) ggggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAAt < 1:1131582/76‑1 (MQ=255) ggggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCC‑CCCCCaaaa > 1:836657/1‑75 (MQ=255) ttACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCCCAAATTTGa > 2:1622375/1‑76 (MQ=255) ttACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGAATCGGCAGGGGGGCCC‑CCCCCAAAattta > 2:1624362/1‑74 (MQ=255) accaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACg < 1:179568/76‑1 (MQ=255) accaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCCAAATTTAAc > 1:1101712/1‑76 (MQ=255) ccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCC‑CCCCCAAATTTGACGATCg > 2:1040296/1‑76 (MQ=255) aGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCC‑CCCACAAAATTTACGATCGccg > 2:1371986/1‑74 (MQ=255) gggTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTgg < 2:407292/76‑1 (MQ=255) gTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGgag < 2:529733/76‑1 (MQ=255) gTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCCCAAAATTGACGATCGCTGgag > 1:1527345/1‑76 (MQ=255) aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATg < 1:1308710/76‑1 (MQ=255) ttgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGa < 1:1601807/76‑1 (MQ=255) cttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCCCAAAATTTACGATCGCCGGGGGTGACGc > 2:975377/1‑76 (MQ=255) ttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCCCAAAATTTACGATCCCCGGGGGTGACggc > 2:1073394/1‑74 (MQ=255) cATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATc < 1:1486221/76‑1 (MQ=255) aTCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATcg < 2:871826/76‑1 (MQ=255) tCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATcgc < 1:1666029/76‑1 (MQ=255) ccTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATcgcg < 1:519442/76‑1 (MQ=255) aCGGCTGTTATTGGTGGCATCGGCAGGGGGGGCT‑CCCACAAAATTTACGATCGCCGGGGGAGACGCCATCGCGAAg > 2:976114/1‑76 (MQ=255) aCGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCCCCAAATTTAAGAACGCCTGGGATGGAGGCATCGCGGAg > 1:1671276/1‑76 (MQ=255) cGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGc < 1:1344169/76‑1 (MQ=255) cGGCTGTTATTGGTGGCATCGGCAGGGGGGCCT‑CCCCCAAATTTGACCATCGCTGGAGACGACGCCCACGCGCAGc > 1:485654/1‑76 (MQ=255) ggCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcc < 2:569952/76‑1 (MQ=255) ggCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcc < 2:574233/76‑1 (MQ=255) tATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccacccac < 1:946257/76‑1 (MQ=255) gtggCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTg < 2:1367700/76‑1 (MQ=255) ggCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCg > 1:978857/1‑76 (MQ=255) ggCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCg > 1:307801/1‑76 (MQ=255) cAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGt > 2:1367855/1‑76 (MQ=255) aGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa < 1:357650/76‑1 (MQ=255) ggggTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTaa > 1:1048303/1‑76 (MQ=255) gggCGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAt > 1:8812/1‑76 (MQ=255) gTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTg > 1:271082/1‑76 (MQ=255) gTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTg > 1:1384642/1‑76 (MQ=255) | GAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCT‑CCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTG > REL606/3983625‑3983774 |
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: ‑ |