Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,437,406 | 0 | G | T | 36.6% | 83.2 / 21.3 | 41 | D435E (GAC→GAA) | ypdD | fused predicted PTS enzymes: Hpr component/enzyme I component/enzyme IIA component |
Reads supporting (aligned to +/- strand): ref base G (15/11); new base T (15/0); total (30/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.99e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCTTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCAT > REL606/2437331‑2437480 | ctctTCCTGGGGAATATTCAGATAAGGAATGTTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTg < 1:4066367/76‑1 (MQ=255) ctctTCCGGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTg < 2:2779256/76‑1 (MQ=255) tctTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGt < 1:1481513/76‑1 (MQ=255) ctTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTc > 2:4813868/1‑76 (MQ=255) ctTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTc > 2:3982728/1‑76 (MQ=255) ttCGTGGGGAATATTCAGATAAGGAATGTTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcg < 2:5795214/76‑1 (MQ=255) ttCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTct > 1:3296709/1‑75 (MQ=255) tCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgc > 1:1196344/1‑76 (MQ=255) ccTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgcc > 2:1034977/1‑76 (MQ=255) ccTGGGGAATATTCAGATAAGGAATGCTTTTACCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgcc < 2:1350318/76‑1 (MQ=255) cgggggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgccc < 1:1669705/74‑1 (MQ=255) cgggggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgccc < 2:2869702/74‑1 (MQ=255) cTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgccc > 1:99377/1‑76 (MQ=255) cTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgccc < 2:3258756/76‑1 (MQ=255) cTGGGGAATATTCAGATAAGGAATCCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgccc < 2:5169589/76‑1 (MQ=255) cTGGGAAATATTCAGATAAGGATTGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgccc < 2:4111798/76‑1 (MQ=255) tGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTcgcccc > 2:5854190/1‑75 (MQ=255) tGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgcccg < 2:3387610/76‑1 (MQ=255) ggggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgccccc > 1:5009456/1‑74 (MQ=255) ggggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgccccc > 1:3694221/1‑74 (MQ=255) ggaaaTATTCAGATAGGGAATGCTTTTTCCCCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgcccgcc < 2:3741016/73‑1 (MQ=255) ggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTcgccccccc > 1:5479876/1‑75 (MQ=255) ggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTTCccgccg > 1:2953967/1‑76 (MQ=255) ggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCcccccg > 1:847707/1‑76 (MQ=255) gAATATTCAGATAAGGAATGCTTTGATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTGTGCCccccgc > 1:3126874/1‑76 (MQ=255) aaTATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTGTGCCCgccgcc > 1:2680758/1‑76 (MQ=255) atatTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCccccccccc > 2:2045653/1‑75 (MQ=255) atatTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCgcccccc > 1:5739232/1‑75 (MQ=255) tatTCATATAAGGATTGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAg < 1:3143642/76‑1 (MQ=255) tatTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCGCCAg > 2:3126287/1‑76 (MQ=255) atTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTTgccgccgccccg > 1:4217752/1‑73 (MQ=255) atTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCCCCAGt > 1:3303881/1‑76 (MQ=255) atTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCgccccccgc > 1:3518052/1‑73 (MQ=255) atTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCCCCCCCact > 1:5465524/1‑74 (MQ=255) ttCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTTGCCCGCCGCCCGTa > 2:2880166/1‑76 (MQ=255) cAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGTCCCCCGCCCGTAGc > 1:308252/1‑76 (MQ=255) aGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCAGTAGCa > 2:5332685/1‑76 (MQ=255) aGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCCgcagca > 1:5216795/1‑76 (MQ=255) aGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTTCCCCCCCGCCCgcagca > 2:6201639/1‑76 (MQ=255) gATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCCGCCGCCGCCAGTACCAc > 2:3358927/1‑76 (MQ=255) aTAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCGCCCGTAGCAcc > 2:4758633/1‑76 (MQ=255) tAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCCCCCGCAGCACct > 2:5490314/1‑76 (MQ=255) tAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTTGCCCGCCGCCCGCAGGAccc > 1:5543286/1‑75 (MQ=255) aaGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTTCCCCCCGCCCGCAGCAcccc > 1:1642731/1‑74 (MQ=255) aaGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCGCCCGTAGCAcccc > 1:441523/1‑74 (MQ=255) aaGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCCCCCGTAGCAcccc > 1:5118335/1‑74 (MQ=255) aaGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCCCCAGTAGCAcccg > 2:196951/1‑74 (MQ=255) tGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCAGCAGCCACCGCTggga > 1:1921307/1‑74 (MQ=255) gCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCGCCCGGAGCACCCGCTgggtg > 1:5199326/1‑73 (MQ=255) gCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCCGTAGCACCCGCTGGTA‑g > 2:5760816/1‑76 (MQ=255) ttttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCTCCCCCCGCCAGTAGAACCCGCTTGGA‑ggg > 2:3689302/1‑75 (MQ=255) tttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCCGCAGGCCCCGCTTGGA‑GGC‑t > 2:1939094/1‑76 (MQ=255) tttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCCCCCCCAGTAGAACCTGCTTGTA‑GGcg > 2:4643877/1‑75 (MQ=255) tcccgccgcCGTTTTCCAGGGGGCGGAAGATAATCGGTTGGCCGCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTc < 1:3027371/73‑1 (MQ=255) tcccgccgcCGATGCCCATCGGGCGAAAGATAATCGGCTTGTCCCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTc < 1:5809420/73‑1 (MQ=255) tATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCGCCCCCGTAGCACCTGCCGGGA‑GGC‑ttt > 2:2880194/1‑75 (MQ=255) aTCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCCGGAGCACCTGCTGGGA‑GGG‑ttta > 2:4293436/1‑74 (MQ=255) tccccgccgATGTCCATCGTGCGAAAGATAAGCGGCTGGCCCCCCCCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCaa < 1:53346/73‑1 (MQ=255) gccgccGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCCCCAGCAGGAACCGCTGGGC‑GGC‑GTCaaaa > 2:5892401/1‑75 (MQ=255) gccgccGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGGCCGCCGCCCGGAGCACCCGCTGGGA‑GGC‑TCTAAAt > 2:3901026/1‑76 (MQ=255) cccgccgATCCCCTGCGTGCGAAAGATAATCGGTTTGCCCCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAAt < 2:5417150/75‑1 (MQ=255) ccgccgATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCCGTAGCACCCGCTGGGA‑GGG‑TTCAAAtt > 2:3244050/1‑76 (MQ=255) ccgccgATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCCGAAGCACCTGCTGGTG‑GGG‑TTCaaaat > 2:3735280/1‑74 (MQ=255) cgATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCCGGAGCACCTGCTGGTA‑GGG‑TTCAAATTGCTc > 1:768770/1‑76 (MQ=255) gATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCCCTAGCACCTGCCGGGA‑GGC‑TTTAAATTACTcc > 1:5621724/1‑76 (MQ=255) ttGTCCATCGTGCGGAAAATAATCGGTTTGCCGCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCt < 1:6145165/75‑1 (MQ=255) tggTCCATCCGCCGGAAGATAACCGGCTTGCCCCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCt < 2:1403401/74‑1 (MQ=255) tGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCAGTAGCACCCGCTGGTA‑GGG‑GTCAAAATTCTcccg > 2:3310000/1‑74 (MQ=255) tGTCCATCGTGCGGAAGATAATCGGCTTGTTGCCCGCCCCCCGGAGCACCCGCTGGTA‑GGCTTTCAATTTGCTCCt > 2:2391103/1‑76 (MQ=255) gTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCCGTAGCACCCGCTGGTA‑GGG‑GTCAAAATTCTCCCGc > 2:6044732/1‑76 (MQ=255) gTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCCGTAGGACCTGCTGGGA‑GGG‑TTCAAAATGCTCCtcc > 1:4807483/1‑74 (MQ=255) tCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCCGTAGCACCTGCTGGGA‑GGG‑CTCAAAATGCTCCtccc > 1:2507905/1‑73 (MQ=255) cgcGTGGGAAAGATAATCGGTTTGCCCCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGt < 2:1482323/74‑1 (MQ=255) aTCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCCGTAGCACCTGCTGGTAGGGT‑TTCAAATT‑CTCCTGCTCGt > 1:3372906/1‑76 (MQ=255) cGTGGGGAAGATAATCGGCTTGTCCCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTcc < 2:3078074/76‑1 (MQ=255) cGTGCGGAAGATAATCGGCTTTTCGCCCGCCCCCAGTAGCACCCGCTGGGA‑GGC‑TTCAAATTGCTCCTCCTCGccc > 2:2130766/1‑76 (MQ=255) gTGCGGAAGATAATCGGCTTTTTGCCCCCCGCCCGTAGCACCCGCTGGGA‑GGG‑TTCAAAATTCTCCTGCtcctccc > 2:4580981/1‑75 (MQ=255) gggAAGATAATCGTTTGGCCCCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGgcg < 2:1153223/75‑1 (MQ=255) ggaaaGATAATGGCCTTGCCCCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGgcg < 1:3647789/73‑1 (MQ=255) cGGAAGATAATCGGCTTTTGGCCCGCCGCCCGTAGGAACTGCTCGTG‑GGG‑TTTAAATTTGTCCTGCTCGGTCGgcg > 2:3636563/1‑76 (MQ=255) ggAAGATAATCGGCTTTTCGCCCGCCGCCCGTAGCACCCGCTGGGG‑GGG‑TTTAAAATTTTCCTGCTCCTCCGgggc > 2:4887911/1‑76 (MQ=255) aGAAAAGCGGCTGGCCCCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCt < 1:6126590/76‑1 (MQ=255) actCGTCTTGCCCCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTcg < 1:3011325/74‑1 (MQ=255) tCGGGTTGGCCCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTcgcg < 1:1556542/76‑1 (MQ=255) tCGGCTTGCCCCCCCCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTcgcg < 2:3303881/76‑1 (MQ=255) cTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATcc > 1:6208401/1‑76 (MQ=255) cTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATcc > 1:4266419/1‑76 (MQ=255) cTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATcc > 1:3995494/1‑76 (MQ=255) cTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCCCGAGcc > 2:1964754/1‑76 (MQ=255) ttGTCGCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCa > 1:843068/1‑76 (MQ=255) ttGTCGCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCa > 1:5502008/1‑76 (MQ=255) ttGTCGCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCa > 2:3072659/1‑76 (MQ=255) tGTCGCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGGTCCat > 1:2418396/1‑76 (MQ=255) tGTCGCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCat > 1:6167165/1‑76 (MQ=255) tGTCGCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCat > 2:5773075/1‑76 (MQ=255) tGTCGCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCat > 1:358690/1‑76 (MQ=255) | CTCTTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTA‑GGC‑TTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCAT > 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: ‑ |