Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,840,640 | 0 | A | G | 26.7% | 140.5 / 16.1 | 60 | I373V (ATT→GTT) | ascF | beta‑glucoside specific PTS enzyme IIBC component |
Reads supporting (aligned to +/- strand): ref base A (6/38); new base G (16/0); total (22/38) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.99e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCT > NC_000913/2840592‑2840689 | cTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAt < 1:1101305/50‑1 (MQ=255) tGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGtt > 1:6223142/1‑50 (MQ=255) tGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGtt > 1:442737/1‑50 (MQ=255) cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAttt < 1:3445969/49‑1 (MQ=255) cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAttt < 1:742752/49‑1 (MQ=255) cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAttt < 1:4481730/49‑1 (MQ=255) cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAttt < 1:3007698/49‑1 (MQ=255) cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAttt < 1:1692274/49‑1 (MQ=255) cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTc < 1:4624638/50‑1 (MQ=255) cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTc < 1:2872431/50‑1 (MQ=255) tCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGttttc > 1:5720336/1‑47 (MQ=39) tCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCg < 1:2121140/50‑1 (MQ=255) aGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGaa > 1:1178298/1‑50 (MQ=39) aGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGa > 1:3590875/1‑49 (MQ=38) aGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGa < 1:5330750/49‑1 (MQ=255) gACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGaaa > 1:2540319/1‑49 (MQ=38) ggCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGGAcccg > 1:641689/1‑48 (MQ=255) gcgcTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGc < 1:4813691/50‑1 (MQ=255) cgcTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGAACCGGCg > 1:3130177/1‑50 (MQ=255) cgcTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTCCCAACCGGCg > 1:6673120/1‑50 (MQ=255) cgcTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCg < 1:1343210/50‑1 (MQ=255) gctggctgCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGAACCcgcgg > 1:1348234/1‑49 (MQ=255) gctggctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGt < 1:2051621/50‑1 (MQ=255) gctggctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGt < 1:1138649/50‑1 (MQ=255) gctggctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGt < 1:2097451/50‑1 (MQ=255) gctggctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGt < 1:778779/50‑1 (MQ=255) ctggctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGtt < 1:1537177/50‑1 (MQ=255) ggctgCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGAAACCGCGGta > 1:318647/1‑49 (MQ=25) gctgCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGGAACGGCGTtata > 1:5933170/1‑50 (MQ=255) tgCGGCGGCATCAGCCATTATGGCGGGGGGTTCCCCAACTGCGTTATTtg > 1:1416148/1‑50 (MQ=255) cggcggCATCAGCCATTATGGCGGGGGTTTTCCAACCGGCGGCTTAtggc > 1:236646/1‑50 (MQ=255) cggcggCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTTTAtggc < 1:1737056/50‑1 (MQ=255) cggcggCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggc < 1:4954829/50‑1 (MQ=255) cggcggCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtg < 1:5587095/48‑1 (MQ=255) cggcggCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtg < 1:1263112/48‑1 (MQ=255) ggcggcATCAGCCATTATGGCGGGGGGTTTCGGACCGGGGGTTTAtgggg > 1:4871568/1‑48 (MQ=255) ggcggcATCAGCCATTATGGCGGGGATTTTCGAACCGGCGGTATAtggcg > 1:919258/1‑50 (MQ=38) ggcggcATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggcg < 1:1781321/50‑1 (MQ=255) gcggcATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggcgt < 1:2804132/50‑1 (MQ=255) cATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggcgtggcg < 1:5049714/50‑1 (MQ=255) aTCAGCCATTATGGCGGGGGGTTTCGGACCCGCGGTTTATGGCGTGGCGa > 1:4753656/1‑50 (MQ=255) aTCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGa < 1:3942698/50‑1 (MQ=255) tCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGAt < 1:4734979/50‑1 (MQ=255) aGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATcc < 1:1250585/50‑1 (MQ=255) aGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATcc < 1:2442026/50‑1 (MQ=255) aGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATcc < 1:840378/50‑1 (MQ=255) ccATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGc < 1:4741565/50‑1 (MQ=255) tGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaaa < 1:3013396/50‑1 (MQ=255) tGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaaa > 1:4603795/1‑50 (MQ=255) tGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaaa > 1:4762595/1‑50 (MQ=255) tGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaaa < 1:2101041/50‑1 (MQ=255) tGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaaa > 1:1940163/1‑50 (MQ=255) tGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaa < 1:5158075/49‑1 (MQ=255) tGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaa < 1:2760646/49‑1 (MQ=255) gCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACg > 1:3156817/1‑50 (MQ=255) ggggATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTc < 1:6390198/50‑1 (MQ=255) gggATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTcc > 1:888308/1‑50 (MQ=255) ggATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCg < 1:5124294/50‑1 (MQ=255) ggATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCg < 1:1005512/50‑1 (MQ=255) gATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGc < 1:4937278/50‑1 (MQ=255) aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt > 1:6669204/1‑50 (MQ=255) aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt < 1:2797897/50‑1 (MQ=255) aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt < 1:2799924/50‑1 (MQ=255) | CTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCT > NC_000913/2840592‑2840689 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |