Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,840,640 | 0 | A | G | 21.3% | 193.6 / 17.7 | 75 | I373V (ATT→GTT) | ascF | beta‑glucoside specific PTS enzyme IIBC component |
Reads supporting (aligned to +/- strand): ref base A (6/53); new base G (16/0); total (22/53) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.73e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCT > NC_000913/2840591‑2840689 | aCTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑a < 1:3644895/50‑1 (MQ=255) cTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑At < 1:5414897/50‑1 (MQ=255) cTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑At < 1:6129793/50‑1 (MQ=255) cTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑At < 1:3785668/50‑1 (MQ=255) tGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑Gtt > 1:3450140/1‑50 (MQ=255) cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑Attt < 1:6206222/49‑1 (MQ=255) cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑Att < 1:4109163/48‑1 (MQ=255) cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑Att < 1:5270294/48‑1 (MQ=255) cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑Att < 1:3765976/48‑1 (MQ=255) cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑ATTTc < 1:5334505/50‑1 (MQ=255) gTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGTG‑Agtttc > 1:1220311/1‑48 (MQ=38) gTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑ATTTcc < 1:1973544/50‑1 (MQ=255) gTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑ATTTcc < 1:850478/50‑1 (MQ=255) gTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑ATTTcc < 1:5040693/50‑1 (MQ=255) gTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑ATTTcc < 1:323032/50‑1 (MQ=255) gTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑ATTTcc < 1:1305512/50‑1 (MQ=255) tCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑ATTTCCg < 1:2984383/50‑1 (MQ=255) cAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑GTTTCCg > 1:2608539/1‑49 (MQ=255) aGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑GTTTTCgca > 1:6143738/1‑48 (MQ=38) aGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑GTTTTCgc > 1:1495234/1‑48 (MQ=38) aGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑AGTTTCgca > 1:3755433/1‑48 (MQ=38) gACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑ATTTCCGAAc < 1:2546261/50‑1 (MQ=255) gACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑ATTTCCGAAc < 1:5506347/50‑1 (MQ=255) gACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑ATTTCCGAAc < 1:7432630/50‑1 (MQ=255) aCGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑GTTTTCGAAcc > 1:3514745/1‑50 (MQ=39) aCGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑ATTTCCGAAcc < 1:3766180/50‑1 (MQ=255) aCGGCGCTCGCTGCGGCGGCATCGGCCATTATGGCGGGG‑ATTTCCGAAcc < 1:6158750/50‑1 (MQ=38) cGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑GTTTTCGCAccc > 1:4978730/1‑49 (MQ=255) cGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑GTTTCCGCAccc > 1:5968394/1‑49 (MQ=255) cGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑ATTTCCGAACCg < 1:5994525/50‑1 (MQ=255) ggCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAATTTCGGAAcgg > 1:2159288/1‑48 (MQ=37) cgcTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCg < 1:4819128/50‑1 (MQ=255) gctggctgCGGCGGCATCAGCCATTATGGCGGGG‑GTTTTCGAACCGGCgg > 1:6727597/1‑49 (MQ=255) gctggctgCGGCGGCATCAGCCATTATGGCGGGG‑GTTTTCGAACCGGCgg > 1:4079482/1‑49 (MQ=255) gctggctgCGGCGGCATCAGCCATTATGGCGGGG‑GGTTTCCAACCcgcgg > 1:6525648/1‑49 (MQ=255) gctggctgCGGCGGCATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGt < 1:6440590/50‑1 (MQ=255) gctggctgCGGCGGCATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGt < 1:3931530/50‑1 (MQ=255) ggctgCGGCGGCATCAGCCATTATGGCGGGG‑GTTTCCGAACCCGCGGta > 1:2799051/1‑49 (MQ=38) gctgCGGCGGCATCAGCCATTATGGCGGGG‑GTTTTCCGACCGGCGTtata > 1:5931046/1‑50 (MQ=255) gctgCGGCGGCATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTtata < 1:2946972/50‑1 (MQ=255) gctgCGGCGGCATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTtat < 1:6789954/49‑1 (MQ=255) gctgCGGCGGCATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTta < 1:6884498/48‑1 (MQ=255) ctgCGGCGGCATCAGCCATTATGGCGGGG‑GTTTTCCAACCCGCGTTttat > 1:5228994/1‑50 (MQ=255) ctgCGGCGGCATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTatat < 1:6175956/50‑1 (MQ=255) tgCGGCGGCATCAGCCATTATGGCGGGG‑GTTTTCGAAACGGCGGTttatt > 1:2763173/1‑49 (MQ=255) cggcggCATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATAtgg < 1:3992407/49‑1 (MQ=255) cggcggCATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATAtgg < 1:147284/49‑1 (MQ=255) ggcggcATCAGCCATTATGGCGGGG‑GTTTTCGAACCCGCGGTTTAtggcg > 1:6621193/1‑50 (MQ=255) ggcggcATCAGCCATTATGGCGGGG‑GTTTTCCGACCGGGGGTTTTtggcg > 1:6938705/1‑50 (MQ=255) ggcggcATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATAtggcg < 1:7295645/50‑1 (MQ=255) gcggcATCAGCCATTATGGCGGGG‑ATTTTCGCAACGGCGGTATAttgcgg > 1:6079926/1‑49 (MQ=25) cggcATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATAtggcgtg < 1:1857643/50‑1 (MQ=255) cggcATCAGCCATTATGGCGGGG‑AGTTTCGAAACGGCGGTATAtggcggg > 1:3493635/1‑48 (MQ=255) ggcATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATAtggcgtgg < 1:1772249/50‑1 (MQ=255) gcATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATAtggcgtggc < 1:6418500/50‑1 (MQ=255) gcATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATAtggcgtggc < 1:3433263/50‑1 (MQ=255) cATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATAtggcgtggcg < 1:1977181/50‑1 (MQ=255) cATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATAtggcgtggcg < 1:6190340/50‑1 (MQ=255) tCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGAt < 1:7273943/50‑1 (MQ=255) aGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATcc < 1:2901980/50‑1 (MQ=255) aGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATcc < 1:4001160/50‑1 (MQ=255) aGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATc < 1:2577869/49‑1 (MQ=255) aGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATc < 1:2555194/49‑1 (MQ=255) aGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATc < 1:2094978/49‑1 (MQ=255) ccATTATGGCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGc < 1:1272372/50‑1 (MQ=255) ccATTATGGCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGc < 1:4415107/50‑1 (MQ=255) aTTATGGCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCt < 1:3717610/50‑1 (MQ=255) ttATGGCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTg < 1:955974/50‑1 (MQ=255) tGGCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaa < 1:3404068/49‑1 (MQ=255) ggCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAAc > 1:6058497/1‑50 (MQ=255) gCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAAc < 1:5218352/49‑1 (MQ=255) gCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAAc < 1:1295146/49‑1 (MQ=255) gggg‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTc < 1:1161426/50‑1 (MQ=255) ggg‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTcc < 1:3760122/50‑1 (MQ=255) ggg‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTcc < 1:4572383/50‑1 (MQ=255) g‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGc < 1:169611/50‑1 (MQ=255) aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt > 1:6408167/1‑50 (MQ=255) aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt > 1:2447789/1‑50 (MQ=255) aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt < 1:2658921/50‑1 (MQ=255) aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt > 1:2963063/1‑50 (MQ=255) aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt < 1:3327128/50‑1 (MQ=255) aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt < 1:819463/50‑1 (MQ=255) | ACTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGG‑ATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCT > NC_000913/2840591‑2840689 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |