| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 2,103,887 | 2 | . | A | 37.5% | 9.1 / 27.8 | 16 | coding (154/216 nt) | ECB_01992 | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base . (10/0); new base A (6/0); total (16/0) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AAACATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑CAGCCAGCCAGCCAGCCAGCCAGCCAGCC‑‑‑‑AGGTGTGAATGTGGATTTTGTACTACGAAAT > REL606/2103813‑2103947 | aaaCATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTca‑‑‑‑c < 2:3381042/76‑1 (MQ=255)aaaCATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTca‑‑‑‑c < 2:2718599/76‑1 (MQ=255) aCATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cag < 1:6917034/76‑1 (MQ=255) aCATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cag < 1:5022721/76‑1 (MQ=255) cATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagc > 1:3802172/1‑76 (MQ=255) aTGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagcc > 1:1056300/1‑76 (MQ=255) aTGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagcc < 2:3122877/76‑1 (MQ=255) aTGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagcc < 2:1101023/76‑1 (MQ=255) aTGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagcc < 1:4039942/76‑1 (MQ=255) tGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagcca < 1:617958/76‑1 (MQ=255) tGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagcca < 2:1400059/76‑1 (MQ=255) ggCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccag < 2:917750/76‑1 (MQ=255) ggCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccag < 1:1290092/76‑1 (MQ=255) ggCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccag > 1:3063613/1‑76 (MQ=255) ggCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccag < 2:2936289/76‑1 (MQ=255) gCTATACAATGAAAACTCCTTCGTTATTGACTTCGAATATAAAGGGGAGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagcgagc > 2:6173700/1‑76 (MQ=255) cTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagcc < 2:6606421/76‑1 (MQ=255) cTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagcc < 1:3759868/76‑1 (MQ=255) tataCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagcca < 1:2491446/76‑1 (MQ=255) tataCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagcca < 1:445194/76‑1 (MQ=255) tataCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagcca > 2:4780293/1‑76 (MQ=255) tataCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑aagccagcca > 1:3978548/1‑76 (MQ=255) ataCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccag > 1:238628/1‑76 (MQ=255) ataCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccag > 1:5076532/1‑76 (MQ=255) ataCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccag < 2:2840041/76‑1 (MQ=255) aCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagcc < 1:4189634/76‑1 (MQ=255) cAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagcca > 1:4254286/1‑76 (MQ=255) cAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagcca > 1:4696485/1‑76 (MQ=255) cAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagcca > 2:939647/1‑76 (MQ=255) cAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagcca < 1:5580590/76‑1 (MQ=255) aaTGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccag > 1:2926125/1‑76 (MQ=255) aaTGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccag < 1:5282305/76‑1 (MQ=255) aaTGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccag > 1:5408036/1‑76 (MQ=255) aaTGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccag > 2:5715158/1‑76 (MQ=255) tGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccagcc > 1:2143203/1‑76 (MQ=255) gAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccagcca > 2:6701680/1‑76 (MQ=255) gAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccagcca > 1:5172296/1‑76 (MQ=255) gAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccagcca < 1:3264666/76‑1 (MQ=255) aaaaCTGTTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGTTTTCA‑‑‑‑cagccagccagccagccag < 2:423213/76‑1 (MQ=255) aaaaCTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccagccag > 1:3053637/1‑76 (MQ=255) aaaaCTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccagccag > 2:1451935/1‑76 (MQ=255) aaaaCTGCTCCGTTATTGACTTCGAATATAAAGGGGTAAAA‑TTTGCA‑TTAGATTCCA‑‑‑‑cggccagccacccagccag < 2:2828223/76‑1 (MQ=255) aaaCTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccagccagc > 2:4765320/1‑76 (MQ=255) accTGCTTGTTAATTGACTTCGAATATAAAGGGGTGAAA‑TTTGAT‑TAAGTTTTCA‑‑‑‑cagccagccagccagccagcc < 2:1907798/74‑1 (MQ=255) aaCTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccagccagcc > 1:4260983/1‑76 (MQ=255) aCTGTTTCGTTATTGACTTCGAATA‑AAAGGGGGGAAATTTTGCA‑TTAGTTTTCA‑‑‑‑cagccagccagccagccagcca < 2:364799/76‑1 (MQ=255) cTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccagccagccag > 1:333838/1‑76 (MQ=255) tGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccagccagccagc > 1:1623954/1‑76 (MQ=255) gCTTGGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTGGCA‑TTAAATTTCA‑‑‑‑cagccagccagccagccagccagcc < 1:1348915/76‑1 (MQ=255) cttcttTATTGCTTTCAAATAAAAGGGGGTAAAT‑TTTGCA‑TTATTTTCCA‑‑‑‑cacccagccacccagccagccagcca < 2:2996974/76‑1 (MQ=255) cTTCGTTATTGACTTCGAATAAAAGGGGGTAAAA‑TTTGCA‑TTAGTTTTCA‑‑‑‑cagccagccagccagccagccagcca < 1:5670972/76‑1 (MQ=255) ttCTTTATTGACTTCGAATATAAAGGGGTAAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccagccagccagccag < 2:3215311/76‑1 (MQ=255) ttCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccagccagccagccag > 2:2043067/1‑76 (MQ=255) ttCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccagccagccagccag > 1:3566386/1‑76 (MQ=255) tCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑cagccagccagccagccagccagccagc > 2:2574522/1‑76 (MQ=255) cGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑Cagccagccagccagccagccagccagcc > 2:5540316/1‑76 (MQ=255) cGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑Cagccagccagccagccagccagccagcc > 1:3841428/1‑76 (MQ=255) cGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑Cagccagccagccagccagccagccagcc > 2:5519026/1‑76 (MQ=255) gTTATTGATTTCAAATAAAAAGGGGTAAAT‑TTTGCA‑TAAGATTTCC‑‑‑‑CAgccagccagccagccagccagccagcc‑‑‑‑a < 1:1129191/76‑1 (MQ=25) gTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑CAgccagccagccagccagccagccagcc‑‑‑‑a > 2:2931957/1‑76 (MQ=255) gTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑CAgccagccagccagccagccagccagcc‑‑‑‑a > 1:5168522/1‑76 (MQ=255) ttATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑CAGccagccagccagccagccagccagcc‑‑‑‑ag > 1:3327344/1‑76 (MQ=255) ttATTGACTTCGAATAAAAAGGGGGTAAA‑TTTGGATTTATTTTCCC‑‑‑‑CAgcagcccagccagccagccagccagcc‑‑‑‑a < 2:2787093/76‑1 (MQ=12) tATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑CAGccagccagccagccagccaggcagcc‑‑‑‑agc > 2:5509259/1‑75 (MQ=255) ttGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑CAGccagccagccagccagccagccagcc‑‑‑‑agcca > 1:2618109/1‑73 (MQ=255) ttGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑CAGccagccagccagccagccagccagcc‑‑‑‑agcca > 1:2537690/1‑73 (MQ=255) aGGGGTGAAA‑TTTGCA‑TTAGATTTCACAGCCAGCCAGCCAGCCAGCCAGCCAGCCAGCC‑‑‑‑AGGTGTGAATGTGGAtt > 2:5687850/1‑76 (MQ=255) ggggTGAAA‑TTTGCA‑TTAGATTTCACAGCCAGCCAGCCAGCCAGCCAGCCAGCCAGCC‑‑‑‑AGGTGTGAATGTGGAttt > 1:5852431/1‑76 (MQ=255) gggTGAAA‑TTTGCA‑TTAGATTTCACAGCCAGCCAGCCAGCCAGCCAGCCAGCCAGCC‑‑‑‑AGGTGTGAATGTGGAtttt > 2:2439951/1‑76 (MQ=255) gggTGAAA‑TTTGCA‑TTAGATTTCACAGCCAGCCAGCCAGCCAGCCAGCCAGCCAGCC‑‑‑‑AGGTGTGAATGTGGAtttt > 1:6061226/1‑76 (MQ=255) gggTGAAA‑TTTGCA‑TTAGATTTCACAGCCAGCCAGCCAGCCAGCCAGCCAGCCAGCC‑‑‑‑AGGTGTGAATGTGGAtttt > 1:1606351/1‑76 (MQ=255) tGCA‑TTAGATTTCACAGCCAGCCAGCCAGCCAGCCAGCCAGCCAGCC‑‑‑‑AGGTGTGAATGTGGATTTTGTACTACGaa > 1:844162/1‑76 (MQ=37) cA‑TTAGATTTCA‑‑‑‑CAGCCAGCCAGCCAGCCAGCCAGCCAGCCAGCAAGGTGTGAATGTGGATTTTGTACTACGAaat > 1:634727/1‑76 (MQ=37) | AAACATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAA‑TTTGCA‑TTAGATTTCA‑‑‑‑CAGCCAGCCAGCCAGCCAGCCAGCCAGCC‑‑‑‑AGGTGTGAATGTGGATTTTGTACTACGAAAT > REL606/2103813‑2103947 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |