Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 543,062 | T→A | 11.6% | intergenic (+29/+199) | glxK → / ← allE | glycerate 2‑kinase 2/(S)‑ureidoglycine aminohydrolase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 543,062 | 0 | T | A | 11.6% | 81.5 / 3.3 | 36 | intergenic (+29/+199) | glxK/allE | glycerate 2‑kinase 2/(S)‑ureidoglycine aminohydrolase |
Reads supporting (aligned to +/- strand): ref base T (11/19); new base A (2/2); total (13/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.27e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-01 |
CTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAA > NC_000913/542924‑543196 | cTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑Cg > 1:216239/1‑143 (MQ=255) cTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGTAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATaa < 2:26351/143‑1 (MQ=255) cACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATATAAAT‑CGCAtccgctt < 2:376267/143‑7 (MQ=255) aCCTGTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTGTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTg < 2:62393/143‑1 (MQ=255) cGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGCATCCGGCAt < 1:197192/143‑1 (MQ=255) aaaCCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGGATGGGTa > 1:48271/1‑143 (MQ=255) aaaCCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTa > 1:618156/1‑143 (MQ=255) ctTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtcctt > 1:125127/1‑143 (MQ=255) gaaCAGCGCGGGAAATATTGCCTGCGCCATTAAGATAGGTGAGGGAATTAAAAACTAACCCTTACCTTTGAAGCGGAGGCGATTTATAT‑CGCAGAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttcct < 2:251945/142‑1 (MQ=255) tcgtgAAATATCGCCTGCTCCATGAAAATAGGTCCGGGAATTATAAACTAACCCTTACCTTTACAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGCGTACCATATCCTTCCTGAAAATc < 2:159905/139‑1 (MQ=255) gcgcgAAATATTGCCTGCGCCATGAAAATAGGTCAGGCAACTAAAAACTAACCCTTACCTTTAAAGCGGATCCGATTTATAT‑CGCATAAGAGTGCACTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATc < 2:285296/143‑1 (MQ=255) gcgcgAAATATTGCCTGCGCAATTAAATTAGGTCAGTGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATc < 2:71842/143‑1 (MQ=255) cgAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATcgcg > 1:170856/1‑143 (MQ=255) cttGCCTGCTCCATTAAAACAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAAtt < 1:602339/142‑1 (MQ=255) ttGCCTGCGACATTAATATAGGACAGGGAATTAAAAACTAACCCTTAGCTTTAAAGCGGATGCGATTTATAG‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTc < 1:102970/143‑1 (MQ=255) gcgcCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTGAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGTTGCGGCATGAGTACCATATCCTTCCTGAAAAACAAGTAAATTCtatata > 1:201250/1‑143 (MQ=255) ctaTTTAAATATGTCAGTGAATTAAAAACGAACCCTTACCATAGAATTGGATTCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTg < 2:390768/141‑1 (MQ=255) aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTGTATATTGCa > 1:495093/1‑143 (MQ=255) gataagcggTCAGGGGATTAAAAACTGACCCTTAGCTTGACAGCGGATGCGATTTATAT‑CGCATGAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAGCCTTCCTGGATATCGCGCAAATTCTATATATTGCaga < 2:12844/136‑1 (MQ=255) tAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCCGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCa < 1:247762/143‑1 (MQ=255) gtcgtgcaaTTAAGAGCTACCCCTTACCTTTAAAGCGGAGTCTATTTATCT‑CGCATAAGAGTGCAGTACGCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGt < 2:353447/136‑1 (MQ=38) tCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGGTGCGAACAAAATACACA‑AAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTa < 1:149235/143‑1 (MQ=255) tCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGCAAAAAGAGCATATTCTTTATCTTGCAGCGCTCATGTa > 1:444597/1‑143 (MQ=255) ggAATTAAATACTAACCCTTAGCGTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGcc < 2:127441/143‑1 (MQ=255) ggAATTAAAAGCTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGcc < 2:127444/143‑1 (MQ=255) gagttAAAAACTAACCCTTACCGTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCt < 2:313092/140‑1 (MQ=255) gacttAACACCTAACCATTACCTTTAAAGCGGATGCGATGTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTGCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCt < 2:190248/140‑1 (MQ=38) aaaaCTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaag > 1:345774/1‑143 (MQ=255) ccTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTCGGGCTGATgagcggagcgga > 2:166708/1‑141 (MQ=255) tttAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGCGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGcc > 2:261408/1‑142 (MQ=255) ttAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCaa < 1:265668/143‑1 (MQ=255) actcttaTGCGATATAAAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGt > 1:20816/7‑143 (MQ=255) gcgattGCGAGTTACAG‑CGCATAAGAGTGCAGTACTCCTGCCGGATGCGGCATGAGTACCATATCCTCCCTGAAACTCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTa < 2:221978/138‑1 (MQ=37) taTGCGATATAAAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaa > 1:390331/2‑143 (MQ=255) gATGCGATTTACAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaa < 1:39768/143‑1 (MQ=255) gCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaa < 1:349902/143‑1 (MQ=255) | CTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATAT‑CGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAA > NC_000913/542924‑543196 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |