Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 2,738,963 | A→T | 12.5% | D219V (GAT→GTT) | murI → | glutamate racemase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 2,738,963 | 0 | A | T | 12.5% | 80.9 / 5.1 | 32 | D219V (GAT→GTT) | murI | glutamate racemase |
Reads supporting (aligned to +/- strand): ref base A (18/10); new base T (2/2); total (20/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.20e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01 |
AACAATATGGAACGATTGTGTTAGGATGCACTCATTTTCCGTTTTTCAAAAACAGTTTTGAAAAACTGTTTGGGATAAAAGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCAAAGTACAAAAGATTGTTTGAAATATTAGATGAAACTCAAAGGAGTCATGTTGGACA > BS168_reference_genome/2738825‑2739101 | aaCAATATGGAACGATTGTGTTAGGATGCACTCATTTTCCGTTTTTCAAAAACAGGTTTGAAAAACTGTTTGGGATAAAAGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGa < 8:86717/139‑1 (MQ=255) gAACGATTGTGTTAGGATGCACTCATTTTCCGTTTTTCAAAAACAGTTTTGAAAAACTGTTTGGGATAAAAGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAAc > 5:203846/1‑137 (MQ=255) gAACGATTGTGTTAGGATGCACTCATTTTCCGTTTTTCAAAAACAGTTTTGAAAAACTGTTTGGGATAAAAGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAAc > 4:41843/1‑137 (MQ=255) ttGTGTTAGGATGCACTCATTTTCCGTTTTTCAAAAACAGTTTTGAAAAACTGTTTGGGATAAAAGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGTTCGTAACCAGCTTgg > 3:96460/1‑139 (MQ=255) ttGTGTTAGGATGCACTCATTTTCCGTTTTTCAAAAACAGTTTTGAAAAACTGTTTGGGATAAAAGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGTTCGTAACCAGCTTgg > 3:99907/1‑139 (MQ=255) ttGTGTTAGGATGCACTCATTTTCCGTTTTTCAAAAACAGTTTTGAAAAACTGTTTGGGATAAAAGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTgg > 5:169493/1‑139 (MQ=255) aTGCACTCATTTTCCGTTTTTCAAAAACAGTTTTGAAAAACTGTTTGGGATAAAAGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGTTCGTAACCAGCTTGGTAAAGGATCg < 4:99907/139‑1 (MQ=255) aTGCACTCATTTTCCGTTTTTCAAAAACAGTTTTGAAAAACTGTTTGGGATAAAAGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGTTCGTAACCAGCTTGGTAAAGGATCg < 4:96460/139‑1 (MQ=255) cACTCATTTTCCGTTTTTCAAAAACAGTTTTGAAAAACTGTTTGGGATAAAAGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCggg < 1:86209/138‑1 (MQ=255) cACTCATTTTCCGTTTTTCAAAAACAGTTTTGAAAAACTGTTTGGGATAAAAGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCggg > 2:86209/1‑138 (MQ=255) cACTCATTTTCCGTTTTTCAAAAACAGTTTTGAAAAACTGTTTGGGATAAAAGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGc > 6:104981/1‑139 (MQ=255) aTTTTCCGTTTTTCAAAAACAGTTTTGAAAAACTGTTTGGGATAAAAGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCaa > 2:108205/1‑138 (MQ=255) ggATAAACGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCaag < 6:130642/138‑1 (MQ=255) ggATAAAAGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCaaga < 3:41843/139‑1 (MQ=255) gATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCa < 6:203846/138‑1 (MQ=255) ttCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTc > 8:88330/1‑96 (MQ=255) ttCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTc < 7:88330/96‑1 (MQ=255) aGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCAAAGTACAAAAGa < 5:65233/139‑1 (MQ=255) tgtCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGCCCAAGAAGTCATATCAAAGTACAAAAGATTGttt > 5:202173/1‑139 (MQ=255) tgtCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCAAAGTACAAAAGATTGttt > 2:32718/1‑139 (MQ=255) tgtCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCAAAGTACAAAAGATTGttt > 4:83787/1‑139 (MQ=255) tgtCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCAAAGTACAAAAGATTGttt > 6:53505/1‑139 (MQ=255) gCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCAAAGTACAAAAGATTGTTTGAAatat > 4:106677/1‑136 (MQ=255) gCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCAAAGTACAAAAGATTGTTTGAAatat > 6:216857/1‑136 (MQ=255) cAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCaaa > 4:123818/1‑112 (MQ=255) cAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCaaa < 3:123818/112‑1 (MQ=255) ttGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCAAAGTACAAAAGATTGTTTGAAATATTAGATGAAACt > 8:173774/1‑138 (MQ=255) ttGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCAAAGTACAAAAGATTGTTTGAAATATTAGATGAAACt > 3:219207/1‑138 (MQ=255) ttCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCaag > 7:164769/1‑81 (MQ=255) ttCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCaag < 8:164769/81‑1 (MQ=255) ttGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCAAAGTACAAAAGATTGTTTGAAATATTAGATGaa < 4:86158/122‑1 (MQ=255) ttGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCAAAGTACAAAAGATTGTTTGAAATATTAGATGaa > 3:86158/1‑122 (MQ=255) gCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCAAAGTACAAAAGATTGTTTGAAATATTAGATGAAACTCAAAGGAGTCATGt < 5:53505/137‑1 (MQ=255) aTCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCAAAGTACAAAAGATTGTTTGAAATATTAGATGAAACTCAAAGGAGTCATGTTGGACa < 6:169493/139‑1 (MQ=255) | AACAATATGGAACGATTGTGTTAGGATGCACTCATTTTCCGTTTTTCAAAAACAGTTTTGAAAAACTGTTTGGGATAAAAGTTGATATGATTTCAGGCAGTGTCGGTACTGCAAAACAGTTGAAAAAAGTTCTTGCAGATCGTAACCAGCTTGGTAAAGGATCGGGCTCAATCACCTTTTTCAATTCTGGCCATAAAATTGTGGACCAAGAAGTCATATCAAAGTACAAAAGATTGTTTGAAATATTAGATGAAACTCAAAGGAGTCATGTTGGACA > BS168_reference_genome/2738825‑2739101 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |