Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 476,746 | 0 | T | G | 36.5% | 166.6 / 18.7 | 85 | Q6H (CAA→CAC) | fsr | predicted fosmidomycin efflux system |
Reads supporting (aligned to +/- strand): ref base T (17/37); new base G (31/0); total (48/37) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.34e-11 | |||||||||||
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 138 total aligned reads displayed.
TGATAGCACCTAAAATACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAC > REL606/476673‑476821 | tGATAGCACCTAAAATACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtg > 2:4698315/1‑76 (MQ=255) tGATAGCACCTAAAATACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtg > 2:5466864/1‑76 (MQ=255) tGATAGCACCTAAAATACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtg > 1:4414734/1‑76 (MQ=255) tGATAGCACCTAAAATACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtg > 1:6384917/1‑76 (MQ=255) gATAGCACCTAAAATACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtgt > 2:7074833/1‑76 (MQ=255) gATAGCACCTAAAATACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtgt > 2:779547/1‑76 (MQ=255) gATAGCACCTAAAATACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtgt > 1:5139708/1‑76 (MQ=255) gATAGCACCTAAAATACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtgt > 1:6498427/1‑76 (MQ=255) aaGCCCCCAAAATTCCAAAACGATGTCGGGGCGTGGGTTGTCGACGCCCCCCCCCCCGCCACAGGCTGGGgttgtt < 1:4552148/75‑1 (MQ=255) aCCTAAAATACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCAcac > 2:6418369/1‑74 (MQ=255) aCCTAAAATACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTc > 2:467761/1‑76 (MQ=255) cTAAAATACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGGGTTCACTCAt > 2:1129946/1‑76 (MQ=255) tAAAATACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCCCTCAtt > 1:4974463/1‑76 (MQ=255) aaaaTACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGGTTACTCATTg > 2:7233023/1‑76 (MQ=255) aaaTACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGc > 2:5993780/1‑76 (MQ=255) aCC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGGCATa > 1:5491104/1‑76 (MQ=255) aCC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCCTa > 1:6682457/1‑76 (MQ=255) c‑AAGCGTTGTTGGGGCCTGGGTCGTCCCCCCCCCCCCCCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAgg < 1:1128308/76‑1 (MQ=255) c‑AAGCGTTGTGCGGCCTTTGGTCGTTGCCCCCCCCCCCCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAgg < 2:7376893/76‑1 (MQ=255) aaaCGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCTCTCATTGCCCTTGGa > 1:4828020/1‑76 (MQ=255) aaaCGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTTATTGCCATAGGa > 1:4925403/1‑76 (MQ=255) aaaCGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCCTTGCCATAggg > 1:34371/1‑75 (MQ=255) aaaCGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGa > 2:5574568/1‑76 (MQ=255) aaaCGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGGTTACTCATTGCCCTAGGa > 1:6317681/1‑76 (MQ=255) aCGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTTACTCCTTGCCATAGGaaa > 1:5256804/1‑76 (MQ=255) aCGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGgaa > 1:6269899/1‑76 (MQ=255) gATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGGAAGt > 1:1486753/1‑76 (MQ=255) aTGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGGTTACTCATTGCCATAGGAAAggc > 1:3372215/1‑74 (MQ=255) gTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCAc > 2:2053093/1‑76 (MQ=255) tCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCCTTGGCATAGGAAAAGCACtt > 2:2166840/1‑76 (MQ=255) gggCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCCTTGGAAAGATACtttt > 2:7282269/1‑76 (MQ=255) gggCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACtttt > 1:2925348/1‑76 (MQ=255) ggCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACttttt > 2:1002622/1‑76 (MQ=255) gCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAATCACtttttt > 2:2069135/1‑75 (MQ=255) ccTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCCTTGCCATAGGGAAGTCACTTTTTCa > 1:2544018/1‑76 (MQ=255) cTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTTACTCATTGGCCCTGGTAAGTTACttttttcg > 2:6234156/1‑73 (MQ=255) cTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAg > 1:3240606/1‑76 (MQ=255) ttGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGGCACTTTTTCAgg > 1:5601119/1‑76 (MQ=255) ttGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAgg > 1:6634879/1‑76 (MQ=255) ggTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCCAAGGAAAGTCACTTTTTCAGGGt > 1:3380193/1‑76 (MQ=255) ggTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTTAGGGt > 2:4802603/1‑76 (MQ=255) ggTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCCGGGt > 1:1565680/1‑76 (MQ=255) ggTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAAGAAAGTCACTTTTTCAgggg > 2:2409501/1‑75 (MQ=255) gTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTTACTCATTTCCCTAGGGAAGTCACTTTTTTAggggt > 1:4301952/1‑74 (MQ=255) cGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGc > 1:3017558/1‑76 (MQ=255) cGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTTACTCATTGCCATAGGAAAGTCACTTTTTTAGGGGTGc > 2:4062280/1‑76 (MQ=255) gACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACCTTTTTAGGGGTGCGATg > 1:2664297/1‑76 (MQ=255) gACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGGCACTTTTTCAGGGGTGCGATg > 2:1419464/1‑76 (MQ=255) gACGCAGCCGCGCCCGCCACAGGCTGGGGATGTTCACACATTGCCCTAGGAAAGGCACTTTTTCAGGGTTGCGATg > 1:7106460/1‑76 (MQ=255) cGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTa > 1:6178271/1‑76 (MQ=255) gCAGCCGCGCCCGCCACAGGCTGGGGGTGTTTACTCATTGCCATAGGAAAGTCACTTTTTTAGGGTTGCGATGTaa > 2:4608764/1‑76 (MQ=255) cAGCCGCGCCCGCCACAGGCTGGGGGTGTTCCCCCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTaaa > 1:4534483/1‑76 (MQ=255) cAGCCGCGCACGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTaaa > 2:4184298/1‑76 (MQ=255) cgCCCCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAg < 2:4937905/75‑1 (MQ=255) aGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAg > 1:1439643/1‑76 (MQ=255) ccccccgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATg < 2:5129202/73‑1 (MQ=255) cgcccgccACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGa > 2:5027097/1‑76 (MQ=255) cccccgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGa < 2:1006697/74‑1 (MQ=255) cccccgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGa < 1:4896182/74‑1 (MQ=255) ccccgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGAt < 2:26837/75‑1 (MQ=255) ccccgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGAt < 2:3019207/75‑1 (MQ=255) ccccgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGAt < 2:3413841/75‑1 (MQ=255) cgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCtt < 1:2797414/76‑1 (MQ=255) cgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCtt < 1:4491962/76‑1 (MQ=255) cgccACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCtt > 2:3193338/1‑76 (MQ=255) gccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTa < 1:2897613/76‑1 (MQ=255) ccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAt < 1:2647261/76‑1 (MQ=255) ccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAt < 2:6423772/76‑1 (MQ=255) ccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAt < 2:7050560/76‑1 (MQ=255) ccACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCCCTTTTTCAGGGGTGCGATGTAAAGAATGATCTTAt > 2:321579/1‑76 (MQ=255) cacaGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAtt < 2:17703/76‑1 (MQ=255) acaGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAttt < 2:4064088/76‑1 (MQ=255) acaGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAttt < 2:6497704/76‑1 (MQ=255) acaGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTTAGGGTTGCGGTGTAAAGAATGATCTTAttt > 2:4207522/1‑76 (MQ=255) ggCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTtgtg < 1:6235118/76‑1 (MQ=255) gCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGa < 2:912617/76‑1 (MQ=255) gCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGa > 1:3229722/1‑76 (MQ=255) gCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGa < 1:3951949/76‑1 (MQ=255) cTGGGGTTGTTCACTCATTGCCCTAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGat > 2:4712041/1‑76 (MQ=255) cTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGat < 1:1681266/76‑1 (MQ=255) cTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGat > 1:4097348/1‑76 (MQ=255) tGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGatt > 1:1411302/1‑76 (MQ=255) tGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGatt < 1:2578669/76‑1 (MQ=255) tGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGatt < 2:276443/76‑1 (MQ=255) ggggTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGatta > 2:3676144/1‑76 (MQ=255) gtgttgttCTCTCATTGCCATATGAAAGTCGCTATTTCAGGGTTGCGATGAAAAGAATGATGTTATTTGTGtttat < 2:3145916/74‑1 (MQ=255) gggTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGattat < 2:694694/76‑1 (MQ=255) gggTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGattat < 2:3145915/76‑1 (MQ=255) ggTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGattatt < 1:5441428/76‑1 (MQ=255) ggTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGattatt > 2:6997132/1‑76 (MQ=255) ggTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGattatt < 2:2793863/76‑1 (MQ=255) gttgttCACTCCTTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGAttatta > 2:609600/1‑76 (MQ=255) gttgttCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGAttatta > 2:2207995/1‑76 (MQ=255) gttgttCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGAttatta < 2:2174883/76‑1 (MQ=255) gttgttCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGAttatta < 1:551585/76‑1 (MQ=255) gttgttCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGAttatta < 1:4559928/76‑1 (MQ=255) gttgttCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGAttatta < 2:3428300/76‑1 (MQ=255) gttgtgCACTCATTGCTATACGAAAGTGACGTTGTCAGGGTTGCGATTTAAATAATGATCTTATTTGTGAttatta < 2:2174886/76‑1 (MQ=255) ttgttCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAc > 2:4057207/1‑76 (MQ=255) ttgttCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAc < 2:5741420/76‑1 (MQ=255) | TGATAGCACCTAAAATACC‑AAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAC > REL606/476673‑476821 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |