Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,180,249 | 0 | T | G | 41.7% | 197.7 / 43.4 | 120 | G24G (GGT→GGG) | yjaH | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (18/52); new base G (50/0); total (68/52) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.70e-19 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
Only 100 of 223 total aligned reads displayed.
GGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAG > REL606/4180174‑4180324 | cAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTgc < 1:613061/76‑1 (MQ=255) cAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTgc < 2:3155785/76‑1 (MQ=255) aGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTgcg < 1:1617128/76‑1 (MQ=255) aGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGgcg > 2:4154786/1‑76 (MQ=255) aCAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑tt < 1:6367070/76‑1 (MQ=255) cAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCCG‑ttt > 2:2537566/1‑76 (MQ=255) cAGTGTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑ttt < 2:5488702/76‑1 (MQ=255) aGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTtc < 1:5286038/76‑1 (MQ=255) aGTTTTAACGAAGGGGTGGTGTCACCCCCTTGGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTtc < 2:3981947/76‑1 (MQ=255) gTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTtct < 2:5347484/76‑1 (MQ=255) tttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCa < 1:1521393/76‑1 (MQ=255) ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGGTT‑TTtttccc > 1:1032410/1‑74 (MQ=255) ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTtcttcc > 2:1094233/1‑73 (MQ=255) tAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑t < 2:1627425/76‑1 (MQ=255) tAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑t < 2:2275424/76‑1 (MQ=255) tAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑t < 1:6647876/76‑1 (MQ=255) cGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTC‑TTTTTCAC‑TGcc > 2:4028814/1‑76 (MQ=255) cGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGGCG‑TTTTTCAC‑TGcc > 1:4440286/1‑76 (MQ=255) gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCt < 2:2396291/76‑1 (MQ=255) gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TGccc > 1:5239265/1‑75 (MQ=255) aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTg < 2:3342528/76‑1 (MQ=255) aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGGTG‑TTTTTTAC‑TTCCTg > 1:4326911/1‑76 (MQ=255) aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTG‑TTTCTTAC‑TGCCTg > 2:430533/1‑76 (MQ=255) aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑GTTCTCAC‑TGCCTg > 2:836361/1‑76 (MQ=255) aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGAGCTG‑TTTCTTAC‑TCCCtt > 1:2589582/1‑75 (MQ=255) aGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑GTTCTCCC‑GGcccgc > 2:3517760/1‑73 (MQ=255) ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTa < 1:6469765/76‑1 (MQ=255) ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCCG‑GTTTTCAC‑TGCCCGTa > 1:2900335/1‑76 (MQ=255) gggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAg < 2:3579681/76‑1 (MQ=255) gggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAg < 2:759284/76‑1 (MQ=255) gggTGGTATCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAg < 2:3637217/76‑1 (MQ=255) gtggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAGcc < 1:6047387/76‑1 (MQ=255) gtggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAGcc < 2:7119878/76‑1 (MQ=255) gtggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAGcc < 2:5726181/76‑1 (MQ=255) gtggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑GTTCTCAC‑CCCCCGTAGcc > 1:5421613/1‑76 (MQ=255) tggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TTCCTGTTGccc > 1:1231950/1‑75 (MQ=255) tggtTTCACCCCTTTTGACTTTCTGGCGTCGATCATTGATGCTGGCAGGGACTGCG‑TTTCTCCC‑TCCCCCTAGCca > 2:6724665/1‑76 (MQ=255) tgtTTCACCCTTTTTGTCTTTCTGGCGTCGTTCATTGATGCTGGCAGGTGCGCTG‑TTTCTCCC‑TGCCTTTAGCcac < 2:1328613/75‑1 (MQ=255) ggtTTCCCCCTTTTTGTCTTTCTGGCTTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCCC‑TGCCTGTAGCcac < 2:568403/76‑1 (MQ=255) gtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAGCcaca < 1:3830517/76‑1 (MQ=255) tttCCCCCTTTTTTTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGTTT‑TTTCTCCC‑TGCCTTTAGCCACaa < 2:4327232/76‑1 (MQ=255) ttCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGTCG‑CTTCTTAC‑TGCCCCTTGCcaccac > 2:392755/1‑76 (MQ=255) tccccctTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCGGGTGCGCTG‑TTTCCCCC‑TGCTTTTAGCCACAAct < 2:3992292/73‑1 (MQ=255) tCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTG‑TTTCTCAC‑TTCCTCTAGCCACaaac > 1:6997732/1‑74 (MQ=255) tCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAAct > 2:308543/1‑76 (MQ=255) aCCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TCCCTGTAGCCACCAccct > 1:2812216/1‑76 (MQ=255) cccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TGCCTGTAGGCACAACTCTTc > 1:6847869/1‑76 (MQ=255) ctttttttttttGTGGCGTCGATCATTTATGCTGGCGGGGGCGCTT‑TTTCTCCC‑TGCCTTTAGCCACAACTCTTCa < 2:1583099/74‑1 (MQ=255) cTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑CGCCCGTTGCCACAActcctccc > 2:1040001/1‑74 (MQ=255) ttttttttttCGGGGGTCGATATGTGAGGCTGGGAGGTTCGTTT‑TTCCACAT‑TCCCTTTAGCCACAACTCTTCACt < 2:1527041/76‑1 (MQ=255) ttttGTTTTTCGGGCGTCGATCATGAAGGCGGGCAGGTGCGTTT‑TTTCTCCC‑TGCCTTTAGCCACAACTCTTCACt < 2:4377352/76‑1 (MQ=255) ttttGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTG‑TTTCTCAC‑TGCCTGTAGCCACAAATCTTTACt > 1:3586814/1‑76 (MQ=255) ttttGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TGCCTGTAGCCCCAAATCTTCACt > 1:4482017/1‑76 (MQ=255) tttGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACtt < 1:117008/76‑1 (MQ=255) tttGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTTTCAC‑TGCCTGTAGCCACAACCCTTCTCtt > 1:7154522/1‑76 (MQ=255) tttGTCTTGTTGGCGTCGATTATGGATGTTGGAGGGTGTGCTG‑TTTCTCCC‑GGCTGTTAGCCACAACTCTTCACtt < 2:6292930/76‑1 (MQ=255) ttGTCTTTCGGGCGTCGATCATTGATGCTGGCGGGTGCGCTG‑TTTCTCCC‑TGCCTTTAGCCACAACTCTTCACTtc < 2:6001464/76‑1 (MQ=255) ttttttCGGGCGTCGATCATGGATGCTGGCGGGGCCGCTT‑TTTCCCCC‑TGCCTTTAGCCACAACTCTTCACTtcct < 2:776381/73‑1 (MQ=255) ttCTTTCGGGGCTCGATTATTGATGCGGGCGGGGCCGCTT‑TTTCTCCC‑TGCCTTTAGCCACAACTCTTCACTtcct < 2:3130253/75‑1 (MQ=255) gTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACTtcct < 1:4216952/76‑1 (MQ=255) gTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCCG‑TTTCTCACTTCCCTGT‑TCCACCAATCTTCAcctcct > 2:1686606/1‑76 (MQ=255) gTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTTACTtcct > 2:3771669/1‑76 (MQ=255) tCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTG‑TTTCTCAC‑TGCCTGTAGCCACCACTCTTCACTtcctc > 2:2728959/1‑76 (MQ=255) tCTTGCGGGGGTCGATAATTGATGCTGGCAGGGGCGTTT‑TTCCTCAC‑TCCTTTTAGCCACAACTCTTCACTtccta < 2:2015889/76‑2 (MQ=255) cTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACTtcctcc < 1:3716759/76‑1 (MQ=255) cTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTG‑TTTTTCAC‑TGCCTGTAGCCACAAATCTTCACTtcctcc > 1:7546255/1‑76 (MQ=255) cTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTG‑TTTCTCCC‑TGCCTGTAGCCCCCACTCTTCACTtcctcc > 1:2037915/1‑76 (MQ=255) cTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACTttctcc > 2:2922362/1‑76 (MQ=255) cTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCCG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACTtcctcc > 2:7506639/1‑76 (MQ=255) tttCTGGCGTGGATCATTGATGTTGGCAGGGGCGCTG‑TTTCTCAC‑TCCCTGTAGCCCCAACTCTTCACTTCCTccc < 1:913498/76‑1 (MQ=255) tttCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACTTCCTccc > 1:7420558/1‑76 (MQ=255) tttCTGGCGTCGATCATTGATGCTGGCAGGGGCGCGG‑TTTCCCAC‑TGGCTGTAGCCCCAACTCCTTACTTcccccc > 1:3910594/1‑76 (MQ=255) ttCTGGTGTCGATTTTGGATGTGGCCGGGTGCTTTG‑TTCCTCTC‑TGCCTGTAGCCACAACTCTTCACTTCCTCCCt < 1:2990893/76‑1 (MQ=255) ttCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACTTCCTCCCt > 2:7562125/1‑76 (MQ=255) tCTGGGGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTCCTCAC‑TGCCTGTAGCCACAACTCTTCACTTCCTCCCtt < 1:2627605/76‑1 (MQ=255) tCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCGC‑TGCCCGTAGCCACAACTCTTTACTTCCTcccct > 2:3834497/1‑74 (MQ=255) tCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TGCCTGTTGGCGCCACTCTTCACTTCCTCCCtt > 2:1812599/1‑76 (MQ=255) cTGGCGTCGATCATTGATGCTGGCAGGGGCGGGGCTTTCT‑AC‑TGTCTCTAGCCACCACTCTTTATTTCCTTCCttt > 1:4086362/1‑76 (MQ=255) cTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCCC‑TGCCCGTAGCCACCACTCTTCACTTCCTCCCttt > 1:5895272/1‑76 (MQ=255) cTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TGCCTGTAGCCCCCACTCTTCACTTCCTCCCttt > 1:6429436/1‑76 (MQ=255) cTGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCCCTTCCTCCCttt > 1:1425089/1‑76 (MQ=255) cTGGCGTCGATCATTGATGCGGGCGGGTGCGTTT‑TTTCTCCC‑TGCTTTTAGCCACAACTCTTCACTTCCTCCCttt < 2:4931141/76‑1 (MQ=255) tGGCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTTTCAC‑TGCCTGTAGCCACAACTCTTCACTTCCTCCCtttt > 1:5051287/1‑75 (MQ=255) tGGCGTCGATCATTGATGCTGGCAGGGGCGCCG‑TTTCTCAC‑TTTCTGTGGCCCCAAGTCTTCACCTTCTCCCTTTa > 2:5047302/1‑76 (MQ=255) gCGTCGATCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTAcc > 1:1860641/1‑76 (MQ=255) cGATCATTGATGCTGGCAGGGGCGCTG‑TTTTTCAC‑TGCCTGTAGCCACAACTCTTCCCTTCCTCCCTTTATCGccc > 2:6214319/1‑75 (MQ=255) aTCATTGATGCTGGCAGGGGGGCTG‑TTTCTCAC‑TGCCCGTAGGCACAACTCTTTACTTCCTCCCTTTTCCGcccgg > 2:6253399/1‑73 (MQ=255) aTCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCCCTTCCTCCCTTTACCGCCAGt > 1:5550424/1‑76 (MQ=255) aTCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCCGt > 1:3424782/1‑76 (MQ=255) aTCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAAATCTTCACTTCCTCCCTTTACCGcccgg > 2:4733036/1‑73 (MQ=255) tCATTGATGCTGGCAGGGGCGCTG‑TTTCTCAC‑CGCCTGTAGCCACAACTCTTCACTTCCTCCCCTTACCCCCaccg > 2:1242470/1‑73 (MQ=255) tGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAGGCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGAtt > 2:6018862/1‑76 (MQ=255) cTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTg > 1:7033445/1‑76 (MQ=255) tGGCAGGCGCGCTG‑TTTCTCAC‑CGGCTGCGACCACAACCCTTTTCTTACTCCCATTACCGCCCGTGGATTTGCTGa > 2:6633124/1‑76 (MQ=255) ggCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGaa > 1:3308300/1‑76 (MQ=255) cgggTGCTCTT‑TTTCTCCC‑TGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGa < 1:4744858/74‑1 (MQ=255) cAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGa > 2:5182254/1‑76 (MQ=255) aGGTGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAc > 2:2626695/1‑76 (MQ=255) ggTGCGTTT‑TTTCTCCC‑TGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAcc < 1:1694136/76‑1 (MQ=255) tGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAg > 1:1782789/1‑76 (MQ=255) | GGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTG‑TTTCTCAC‑TGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAG > REL606/4180174‑4180324 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |