Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,335,493 | G→A | 14.9% | P336L (CCC→CTC) | PP_2052 ← | bifunctional sugar‑phosphatase/mannitol‑1‑phosphate 5‑dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,335,493 | 0 | G | A | 14.9% | 55.7 / 2.3 | 26 | P336L (CCC→CTC) | PP_2052 | bifunctional sugar‑phosphatase/mannitol‑1‑phosphate 5‑dehydrogenase |
Reads supporting (aligned to +/- strand): ref base G (10/12); new base A (2/2); total (12/14) | |||||||||||
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.16e-01 |
AACTCGGCTTGCACATGACGCCGCACGAACGCTGCACCGCCAACCTTGTTCAGCACGATCAGCAGGGTCAGCTCGCGGCCTCGGCGCTCGAAGCGTTGCAGCAGGCCTTGGGCGATCACCCTGGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGAACGACGTGGCACCGTAGCGCACGCTGTAGCGGCCGAAGGCGCTGACGCT > NC_002947/2335346‑2335641 | aaCTCGGCTTGCACATGACGCCGCACGAACGCTGCACCGCCAACCTTGTTCAGCACGATCAGCAGGGTCAGCTCGCGGCCTCGGCGCTCGAAGCGTTGCAGCAGGCCTTGGGCGATCACCCTGGCCTGGCTTCTGATGGCCTGCTCGGGc > 2:192037/1‑150 (MQ=255) aTGACGCCGCACGAACGCTGCACCGCCAACCTTGTTCAGCACGATCAGCAGGGTCAGCTCGCGGCCTCGGCGCTCGAAGCGTTGCAGCAGGCCTTGGGCGATCACCCTGGCCTGTCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCgat < 2:627/150‑1 (MQ=255) aTGACGCCGCACGAACGCTGCACCGCCAACCTTGTTCAGCACGATCAGCAGGGTCAGCTCGCGGCCTCGGCGCTCGAAGCGTTGCAGCAGGCCTTGGGCGATCACCCTGGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCgat < 1:309396/150‑1 (MQ=255) gACGCCGCACGAACGCTGCACCGCCGACCTTGTTCAGCACGATCAGCAGGGTCAGCTCGCGGCCTCGGCGCTCGAAGCGTTGCAGCAGGCCTTGGGCGATCACCCTAGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCgatga < 2:400063/150‑1 (MQ=255) gcACGAACGCTGCACCGCCAACCTTGTTCAGCACGATCAGCAGGGTCAGCTCGCGGCCTCGGCGCTCGAAGCGTTGCAGCAGGCCTTGGGCGATCACCCTGGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCGATGATTTCgg < 1:192037/150‑1 (MQ=255) cAACCTTGTTCAGCACGATCAGCAGGGTCAGCTCGCGGCCTCGGCGCTCGAAGCGTTGCAGCAGGCCTTGGGCGATCACCCTGGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCtga < 2:405866/150‑1 (MQ=255) aGCACGATCAGCAGGGTCAGCTCGCGGCCTCGGCGCTCGAAGCGTTGCAGCAGGCCTTGGGCGATCACCCTGGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGcgtcg > 2:132855/1‑150 (MQ=255) cACGATCAGCAGGGTCAGCTCGCGGCCTCGGCGCTCGAAGCGTTGCAGCAGGCCTTGGGCGATCACCCTGGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGcgtcgtc < 2:63088/150‑1 (MQ=255) aTCAGCAGGGTCAGCTCGCGGCCTCGGCGCTCGAAGCGTTGCAGCAGGCCTTGGGCGATCACCCTGGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGcgtcgtcgtcg > 2:303083/1‑150 (MQ=255) gcagGGTCAGCTCGCGGCCTCGGCGCTCGAAGCGTTGCAGCAGGCCTTGGGCGATCACCCTGGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCa > 2:101885/1‑150 (MQ=255) gcagGGTCAGCTCGCGGCCTCGGCGCTCGAAGCGTTGCAGCAGGCATTGAGCGATCACCCTGGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGGCGATGAGTTCGGCAGTGGTGTACATGATGATGACCGCGTCGTCGTCGTCCa > 1:474124/1‑150 (MQ=255) tCAGCTCGCGGCCTCGGCGCTCGAAGCGTTGCAGCAGGCCTTGGGCGATCACCCTGGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGa > 1:78844/1‑150 (MQ=255) cTCGCGGCCTCGGCGCTCGAAGCGTTGCAGCAGGCCTTGGGCGATCACCCTGGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCAt < 1:303083/150‑1 (MQ=255) cagcagGCCTTGGGCGATCACCCTGACCTGGCTTCTGATGGCCTGCTCGGGCAGGCCCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGaa > 2:62079/1‑150 (MQ=255) agGCCTTGGGCGATCACCCTGGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGAacgac < 1:132855/150‑1 (MQ=255) gCCTTGGGCGATCACACTGGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGAACGACGt < 2:152047/150‑1 (MQ=255) ccTTGGGCGATCACCCTGGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGAACGACGTg > 1:240389/1‑150 (MQ=255) tgagCCTGCCCGAGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGaa > 2:115119/4‑113 (MQ=255) tgagCCTGCCCGAGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGaa < 1:115119/110‑1 (MQ=255) gagCCTGCCCGAGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGAACGACGTgg < 2:158718/118‑1 (MQ=255) gagCCTGCCCGAGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGAACGACGTgg > 1:158718/3‑120 (MQ=255) ccTGCTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGAACGACGTGGCACCGTAGCGCACGCTGTAGCGGCCGAAGgcgc < 1:101885/150‑1 (MQ=255) cTGCTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGAACGACGTGGCACCGTAGCGCACGCTGTAGCGGCCGAAGGCGCt < 1:58108/150‑1 (MQ=255) tGCTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGAACGACGTGGCACCGTAGCGCACGCTGCAGCGGCCGAAGGCGCTg > 2:320194/1‑150 (MQ=255) cTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGAACGACGTGGCACCGTAgc < 2:136447/122‑1 (MQ=255) cTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGAACGACGTGGCACCGTAgc > 1:136447/1‑122 (MQ=255) gggCAGGCTCAGGCCGATGATCTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGAACGACGTGGCACCGTAGCGCACGCTGTAGCGGCCGAAGGCGCTGACGCt > 2:365518/1‑150 (MQ=255) | AACTCGGCTTGCACATGACGCCGCACGAACGCTGCACCGCCAACCTTGTTCAGCACGATCAGCAGGGTCAGCTCGCGGCCTCGGCGCTCGAAGCGTTGCAGCAGGCCTTGGGCGATCACCCTGGCCTGGCTTCTGATGGCCTGCTCGGGCAGGCTCAGGCCGATGATTTCGGCAGTGGTGTACATGCTGATGACCGCGTCGTCGTCGTCCATGTCGATCATGCGCACGTTGTCGATGGTCTGGTCGAACGACGTGGCACCGTAGCGCACGCTGTAGCGGCCGAAGGCGCTGACGCT > NC_002947/2335346‑2335641 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |