Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 124,241 | C→T | 100% | Q409* (CAG→TAG) | aceE → | pyruvate dehydrogenase E1 component |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 124,241 | 0 | C | T | 98.3% | 38.4 / ‑6.7 | 18 | Q409* (CAG→TAG) | aceE | pyruvate dehydrogenase E1 component |
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (6/11); minor base G (1/0); total (7/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.89e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GTCACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTAT > NC_000913/124107‑124369 | gtCAGGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGGAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAg < 2:381038/143‑1 (MQ=255) tCCGAAGAAAATGTACGCTGCATGCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGGGCAGACCATTAAAGGGTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGa < 2:358913/143‑1 (MQ=255) aTCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACgg > 2:511643/1‑143 (MQ=255) tACGCTGCATTCAAGAAAGCGCAGGAAACGAAAGGCAAGGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGCAGGTAAAAACATCGCGCAGTAGGTTAAGGTAGTGAACATGGTCGgttt > 1:377016/1‑141 (MQ=255) tGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTc > 2:494834/1‑143 (MQ=255) ttCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGGGGCTGAAGATAAAAACATCGCGCGCTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCacat > 2:374104/1‑143 (MQ=255) tacaGCGACATGAATCCTTGCTCCTGCCAGTAAAGCTTACTGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGgtgt < 2:4631/140‑1 (MQ=255) cAAAGCGACAGTAATCCGTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGgtgt < 2:157941/143‑1 (MQ=255) aaaGCGACAGTACTCCTTGCTCATACCATTAAAGGTTACGGCATGGCCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTc < 1:110361/143‑1 (MQ=255) cGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGAt < 1:536654/143‑1 (MQ=255) aaTCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGTGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGata < 1:127316/143‑1 (MQ=255) tCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATc > 2:524680/1‑143 (MQ=255) cTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACGAGGGTAAGAAAATGAACATGGACGGTGTGCGTCATATGCGCGCCCGCTTCAATGTGCCGGGGTCTGATGCAGACATCGa > 1:145629/1‑143 (MQ=255) tGCTCATGCCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAATCATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGaaa < 2:26784/143‑1 (MQ=255) aTACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTg < 1:511643/143‑1 (MQ=255) gggCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAg < 1:436319/143‑1 (MQ=255) gtgggaCAGCGCGCGCTAGGTTAGGAAGAGGGACACGGACGGTGTGCCTCATCTCCGCGACCGTGCCAATGTGCCGGTGTCTGCTGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGCGCATACCt < 2:287647/138‑1 (MQ=255) aaaaCATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTGCCCGGAAGGTTCTGAAGAGCAGACCTAt > 2:413304/1‑143 (MQ=255) | GTCACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTAT > NC_000913/124107‑124369 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |