Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 275,418 | 0 | A | C | 28.1% | 58.7 / 11.4 | 32 | F528V (TTT→GTT) | yagW | predicted receptor |
Reads supporting (aligned to +/- strand): ref base A (18/5); new base C (0/9); total (18/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.14e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTT > REL606/275343‑275490 | tGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCaaa > 2:900924/1‑76 (MQ=255) tagcagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCa > 2:213269/2‑76 (MQ=255) cagcagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCa < 1:1173795/76‑1 (MQ=255) agcagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAg < 1:1157432/76‑1 (MQ=255) cagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGcc < 1:791123/76‑1 (MQ=255) gggCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGt > 2:930877/1‑76 (MQ=255) ccTTAGTTAATGTTACGCCCCGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTat < 2:1129736/76‑1 (MQ=255) acatAGTAAGTGTTACGCAACGGCGCCGGAGCATGATTCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTat < 2:1129747/73‑1 (MQ=255) tAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTatcat > 2:455538/1‑76 (MQ=255) ttAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTa > 1:671392/1‑76 (MQ=255) tattGTTACGCCCCGTCGCCTGATCATGATTCTTTCTTGAAGCGTTGTTTTCGCCAACCCAGCCGTTATCATCTAc < 2:352300/73‑1 (MQ=255) tAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTAc > 1:26569/1‑76 (MQ=255) aaTGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTAcc > 1:161767/1‑76 (MQ=255) cGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCa < 2:552465/76‑1 (MQ=255) cGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCa < 1:824494/76‑1 (MQ=255) cACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCg > 1:740500/1‑76 (MQ=255) cACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCg > 2:15456/1‑76 (MQ=255) gTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGaaa > 2:164305/1‑76 (MQ=255) ccTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGc > 1:432393/1‑76 (MQ=255) ccTAAACATGAATCTCCCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGc < 1:998116/76‑1 (MQ=255) gAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGtt < 2:1206900/76‑1 (MQ=255) gAACATGAATCTCCCCGGAAGCGCTGACTTCCCAAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGtt < 1:963669/76‑1 (MQ=255) aTGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATcc > 1:948551/1‑76 (MQ=255) gAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCAt > 2:30450/1‑76 (MQ=255) aaTCTCCCCTGAAGCGCTGCCTTCCCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCAtt < 1:1059160/76‑1 (MQ=255) ctcCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAAt > 2:500640/1‑76 (MQ=255) cccTGAAGCGCTGCCTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATc < 1:1188562/75‑1 (MQ=255) ccTGAAGCGCTGCCTTCGCAAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCg < 1:1032158/76‑1 (MQ=255) cTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGa > 1:240193/1‑76 (MQ=255) tGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaa > 2:345853/1‑76 (MQ=255) aaGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaaaa > 2:446317/1‑76 (MQ=255) ggcgcTGTTTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAAt < 1:864095/75‑1 (MQ=255) gcTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTg > 2:564042/1‑76 (MQ=255) gcTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTg > 1:692020/1‑76 (MQ=255) cTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGa > 1:303808/1‑76 (MQ=255) ttCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACtgtg > 2:319028/1‑76 (MQ=255) cAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCtt < 2:1196640/76‑1 (MQ=255) cAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCtt < 2:1196647/76‑1 (MQ=255) cAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCtt > 1:993287/1‑76 (MQ=255) | TGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTT > REL606/275343‑275490 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |