Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 275,418 | 0 | A | C | 37.8% | 42.5 / 19.5 | 37 | F528V (TTT→GTT) | yagW | predicted receptor |
Reads supporting (aligned to +/- strand): ref base A (22/1); new base C (0/14); total (22/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.46e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTC > REL606/275343‑275488 | tGTCAGCAGGGCCTTAGTTAATGTTACGCCCCGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTAGCCaaa > 1:218527/1‑76 (MQ=255) tGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCaaa > 1:146633/1‑76 (MQ=255) tGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCaaa > 1:218559/1‑76 (MQ=255) cagcagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCa < 1:832895/76‑1 (MQ=255) cagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGcc > 2:109972/1‑76 (MQ=255) agGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCg > 2:104886/1‑76 (MQ=255) ccTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTat > 2:434799/1‑76 (MQ=255) aGTTAATGTTACGCCCCGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTatcatc < 2:12127/76‑1 (MQ=255) aGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTatcatc > 1:800358/1‑76 (MQ=255) aaTGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTAcc > 1:415552/1‑76 (MQ=255) tGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGt > 2:595865/1‑76 (MQ=255) ttACGCCCCGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTAc < 1:560822/76‑1 (MQ=255) ttACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTAc > 2:267712/1‑76 (MQ=255) ttACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTAc > 2:1311662/1‑76 (MQ=255) cGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCa < 1:217275/76‑1 (MQ=255) gCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAg > 1:770286/1‑76 (MQ=255) tCGCCTGAACATGAACCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAAt < 1:1269237/76‑1 (MQ=255) ccccTGAAAATAAATCTCCCCTGAAGCGCTGCCTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATg < 1:875426/74‑1 (MQ=255) cGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATg < 1:667169/76‑1 (MQ=255) cGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATg > 1:82555/1‑76 (MQ=255) ccTGAACATGAATCTCCCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGc < 2:335898/76‑1 (MQ=255) gAACATGAATCTCCCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGtt < 1:614717/76‑1 (MQ=255) aaCATGAATCCCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTa < 2:730587/76‑1 (MQ=255) cATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATc > 2:293249/1‑76 (MQ=255) aaTCTCTCCTGAAGCGCTGCTTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCAtt < 2:132028/76‑1 (MQ=255) aaTCTCCCCTGAAGCGCTGCCTCCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCAtt < 1:674176/76‑1 (MQ=255) aTCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTg < 2:518207/76‑1 (MQ=255) tctctcCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTgg > 1:769433/1‑76 (MQ=255) tctcccCGGAAGCGCTGACTCCGCAAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTgg < 1:770281/76‑1 (MQ=255) tctcCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGaa > 2:796202/1‑76 (MQ=255) ctcCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAAt > 2:1055871/1‑76 (MQ=255) ccTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCg > 2:1018056/1‑76 (MQ=255) ccTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCg > 2:1237766/1‑76 (MQ=255) ccTAAAGCGCTCACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCg < 1:51342/76‑1 (MQ=255) cTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGa > 2:645830/1‑76 (MQ=255) gAAGCGCTGCCTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaaa < 1:285892/76‑1 (MQ=255) aaGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATACATTGGAATCGaaaa > 1:1230056/1‑76 (MQ=255) cgcTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAAttt < 2:316872/76‑1 (MQ=255) gcTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTg > 2:41582/1‑76 (MQ=255) gACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACt > 2:918405/1‑76 (MQ=255) ttCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACtgtg < 1:785999/76‑1 (MQ=255) ttCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTAACCATTGGAATCGAAAATTTGACtgtg > 1:561039/1‑76 (MQ=255) cGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAGATGGCGTTATCCATTGGAATCGAAAATTTGACTGTggn < 1:37263/76‑2 (MQ=255) gCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTc > 2:947585/1‑76 (MQ=255) | TGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTC > REL606/275343‑275488 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |