| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 275,418 | 0 | A | C | 28.6% | 67.7 / 11.1 | 36 | F528V (TTT→GTT) | yagW | predicted receptor |
| Reads supporting (aligned to +/- strand): ref base A (20/5); new base C (0/10); total (20/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.64e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTA > REL606/275343‑275491 | tGTCAGCAGGGCCTTAGTAAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCaaa < 2:1496729/76‑1 (MQ=255) tCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAAcc < 1:345683/76‑1 (MQ=255) tCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAAcc > 1:205460/1‑76 (MQ=255) cagcagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAGCCa < 1:731766/76‑1 (MQ=255) agcagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAg < 2:731734/76‑1 (MQ=255) agGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCg > 2:55651/1‑76 (MQ=255) gcccTTAGTTAATGTTCCGCCACGTCGCCTGAACATAAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGtt < 2:803437/74‑1 (MQ=255) cTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTatc > 1:1636676/1‑76 (MQ=255) cTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTatc > 2:701890/1‑76 (MQ=255) ttAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTatca > 1:1913807/1‑76 (MQ=255) tAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTAc > 2:958470/1‑76 (MQ=255) aaTGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTAcc < 2:1805448/76‑1 (MQ=255) ttccGCCACGTCCCCTAAACAGGAATCTCCCCTGACGCGCTGACTTCCCCAACCCAGCCGTTATCATCTACCGTAc < 2:1240622/73‑1 (MQ=255) cgccgcgccgCCTAAACATAATCCTCTCATGACGCGCTGCCTTCGCCAAACCAGCCGTTATCATCTACCGTACGCa < 1:660516/76‑1 (MQ=255) gCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAg > 2:58115/1‑76 (MQ=255) ccACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGc > 2:1442181/1‑76 (MQ=255) ccGTCGCCTAAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGa < 1:1047338/75‑1 (MQ=255) cGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGan < 1:1654563/76‑2 (MQ=255) gCCTAAACATAAATCTCTCCTGAAGCGCTGCCTCCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATgg < 1:416487/76‑1 (MQ=255) ccTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGc > 2:1517810/1‑76 (MQ=255) cTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCg > 2:1561749/1‑76 (MQ=255) cTGAAAATGAATCTCTCCTGAAGCGCTGCCTCCCCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCg < 1:694230/76‑1 (MQ=255) gAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGtt > 1:1232521/1‑76 (MQ=255) gAACATGAATCTCCCCTGAAGCGCACACTCCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGtt < 1:338211/76‑1 (MQ=255) aTGAATCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATcc < 1:267393/76‑1 (MQ=255) gAATCTCTCCTGAGGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCAt < 2:179653/76‑1 (MQ=255) gAATCTCTCCGGAAGCGCTGCCTCCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCAt < 1:1379258/76‑1 (MQ=255) gAATCTCCCCGGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCAt < 2:526395/76‑1 (MQ=255) aaTCTCTCCTGAAGCGCTGCCTTCCCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCAtt < 1:2020881/76‑1 (MQ=255) aTCTCTCCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTg < 1:1824078/76‑1 (MQ=255) ctctcCTGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGa < 1:1361282/76‑1 (MQ=255) ctcCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAAt > 2:1128927/1‑76 (MQ=255) ctcCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAAt > 1:659486/1‑76 (MQ=255) ccccTGAAGCGCTGCCTTCCCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAAt < 2:1277149/74‑1 (MQ=255) ccTGAAGCGCTGCCTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCg < 1:1399576/76‑1 (MQ=255) tGAAGCGCTGACTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaa < 2:2003878/76‑1 (MQ=255) tGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaa > 1:1043818/1‑76 (MQ=255) tGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaa > 2:204504/1‑76 (MQ=255) tGAAGCGCTGACTTCGCAAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaa < 1:237011/76‑1 (MQ=255) gAAGCGCTGCCTCCCCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaaa < 1:278055/76‑1 (MQ=255) gAAGCGCGGACTCCCCCACCCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaaa < 2:1028744/76‑1 (MQ=255) aaGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaaaa > 2:320984/1‑76 (MQ=255) aGCGCTGCCTTCCCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAAt < 1:1152671/76‑1 (MQ=255) tcaCTTCGCCAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGAc < 2:1884682/74‑1 (MQ=255) tGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGAc > 2:432802/1‑76 (MQ=255) tGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGAc > 1:872732/1‑76 (MQ=255) gACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACt > 2:172292/1‑76 (MQ=255) gACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACt > 1:77971/1‑76 (MQ=255) cAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCtt < 1:1183959/76‑1 (MQ=255) aaaCCAGCCGTTATCATCTACCGTATGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTa > 2:137124/1‑76 (MQ=255) aaaCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTa > 1:1845114/1‑76 (MQ=255) | TGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTA > REL606/275343‑275491 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |