Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,549,418 | 0 | T | G | 34.7% | 146.5 / 14.3 | 72 | T6P (ACC→CCC) | yjiO | multidrug efflux system protein |
Reads supporting (aligned to +/- strand): ref base T (13/34); new base G (25/0); total (38/34) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.85e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGG > REL606/4549368‑4549468 | gTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGt < 1:3319502/51‑1 (MQ=255) cATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑aa < 1:2174217/51‑1 (MQ=255) cATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑aa < 1:4614902/51‑1 (MQ=255) aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑aaa < 1:2696128/51‑1 (MQ=255) tACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑aaaa < 1:5563275/51‑1 (MQ=255) aCAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑aaaaa < 1:6385968/51‑1 (MQ=255) aCAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaa > 1:5466814/1‑51 (MQ=255) cAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaa > 1:1720385/1‑51 (MQ=255) aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑aaaaaaa > 1:5093635/1‑51 (MQ=255) aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:3443273/1‑51 (MQ=255) aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:2962714/1‑51 (MQ=255) aaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAAc < 1:1798486/51‑1 (MQ=255) aaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAAc < 1:84668/51‑1 (MQ=255) aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGt < 1:5492839/51‑1 (MQ=255) tCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:2488217/1‑47 (MQ=255) tCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:7267912/1‑47 (MQ=255) tCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:2856615/1‑47 (MQ=255) cAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTgg < 1:3863358/51‑1 (MQ=255) aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:6832533/1‑44 (MQ=255) aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:1710024/1‑44 (MQ=255) cGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCAt < 1:3105810/51‑1 (MQ=255) gCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCATa < 1:1004771/51‑1 (MQ=255) ccATCTGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCATaa < 1:3881675/51‑1 (MQ=255) ccATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:1538827/1‑41 (MQ=255) ccATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:2516832/1‑41 (MQ=255) cATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCATaac < 1:6775235/51‑1 (MQ=255) cATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:91211/1‑40 (MQ=255) cATCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:672179/1‑40 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:1216695/1‑39 (MQ=255) aTCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:6738319/1‑39 (MQ=255) tCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:1064637/1‑38 (MQ=255) tCGGGAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAAAGTGGCATaacaa > 1:6361609/1‑51 (MQ=38) cGGGAAAAACAGCGTGGCGGCATGGTGGGT‑AAAAAAACGTGGCATaacaaa > 1:911138/1‑51 (MQ=255) cGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:184220/1‑37 (MQ=255) cGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:1115612/1‑37 (MQ=255) cGGGAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaa > 1:3422124/1‑36 (MQ=255) gggAAAAACAGCGTGGCGGCATGGCGGGT‑aaaaaaa > 1:897518/1‑36 (MQ=255) ggAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCATaacaaaca < 1:4766572/51‑1 (MQ=255) ggAAAAACAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:6443671/1‑35 (MQ=255) ggAAAAACAGCGTGGCGGCATGGCGGGG‑AAAAAAAAGTGGCATaacaaaca > 1:5054384/1‑51 (MQ=38) gAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCATAacaaacaa < 1:459518/51‑1 (MQ=255) aaaaaCAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCATAACAaacaac < 1:3851069/51‑1 (MQ=255) aaaaCAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCATAACAAACAACt < 1:1514771/51‑1 (MQ=255) aaaaCAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCATAACAAACAACt < 1:6689919/51‑1 (MQ=255) aaaaCAGCGTGGCGGCATGGCGGGG‑AAAAAAAAGTGGCATaacaaa > 1:2819037/1‑46 (MQ=38) aaaCAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTc < 1:461398/51‑1 (MQ=255) aaaCAGCGTGGCGGCATGGCGGGG‑aaaaaaa > 1:1420548/1‑31 (MQ=255) aaCAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTcc < 1:5258305/51‑1 (MQ=255) cAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTCCtt < 1:3913190/51‑1 (MQ=255) aGCGTGGCGGCATGGCGGGG‑AAAAAAAAGTGGCATAacaaacaa > 1:3420003/1‑44 (MQ=38) gTGGCGGCATGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTcc < 1:4058042/51‑1 (MQ=255) gcggcATGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGga < 1:1316278/51‑1 (MQ=255) ggcATGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGgaga < 1:820304/51‑1 (MQ=255) gcATGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGaa > 1:1143524/1‑51 (MQ=255) gcATGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGaa < 1:6777126/51‑1 (MQ=255) gcATGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGaa > 1:1600901/1‑51 (MQ=255) gcATGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGaa > 1:5625090/1‑51 (MQ=255) cATGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAAc < 1:5745538/51‑1 (MQ=255) tGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAAccc > 1:2157481/1‑51 (MQ=255) tGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAAccc < 1:1832815/51‑1 (MQ=255) tGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAAccc > 1:4292763/1‑51 (MQ=255) ggCGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAAcccc < 1:1342963/51‑1 (MQ=255) cGGGTAAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCa < 1:4885758/51‑1 (MQ=255) cGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCaa > 1:642978/1‑51 (MQ=255) cGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCaa < 1:3784268/51‑1 (MQ=255) cGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCaa < 1:294733/51‑1 (MQ=255) cGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCaa > 1:304765/1‑51 (MQ=255) cGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCa > 1:6744164/1‑50 (MQ=255) gggT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAAc < 1:478468/51‑1 (MQ=255) gggT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAAc < 1:3210508/51‑1 (MQ=255) gggT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAAc > 1:3934591/1‑51 (MQ=255) gggT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAAc > 1:3102558/1‑51 (MQ=255) ggT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACg < 1:4168691/51‑1 (MQ=255) t‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACggg < 1:3527772/51‑1 (MQ=255) t‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACggg > 1:1870632/1‑51 (MQ=255) t‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACggg > 1:2772018/1‑51 (MQ=255) | GTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGT‑AAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGG > REL606/4549368‑4549468 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |