Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,958,932 | 0 | G | C | 37.5% | 20.6 / 11.5 | 16 | A112G (GCT→GGT) | cutC | protein CutC |
Reads supporting (aligned to +/- strand): ref base G (10/0); new base C (0/6); total (10/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.25e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTATATAAAGGGTTAGCGCACATATCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCC > NC_000913/1958869‑1958990 | gtATATAAAGGGTTAGCGCACATATCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGgcagg > 1:2933397/1‑64 (MQ=255) taAAGGGTTAGCGCACATATCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGGAGCCa > 1:1019417/1‑65 (MQ=255) taAAGGGTTAGCGCACATATCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGGAGCCa > 1:1567175/1‑65 (MQ=255) taAAGGGTTAGCGCACATATCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCa > 1:2635748/1‑65 (MQ=255) aaGGGTTAGCGCACATATCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCatt > 1:1903195/1‑65 (MQ=255) gggTTAGCGCACATATCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCattat > 1:2711570/1‑65 (MQ=255) ggTTAGCGCACATATCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCGGCCattatt > 1:419295/1‑65 (MQ=255) ggTTAGCGCACATATCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCattatt > 1:821520/1‑65 (MQ=255) ggTTAGCGCACATATCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCactatt > 1:1203429/1‑65 (MQ=255) ggTTAGCGCACATATCGAAGGCGCGATGAAAAGCCACTGCCAGCGGACCGGCAGCAGCCcttatt > 1:3814804/1‑65 (MQ=255) cgcACATATCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGGAGCCCTTATTTTTTcc > 1:1580444/1‑65 (MQ=255) cgcACATATCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGGAGCCATTATTTTTTcc > 1:2668727/1‑65 (MQ=255) cgcACATATCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTcc > 1:3935549/1‑65 (MQ=255) tatCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTTTTTTTTCCATTCGt > 1:1520904/1‑65 (MQ=255) tatCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGt > 1:904138/1‑65 (MQ=255) atCGAGGGCGCGATAAAAAGCCACGGCCAGCGGCCCGCCAGCAGCCATTATTTTTTCCATTCGTg < 1:448734/65‑1 (MQ=255) atCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTg > 1:18428/1‑65 (MQ=255) atCAAAGGCGGGATAAAAAGCCACTCCCAGCGGCCCGGCCCCAGCCATTATTTTTTCCATTCGTg < 1:2087033/65‑1 (MQ=255) ggaGGCGCAATAAAACGTCTCCGCCGCCGGACCGGCACCAGCCATTATTTTTTCCATTCGTGGCa < 1:3790875/63‑1 (MQ=255) aaaGGCGCGATGAAAAGCCCCGCCGAGCGGCCCGGCGCCAGCCATTATTTTTTCCATTCGTGGCa < 1:1067314/64‑1 (MQ=255) aaGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCGGCACTTCTTTTTTCCATTTGTgggcc > 1:1208412/1‑62 (MQ=255) aaGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCCGCCATTATTTTTTCCATTCGTGGCaa > 1:763755/1‑64 (MQ=255) aaGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTCTTTTTTCCATTTGTGGAat > 1:2457859/1‑65 (MQ=255) aaGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCCTTCCTGGCat > 1:2522348/1‑65 (MQ=255) aaGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCat > 1:2944652/1‑65 (MQ=255) aaGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCat > 1:3186365/1‑65 (MQ=255) aaGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGGGGGat > 1:402643/1‑65 (MQ=255) aaGGCGCGAAAAAAAGTCCCGCCCGGCGGACCGGCACCAGCCATTATTTTTTCCATTCGTGGCat < 1:1953676/65‑1 (MQ=255) ggCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGTAGCCCTTCTTTTTTTCATTCCTGGCatat > 1:3896785/1‑65 (MQ=255) ggCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCCTTATTTTTTTCCTTCGTGGCatat > 1:2011916/1‑65 (MQ=255) ggCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCatat > 1:199231/1‑65 (MQ=255) ggCGCGATAAAAAGCCCCGCCCAGCGGCCCGGCCCCAGCCATTATTTTTTCCATTCGTGGCatat < 1:1799969/65‑1 (MQ=255) ggCGCGAAAAAAACTCCCCCCCAGCGGCCCGCCCCCAGCCATTATTTTTTCCATTCGTGGCatat < 1:3528590/65‑1 (MQ=255) gcggcgATGAAAAGACGCCCCCCGCCGGCCGGCCCCAGCCATTATTTTTTCCATTCGTGGCatat < 1:3137200/62‑1 (MQ=255) cgATGAAAAGTCTCTGCCAGCGACCCGGCCCCAGCCATTATTTTTTCCATTCGTGGCATATCGAc < 1:438828/65‑1 (MQ=255) aaaaTTCCCTCCCCGCCGCCCGGCCGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTcc < 1:829266/65‑1 (MQ=255) aaaaGCCCCTCCCGGGGGCCCGGCCCCAGCCATTATTTTTTCCATTCGTCGCATATCGACATTcc < 1:3115181/65‑1 (MQ=255) ctgTCACTGCCGGCGGCCGGCCCCCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCg < 1:169683/63‑1 (MQ=255) aGTCACTGCCAGCGGACCGGCAGCAGCCCTTCTTTTTTCCAATGGTGGCGTTTCTGCATTCCCGt > 1:1380806/1‑65 (MQ=255) ccAGCGGACCGGCAGCAGCCCTTATTTTTTCCATTTGTGGCATATCGACATTCCCGTCAACACCg > 1:803745/1‑65 (MQ=255) ccGGCACCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGcc < 1:516720/65‑1 (MQ=255) ccGGCACCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGcc < 1:1004653/65‑1 (MQ=255) | GTATATAAAGGGTTAGCGCACATATCGAAGGCGCGATGAAAAGTCACTGCCAGCGGACCGGCAGCAGCCATTATTTTTTCCATTCGTGGCATATCGACATTCCCGTCAACATCGAGAACGCC > NC_000913/1958869‑1958990 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |