Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,736,356 | C→T | M274I (ATG→ATA) | bax ← | putative glycoside hydrolase Bax |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,736,356 | 0 | C | T | 83.3% | 65.4 / 7.3 | 30 | M274I (ATG→ATA) | bax | putative glycoside hydrolase Bax |
Reads supporting (aligned to +/- strand): ref base C (0/5); new base T (14/11); total (14/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.47e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01 |
GTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGT > NC_000913/3736279‑3736430 | gTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATtatata > 1:668721/1‑81 (MQ=255) cTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATgcg < 1:668696/81‑1 (MQ=255) aaCTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGcgcc < 2:668697/80‑1 (MQ=255) aCTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGcgccgc < 1:668698/81‑1 (MQ=255) ccTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGAtt < 2:668699/81‑1 (MQ=255) agaCAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAAc < 1:668700/80‑1 (MQ=255) agaCAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAAc < 1:668701/80‑1 (MQ=255) agaCAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAAc < 2:668702/80‑1 (MQ=255) gTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGtt < 1:668703/81‑1 (MQ=255) aaaGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATcc < 1:668704/81‑1 (MQ=255) aaaGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATcc < 2:668705/81‑1 (MQ=255) aaGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATCCt > 1:668739/1‑81 (MQ=255) gCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATCCTgg > 1:668737/1‑81 (MQ=255) cTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATCCTGGTa > 2:668738/1‑81 (MQ=255) ttATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATCCTGGTa > 2:668732/1‑80 (MQ=255) tATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACa > 1:668731/1‑81 (MQ=255) aTCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACAt > 1:668726/1‑81 (MQ=255) tCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACAtt > 1:668730/1‑81 (MQ=255) ccGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTg > 1:668755/1‑81 (MQ=255) aCTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCg > 2:668734/1‑80 (MQ=255) aCAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAg < 2:668706/81‑1 (MQ=255) gTCGGATAAGGCTCAAGGAAAAGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAg < 1:668707/81‑1 (MQ=255) ggATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAgttg > 2:668736/1‑81 (MQ=255) ggATAAGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAgtt < 2:668708/80‑1 (MQ=255) aTAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAgttgtt < 2:668709/81‑1 (MQ=255) aaGGCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTta > 2:668742/1‑80 (MQ=255) aGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTtata < 1:668710/81‑1 (MQ=255) gCTCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAAc > 2:668729/1‑81 (MQ=255) tCAAGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAActc > 2:668752/1‑81 (MQ=255) aaGGAAATGCAATTATATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCtt > 1:668745/1‑81 (MQ=255) cAATTATATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGt < 1:668711/81‑1 (MQ=255) | GTAGCTCAACTGCCTGAGACAAAGTAAAAAGCCTTATCCGACTGACAAGTCGGATAAGGCTCAAGGAAATGCAATTACATATGCGCCGCGATTAACCGTTGGTTATCCTGGTACATTGCGAACAGGTAGTTGTTATAACTCTTCCCTTTGGT > NC_000913/3736279‑3736430 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |