Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,586,707 | G→A | 63.6% | intergenic (+269/‑115) | cysK → / → ptsH | cysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,586,707 | 0 | G | A | 63.6% | 42.4 / 44.9 | 74 | intergenic (+269/‑115) | cysK/ptsH | cysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr) |
Reads supporting (aligned to +/- strand): ref base G (12/15); new base A (25/22); total (37/37) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.30e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01 |
ATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCC > minE/1586650‑1586771 | aTGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTg > 1:416656/1‑69 (MQ=255) gTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTgcg > 1:3108790/1‑69 (MQ=255) gTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTgcg > 1:387580/1‑69 (MQ=255) tGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTgcgtcgcg < 1:2083597/68‑1 (MQ=255) gACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTgcgtcgcg < 1:103124/67‑1 (MQ=255) gTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCt > 1:1459101/1‑67 (MQ=255) gTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTg > 1:358674/1‑68 (MQ=255) gTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTg > 1:3074572/1‑68 (MQ=255) gTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTg > 1:1514797/1‑68 (MQ=255) gTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTg > 1:3246262/1‑68 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:328697/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:377120/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:1922685/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:3428910/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:3337771/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:2067389/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:2692853/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:3197262/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:3136994/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:2256393/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:2828849/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:2472048/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:389225/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:656873/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:577124/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:1284943/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:1342821/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:1108709/59‑1 (MQ=255) aCGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:1101493/59‑1 (MQ=255) cGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:969566/58‑1 (MQ=255) aaaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTaaa < 1:424200/68‑1 (MQ=255) aaaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTaaa < 1:1100082/68‑1 (MQ=255) aaaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTaaa < 1:423056/68‑1 (MQ=255) aaaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTaaa < 1:2212781/68‑1 (MQ=255) aaaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTaaa < 1:864842/68‑1 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCtt > 1:499536/1‑47 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:55511/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:1055180/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:3062674/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:3033407/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:871688/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:881522/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:1236706/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:1679910/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:1613176/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:1525762/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:1442605/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:2009972/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:2069607/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:212808/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:1290268/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:1216071/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:2272378/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:2367004/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:1199777/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:1156537/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:2410741/1‑69 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaa > 1:11218/1‑66 (MQ=255) aaCGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATCACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:982544/1‑69 (MQ=255) aaCGTCAGAGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTGCGTCGCGTCTGTTaaaaa > 1:1416272/1‑69 (MQ=255) aGCGGTCAACACCCGCCAGCAATGGACTGTATTgcg > 1:1900253/1‑36 (MQ=255) aGCGGTCAACACCCGCCAGCAATGGACTGTATTgcg > 1:849240/1‑36 (MQ=255) aGCGGTCAACACCCGCCAGCAATGGACTGTATTgcg > 1:318570/1‑36 (MQ=255) aGCGGTCAACACCCGCCAGCAATGGACTGTATTgcg > 1:2938310/1‑36 (MQ=255) cGGTCAACACCCGCCAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:209580/47‑1 (MQ=255) tCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAAc < 1:1967341/69‑1 (MQ=255) tCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAAc < 1:1465320/69‑1 (MQ=255) tCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAAc < 1:2188786/69‑1 (MQ=255) tCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAAc < 1:2258622/69‑1 (MQ=255) tCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAAc < 1:678065/69‑1 (MQ=255) tCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAAc < 1:2402726/69‑1 (MQ=255) tCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAAc < 1:2550003/69‑1 (MQ=255) aaCACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACaa < 1:3338921/69‑1 (MQ=255) cccGACAGCAATGGACTGTATTACGCTCTTCGTgcgtc < 1:1417323/38‑1 (MQ=38) aaTGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAg > 1:2995773/1‑68 (MQ=255) tATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGcc > 1:2557817/1‑69 (MQ=255) tATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGcc > 1:2831914/1‑69 (MQ=255) tATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGcc > 1:1115758/1‑69 (MQ=255) | ATGTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCC > minE/1586650‑1586771 |
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: ‑ |