Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,486,554 | C→G | 8.2% | I422M (ATC→ATG) | purD → | phosphoribosylamine‑‑glycine ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,486,554 | 0 | C | G | 8.2% | 153.1 / 5.1 | 49 | I422M (ATC→ATG) | purD | phosphoribosylamine‑‑glycine ligase |
Reads supporting (aligned to +/- strand): ref base C (21/24); new base G (2/2); total (23/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01 |
GCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGA > NC_002947/5486406‑5486685 | gCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACTACAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCg > 1:513001/1‑150 (MQ=255) accaACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGAAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAgcgc > 1:39975/1‑150 (MQ=255) caACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTCGCGGGCGATGGCCcgata > 2:465945/1‑141 (MQ=255) caACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTCGCGGGCGATGGCCcgata < 1:465945/144‑4 (MQ=255) aaCGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGt < 1:1371683/150‑1 (MQ=255) aaCGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGt > 2:716252/1‑150 (MQ=255) gggcgCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCAc < 2:629377/150‑1 (MQ=255) ggcgCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCAcc > 2:497683/1‑150 (MQ=255) gcgcGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACca < 2:1382409/150‑1 (MQ=255) cTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCagtaa < 1:11431/150‑1 (MQ=255) gtgtgCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCagtaagt > 2:673719/1‑150 (MQ=255) gtgtgCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCagtaa < 1:1269991/148‑1 (MQ=255) gtgtgCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCagtaa > 2:1269991/1‑148 (MQ=255) tgCCACCGCCATGGGCAGCACCGTTGCCGACGTGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGa < 2:569256/123‑1 (MQ=255) tgCCACCGCCATGGGCAGCACCGTTGCCGACGTGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGa > 1:569256/1‑123 (MQ=255) ccACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGtactac < 1:814365/150‑1 (MQ=255) aCCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCg < 1:1376407/150‑1 (MQ=255) ccGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCgg > 2:3222/1‑150 (MQ=255) ccATGGGCAGAACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTAt < 2:295987/150‑1 (MQ=255) cagcaCCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTa > 2:1274889/1‑149 (MQ=255) cagcaCCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTa < 1:1274889/149‑1 (MQ=255) gCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGc < 1:497683/150‑1 (MQ=255) cGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCt < 2:296634/150‑1 (MQ=255) gACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTg > 2:848004/1‑150 (MQ=255) cagcagcagGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGc > 2:877447/1‑126 (MQ=255) cagcagcagGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGc < 1:877447/126‑1 (MQ=255) agcagcagGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTg < 1:243693/143‑1 (MQ=255) agcagcagGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTg > 2:243693/1‑143 (MQ=255) agcagcagGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATc < 2:1398496/150‑1 (MQ=255) gcagGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGc > 1:740800/1‑150 (MQ=255) cTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCtcgtc > 1:896866/1‑150 (MQ=255) ccTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTCTGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAggg > 2:1029023/1‑150 (MQ=255) aaGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTc < 2:560662/114‑1 (MQ=255) aaGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTc > 1:560662/1‑114 (MQ=255) aaGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCGATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGaa < 1:1373212/150‑1 (MQ=255) ggAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATg < 2:39975/150‑1 (MQ=255) gCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGgccttgc < 1:716252/150‑1 (MQ=255) cTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCg < 2:369169/93‑1 (MQ=255) cTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCg > 1:369169/1‑93 (MQ=255) cTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGgccttgcc > 1:896820/1‑150 (MQ=255) cTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGgccttgcc < 1:258539/150‑1 (MQ=255) gCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACtt > 2:721580/1‑150 (MQ=255) aCCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGcgcgccg > 1:120361/1‑150 (MQ=255) gACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGt < 1:3222/150‑1 (MQ=255) cATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGa < 2:740800/150‑1 (MQ=255) ggCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGgc > 2:766793/1‑113 (MQ=255) ggCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGgc < 1:766793/113‑1 (MQ=255) tcgcGGGCGATGGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGc < 2:612382/83‑1 (MQ=255) tcgcGGGCGATGGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGc > 1:612382/4‑86 (MQ=255) | GCAAGGTAACCACCAACGGCGGGCGCGTGCTGTGTGCCACCGCCATGGGCAGCACCGTTGCCGACGCGCAGCAGCAGGCTTATCGCCTGGCCAAGGAAGTCAGCTGGAACGGCAGCTTCTACCGTTCCGACATCGGCTATCGGGCCATCGCCCGCGAGCGCGGTGAGCACCAGCAGTAAGTACTACCGGTATGCCGTAGCGGTTTCAGCCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGA > NC_002947/5486406‑5486685 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |