Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,750,507 | C→A | 100% | S236S (TCG→TCT) | yfiF ← | predicted methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,750,507 | 0 | C | A | 75.0% | 36.7 / 13.3 | 32 | S236S (TCG→TCT) | yfiF | predicted methyltransferase |
Reads supporting (aligned to +/- strand): ref base C (5/3); new base A (24/0); total (29/3) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.13e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCAT > minE/1750443‑1750572 | tgttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGAc > 1:2590180/1‑68 (MQ=255) tgttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGAc > 1:818905/1‑68 (MQ=255) tgttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGAc > 1:2953541/1‑68 (MQ=255) tgttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGAc > 1:2947002/1‑68 (MQ=255) tgttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGAc > 1:2042734/1‑68 (MQ=255) tgttgtCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGAc > 1:116104/1‑68 (MQ=255) cAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCCGACTCCAGCAGTg > 1:462531/1‑69 (MQ=255) cAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCCGACTCCAGCAGTg > 1:1115469/1‑69 (MQ=255) ggCTGAACGTGCTCTGCTCCGCCTTCTGCGGTACGGATAGCCGCCCCCGACTCCAGCAGTTCCGCATcc < 1:190356/69‑1 (MQ=255) ggCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCCGACTCCAGCAGTGCCGCATcc < 1:2402643/69‑1 (MQ=255) cGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCCGACTCCAGCAGTGCCGCATCCTGcacc > 1:2968362/1‑68 (MQ=255) cGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCCGACTCCAGCAGTGCCGCATCCTGcacc > 1:1481297/1‑68 (MQ=255) tGCGCCGCCTTCTGCGGTCCGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATGCTGcaccaca > 1:36960/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGa > 1:1771840/1‑36 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:3177464/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:748464/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:719746/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:180634/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:417064/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:3395566/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:3381612/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:3313638/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:3288344/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:2197508/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:3113243/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:2984996/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:1766808/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:144891/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:2574432/1‑63 (MQ=255) tGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCAGACTCCAGCAGTGCTGCATCCTGcaccaca > 1:2486030/1‑63 (MQ=255) ccgccTTCTGCGGTACGGATAGCCGCCCCCGACTCCAGCAGTGCCGCATCCTTCACCACGACACCttt < 1:2484136/68‑1 (MQ=255) aCGGATAGCCGCCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGa > 1:1335709/1‑62 (MQ=255) ccccGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCATCTGCGcat < 1:2521435/69‑1 (MQ=255) | TGTTGTCGCCAGTAATCGGCTGAACGTGCTCTGCGCCGCCTTCTGCGGTACGGATAGCCGCCCCCGACTCCAGCAGTGCCGCATCCTGCACCACGACACCTTTCACGCCGAAGTGCGCGCAGCTGCGCAT > minE/1750443‑1750572 |
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: ‑ |