Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,858,656 | T→C | 100% | G310G (GGT→GGC) | amiB → | N‑acetylmuramoyl‑l‑alanine amidase II |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,858,656 | 0 | T | C | 82.9% | 33.0 / ‑0.1 | 23 | G310G (GGT→GGC) | amiB | N‑acetylmuramoyl‑l‑alanine amidase II |
Reads supporting (aligned to +/- strand): ref base T (1/3); new base C (0/19); total (1/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.74e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCT > minE/2858593‑2858696 | tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:2821696/69‑1 (MQ=255) tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:572866/69‑1 (MQ=255) tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:523807/69‑1 (MQ=255) tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:696578/69‑1 (MQ=255) tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:719865/69‑1 (MQ=255) tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:3280250/69‑1 (MQ=255) tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:3118653/69‑1 (MQ=255) tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:3038242/69‑1 (MQ=255) tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:3015509/69‑1 (MQ=255) tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:2465831/69‑1 (MQ=255) tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:730284/69‑1 (MQ=255) tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:201095/69‑1 (MQ=255) tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:738665/69‑1 (MQ=255) tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:1611630/69‑1 (MQ=255) tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:1228625/69‑1 (MQ=255) tCGCGCCAACAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:118169/69‑1 (MQ=255) aaCAGCGAAATGGCAAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:2068909/62‑1 (MQ=255) ggctggAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:860674/43‑1 (MQ=37) ggctggAGCAGCACGAGAAACAGTCGGAGCTACTGGGCGGAGc < 1:1663768/43‑1 (MQ=37) tggAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCt < 1:2227231/68‑1 (MQ=255) tggAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCt < 1:1192971/68‑1 (MQ=255) tggAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGGGATGTGCTGGCGAACAGTCAGTCt < 1:399418/68‑1 (MQ=255) gctgctGGGTGGGGCTGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCt > 1:368973/1‑50 (MQ=255) | TCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGAACAGTCAGTCTGACCCCT > minE/2858593‑2858696 |
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: ‑ |