Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 125,387 | C→T | 100% | Q791* (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 | 125,387 | 0 | C | T | 98.4% | 72.3 / ‑6.7 | 30 | Q791* (CAG→TAG) | aceE | pyruvate dehydrogenase E1 component |
Reads supporting (aligned to +/- strand): ref base C (2/0); new base T (15/12); total (18/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.98e-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: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
ATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAG‑GATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAG > NC_000913/125249‑125524 | aTCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGCTGGTCAG‑GATTGTGAACGCTGGAGCATGCTGCACCCGCTGGAAACTCGGCGCGTTCCGTATAGCGCTGAggg > 1:563000/1‑142 (MQ=255) ccTGGCGAAAGAGTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGAAAAAT‑ACAAACTGTTATTGTGATCGCTCGAACATTCTTCACCCGCTGTAAACTCCGCGCGTTCCGTATATCGCTTAGGtga < 1:814471/143‑1 (MQ=255) cTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAG‑GATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGGTCCGTATATCGCTCAGGtgat > 1:720842/1‑143 (MQ=255) tGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAG‑GATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgatg > 1:410791/1‑143 (MQ=255) aGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAG‑GATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacgacg > 1:743577/1‑143 (MQ=255) gCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAG‑GCTTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGGTGAACGACGCTCCGGCa > 2:141557/1‑143 (MQ=255) gCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAG‑GATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCa > 1:359396/1‑143 (MQ=255) taagTTCGGGCGTTTATCACAGGACCGCCTTCACCGAGCTGGCGCGTGATGGTCAG‑GATTGTGAACGCTGGAACATGCTGCTCCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTg < 2:319601/140‑1 (MQ=255) ggTTCTGACTTTTATATCGTGATCTCCTTCACCGAGCTGGCGCGGGATGGTCAG‑GATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGc < 1:624209/143‑1 (MQ=255) tttATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAG‑GATTGTGAACGCTGGAACATGCTGCACCCGCGGGAAACTCCGCGCGTTCCGTATGTCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGAc > 1:70417/1‑143 (MQ=255) cGTGGCCTCCTTCACCGAGCTGGCGCGTGATGGTCAG‑GACTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGa < 2:659110/143‑1 (MQ=255) ctTCACCGAGCTGGCGCGTGATGGTCAG‑GATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCg < 1:666383/143‑1 (MQ=255) tCAG‑GATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTCGGTGAAGAACGCCGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGACCAGGTCCGTACTTACGTAc > 2:160988/1‑143 (MQ=255) tCAG‑GATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACCATATGAAACTGTTCGCTGAGCCGGTCCGTACCTACGTAc > 2:180493/1‑143 (MQ=255) cAG‑GATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTAcc > 1:662706/1‑143 (MQ=255) aG‑GATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCg < 1:270863/143‑1 (MQ=255) gATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGATGGGCAGGTCCGCACTTACGTACCGGc > 2:128394/1‑143 (MQ=255) gccgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg < 1:676214/140‑1 (MQ=255) cccGCCGGAACCTCCGCGCGTTCCGTAGTTCGGTTAGGTGGTGAACGACGCTCCGGCAGTGGCATCTACCGACTGTATGAAACTGTTCGCTGAGCAGGTCCGTACGTACGTACCGGCTGACGACTACCGCGTACTGGGTACTg < 2:472251/143‑1 (MQ=255) ccGCGGGAAACTCCGCGCGTTCCGTATATCGCTCGGGGGGTGAACGACGCGCGGGCCGTGTCGTCTACCGCCCGTATGACACTGTTCCCTGAGCAGATCCAAAACAACGTACACACTGACCACCACCGCATACCAGGTACTgc > 2:642619/1‑142 (MQ=255) cGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGAt > 1:220351/1‑143 (MQ=255) cGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTCTTCGCCGAGCAGGTCCGTACTTACCTACCGGCTGACCACTACAGCGGACCGGGTACTgcg > 2:768248/1‑141 (MQ=255) gCTGGAACCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATg < 1:609371/143‑1 (MQ=255) gCTGGAACCTCCGCGCGTGCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAGCTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATg < 2:806500/143‑1 (MQ=255) aaaCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGCTGGCTTc > 2:694396/1‑143 (MQ=255) tCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTCCCGCGTACTGGGTACTGATGGCTTCGGTc > 1:631889/1‑143 (MQ=255) tCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTc > 2:78141/1‑143 (MQ=255) atatCGCTTAGGTCATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGt < 1:506735/143‑1 (MQ=255) tatCGCTTAGGTGAGGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTg < 1:172678/143‑1 (MQ=255) tCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGagag > 1:573109/1‑143 (MQ=255) | ATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAG‑GATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAG > NC_000913/125249‑125524 |
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: ‑ |