Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 387,343 | G→C | 7.8% | E160D (GAG→GAC) | glyA‑I → | serine hydroxymethyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 387,343 | 0 | G | C | 7.8% | 155.8 / 5.3 | 51 | E160D (GAG→GAC) | glyA‑I | serine hydroxymethyltransferase |
Reads supporting (aligned to +/- strand): ref base G (24/23); new base C (2/2); total (26/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01 |
CCGGTGACACCATCCTGGGCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGGCGCACTGCTGTT > NC_002947/387198‑387459 | ccGGTGACACCATCCTGGGCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGcc < 2:986439/150‑1 (MQ=255) cGGTGACACCATCCTGGGCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGACCGcca > 2:407305/1‑149 (MQ=255) tGACACCATCCTGGGCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGc < 2:173203/150‑1 (MQ=255) acCATCCTGGGCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTc > 1:22099/1‑150 (MQ=255) cATCCTGGGCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGa > 1:678327/1‑150 (MQ=255) tCCTGGGCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTGCGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGc > 1:356244/1‑150 (MQ=255) ggCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAg > 2:634438/1‑150 (MQ=255) ggCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGACCGCCAGGCGCTCGacttcgtc < 1:407305/150‑8 (MQ=38) cATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGcc > 2:477916/1‑150 (MQ=255) aGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGACCGCCAGGCGCTCGa < 2:1062363/137‑1 (MQ=255) aGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGACCGCCAGGCGCTCGa > 1:1062363/1‑137 (MQ=255) gCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGaaaa < 1:480034/150‑1 (MQ=255) ccACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGt > 1:578354/1‑150 (MQ=255) ggcggTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGcc < 2:985899/150‑1 (MQ=255) ggccgTCGCCTGACCCACTGTGGCAAGGGGGCTGCCTCTGGCAAACTGTACAACGCCGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTCTCTGTCGAGCGCAAGCCCAAAATGCTCGTTGcc < 2:82321/150‑1 (MQ=17) ggTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGt < 1:471819/150‑1 (MQ=255) ggTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGt < 2:578354/150‑1 (MQ=255) ggCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTcc < 1:634438/150‑1 (MQ=255) ccAAAGTTTCGTCATCGGTCAAACTGTACACCGCACGCCAGTACGGCATCAATACCAACTCCCGCCTTATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCTGCCTCCTCCAAGa < 2:264537/150‑1 (MQ=21) cAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGttt < 2:229508/132‑1 (MQ=255) cAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGttt > 1:229508/1‑132 (MQ=255) cAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGAc > 1:1093802/1‑150 (MQ=255) aGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCGGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGAccc > 2:876051/1‑150 (MQ=255) gtgtCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGGTCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTc < 1:335498/150‑1 (MQ=255) gtgtCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTc > 1:280074/1‑150 (MQ=255) tgtCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCg > 1:334110/1‑150 (MQ=255) tgtCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCg > 1:453426/1‑150 (MQ=255) tcgtaCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGAc > 2:720510/1‑150 (MQ=255) cGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCAc < 1:876051/150‑1 (MQ=255) aaCTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTcc < 1:55551/150‑1 (MQ=255) cTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCcgc > 1:1206627/1‑150 (MQ=255) cTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTAATCCAAGACCCTCGACTTCCCACGTTTCcgc < 2:678327/150‑1 (MQ=255) tGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCcgcg > 2:478074/1‑150 (MQ=255) aaCGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCg > 2:914720/1‑108 (MQ=255) aaCGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCg < 1:914720/108‑1 (MQ=255) cAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATCATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAg > 2:1108743/1‑150 (MQ=255) gTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGt < 2:1093802/150‑1 (MQ=255) ggCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGgc > 2:113037/1‑150 (MQ=255) ggCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGgc < 2:152746/150‑1 (MQ=255) cATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGGCATCGCCGACAAGGTGGgcgc > 2:659056/1‑150 (MQ=255) tCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCATGGCGGTCGAGCACAACCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAAGTGGGCGCAc < 1:139894/150‑1 (MQ=255) gACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAACATGATCGTCGACGGTTTCTCGGCTGGTGCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGGCGCActg < 1:478074/150‑1 (MQ=25) acacCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGGCGCActgc > 2:506929/1‑150 (MQ=255) acacCAACACCGGACTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGGCGCActgc < 1:28079/150‑1 (MQ=255) cAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGGCGCACTGCTGtt > 2:340002/1‑150 (MQ=255) acCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGa > 2:475646/1‑98 (MQ=255) acCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGa < 1:475646/98‑1 (MQ=255) cGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGa > 1:814127/1‑89 (MQ=255) cGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGa < 2:814127/89‑1 (MQ=255) tCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTccc > 2:19091/1‑94 (MQ=255) tCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTccc < 1:19091/94‑1 (MQ=255) | CCGGTGACACCATCCTGGGCATGAGCCTGGCCCACGGCGGTCACCTGACCCACGGCGCCAAGGTGTCGTCCTCGGGCAAACTGTACAACGCAGTCCAGTACGGCATCGACACCAACACCGGCCTGATCGACTACGACGAAGTCGAGCGCCTGGCGGTCGAGCACAAGCCGAAAATGATCGTTGCCGGTTTCTCGGCCTACTCCAAGACCCTCGACTTCCCACGTTTCCGCGCCATCGCCGACAAGGTGGGCGCACTGCTGTT > NC_002947/387198‑387459 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |