Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BX571857 | 10,884 | G→A | 33.3% | intergenic (‑299/‑9) | SAS0007 ← / → SAS0008 | conserved hypothetical protein/putative histidine ammonia‑lyase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BX571857 | 10,884 | 0 | G | A | 33.3% | 46.9 / 35.0 | 33 | intergenic (‑299/‑9) | SAS0007/SAS0008 | conserved hypothetical protein/putative histidine ammonia‑lyase |
Reads supporting (aligned to +/- strand): ref base G (11/11); new base A (6/5); total (17/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 = 1.00e+00 |
ATAATAGGTTATAAATGTACGTAAAACCATTGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGAT > BX571857/10791‑10979 | ataataGGTTATAAATGTACGTAAAACCATTGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGATaa < 1:87841/97‑1 (MQ=255) taataGGTTATAAATGTACGTAAAACCATTGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGATAAc < 1:56235/97‑1 (MQ=255) aataGGTTATAAATGTACGTAAAACCATTGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGATACTTAATCTAAACAAGGGGGTaaaa > 2:163362/1‑95 (MQ=255) taGGTTATAAATGTACGTAAAACCATTGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACata < 1:158174/97‑1 (MQ=255) tataAATGTACGTAAAACCATTGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGa < 1:54329/97‑1 (MQ=255) tataAATGTACGTAAAACCATTGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGa < 1:41880/97‑1 (MQ=255) taAATGTACGTAAAACCATTGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGATAACATATATGACt > 1:36683/1‑97 (MQ=255) aaaCCATTGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGAt > 1:128099/1‑97 (MQ=255) aaaCCATTGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGAt > 2:218288/1‑97 (MQ=255) aaaCCATTGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGAt > 1:220697/1‑97 (MQ=255) aaaCCATTGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGAt > 2:127427/1‑97 (MQ=255) aaaCCATTGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGATAACATATATGACTTTATATTTAGAt > 2:149958/1‑97 (MQ=255) aaaCCATTGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGATAACATATATGACTTTATATTTAGAt < 1:176818/97‑1 (MQ=255) ttGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGATAACATATATGACTTTATATTTAGATGGTGaa > 2:259578/1‑97 (MQ=255) tAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACa > 1:82098/1‑97 (MQ=255) aaaTGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAAtt < 1:266537/97‑1 (MQ=255) aaaTGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAAtt < 1:266579/97‑1 (MQ=255) aaaTGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAAtt < 1:22543/97‑1 (MQ=255) aaTGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTg > 1:168375/1‑97 (MQ=255) aaaCGTATACTTCAAGAAGGATGGGTTACTTAAGATAAACAAGGGGATAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGat < 1:145733/97‑1 (MQ=255) cGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATAtt > 2:67961/1‑97 (MQ=255) taCTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAAt > 1:71378/1‑97 (MQ=255) cTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCa < 1:140405/97‑1 (MQ=255) tCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCAtt < 1:344523/97‑1 (MQ=255) tCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCAtt < 1:11833/97‑1 (MQ=255) aagaagGATGGGTTACTTAATATAAACAAGGGGATAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCAtttt > 2:304703/1‑97 (MQ=255) agaagGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCAttttt < 1:118984/97‑1 (MQ=255) ggTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTacaacaaca > 1:144314/1‑97 (MQ=255) ggTTACTTAATATAAACAAGGGGATAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTacaacaaca > 1:234813/1‑97 (MQ=255) gTTACTTAATATAAACCAGTGGATAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTacaacaacgg < 2:89494/97‑3 (MQ=255) gTTACTTAATATAAACAAGGGGATAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTAcaacaacaa > 1:372691/1‑97 (MQ=255) aCTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCa > 2:319650/1‑97 (MQ=255) aaaCAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGaa < 1:210184/97‑1 (MQ=255) ggTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATtgat < 2:248371/97‑1 (MQ=255) | ATAATAGGTTATAAATGTACGTAAAACCATTGTTTCAATAAAAATGAAAACGTATACTTCAAGAAGGATGGGTTACTTAATATAAACAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGAT > BX571857/10791‑10979 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |