Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,907,454 | C→G | 20.0% | R127P (CGG→CCG) | PP_4316 ← | D‑isomer specific 2‑hydroxyacid dehydrogenase family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,907,454 | 0 | C | G | 20.0% | 42.0 / 6.9 | 20 | R127P (CGG→CCG) | PP_4316 | D‑isomer specific 2‑hydroxyacid dehydrogenase family protein |
Reads supporting (aligned to +/- strand): ref base C (8/8); new base G (2/2); total (10/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01 |
GTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCT > NC_002947/4907313‑4907598 | gtgCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTg > 1:20294/1‑150 (MQ=255) gTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGcgcg > 1:237404/1‑150 (MQ=255) cTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCggggg < 2:302854/150‑1 (MQ=255) cTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCggggg < 2:20294/150‑1 (MQ=255) ccAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGGTGTCGCGCAGCAGCGCca > 2:205186/1‑150 (MQ=255) gCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGAt > 1:117982/1‑121 (MQ=255) gCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGAt < 2:117982/121‑1 (MQ=255) ttGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATTTCGCGCAGCAGCGCCAGCAGCAGGGCCAgg < 1:355180/150‑1 (MQ=255) cGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGgtgt > 2:72343/1‑150 (MQ=255) aTCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGgtgtcc < 2:155324/150‑3 (MQ=255) cTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCGGCGTGCTGGCATCGGCGCggggg > 1:266489/1‑98 (MQ=255) cTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCGGCGTGCTGGCATCGGCGCggggg > 1:266518/1‑98 (MQ=255) cTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCGGCGTGCTGGCATCGGCGCggggg < 2:266518/98‑1 (MQ=255) cTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCGGCGTGCTGGCATCGGCGCggggg < 2:266489/98‑1 (MQ=255) cAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCAtt < 1:162425/150‑1 (MQ=255) aGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGcccc < 2:331378/140‑1 (MQ=255) aGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGcccc > 1:331378/1‑140 (MQ=255) ggCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTgg < 1:72343/150‑1 (MQ=255) aCCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGTCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAgcggcggcg > 1:69901/1‑150 (MQ=255) ccATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGACCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTc > 1:51152/1‑150 (MQ=255) cgcgccgcgTGCTGGCATCGGCGCGGGGGATGTCGCGACGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCt > 1:498880/1‑150 (MQ=255) | GTGCGGCTGTGGTAGCTGACGTGCATGTCGAAGCCCAAGTGGGCCCGCTTGGCGATCGCCAGCCCCACTGCGCCGAGGCCGAGAATGCCCAGGCGCTTGCCACTCACCGAGGGGCTGATAACCCGGTTCCATTCACCGCGCCGCGTGCTGGCATCGGCGCGGGGGATGTCGCGCAGCAGCGCCAGCAGCAGGGCCAGGGTGTGGTCGGCGACGGCCGCGGCATTGGCCCCGGCGCCATTGGTGACGGTTATGCCGCGGGCAGCGGCGGCGGCCAGGTCGACCTGCT > NC_002947/4907313‑4907598 |
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: ‑ |