Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 4,407,692 | C→T | 46.3% | R325C (CGT→TGT) | hflK → | modulator for HflB protease specific for phage lambda cII repressor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 4,407,692 | 0 | C | T | 46.3% | 8.4 / 54.3 | 41 | R325C (CGT→TGT) | hflK | modulator for HflB protease specific for phage lambda cII repressor |
Reads supporting (aligned to +/- strand): ref base C (12/10); new base T (9/10); total (21/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.58e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGG > W3110S.gb/4407627‑4407759 | gTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTg < 1:4491359/71‑1 (MQ=255) gTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTg < 1:4616422/71‑1 (MQ=255) aaGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGta > 1:4916694/1‑70 (MQ=255) tGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATc < 1:1985262/71‑1 (MQ=255) ggCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCt < 1:4220870/63‑1 (MQ=255) ggCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATc > 1:367083/1‑70 (MQ=255) tGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATg > 1:2376696/1‑70 (MQ=255) tGCTAAACTTCTGCCGAAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATg > 1:2355452/1‑70 (MQ=255) gCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATgg > 1:2582902/1‑70 (MQ=255) aCTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGaaaa < 1:4197954/69‑1 (MQ=255) aCTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAg < 1:581620/71‑1 (MQ=255) tCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGtat < 1:1858003/49‑1 (MQ=255) cTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGt > 1:2302756/1‑42 (MQ=255) cTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGa > 1:1577737/1‑58 (MQ=255) cTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAgtg < 1:4967455/69‑1 (MQ=255) gAATATAAAGCCGCGCCTGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGaaaa < 1:459966/59‑1 (MQ=255) gAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTaa > 1:2012168/1‑71 (MQ=255) atatAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAaca < 1:4168900/71‑1 (MQ=255) tataAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAacac < 1:459528/71‑1 (MQ=255) taAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACAccc > 1:1117897/1‑71 (MQ=255) aaGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGtt > 1:3334778/1‑58 (MQ=255) aGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGc > 1:357152/1‑70 (MQ=255) ccGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGAt > 1:1663465/1‑43 (MQ=255) cGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTT‑GGTAACACCCGCaa < 1:4461363/63‑1 (MQ=255) gAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAg > 1:2383698/1‑64 (MQ=255) gAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGtt < 1:5071264/47‑1 (MQ=255) ttACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAg < 1:66334/39‑1 (MQ=255) ttACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTaa < 1:3146983/49‑1 (MQ=255) ttACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTaa < 1:4866725/49‑1 (MQ=255) ttACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAAc < 1:5033972/71‑1 (MQ=255) tACTCGCGAGTGTCTGTATATCGAGACGATGGaaaa > 1:4973204/1‑36 (MQ=255) tACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAacacac > 1:1078029/1‑51 (MQ=255) tACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAAc > 1:4288539/1‑70 (MQ=255) tACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAAc > 1:2580325/1‑70 (MQ=255) tACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAAc > 1:2505062/1‑70 (MQ=255) tACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAAc > 1:1153786/1‑70 (MQ=255) tACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTgg > 1:272634/1‑65 (MQ=255) aCTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAaca > 1:3017625/1‑49 (MQ=255) gAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAgg < 1:2458194/71‑1 (MQ=255) gAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAgg < 1:2452348/71‑1 (MQ=255) aGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGt < 1:3211839/38‑1 (MQ=255) | GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGG > W3110S.gb/4407627‑4407759 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |