Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 2,907,559 | C→A | 10.5% | L179L (CTG→CTT) | sdhA ← | succinate dehydrogenase (flavoprotein subunit) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 2,907,559 | 0 | C | A | 10.5% | 87.8 / 2.9 | 35 | L179L (CTG→CTT) | sdhA | succinate dehydrogenase (flavoprotein subunit) |
Reads supporting (aligned to +/- strand): ref base C (17/14); new base A (2/2); total (19/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01 |
TCCCTGCTGATAAACGATTGACGCAGCAGATCCTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCTAAAACAGCGCCAAGGAATTCCCAGCCTTCGTATTTTGTGACCAGACCCGCAACTTCGTATCTGCGCACTTGCTCGTCCAGTGCGTAA > BS168_reference_genome/2907424‑2907688 | tCCCTGCTGATAAACGATTGACGCAGCAGATCCTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAgg > 4:83490/1‑139 (MQ=255) gctgATAAACGATTGACGCAGCAGATCCTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCTTAGCCAAAATGACAGCTTCTGAACGGAAGGAATCAATTTGCATGTTTGTCAGGTTTtg > 7:242953/1‑139 (MQ=255) gctgATAAACGATTGACGCAGCAGATCCTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTtg > 2:129426/1‑139 (MQ=255) aaCGATTGACGCAGCAGATCCTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGAGTCTGAACGGAAGGATTAAATTTGCATGTTTGTCAGGTTTTGTGgaac > 6:38079/1‑138 (MQ=255) cGATTGACGCAGCAGATCCTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTAAGGTTTTGTGcaac < 3:158832/136‑1 (MQ=255) cGATTGACGCAGCAGATCCTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTAAGGTTTTGTGcaac > 8:190313/1‑136 (MQ=255) cGATTGACGCAGCAGATCCTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTAAGGTTTTGTGcaac > 4:158832/1‑136 (MQ=255) cGATTGACGCAGCAGATCCTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTAAGGTTTTGTGcaac < 7:190313/136‑1 (MQ=255) aTTGACGCAGCAGATCCTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCTATTTGCATGTTTGTCAGGTTTTGTGCAGCTATTc > 8:122564/1‑139 (MQ=255) gcagcaGATCCTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGc < 7:58584/139‑1 (MQ=255) gcagcaGATCCTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGc < 7:112894/139‑1 (MQ=255) cgATCCGTTATTAATCAGGGAGTTCGTGGATGTTCCAAAGATGATTCCCGGGCCGCCCGTCGCCATAATGACCGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAgg < 2:103351/136‑1 (MQ=255) agATCCTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGtt < 8:243156/139‑1 (MQ=255) ccTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGt > 8:261636/1‑139 (MQ=255) ccTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGt > 8:33607/1‑139 (MQ=255) tATTAATCATGGAGTTCTTTGATTTTCCAAAGAGGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGtca < 4:105785/137‑1 (MQ=255) tATTAATCATGGAGTGCGTTGATTTTCAAAAGATGATTCCCGGTCCGCCCGGCGCCATAATGACAGCGTCGGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGtca < 6:85158/137‑1 (MQ=255) aTCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTtg < 5:58263/104‑1 (MQ=255) aTCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTtg > 6:58263/1‑104 (MQ=255) aTCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCt > 6:203287/1‑138 (MQ=255) aTCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCt > 5:252173/1‑138 (MQ=255) ggAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTAGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCTAAAACa > 1:157742/1‑139 (MQ=255) ggAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATCACAGCGTCTGAACGGAAGGATTAAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCTGGTTCTGTCATCATCTaaaa > 5:290588/1‑137 (MQ=255) cGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCTAAAACAGCGcc > 7:131765/1‑138 (MQ=255) cGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCTAAAACAGCGcc > 8:132690/1‑138 (MQ=255) cGTTGATTTTCCAAAGATGATTCCCGGTAAGCACGTCGCCATAATGACAGCGTCTGAACTGAAGTATTCAATTTGAATGTTTGTCATGTTTTGTTCAACAATTCCGCGGCAGGTACTTTCATCATATAAAACAGCGcc > 3:88140/1‑138 (MQ=255) gATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCTAAAACAGCGCCAAgg < 1:129426/138‑1 (MQ=255) ccATAATGACTGCGTCGGACCGGAAGGATTAAATTTGAATGTTTGGCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCTAAAACAGCGCCAAGGAATTCCCAGCCTTCGTATTTTGTGACCAGACCCGc < 5:38079/139‑1 (MQ=255) ccATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCTAAAACAGCGCCAAGGAATTCCCAGCCTTCGTATTTTGTGACCAGACCCGc < 7:132690/139‑1 (MQ=255) aTGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCTAAAACAGCGCCAAGGAATTCCCAGCCTTCGTATTTTGTGACCAGACCCGCAACtt < 3:83490/139‑1 (MQ=255) aaCGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCTAAAACAGCGCCAAGGAATTCCCAGCCTTCGTATTTTGTGACCAGACCCGCAACTTCGTATCTGCGCa < 6:252173/138‑1 (MQ=255) aTTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCTAAAACAGCGCCAAGGAATTCCCAGCCTTCGTATTTTGTGACCAGACCCGCAACTTCGTATCTGCGCACTTGCTCGTc > 6:72414/1‑139 (MQ=255) cAATTTGCATGTTTGGCCGGCTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCTAAAACAGCGCCAAGGAATTCCCAGCCTTCGTATTTTGTGACCAGACCCGCAACTTCGTATCTGCGCACTTGCTCGTCCa < 5:87438/138‑1 (MQ=255) cATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCTAAAACAGCGCCAAGGAATTCCCAGCCTTCGTATTTTGTGACCAGACCCGCAACTTCGTATCTGCGCACTTGCTCGTCCAGTGCGTaa < 5:203287/139‑1 (MQ=255) cATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCTAAAACAGCGCCAAGGAATTCCCAGCCTTCGTATTTTGTGACCAGACCCGCAACTTCGTATCTGCGCACTTGCTCGTCCAGTGCGTa > 4:128492/1‑138 (MQ=255) | TCCCTGCTGATAAACGATTGACGCAGCAGATCCTGTATTAATCATGGAGTTCGTTGATTTTCCAAAGATGATTCCCGGTCCGCCCGTCGCCATAATGACAGCGTCTGAACGGAAGGATTCAATTTGCATGTTTGTCAGGTTTTGTGCAACAATTCCGCGGCAGGTTCTGTCATCATCTAAAACAGCGCCAAGGAATTCCCAGCCTTCGTATTTTGTGACCAGACCCGCAACTTCGTATCTGCGCACTTGCTCGTCCAGTGCGTAA > BS168_reference_genome/2907424‑2907688 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |