Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 1,191,688 | C→T | 22.2% | T370T (ACC→ACT) | icd → | e14 prophage; isocitrate dehydrogenase, specific for NADP+ |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 1,191,688 | 0 | C | T | 22.2% | 1027.3 / 288.9 | 443 | T370T (ACC→ACT) | icd | e14 prophage; isocitrate dehydrogenase, specific for NADP+ |
Reads supporting (aligned to +/- strand): ref base C (200/144); new base T (48/50); total (249/194) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.33e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01 |
Only 100 of 457 total aligned reads displayed.
ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACT > CP009273/1191624‑1191755 | aCAAAGTAAATCCTGGCTCTATTATTCTTTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGaa > 1:27890539/1‑68 (MQ=255) aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGaa > 1:10660043/1‑68 (MQ=255) aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGaa > 1:21781125/1‑68 (MQ=255) aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGaa > 1:27634177/1‑68 (MQ=255) aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATAATGCTGCGCCATATGGGTTGGACTGaa > 1:30522597/1‑68 (MQ=37) aaGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCg < 1:3045746/68‑1 (MQ=255) gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg > 1:34369573/1‑66 (MQ=255) gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg > 1:26562594/1‑66 (MQ=255) gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg > 1:9247988/1‑66 (MQ=255) gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg > 1:9123477/1‑66 (MQ=255) gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg > 1:22593318/1‑66 (MQ=255) aaTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGa < 1:14912263/68‑1 (MQ=255) aTCCTGTATCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGAc < 1:17593443/68‑1 (MQ=255) aTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGAc < 1:36800342/68‑1 (MQ=255) tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGC‑TGa > 1:1319288/1‑66 (MQ=255) tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGC‑TGa > 1:32923301/1‑66 (MQ=255) tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACt > 1:37421567/1‑68 (MQ=255) tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACt < 1:5310973/68‑1 (MQ=255) tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACt > 1:7746350/1‑68 (MQ=255) tCCTGGCTCAATTATTCTCTCCGCTGAGATGATGCAGCGCCATATGGTTTGGACTGAAGCGGC‑TGAcc > 1:4285871/1‑67 (MQ=21) ccTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGC‑TGa > 1:35465346/1‑65 (MQ=255) ccTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACtt > 1:1609980/1‑68 (MQ=255) cTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGC‑TGAcctg > 1:7244260/1‑65 (MQ=21) cTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTa > 1:30320740/1‑68 (MQ=255) ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTg > 1:37726194/1‑68 (MQ=255) ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTg > 1:17055369/1‑68 (MQ=255) tctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGC‑TGACCTGATTg < 1:23005192/67‑1 (MQ=17) tctATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGC‑TGACCTGATTGt < 1:24712715/68‑1 (MQ=17) tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGt < 1:11889752/68‑1 (MQ=255) tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGGTGGACCGAAGCGGC‑TGACTTAATTGt < 1:37990376/68‑1 (MQ=255) tctATTATTCTCACCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGCAGCGGC‑TGACTTAATTGt < 1:8258136/68‑1 (MQ=255) ttattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTaaa < 1:10952025/68‑1 (MQ=255) ttattCTCGCCGCTGAGATGATGCTGCGCCACATGGGTTCGACCGAAGCGGA‑TGAATTAATTGTTaaa > 1:27156090/1‑68 (MQ=255) attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGC‑TGACCTGATTGTTAAAgg < 1:13553184/68‑1 (MQ=14) attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGC‑TGACCTGATTGTTAAAgg < 1:2359402/68‑1 (MQ=14) attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGC‑TGACCTGATTGTTAAAgg < 1:27367304/68‑1 (MQ=14) attCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGC‑TGACCTGATTGTTAAAgg < 1:29530465/68‑1 (MQ=14) attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAgg < 1:14151463/68‑1 (MQ=255) attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAgg < 1:2133593/68‑1 (MQ=255) attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAgg < 1:39174069/68‑1 (MQ=255) attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAgg < 1:13497185/68‑1 (MQ=255) attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAgg < 1:3944411/68‑1 (MQ=255) attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAgg < 1:34953933/68‑1 (MQ=255) attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAgg < 1:31366841/68‑1 (MQ=255) attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAgg < 1:16843348/68‑1 (MQ=255) attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAgg < 1:5398667/68‑1 (MQ=255) attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAgg < 1:11441118/68‑1 (MQ=255) attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAgg < 1:5849983/68‑1 (MQ=255) attCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAgg < 1:6590777/68‑1 (MQ=255) ttCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGt < 1:29984540/68‑1 (MQ=255) tctctcCGCTGAGGTGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTa > 1:38299685/1‑68 (MQ=255) tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGC‑TGACCTGATTGTTAAAGGTa > 1:28915819/1‑68 (MQ=11) tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGC‑TGACCTGATTGTTAAAGGTa < 1:36034261/68‑1 (MQ=11) tctctcCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGC‑TGACCTGATTGTTAAAGGTa > 1:39546708/1‑68 (MQ=11) tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTa > 1:10021107/1‑68 (MQ=255) tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTa > 1:25413835/1‑68 (MQ=255) tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTa > 1:24167966/1‑68 (MQ=255) tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTa > 1:4865150/1‑68 (MQ=255) tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTa > 1:1471927/1‑68 (MQ=255) tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTa > 1:9651402/1‑68 (MQ=255) tctctcCGCTGAGATGATGCTGCGACACATGGGTTGGACCGAAGCGTC‑TGACTTAATTGTTAAAGGTa > 1:11123273/1‑68 (MQ=25) tctcCGCTGAGATGGTGCTGCGCCACATGGGTGGGACCGAAGCGGC‑TGGGTTAATTGTTAAAGGTATg > 1:28559491/1‑68 (MQ=21) tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATg > 1:19372987/1‑68 (MQ=255) tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATg > 1:37814002/1‑68 (MQ=255) tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATg > 1:15114413/1‑68 (MQ=255) gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAgg < 1:33552556/68‑1 (MQ=255) gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAg > 1:1901114/1‑67 (MQ=255) gatgatGCTGCGCCACATGGTTTGGA‑CGAAACGGCGTGCCTTAA‑TGTTAAAGGTATTGAAGGCCCaa > 1:1631780/1‑67 (MQ=9) gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCaa > 1:17960990/1‑68 (MQ=255) ctaatgCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCaa > 1:20203228/4‑67 (MQ=255) tgctgcGCCACATGGGTTGGCCCGCAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc > 1:30858457/1‑68 (MQ=255) tgctgcGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc > 1:6211303/1‑68 (MQ=255) tgctgcGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc > 1:38527340/1‑68 (MQ=255) gctgcGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc > 1:23900226/1‑67 (MQ=255) gctgcGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc > 1:15723313/1‑67 (MQ=255) gctgcGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc > 1:22743775/1‑67 (MQ=255) gctgcGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc > 1:569288/1‑67 (MQ=255) gctgcGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTCTTAAAGGTATGGAACGCGCAATCAAc > 1:33199377/1‑67 (MQ=255) ctgcGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc > 1:28191067/1‑68 (MQ=255) ctgcGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc > 1:6975861/1‑68 (MQ=255) ctgcGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgc > 1:33816752/1‑68 (MQ=255) tgcGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcgcg > 1:12392528/1‑68 (MQ=255) gcCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaa > 1:14039259/1‑68 (MQ=255) cacaTGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa < 1:10470622/67‑1 (MQ=255) cacaTGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc < 1:25938527/68‑1 (MQ=255) cacaTGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAc < 1:31209339/68‑1 (MQ=255) acaTGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAAcc < 1:7490749/68‑1 (MQ=255) gggTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTa > 1:14416158/1‑67 (MQ=255) ggTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa > 1:18007505/1‑67 (MQ=255) gTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc < 1:21026645/67‑1 (MQ=255) gTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc < 1:23280836/67‑1 (MQ=255) gTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc < 1:34605464/67‑1 (MQ=255) gTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc < 1:3196648/67‑1 (MQ=255) ttGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt > 1:36505013/1‑68 (MQ=255) ttGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt > 1:35758967/1‑68 (MQ=255) ttGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt > 1:1270190/1‑68 (MQ=255) tGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt < 1:12286779/67‑1 (MQ=255) aCCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg > 1:19737133/1‑67 (MQ=255) aCCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg > 1:38897849/1‑67 (MQ=255) aCCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg > 1:8748822/1‑67 (MQ=255) | ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGC‑TGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACT > CP009273/1191624‑1191755 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |