Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,363,720 | C→T | 13.8% | A32T (GCC→ACC) | glmM ← | phosphoglucosamine mutase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,363,720 | 0 | C | T | 13.8% | 62.9 / 4.1 | 29 | A32T (GCC→ACC) | glmM | phosphoglucosamine mutase |
Reads supporting (aligned to +/- strand): ref base C (15/10); new base T (2/2); total (17/12) | |||||||||||
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.38e-01 |
ATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATG > NC_002947/5363596‑5363860 | aTCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAg > 2:70711/1‑150 (MQ=255) tCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTc < 1:249766/146‑1 (MQ=255) tCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTc > 2:249766/1‑146 (MQ=255) tCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGc < 1:70711/150‑1 (MQ=255) tCGGCACCTGCCGCGGACAGGCCGACTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGc < 2:589804/150‑1 (MQ=255) tGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCgg < 2:81584/150‑1 (MQ=255) aGGCCGGCTTCGAGTGCGGACTCGAACATGTAGACGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCg > 1:557653/1‑150 (MQ=255) ccGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTTCGGAAGGCCATGCCGGCCGCCCAGCCAAGCCTCCGCATGAAGTCAGGCGTGATCGGGTATTCGCCGAcg > 1:409241/1‑150 (MQ=255) cGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGAcgc < 1:209434/150‑1 (MQ=255) gAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACg > 1:254141/1‑150 (MQ=255) aCGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCt < 2:114187/150‑1 (MQ=255) gcgTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAGAGTATTTTCTGCTCa > 1:197093/1‑150 (MQ=255) ttGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCgg > 1:24947/1‑121 (MQ=255) ttGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCgg < 2:24947/121‑1 (MQ=255) ttGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACt > 2:229559/1‑150 (MQ=255) cccACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCa > 1:602362/1‑150 (MQ=255) caGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGTCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCg > 1:576364/1‑109 (MQ=255) caGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGTCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCg > 1:576326/1‑109 (MQ=255) caGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGTCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCg < 2:576326/109‑1 (MQ=255) caGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGTCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCg < 2:576364/109‑1 (MQ=255) caGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGc > 1:551022/1‑98 (MQ=255) caGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGGATTCGCCGACGCGGc < 2:551022/98‑1 (MQ=255) tagccagaGGCTATGGAGGCCGACCAGCTAATCGTCAGCATGGAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGc > 2:499997/8‑150 (MQ=255) tgcGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCa < 1:229559/150‑1 (MQ=255) gAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGa > 1:504842/1‑150 (MQ=255) aaGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGctc < 1:96999/123‑1 (MQ=255) aaGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGctc > 2:96999/1‑123 (MQ=255) ccATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATg > 1:284408/1‑150 (MQ=255) ccATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATg > 1:260786/1‑150 (MQ=255) | ATCACATCGGCACCTGCCGCGGACAGGCCGGCTTCGAGTGCGGACTCGAACATGTAGCCGGAGATACGCGTGTCCTTGCCCACCAGCACGCGGCAATGGCCTTGCTTGCGGAAGGCCATGCCGGCCGCCCAGCCAAGCTTCAGCATGAAGTCAGGCGTGATCGGGTATTCGCCGACGCGGCCACGAATGCCGTCGGTACCAAAGTATTTTCTGCTCATAGGGACTCCAATGCTCTTATTCGGCGTTTTGCACCGCAGCGATCATG > NC_002947/5363596‑5363860 |
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: ‑ |