| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 2,437,406 | 0 | G | T | 36.7% | 42.5 / 14.6 | 30 | D435E (GAC→GAA) | ypdD | fused predicted PTS enzymes: Hpr component/enzyme I component/enzyme IIA component |
| Reads supporting (aligned to +/- strand): ref base G (13/6); new base T (11/0); total (24/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.13e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTCTTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCAT > REL606/2437331‑2437480 | ctctTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGTTTg > 2:7291156/1‑76 (MQ=255)ctctTCCTGGGGAATATTAAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTg < 2:2992943/76‑1 (MQ=255) tctTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGt > 2:3687825/1‑76 (MQ=255) tctTCCTGGGGAATATTCAGAGAAGGAATGCTTTTACCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGt < 1:260536/76‑1 (MQ=255) ttCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcg < 2:1377948/76‑1 (MQ=255) tCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTccc > 1:4692702/1‑74 (MQ=255) ccTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgcc > 1:7016865/1‑76 (MQ=255) cTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgccc > 1:7169123/1‑76 (MQ=255) cTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTctccc > 1:3724421/1‑76 (MQ=255) cTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgccc < 2:1878020/76‑1 (MQ=255) tGGGGAATATTCAGATAAGGAATTCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgcccg < 1:5964834/76‑1 (MQ=255) tGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgcccc > 1:2060034/1‑75 (MQ=255) tGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgcccc > 1:5940516/1‑75 (MQ=255) tGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTcgcccc > 2:264072/1‑75 (MQ=255) tGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgcccg > 1:5025007/1‑76 (MQ=255) ggggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgccccc > 1:2121229/1‑74 (MQ=255) ggggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgccccc > 1:5056518/1‑74 (MQ=255) gggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTctcccccc > 1:2773572/1‑76 (MQ=255) ggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTcgcccgccc > 1:5917251/1‑75 (MQ=255) ggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTcgccccccc > 1:7419888/1‑75 (MQ=255) gaaaTATTCATATAAGAAATGTTTTTATCGCCGCCGATGTCCATCGGGCGGAAGATAATCGGCTTGTCGCccgccg < 2:7929839/74‑1 (MQ=255) gaaaTATCAAGAAAAGGATTTCTTTTACCCCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCccgccg < 2:2741897/74‑1 (MQ=255) aaTATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCcgcgcc > 2:2594906/1‑76 (MQ=255) aaTATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTTGCCCcccgcc > 1:3374767/1‑76 (MQ=255) aaTATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCTCCCgccccc > 2:1994172/1‑76 (MQ=255) atatTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCgcccccc > 1:4106881/1‑75 (MQ=255) atatTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCgcccccc > 1:1085572/1‑75 (MQ=255) atatTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCTCCCgcccccc > 2:1788588/1‑75 (MQ=255) atatTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCgccgccc > 2:8114375/1‑75 (MQ=255) atatTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCACCCCGCGCCa > 1:7333679/1‑76 (MQ=255) atatTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCcccgccc > 1:1291719/1‑75 (MQ=255) tatTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCCCCAg > 1:2134986/1‑76 (MQ=255) tatTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCCCCCCCac > 2:4211715/1‑75 (MQ=255) atTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAGt > 2:5785421/1‑76 (MQ=255) ttCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCCCCACTa > 2:3650789/1‑76 (MQ=255) tCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCCGTA‑g > 1:6365989/1‑76 (MQ=255) cAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCCCCCATA‑Gc > 1:4407128/1‑76 (MQ=255) aGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCCCCAGGA‑ggc > 1:1899136/1‑74 (MQ=255) aGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCGCCCGTA‑gga > 2:819184/1‑74 (MQ=255) gATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCTCCCCCCCCCCGCA‑GCAc > 1:2277831/1‑76 (MQ=255) aTAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCGCCCGTA‑Gccgc > 1:7802284/1‑73 (MQ=255) aTAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCACGA‑GCAcc > 2:4908662/1‑76 (MQ=255) aTAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCCCCCCGCCCCCGCA‑GCAcc > 2:3161974/1‑76 (MQ=255) tAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCACCCCCCGCCAGCA‑GGAccc > 2:4277414/1‑75 (MQ=255) aaTGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCCCCCGGA‑GCAActgcccg > 2:1216507/1‑73 (MQ=255) aaTGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCCGCA‑GCACCCGCTgg > 1:7556870/1‑76 (MQ=255) aaTGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCACGA‑GCACCCGCCgg > 2:1052012/1‑76 (MQ=255) aaTGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGCCTTGTCGCCCGCCGCCCGTA‑ACACCAGCTgg > 2:5936719/1‑76 (MQ=38) aTGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTTCGCCGCCGCCAGTG‑CAACCTGCTggg > 2:8124601/1‑75 (MQ=37) gCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCGCCCGTA‑ACACCTGCTggggg > 1:7565673/1‑73 (MQ=255) gCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCCCCAGTC‑GCACCTCCTGGTA‑g > 2:5144557/1‑76 (MQ=255) gCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCAGGA‑GCACCTGCCGGGA‑g > 2:5750525/1‑76 (MQ=255) tttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCGCCCGTA‑GCACCTCCTTGGA‑GGC‑t > 1:971328/1‑76 (MQ=255) tttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCCCCCCCCGCCAGAA‑ACAACTGCTGGCA‑GGca > 1:5905501/1‑75 (MQ=255) ttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCCCCCGTA‑GCACCCGCTGGGA‑GGcgt > 2:3303698/1‑74 (MQ=255) ttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCACCCCCCGCCACCT‑GCACCTCCCGGTT‑GGt‑tt > 2:5711801/1‑76 (MQ=255) tATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCCGTAGGCACC‑GCCGGGA‑GGC‑ttt > 2:7821435/1‑75 (MQ=255) tATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCCGTA‑GCACCTGCTGGTA‑GGC‑ATc > 1:3800294/1‑76 (MQ=255) aTCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCCATA‑GCACCCGCCGGGA‑GGG‑tttc > 2:7305034/1‑74 (MQ=255) tCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCAGTA‑GCACCTGCTGGGA‑GGC‑TTcca > 2:5950497/1‑74 (MQ=255) ccgccgcCGATGTCCATCGTGCGAAAGATAATCGCCTTGCCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCaa < 1:732818/75‑1 (MQ=38) cccccgccgATGTCCATCGTGCGGAAGATAATCGGCTTGCCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCaa < 1:6720983/73‑1 (MQ=255) cgccgcCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCCGGA‑GCAACTGCTGGGA‑GGG‑TTCcaa > 1:263953/1‑76 (MQ=255) cgccgcCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCCGTA‑GCACCTCCTGGGA‑GGG‑TTCcaa > 1:870673/1‑76 (MQ=255) gccgccGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAGTA‑GCACctgc > 2:3740240/1‑61 (MQ=255) cccgccgATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAGTA‑GCACctgc < 1:3740240/60‑1 (MQ=255) ccgccgATGCCCGTGGTGGGGAAGATAAGCGCTTTGCCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAAtt < 1:4083889/76‑1 (MQ=255) ccgccgATGCCCAGCGGGCGGAAGATAACCGGCTGGCCCCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAAtt < 1:7277886/76‑1 (MQ=255) gccgATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCCCCAGTA‑GGACCCGCTGGGA‑GGG‑TTTCAATTgt > 1:2335437/1‑75 (MQ=255) gccgATGTCCATCGTGCGGAAGATAATCGGCTTGTTGCCCCCCGCCCGTA‑GCACCCGCTGGGT‑GGG‑TTTTAATTGc > 1:1800851/1‑76 (MQ=255) ccgATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCCGTA‑GCACCCGCCGGGA‑GGCATCCAAATTgt > 1:5533516/1‑75 (MQ=255) cagATGCCCATGGGGGGAAAGATAATCGGTTTGCCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCt < 2:2860606/74‑1 (MQ=255) gATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCAGTA‑GAACCCGCTGGGA‑GGC‑CTTAAAATTCTcc > 1:4939908/1‑76 (MQ=255) gTCCATCGTGCGGAAGATAATCGGCTTTTCTCCCCCCGCCCCAA‑GCACCCGCTGGGA‑TGG‑TTCAAAATGATCCTGc > 1:4332801/1‑76 (MQ=255) tCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCAGAA‑GGACCTCCTTGTT‑TGA‑CTTAAATTGGTCCttgc > 2:565848/1‑73 (MQ=255) ccATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCAGTA‑GCACCTGCTGGGAGGGC‑TCCAAATT‑CTCCCGCTc > 1:4219457/1‑76 (MQ=255) ccAGCGGGCGGAAGATAATCGGCTTGCCCCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTc < 2:3987874/76‑1 (MQ=255) cATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCAGGA‑GCACCCGCTGGGA‑GGG‑TTCCAAATGCTCCTGCCCg > 1:2030037/1‑76 (MQ=255) cGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCCGTA‑GCACCTGCTGGGT‑GGC‑GTCCAAATGCTCCCGCTCGTcc > 1:1401979/1‑76 (MQ=255) gTGCGGAAGATAATCGGCTTTTCGCCCGCCCCCCGGA‑GCACCCCCTGGGA‑GGC‑GTCAAATTGCTCCTGCCCCTCCg > 2:2638412/1‑76 (MQ=255) gTGCGGAAGATAATCGGCTTTTCGCCCGCCCCCAGTA‑GCAACTGCTGGTG‑GGG‑GTTTAATTGCTCCCTCTCGTccc > 2:3218245/1‑75 (MQ=255) gggggAAAATAACGGGCTTGCCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCgg < 1:4515728/73‑1 (MQ=255) ggAAGATAATCGGCTTTTTGCCCGCCGCCCGTA‑GGACCCTCTGGGA‑GGG‑TTTCAACTGCTCCTGCTCGTCCGgcgc > 2:8063541/1‑76 (MQ=255) aaGATAATCGGCTTTTTGCCCCCCCCCAGTA‑GGACCTGCCGGTA‑GGC‑CTCAAAATTCTCCTTCTCGTCCGgcgcgg > 2:4685671/1‑75 (MQ=255) tAATCGCCTTGCCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTc < 1:6694701/76‑1 (MQ=255) aaTCGCCTTGCGGCCCTCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTcg < 1:6547082/76‑1 (MQ=255) atgggCTGGTCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTcgc < 2:7955705/73‑1 (MQ=255) ccGGTTTGTGACCCCCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTcgcg < 1:3782187/75‑1 (MQ=255) gCTTGTCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTTTCGCGGTc > 1:5988616/1‑76 (MQ=255) cTTGTCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGGGAGCTGTCCCGCTcc > 2:561066/1‑76 (MQ=255) cTTGTCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATcc > 2:4885761/1‑76 (MQ=255) cTTGTCGCCCGCCGCCAGTA‑GCACCTGCTGGGA‑GGC‑TTCAAATTGGTCCTGCCTGACCGGCCCGCTGTCGCGacac > 2:1755505/1‑73 (MQ=255) ttGTCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCa > 2:5244264/1‑76 (MQ=255) ttGTCGCCCGCCGCCAGTA‑GCACCTGCTGGAA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCa > 2:2510421/1‑76 (MQ=255) tGTCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCat > 2:2289170/1‑76 (MQ=255) tGTCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCat > 1:1680805/1‑76 (MQ=255) tGTCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCCCTGTTGCGATCCat > 2:6597717/1‑76 (MQ=255) | CTCTTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAGTA‑GCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCAT > REL606/2437331‑2437480 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |