Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 275,418 | 0 | A | C | 31.1% | 91.3 / 10.8 | 45 | F528V (TTT→GTT) | yagW | predicted receptor |
Reads supporting (aligned to +/- strand): ref base A (28/3); new base C (0/14); total (28/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.07e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTAT > REL606/275344‑275492 | gTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑c > 1:115538/1‑76 (MQ=255) tCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑cc > 2:805162/1‑76 (MQ=255) agcagGGCCTTAGGTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAg > 2:927528/1‑76 (MQ=255) cagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGcc > 2:393156/1‑76 (MQ=255) cagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGcc < 2:1721343/76‑1 (MQ=255) cagGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGcc < 2:1073876/76‑1 (MQ=255) agGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCg > 2:812076/1‑76 (MQ=255) agGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCg < 2:2220538/76‑1 (MQ=255) ggCCTTAGTTATTGTTGCGCCACGTCGCCTGAACATGAATCTCCCCGGAAGCGCTGACTTCGCCAAC‑CCAGCCGtt < 2:146886/76‑1 (MQ=255) ggCCTTAGTTAATGTTCCCCCCCGTCGCCTGACCATGAATCTCTCCTGAAGCGCTGACTTCGCCAAC‑CCAGCCGtt < 2:887821/76‑1 (MQ=255) cccTTACTAAATGTTACGCCACGCCGCCTTAAAAAGAATCTCTCCTAAAGCCCTGCCTTCGCAAAA‑CCAGCCGTTa < 2:1433406/75‑1 (MQ=255) cTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTatc > 2:1641465/1‑76 (MQ=255) tAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTatcat > 1:940812/1‑76 (MQ=255) ttAATTTTACCCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAC‑CCAGCCGTTATCATCTa < 2:865769/76‑1 (MQ=255) aaTGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTAcc > 1:92522/1‑76 (MQ=255) aTGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCg > 2:1364368/1‑76 (MQ=255) tGTTACGCCACGTGGCCGGAACAAAAATCTCCCCAAAGGCGCTGACTTCGCCAAC‑CCAGCCGTTATCATCTACCGt < 2:551687/76‑1 (MQ=255) tGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGt > 1:1204099/1‑76 (MQ=255) ttACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTAc > 1:47753/1‑76 (MQ=255) ttACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTAc > 2:916927/1‑76 (MQ=255) tACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACg > 2:2176725/1‑76 (MQ=255) tACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACCTCGCCAAA‑CCAGCCATTATCATCTACCGCACg > 2:1805919/1‑76 (MQ=255) aCGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGc > 1:1389080/1‑76 (MQ=255) ccACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGc > 2:153211/1‑76 (MQ=255) gTCGCCTGAACATGAATCTCTCGTGAAGCGCTGACTTCGCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGaaa < 1:529600/76‑1 (MQ=255) gTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGaaa < 1:2095229/76‑1 (MQ=255) gTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGCGaaa > 2:218727/1‑76 (MQ=255) gTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGCGaaa > 2:846001/1‑76 (MQ=255) gTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGCGaaa > 2:1902140/1‑76 (MQ=255) tCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAAt < 2:1438463/76‑1 (MQ=255) gCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATgg > 1:1312705/1‑76 (MQ=255) tGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGt > 1:1006769/1‑76 (MQ=255) aaccATGAACCCCTCCTAAGGCGCTGACTTCGCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGtt < 2:1120643/73‑1 (MQ=255) aaaCATGACGCGCTCCGAAAGCCTTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGtt < 2:1120629/75‑1 (MQ=255) aaCATGAATCTCTCCTGAAGCGCTGCCTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTa < 1:602293/76‑1 (MQ=255) aaCATAAATCTCTCCTGAAGCGCTGACTTCCCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTa < 1:1960498/76‑1 (MQ=255) aCATGAATCTCCCCGGAAGCGCTGACTTCGCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTAt < 2:718960/76‑1 (MQ=255) tGACTCTCTCCGGAAGCGCTGCCTTCGCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCa < 1:1762985/76‑1 (MQ=255) tGAATCCCTCCTAAAGCGCTGACTTCGCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCa < 1:1652724/76‑1 (MQ=255) aaTCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCAtt > 1:999751/1‑76 (MQ=255) aTCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTg > 2:781116/1‑76 (MQ=255) tctcccCTGAAGCGCTGACTTCGCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTgg < 2:694066/76‑1 (MQ=255) ctctcCGGAAGCGCTGACTCCCCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGa < 1:915425/76‑1 (MQ=255) tctcCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGGACGCAGCGAAATGGCGCTATCCATTAGaa > 2:775533/1‑76 (MQ=255) tccccTGAAGCGCTGACTCCGCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGaa < 2:140948/73‑1 (MQ=255) ctcCTGAAGCGCTGACTCCGCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAAt < 2:1266833/76‑1 (MQ=255) ccccTGAAGCGCTGACTCCGCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGaan < 1:953557/74‑2 (MQ=255) ccTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCg > 1:209867/1‑76 (MQ=255) tGAGGCGCTGCCTCCGCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaa < 1:805162/76‑1 (MQ=255) tGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGaa > 1:1502896/1‑76 (MQ=255) ggtGACTTCGCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTg < 2:166235/74‑1 (MQ=255) gcTGACTTCGCCAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTg < 2:126708/76‑1 (MQ=255) cTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGa > 1:742493/1‑76 (MQ=255) gACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACt > 1:635634/1‑76 (MQ=255) cccccAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTgg < 1:225879/73‑1 (MQ=255) ccAAACCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACtgtgngc < 2:1588193/76‑4 (MQ=255) cAAC‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCtt < 2:2126884/76‑1 (MQ=255) aaa‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTa > 2:608347/1‑76 (MQ=255) aaa‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTa > 1:94846/1‑76 (MQ=255) aaa‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTa > 1:1074439/1‑76 (MQ=255) aa‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTAt > 2:1099123/1‑76 (MQ=255) | GTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAA‑CCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTAT > REL606/275344‑275492 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |