| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 275,418 | 0 | A | C | 37.0% | 79.2 / 13.2 | 46 | F528V (TTT→GTT) | yagW | predicted receptor |
| Reads supporting (aligned to +/- strand): ref base A (21/8); new base C (0/17); total (21/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.12e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTA > REL606/275346‑275491 | cagcagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCa > 1:393151/1‑76 (MQ=255) agcagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAg < 1:1671747/76‑1 (MQ=255) agcagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAg < 1:1868406/76‑1 (MQ=255) gcagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGc < 1:225977/76‑1 (MQ=255) agGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCg > 1:1153405/1‑76 (MQ=255) agGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCg > 1:942491/1‑76 (MQ=255) ggCCTTAGTTAATGTTACCCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGtt < 2:687653/76‑1 (MQ=255) aGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTatcatc > 1:1117375/1‑76 (MQ=255) ttAATGTTCCCCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTa < 2:704185/76‑1 (MQ=255) ttAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTa > 1:1616448/1‑76 (MQ=255) gTTCCGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTa < 1:196759/76‑1 (MQ=255) gTTACCCCACGTCCCCTAAACATGAATCTCTCCTGAAGCGCTGGCTTCCCCAACCCAGCCGTTATCATCTACCGTa < 1:1201104/76‑1 (MQ=255) cttACCCCCCGCCCCCTGAACAAGAACCTCTCCGGAAGCGCTGCCTTCGCCAACCCAGCCGTTATCATCTACCGTa < 2:353719/75‑1 (MQ=255) ttACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTAc < 1:605964/76‑1 (MQ=255) cccACGTCGCCTGAACATAAATCTCTCCTGAAGCGCTGTCTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAg < 1:1291890/75‑1 (MQ=255) cACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCg > 1:859758/1‑76 (MQ=255) ccGTCGCCTGAACATGAACCTCTCCTGAAGCGCTGCCTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGa < 1:1818420/75‑1 (MQ=255) aCGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGa > 1:1522358/1‑76 (MQ=255) cGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGaa > 2:11829/1‑76 (MQ=255) gTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGaaa < 1:128795/76‑1 (MQ=255) tCGCCTGAACATCAATCTCTCCTGAAGCTCTGCCTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAAt < 1:23174/76‑1 (MQ=255) ccTGAACATGAATCTCTCCTGAAGCGCTGACTTCCCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGc < 2:661898/76‑1 (MQ=255) gAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGtt < 2:1146618/76‑1 (MQ=255) gAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGtt > 2:1133591/1‑76 (MQ=255) aaaCAGGAATCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGtt < 2:868006/75‑1 (MQ=255) aaaTCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCAt < 1:1611222/75‑1 (MQ=255) aaTCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCAtt > 2:987636/1‑76 (MQ=255) ctctcCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGa < 1:1678063/76‑1 (MQ=255) ctcCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAAt < 1:888848/76‑1 (MQ=255) ctcCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAAt > 2:467971/1‑76 (MQ=255) ctcCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAAt > 2:1393032/1‑76 (MQ=255) ccTGAAGCGCTGCCTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCg < 1:786796/76‑1 (MQ=255) ccTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCg < 1:876655/76‑1 (MQ=255) tGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaa > 1:531631/1‑76 (MQ=255) aaGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaaaa < 1:1610950/76‑1 (MQ=255) gcTGACTCCCCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTg < 1:200364/76‑1 (MQ=255) cTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGa > 1:404701/1‑76 (MQ=255) cTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGa > 1:215167/1‑76 (MQ=255) tGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGAc > 1:449214/1‑76 (MQ=255) aCTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACtg > 2:345161/1‑76 (MQ=255) aCTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACtg > 2:626480/1‑76 (MQ=255) cTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACtgt < 2:456241/76‑1 (MQ=255) ttCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACtgtg < 2:1304286/76‑1 (MQ=255) ttCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACtgtg > 1:341405/1‑76 (MQ=255) tccGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACtgtg < 1:1740195/74‑1 (MQ=255) cGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGt > 1:901086/1‑76 (MQ=255) cAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCtt < 1:1497858/76‑1 (MQ=255) cAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCtt < 2:155882/76‑1 (MQ=255) cAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCtt > 2:1464078/1‑76 (MQ=255) aaaaCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCtt < 2:1739210/75‑1 (MQ=255) aaaCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTa < 2:1665079/76‑1 (MQ=255) aaaCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTa > 1:1719076/1‑76 (MQ=255) | CAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTA > REL606/275346‑275491 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |