Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 125,085 | Δ1 bp | 100% | coding (2069/2664 nt) | aceE → | pyruvate dehydrogenase E1 component |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 125,085 | 0 | T | . | 89.7% | 100.0 / 1.7 | 29 | coding (2069/2664 nt) | aceE | pyruvate dehydrogenase E1 component |
Reads supporting (aligned to +/- strand): ref base T (1/0); major base . (14/12); minor base A (0/1); total (16/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.81e-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. |
CCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTAT > NC_000913/124954‑125220 | ccACATTCAGTCGCTGACTATCCCGACCTGTATCTCTTACGACCCGGCTTACACTTAAGAAATTGCTATCATAATACATAACGGACTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGaa > 1:143322/1‑143 (MQ=255) cacaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTAATACA‑CACTACGCTGAAc > 1:472985/1‑143 (MQ=255) acaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTCCTGTCCTCATACATAACGAACTGGAGCGTATGTACGGTGAAAAAAAAGAGAACGATAACTACTACA‑CACTACGCTCaaag > 1:218680/1‑141 (MQ=255) cTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGAAGGTCTGGAGCGAATGTACGGAGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCCAAACAAAAACAACcaa > 2:425400/1‑142 (MQ=255) tCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCCTGCATGACGGTCTGAAGCCCATGTACGGTGAAAAACAACAGAAAGTGTAATAGTTC‑CAACTAAGATAAACGAAAACTACCACATGccag > 1:122849/1‑141 (MQ=255) tctcTTACGGCCCGGCGTACGCGTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAg < 2:390354/143‑1 (MQ=255) ccGGCTTACCCTTACGAAGTTGCTGTCATAATGCATGACTGTCTGTAGCGTATGAACGGTGAAAAAGAAGAGAACGTCTACTACTACA‑CACTACGCTCAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGaa < 1:139649/143‑1 (MQ=255) cGGCTTAGGCTGACGAAGTTGCTGTCATCATGAATGACGGGCTGGAGCGGATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAg < 1:329025/143‑1 (MQ=255) cGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAg < 1:438572/143‑1 (MQ=255) tacgAAGTTGCTGTCATCATGCATGGCGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTaaa < 2:329181/143‑1 (MQ=255) tacgAAGTTGCTGTCATCATGCATGACGTTTTGGAGGGTATGTACGGTGAAAAACAAGAGAACGTTTACTAATACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTaaa < 1:55287/143‑1 (MQ=255) gAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGt < 1:419642/143‑1 (MQ=255) aaGTTGGTGTCATCATGCATGACGGTCTAGAGCTTATGTACGGTGAAATACAAGGTGACGTTTACTACTACA‑CGCTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTa < 2:370241/143‑1 (MQ=255) ttGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATACGTAAAGatatag > 2:226652/1‑141 (MQ=255) gCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGACGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTAc > 1:102044/1‑143 (MQ=255) tGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGGTGAGGAAGGTATCCGTAAAGGTATCTACaa > 2:461229/1‑143 (MQ=255) gTCAGCATGCATGGTGGTCTTGAGCGTGTGTGTGGTGAAAAACAAGATTATTGTTACTTCTACA‑CACGAAGCTGAACGAGAACTATCACATGGCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACaaa < 2:370046/143‑1 (MQ=255) catgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACt > 1:527878/1‑143 (MQ=255) gAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTCTCCGTAACGGTATCTACAAACTCGAAACTAGTGAAGGTAGc > 2:217424/1‑143 (MQ=255) gAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAggacga > 1:172242/1‑139 (MQ=255) gCGTATGTACGGTGAAAAACAAGAGAGCGTTTACTACTACA‑CGCTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCaa < 2:363170/143‑1 (MQ=255) tccGGTGAAAAACAAGAGAACGTTTACTGCTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATTCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTaaa < 2:3388/141‑1 (MQ=255) tACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGCAGGTGCTGCGGAAGGTATCCGTAAAGGTATCTAAAAACTCGAAAATATTGAAGGTAGAAAAGGTaaa > 1:407803/1‑143 (MQ=255) aCGGTTAAAAACAAAAAAATGTTAACTAATA‑AACAATATACTGAACGAAAACTACCAAATGCCGGCAATGCAGGAAGGTGATGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAg < 1:48177/143‑1 (MQ=255) ggTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAAGTGAAGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGtt > 1:165778/1‑143 (MQ=255) aaCAAGAGAACATCTACTATTA‑TTTATTACGATGAACGAAAACTACCACATGCCGACAATGACGGAAGGTCCTGAGGAAGGTATCCGTACAGGTATCTACAAACTAGAGACTATTGAAGGTAGCACAGGTAAAGTTCAgctcc > 2:223206/1‑141 (MQ=255) aaGAGAACGTTTACTACTATA‑CTCTACGCTGAACGTAAACTGTCTTTTTTCTACAAAACAGGACGGTGCTATATAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCg < 1:548699/143‑1 (MQ=255) agaACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACAACATGCCGACAATGCCGCAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACCCGAAACGATTGCAGCTACCAAAGGTAAAATTCAGCTGATCGGCt > 2:14091/1‑143 (MQ=255) ctactacA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAAGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGCTTCTAt > 1:455781/1‑143 (MQ=255) | CCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTAT > NC_000913/124954‑125220 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |