Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 240,531 | G→C | 23.1% | A276P (GCC→CCC) | PP_0188 → | uroporphyrin‑III C‑methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 240,531 | 0 | G | C | 23.1% | 42.2 / 13.2 | 26 | A276P (GCC→CCC) | PP_0188 | uroporphyrin‑III C‑methyltransferase |
Reads supporting (aligned to +/- strand): ref base G (11/9); new base C (3/3); total (14/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 = 2.58e-01 |
CCGAGTTCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATGCA > NC_002947/240388‑240674 | ccGAGTTCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGCCCggcc > 2:232220/1‑148 (MQ=255) tCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGCCCGGCCAGCAAc < 1:232220/150‑1 (MQ=255) gcgcCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATG‑AAACGTGCGTCCGTTGCTGGCCg < 2:85536/116‑1 (MQ=255) gcgcCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATG‑AAACGTGCGTCCGTTGCTGGCCg > 1:85536/1‑116 (MQ=255) tGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTTTCCCTGAGCCTGa < 2:550310/150‑1 (MQ=255) ccgccgACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGGGCCTgctcg > 1:31724/1‑146 (MQ=255) tCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAATCGTGCGTCCGGTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGCGGGCGGCGCTCAAt > 2:208515/1‑150 (MQ=255) tCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAAt > 1:132244/1‑150 (MQ=255) cAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGt < 2:132244/150‑1 (MQ=255) gtggGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCTGCGCTCAATGGTGATGc < 1:26817/150‑1 (MQ=255) gCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGTGCCTGGCCCTGAGCCTGACTAt < 1:349163/110‑1 (MQ=255) gCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGTGCCTGGCCCTGAGCCTGACTAt > 2:349163/1‑110 (MQ=255) cTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTggg < 1:525767/98‑1 (MQ=255) cTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTggg > 2:525767/1‑98 (MQ=255) aTGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAg < 1:672873/110‑1 (MQ=255) aTGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAg > 2:672873/1‑110 (MQ=255) aCGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCa > 2:390988/1‑150 (MQ=255) cGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCaa < 2:694397/150‑1 (MQ=255) ccGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAAt > 2:635302/1‑150 (MQ=255) gTTGCTGCCCGGCCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATgg < 1:573340/84‑1 (MQ=255) gTTGCTGCCCGGCCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATgg > 2:573340/1‑84 (MQ=255) ttGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATACa > 1:184960/1‑150 (MQ=255) gCTGCCCGGCCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCaa > 2:50855/1‑138 (MQ=255) gCTGCCCGGCCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCaa < 1:50855/138‑1 (MQ=255) tGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATg < 1:800394/144‑1 (MQ=255) tGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATg > 2:800394/1‑144 (MQ=255) | CCGAGTTCGACAGTAACGCCGATGCGCTTGGCGCCCTGACCGCCGACGGTGATGGTGCCAGCAGGCTGTCGCAGTGGTGGGCAGAAATATCCAAGTACTTCCAGATCGACTTCAATGCCGATGAAAACGTGCGTCCGTTGCTGGCCGGGCAGCAACTGAACCAGCTGCGCCTGGCCCTGAGCCTGACTATCGAGCAGGCCCAGTGGGCGGCGCTCAATGGTGATGCCAAGGTTTACACGCAGGCGCTGGACGACGCCCGCAGCGTACTGCTGGCCAACTTCAATGCA > NC_002947/240388‑240674 |
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: ‑ |