Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 417,321 | 0 | T | G | 31.0% | 39.5 / 12.0 | 29 | T27P (ACA→CCA) | ECB_00378 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (8/12); new base G (9/0); total (17/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.29e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCGTCGGGACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTCATACATAATACGGCTAT > REL606/417250‑417396 | ggCGTCGGGACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAATGGCTGgtttt > 1:284746/1‑76 (MQ=255) tCGGGACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTTCAg > 1:198283/1‑76 (MQ=255) cGGGACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAgg < 2:220293/76‑1 (MQ=255) ggACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTTCAGGCa > 1:518556/1‑76 (MQ=255) ggACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTCCAGGct > 1:74862/1‑75 (MQ=255) ggACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTCCAGGcg > 1:518545/1‑75 (MQ=255) ggACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTCCAGGCa > 1:74861/1‑76 (MQ=255) aCTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTCCCGGCAtt > 2:602986/1‑76 (MQ=255) cTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTTTCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTa < 2:125524/76‑1 (MQ=255) aGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTTCAGGCATTATg > 1:414991/1‑76 (MQ=255) ttttGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTCCAGGCATTATGAc > 1:641081/1‑76 (MQ=255) ttGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACct > 1:491411/1‑76 (MQ=255) aCCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGGTTTTTCAGGCATTATGGCCTCTTttatta > 1:371336/1‑76 (MQ=255) cgcATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACtga > 1:80599/1‑76 (MQ=255) cgcATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACtga < 2:453820/76‑1 (MQ=255) gcATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACtgaa < 2:116566/76‑1 (MQ=255) aTTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACtgaatg < 2:387773/76‑1 (MQ=255) tttCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGaaa < 2:203419/76‑1 (MQ=255) gCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTa < 2:750901/76‑1 (MQ=255) ttttATCGACATCGGGTAGAGTGGCTGGTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAggg > 1:628775/1‑76 (MQ=255) tttATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGa < 1:136748/76‑1 (MQ=255) tCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTca < 2:692039/76‑2 (MQ=255) cGACATCGGGTAGAGTGGCTGGTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGa > 1:499516/1‑76 (MQ=255) aTCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGGCGAtctc > 2:354187/1‑76 (MQ=255) cGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTcat > 2:82180/1‑76 (MQ=255) ggTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTcatac > 1:685767/1‑76 (MQ=255) tAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTcatacat > 2:558156/1‑76 (MQ=255) tAGAGTGGCTGGTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTGGGGAGTCGATCTcatacat > 2:517252/1‑76 (MQ=255) agagTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTcatacata < 2:54423/76‑1 (MQ=255) ggCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTCATACATAATACg < 1:552364/76‑1 (MQ=255) tGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTCATACATAATACGGCt > 2:52479/1‑76 (MQ=255) tGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTCATACATAATACGGCt > 2:640745/1‑76 (MQ=255) gTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTCATACATAATACGGCTa < 2:587680/76‑1 (MQ=255) gTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTCATACATAATACGGCTa > 1:651871/1‑76 (MQ=255) tttttCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTCATACATAATACGGTTAt < 1:67106/76‑1 (MQ=255) | GGCGTCGGGACTAGTTTTGAGTTCGACCCGCGCATTTTTCAGCGTTTTATCGACATCGGGTAGAGTGGCTGTTTTTCCAGGCATTATGACCTCTTTGATTAACTGAATGAAATCGTAGGGAGTCGATCTCATACATAATACGGCTAT > REL606/417250‑417396 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |