Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,989,940 | C→T | 10.4% | R276H (CGC→CAC) | hemC ← | hydroxymethylbilane synthase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,989,940 | 0 | C | T | 10.4% | 96.9 / 3.1 | 40 | R276H (CGC→CAC) | hemC | hydroxymethylbilane synthase |
Reads supporting (aligned to +/- strand): ref base C (19/16); new base T (2/2); total (22/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01 |
AGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGT > NC_000913/3989802‑3990074 | agaCGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTc > 2:502836/1‑143 (MQ=255) gTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCCCCGCGaata > 1:893802/1‑143 (MQ=255) ggATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCCCCGCGGCGTTCACCGCGACTAATCTGCg > 1:756712/1‑143 (MQ=255) atccTCATGCCGGGGCGTCTCCGTTCTAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGaa < 2:740175/140‑1 (MQ=255) aTACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCCCCGCGGCGTTCACCGCGAATAATCTGCGaa > 1:315165/1‑143 (MQ=255) tCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGGTATTCAGTAGCACGTCTGCGCGCGAAATCGCCATTTGATCGGCATTTGCCGCCGCACAGGGCCTTTCACAGCGCATAATCTGCGAACCgc > 2:770584/1‑142 (MQ=255) cATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTc < 2:508635/143‑1 (MQ=255) atgagggggCGTCGCCGTCATAGACTTCAGCGAGGATCTCGCGCGCGCCGTGATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCAGCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTcc < 1:140024/138‑1 (MQ=255) aTGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTcc > 1:496792/1‑143 (MQ=255) tGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCg < 1:502836/143‑1 (MQ=255) cGGGGCTTCCCGGTTATAGACTTCAGCGAGGATCCGGCGCGCGCCGTTTGTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGgcg < 2:653818/143‑1 (MQ=255) gtaGACGTCAGCGAGGATCTCGCGCGCGCCCTGCTTCAGTAGCTCTGCTGCCAGCGAAATCCCCATTTGTTCGGCTTCTTGCGGCGCACCGCGGCGTTCACCTCGAATAATCTGCGATCCGTCCGGCGCGCCGACCAGcgcac < 2:240482/142‑1 (MQ=255) ataGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGTGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcac > 2:799350/1‑143 (MQ=255) gACTTCAGGGAGGCGCGCGTGCGCGCCGTTATTCCGTAGCTCTTCTGCCCGCGATTTGCGCATTTGTTCGGGATCTTGCGGCGCCCCGCGGCGTCCCCCGCGAATGGCCTGCGAACCGTCCGGCGCGCCGACCAGcgcacgca < 2:830422/143‑1 (MQ=255) ttCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCCTCTTGCGGCGCACCGCGGCGTTCACCGCGAGTAAGCTGCGAACCGTCCGGCGCGCCGACCAGCGCCCGCAGcc > 1:606595/1‑143 (MQ=255) ttCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGGGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGcc > 1:410211/1‑143 (MQ=255) cGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGTGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAtt > 2:327205/1‑143 (MQ=255) gAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGCCCAGCGCACGCAGCCAGAttt > 2:256151/1‑143 (MQ=255) gATCTCGCGCGCGCCGCTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGc < 1:578249/143‑1 (MQ=255) aTCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGCTCACCGCGAATAATCTGCGCACCGTCCGGCGCGCCGACCAGCGCCCGCAGCCCGAGTTCGcc > 2:811899/1‑143 (MQ=255) tcGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGCATAATCTGCGAACCGTCCGCCGCGCCGACCCGCGCACGCAGCCAGAGTTCGCCATc > 2:301363/1‑143 (MQ=255) gcgcCGTTATTCAGTAGCTCTTCTGCCAGCGAAACCCCCATTTTTTCTCCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAg < 1:884597/143‑1 (MQ=255) cGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATCAGCCCg > 2:494708/1‑143 (MQ=255) tATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCg > 1:544460/1‑143 (MQ=255) tCAGTAGCTCTTCTGCCAGCGAAACCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAg < 1:390118/143‑1 (MQ=255) gTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCCTCAAGAAGGCCCGCGTAGCTa > 2:459090/1‑143 (MQ=255) ctctTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGTGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCaa < 1:327205/143‑1 (MQ=255) ttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTg < 1:726107/143‑1 (MQ=255) ttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAACGCACGCAGCCAGACTTCGCCATCAATAAGCACGGCGTCGCTACCaaata > 2:377372/1‑140 (MQ=255) tctGCCATCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCGCCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTgg < 2:697788/143‑1 (MQ=255) gCCAGCGCAAGCCCCATTTGTGCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCAc < 2:111552/143‑1 (MQ=255) ccAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCCACCAAGTGGCAcc > 1:552469/1‑143 (MQ=255) gCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGa < 1:875684/143‑1 (MQ=255) ccATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCt > 1:681215/1‑143 (MQ=255) tttggTTCGGCATCGTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTc < 2:570685/139‑1 (MQ=255) tttGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGCATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCg < 1:746638/143‑1 (MQ=255) gCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAACCTGCGAACCGTCCGGCGCGCCGAAATGCGCACGACCCCCGATTTCGCCATCCATCACCTCGGCGTAGCTAGCGCTTGCCCCCTCTCCTCCGCCTTCgacgcgcgt > 1:562608/1‑136 (MQ=255) tCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGAGTact > 1:548867/1‑141 (MQ=255) cgcACCGCGGTGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGgcgc < 1:799350/143‑1 (MQ=255) cACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGCCATCCGCCTTCGAGACGCGCCTTCATGGCGCGt > 1:12489/1‑143 (MQ=255) | AGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGT > NC_000913/3989802‑3990074 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |