Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,437,406 | 0 | G | T | 37.5% | 57.0 / 13.2 | 32 | D435E (GAC→GAA) | ypdD | fused predicted PTS enzymes: Hpr component/enzyme I component/enzyme IIA component |
Reads supporting (aligned to +/- strand): ref base G (12/8); new base T (12/0); total (24/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.35e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCTTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCAT > REL606/2437331‑2437480 | ctctTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTg < 1:2592190/76‑1 (MQ=255) ctctTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTg < 2:4054407/76‑1 (MQ=255) ctctTCCTGGGAAATATTCAGATAAGAAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTg < 2:3840491/76‑1 (MQ=255) ctctTCCGGGGAAATATACAGATAAGAAAGGTTTTTACCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTg < 2:3212536/76‑1 (MQ=255) ctctGCCTGGGAAATATTCAGATAAGGATTGCTTTTACCCCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTg < 2:4154997/76‑1 (MQ=255) tctcCCGGGGAAATATTCAAAAAAGGAATCCTTTTACCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGt < 2:1865655/76‑1 (MQ=255) tctTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGt < 2:5029997/76‑1 (MQ=255) ttCCGGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTcg < 2:610450/76‑1 (MQ=255) tCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTccc > 1:2430695/1‑74 (MQ=255) ccgggggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTcgcc < 2:2766632/73‑1 (MQ=255) ccTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTtgcc > 1:3677921/1‑76 (MQ=255) ccTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTcgcc > 2:5734035/1‑76 (MQ=255) ccTGGGGAATATTCAGATAAGAAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTcgcc < 1:5517207/76‑1 (MQ=255) ccTGGGAAATATTAAGATAAGGAATGTTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTcgcc < 2:5114814/76‑1 (MQ=255) cTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTtgccc > 2:2261315/1‑76 (MQ=255) cTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTtgccc > 2:3030250/1‑76 (MQ=255) cTGGGAAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTcgccc < 2:5231738/76‑1 (MQ=255) ggggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTtggccgc > 2:335967/1‑76 (MQ=255) ggggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTtgcccgc > 2:5527849/1‑76 (MQ=255) ggggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTcgtcccc > 2:4351280/1‑74 (MQ=255) gAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTcgcccccccc > 2:4990903/1‑74 (MQ=255) atatTCAGAAAAGGAATGTTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTCGCCCGCCGCCa < 1:619499/76‑1 (MQ=255) atTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTTGCCCcccgccctt > 1:2796238/1‑73 (MQ=255) atTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTCTCCCCCCCCCCGt > 1:4118046/1‑76 (MQ=255) atTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTCGCCCcccgcccgc > 2:316194/1‑73 (MQ=255) cAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTTCGCCGCCCCCAATAGc > 1:1854113/1‑76 (MQ=255) tAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTTCCCCGCCCCCAGTCGGA‑ccc > 1:1489671/1‑75 (MQ=255) tAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCCCCCGTCGCA‑Cct > 2:1166166/1‑76 (MQ=255) aGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTTTCCCCCCCCCCGCAGCA‑Cccgc > 1:1490542/1‑76 (MQ=255) gAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTTGCCCCCCCCCCGCAGCA‑Cccgct‑g > 2:504425/1‑76 (MQ=255) aaTGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCGCCCCGAGCA‑CCCGCC‑gg > 1:5006247/1‑76 (MQ=255) aaTGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCCCCCCCCATAGGC‑Cctgctcg > 2:929392/1‑74 (MQ=255) gCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTTGCCCCCCGGCCGTAGCA‑CCTGCT‑GGGA‑g > 1:3088393/1‑76 (MQ=255) gCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTGGCCCCCCGCCCGTAGGA‑CCCTCT‑GGGA‑g > 2:5922021/1‑76 (MQ=255) gCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTGGCCCCCCGCCCGGAGCA‑CCTGCC‑GGtgg > 2:1521111/1‑74 (MQ=255) ttttATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCCCCCGGAGCA‑CCCTCT‑GGTA‑GGc > 2:3212797/1‑76 (MQ=255) ttttATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCCCCAGTAGCA‑CCTGCT‑GGTA‑ggg > 2:1275720/1‑75 (MQ=255) ttttATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCCCCGCCAGTAGGA‑CCTGCC‑GGGA‑ggg > 1:4717643/1‑75 (MQ=255) ttttATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTCGCCCGCCGCCCGGAGCA‑CCCGCT‑GGTA‑ggg > 1:94796/1‑75 (MQ=255) tttATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCCCCCGTAGCC‑CCCGCC‑GGGA‑gggc > 2:2649086/1‑74 (MQ=255) tttATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCCCCCCCCGTAGCA‑CCTGCT‑GGGA‑GGCt > 1:2955623/1‑76 (MQ=255) ttATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTTGCCCGCCCCCCGAAGCA‑CCCGCC‑GGTT‑gggct > 1:3402534/1‑73 (MQ=255) ttATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCACCAGTAGCA‑CCTGCT‑GGGA‑gggct > 2:3567343/1‑73 (MQ=255) aTCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCCTCCGCCGCCCCTAGTA‑CCCGCT‑GGGA‑GGGTTcc > 2:6298946/1‑75 (MQ=255) tCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTTGCCCGCCGCCAGCAGCA‑CCCGCT‑GGGA‑GGCTTCaa > 2:1471682/1‑76 (MQ=255) ccgccgcCGATGCCCTTGGGCCGAAAG‑ATAATCGTGTTGTCGCCCCCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCaa < 1:775618/75‑1 (MQ=255) cgccgcCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCGCCCGTAGCA‑ACCGCTGGGTG‑GGCTTCaa > 1:3496667/1‑76 (MQ=255) gccgccGATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTTGCCCGCCCCCAAGAGGA‑CCTCCT‑GGGA‑GGCCTCtaaa > 2:5434048/1‑75 (MQ=255) cccgccgGTGCCCATCGTGGGAAAG‑ATAACCGGTTTGCCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAAt < 1:6197414/75‑1 (MQ=255) cccgcagATCTCCATCGTGGGGAAT‑ATAATGGGTTTGTCCCCCCCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAAt < 2:3578498/75‑1 (MQ=255) ccgccgATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCGCCCGTAGCA‑CCTCCT‑GGGAGGGTTTCAAAt > 2:59436/1‑76 (MQ=255) ccgccgATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGGTTCCAAtt > 1:3477966/1‑76 (MQ=255) ccgccgATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTCGCCCGCCGCCCGTAGCA‑CCTGCT‑GGTA‑GGGTTTaaaat > 1:1485715/1‑74 (MQ=255) cgccgTTGCCCATGGGGCGGAAA‑ATAACCGCTTTGCCCCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTg < 2:2298800/76‑1 (MQ=255) ccccgAGGTCCATCGGGGGAAAA‑AAAGCCGTCTTGCCCCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTg < 2:1873538/74‑1 (MQ=255) ccgATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGGA‑GGGTTCTAATTGCt > 1:5194302/1‑76 (MQ=255) ccgATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGGA‑GGGGTCCAAATGCt > 2:648001/1‑76 (MQ=255) ccgATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCGCCAGTAGCA‑CCCGCT‑GGTA‑GGCGTCTAAATGCt > 2:6170393/1‑76 (MQ=255) ccgATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCCCCGCCCCAAGCA‑CCTGCT‑GGTAGGGCTTCAATTTgt > 2:1402207/1‑75 (MQ=255) cgATGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCGCCCGTAACA‑CCCGCT‑GGGA‑GGGCTTAAAATGCTc > 1:281265/1‑76 (MQ=255) cgATGTCCATCGTGCGGAAG‑ATAAACGGCTTTTCGCCCGCCGCCCGCAGCA‑CCCGCT‑GGGA‑GGGTTCAAAATGCTc > 1:5073315/1‑76 (MQ=255) tGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCCCCCGTAGCA‑CCCGCT‑GGGA‑GGGTTCCAATTGCTCCtc > 2:5084956/1‑75 (MQ=255) tGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCCCCGCCCGTAGCA‑CCTGCT‑GGGA‑GGGTTTAAAATTCTTCTg > 2:5872162/1‑76 (MQ=255) tGTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCCCCGCCAGTAGCA‑CCTGCC‑GGTA‑GGGTTCAAATTTCTccct > 2:790484/1‑74 (MQ=255) gTCCATCGTGCGGAAG‑ATAATCGGCTTTTTGCCCCCCGCCCGGAGCA‑CCTGCT‑GGGA‑GGGGTCCAAATGCTCCTGc > 1:5301142/1‑76 (MQ=255) gTCCATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCCCCCCCAGGAGAACCCTCCG‑GGGG‑GGCTTTAAATTGTTCCtc > 1:1670623/1‑75 (MQ=255) gTCCATCGTGCGGAAG‑ATAATCGGCTTGTCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGGA‑GGGTTCCAAATTTTCCTGc > 1:3382868/1‑76 (MQ=255) ccATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCCCCAGGAGCA‑CCTGCT‑GGGA‑GGGGTTAAATTTTTCCTGCTc > 1:5022510/1‑76 (MQ=255) ccATCGGGCGAAAG‑ATATCGGGTTTGCCCCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTc < 1:3273625/76‑1 (MQ=255) cATGGTGGGGAAG‑ATAATGGGCTTGCCCCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCg < 1:2362822/76‑1 (MQ=255) cATCGTGCGGAAG‑ATAATCGGCTTTTCGCCCGCCGCCCGTAGCA‑CCCGCC‑GGGA‑GGCTTCAAAATGCTCCTCCTCg > 1:4196938/1‑76 (MQ=255) aTCGTGCGGAAG‑ATAATCGGCTTGTCGCCCGCCCCCCGTAGCA‑CCTGCT‑GGTC‑GGGTTTAAATTGCTTCCGCTCGt > 1:6336309/1‑76 (MQ=255) cGTGCGGAAG‑ATAATCGGCTTTTCGCCCCCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTTCTCCTGCtcctcc > 1:3721336/1‑76 (MQ=255) gcgCGGAAA‑ATAACCGGTTTGCCGCCCCCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCGTCCg < 1:825725/74‑1 (MQ=255) gTGCGGCAGAATAATCGGTTGGCCCCCCCCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCGTcc < 2:2585918/76‑1 (MQ=255) ggAAG‑ATAATCGTCGTCGCCCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCGTCCGgcgc < 2:201846/76‑1 (MQ=255) tAATCGTCTTGCCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTc < 1:775049/76‑1 (MQ=255) tCGTCTTGCCCCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTcgcg < 2:861028/76‑1 (MQ=255) cGGCTTGTCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGGTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGa > 1:5438521/1‑76 (MQ=255) gCTTGTCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCGTCCGGCCCGCTGTCGCGATc > 1:3024428/1‑76 (MQ=255) cTTGTCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATcc > 1:1040839/1‑76 (MQ=255) cTTGTCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCCATcc > 2:3239771/1‑76 (MQ=255) ttGTCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCa > 1:2140371/1‑76 (MQ=255) tGTCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCat > 2:5683530/1‑76 (MQ=255) tGTCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCat > 2:3097049/1‑76 (MQ=255) tGTCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCat > 2:2928453/1‑76 (MQ=255) tGTCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCat > 2:4706300/1‑76 (MQ=255) | CTCTTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAG‑ATAATCGGCTTGTCGCCCGCCGCCAGTAGCA‑CCTGCT‑GGTA‑GGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCAT > REL606/2437331‑2437480 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |