Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,437,406 | 0 | G | T | 34.2% | 68.4 / 20.6 | 38 | 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/13); new base T (13/0); total (25/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.09e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCTTCCTGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCAT > REL606/2437331‑2437480 | ctctTCCTGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTg < 2:7049215/76‑1 (MQ=255) ctctTCCTGGGGAATATTCAGAAAAGG‑ATGCTTTTTACCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTg < 2:4104498/76‑1 (MQ=255) ctctTCCGGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTg < 2:7158626/76‑1 (MQ=255) ttttCCTGGGAAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGt < 2:5766288/74‑1 (MQ=255) tctTCCTGGGAAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGt < 2:5837804/76‑1 (MQ=255) ttCCTGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTct > 2:7413165/1‑75 (MQ=255) ttCCTGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTcg > 2:4848072/1‑76 (MQ=255) ttCCTGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTcc > 1:6336386/1‑75 (MQ=255) tCCTGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTccc > 1:1218339/1‑74 (MQ=255) ccTGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgcc > 1:5942346/1‑76 (MQ=255) ccTGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTcgcc > 2:1233640/1‑76 (MQ=255) ccTGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTcgcc > 1:73216/1‑76 (MQ=255) ccTGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgcc < 1:3463507/76‑1 (MQ=255) ccTGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgcc > 2:6499019/1‑76 (MQ=255) cTGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgccc > 1:1251232/1‑76 (MQ=255) cTGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTctccc > 1:3156513/1‑76 (MQ=255) cTGGGAAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgccc < 2:957649/76‑1 (MQ=255) tGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgcccc > 2:7361941/1‑75 (MQ=255) tGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgcccg < 1:6140045/76‑1 (MQ=255) tGGGGAATATACAAATAAGGATTGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgcccg < 2:437570/76‑1 (MQ=255) tGGGAAATATTCAGATAAGGATTGC‑TTTTAGCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgcccg < 2:4622153/76‑1 (MQ=255) tGGGAAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgcccg < 1:3126350/76‑1 (MQ=255) tGGAGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgcccg < 2:4237149/76‑1 (MQ=255) ggggAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTccccccc > 1:2239197/1‑74 (MQ=255) gggAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgcccgcc > 1:4464855/1‑76 (MQ=255) aaTATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCcccccccc > 1:5345267/1‑76 (MQ=255) aaTATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTGTGCCCgccgcc > 2:6401660/1‑76 (MQ=255) atatTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCa > 2:5972033/1‑76 (MQ=255) atatTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTTGCCCgcccccc > 2:3994088/1‑75 (MQ=255) atatTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCa < 1:6867997/76‑1 (MQ=255) atTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCGCCCGt > 1:4970298/1‑76 (MQ=255) atTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCGCCAGt > 1:730859/1‑76 (MQ=255) atTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCGGCCCCCAGt > 1:1939307/1‑76 (MQ=255) ttCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCCGTa > 2:56128/1‑76 (MQ=255) ttCAGATAAGGAATGC‑TTTTACCCCCG‑CCGATGTCCATCGTGAGGAAGATAATCGGCTTGTCGCCCGCCGCCAGTa < 1:3721762/76‑1 (MQ=255) tCAGAAAAGGATTTC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAGTAg < 1:1659462/76‑1 (MQ=255) gcAGATAAGGAATGT‑TTTTTCCCCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAGTAg < 1:5702479/75‑1 (MQ=255) cAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCTCCCGCCGCCCGTAga > 2:644291/1‑75 (MQ=255) cAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTTGCCCGCCGCCagcagc > 1:5600713/1‑76 (MQ=255) aGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCCGCCGTAGCa > 2:5360244/1‑76 (MQ=255) aGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCCGTAgaa > 1:3580500/1‑74 (MQ=255) gATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTCGCCGCCCCCCTTAGCAc > 1:5012394/1‑76 (MQ=255) aTAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCCGTAGCAcc > 2:6633769/1‑76 (MQ=255) tAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCCCCCGAAGCAccc > 1:3887218/1‑75 (MQ=255) tAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCAGTAGCAccc > 2:3717555/1‑75 (MQ=255) aaGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCCCCCGAAGAAcccc > 2:4364903/1‑74 (MQ=255) aaGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCCACAGCAcccg > 1:3604148/1‑74 (MQ=255) aGGAGTGC‑TTTTCCCGCCG‑CGGATGTCCATCTTGGGGAAGATAACCGGCTTGCCCCCCGCCGCCAGTAGCACctgc < 1:3461154/76‑1 (MQ=255) aGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCGCCCGTAGCACccgc > 2:2040363/1‑76 (MQ=255) ggAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCCCCCGCCGCCCGGAGCACccgcg > 2:1818847/1‑75 (MQ=255) aTGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTCCCCGCCGCCACTAGCCCCTGCTggg > 2:7595756/1‑75 (MQ=255) aTGC‑TTTTATCGCCG‑CCGATGTCCATCGCGCGGAAGATAATCGGCTTTTTTCCCCCCCCCCGTAGCACCCGCTggg > 1:2954184/1‑75 (MQ=255) ttgttttGTCGCCCCCCGATGCCCATCGTGCGGAAGATAATCGGCTTGCCCCCCGCCGCCAGTAGCACCTGCTGGt < 1:2693569/73‑1 (MQ=255) gC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCCCCCGTAGCACCTGCTGGttg > 2:3112009/1‑74 (MQ=255) tttttACCCCCG‑CCGATGCCCATCGGGCGGAAGATAATCGGCTTGTCCCCCGCCGCCAGTAGCACCTGCTGGTAgg < 1:5170031/75‑1 (MQ=255) tttATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGGCCCCCCGAGGACCCGCTGGCAgggc > 1:352947/1‑74 (MQ=255) tttATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCCGTAGCACCTGCCGGGAGGCt > 1:1376139/1‑76 (MQ=255) aTCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTTCCCGCCGCCAACAGCACCCGCCGGTCGGCGT‑Ca > 2:4005241/1‑76 (MQ=255) aTCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCCGCAGGCCCTGCCGGGAGGGTT‑Ca > 1:5278062/1‑76 (MQ=255) cccccg‑cggTTGCCCATCGTGCGGAAGATAATGGGCTTGTCCCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑Caa < 1:1376995/73‑1 (MQ=255) cgccg‑cCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCGCCCCCAGCACCCGCTTGTGGGGTT‑Ccaa > 2:30543/1‑76 (MQ=255) cccg‑ccgTTGCCCTTCGCGCGAAAAAAAACCGGTTGGCCCCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAAt < 2:5119036/75‑1 (MQ=255) cccg‑ccgATGCCCATCGGGCGAAAGATAAGCGGCTTGCCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAAt < 1:6499019/75‑1 (MQ=255) cg‑ccgATGTCCATCGTGCGGAAGATAATCGGCTTTTTGGCCCCCCCCCGTAGCACCCGCTGGTGGGGTTCCAAAtt > 2:6103861/1‑76 (MQ=255) cg‑ccgATGTCCATCGTGCGGAAGATAATCGGCTTATGGCCCGCCGCCAGAAGCACCTTCTGGTCGGCTT‑CTAAttt > 2:7371671/1‑75 (MQ=255) g‑ccgATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCCGTAGCACCTGCTGGGAGGCCT‑CCAAATGc > 1:3470927/1‑76 (MQ=255) g‑ccgATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCAGTAGCACCCCCTGGGAGGGTT‑TAAAATGc > 1:3640041/1‑76 (MQ=255) cgATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCGCCCGTAGCACCTGCTCGGAGGGTT‑CAAAATGCTc > 1:7496429/1‑76 (MQ=255) cgATGTCCATCGTGCGGAAGATAATCGGCTTTTCCCCCGCCCCCCGTAGAACCTGCTGGGTGGCGT‑CACAATGCTc > 1:241314/1‑76 (MQ=255) cgATGCCCATCGGGCGGAAGAAAATCGGCTTGTCCCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTc < 1:1499290/76‑1 (MQ=255) gATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCCGTAGCACCTGCTGGGAGGCCT‑CCAATTGCTcc > 1:5700168/1‑76 (MQ=255) aTGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCCCCAGGAGCACCTGCTGGTAGGGTT‑CAAATTTCTCCt > 1:4318868/1‑76 (MQ=255) aTGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCCGTAGCACCTGCTGGGAGGGTT‑CAAATTTCTCCt > 1:4093771/1‑76 (MQ=255) aTGCCCATCGGGGGAAAGATAACCGGTTTCCCCCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCt < 1:5197506/76‑1 (MQ=255) gccccAGCGTGGGAAAGATAATGGGTTTGCCCCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTg < 1:5190396/73‑1 (MQ=255) ccATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCCCCCGTAGCCCCTGCCGGGAGGGTT‑CAAAATTCTCCCGCTc > 1:6066499/1‑76 (MQ=255) cATCGGGGGGAAGATAATCGGCTTGCCCCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCg < 1:3782429/76‑1 (MQ=255) aTCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCCGTAGCACCTGCTGGGAGGGTT‑CAAAATGCTCCTGCTCGt > 1:7298979/1‑76 (MQ=255) aTCGGGCGAAAGATAATCGGTTTGTCCCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGt < 2:4008680/76‑1 (MQ=255) gTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCCGTAGCACCCGCTGGGAGGGTT‑TTAATTTCTCCTGCTCGTCCg > 2:6142047/1‑76 (MQ=255) gTGCGAAAAATAGCCGGCTGGCCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGTCCg < 2:1917010/76‑1 (MQ=255) gggggAATATAGTCGTCTTCCCCCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGTCCgg < 2:3766614/73‑1 (MQ=255) gCGAAAAATAACCGCCTTGCCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGTCCGgc < 1:2068015/76‑1 (MQ=255) gaaaGATAATCGTTTTGTCCCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGTCCGgcgc < 1:5606838/74‑1 (MQ=255) caaGATAATCGGCTTTTCGCCCGCCGCCAGTAGCACCCGCTGGGAGGCTT‑CCAATTGCTCCCGCTCGTCCGgcgcg > 1:2131551/2‑76 (MQ=255) gtcTTGCCCCCCCCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGAt < 2:5055155/74‑1 (MQ=255) gCTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGGTc > 2:2606230/1‑76 (MQ=255) gCTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATc > 1:3033109/1‑76 (MQ=255) gCTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCCGCTCGGCCGGCGCGCTGTCGCGATc > 2:4713914/1‑76 (MQ=255) cTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATcc > 1:6338367/1‑76 (MQ=255) cTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATcc > 2:2853639/1‑76 (MQ=255) ttGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATcct > 2:6202591/1‑75 (MQ=255) ttGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCa > 2:6462737/1‑76 (MQ=255) ttGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCa > 1:855324/1‑76 (MQ=255) tGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCat > 1:5732063/1‑76 (MQ=255) | CTCTTCCTGGGGAATATTCAGATAAGGAATGC‑TTTTATCGCCG‑CCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTT‑CAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCAT > REL606/2437331‑2437480 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |