Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,141,212 | T→A | 22.7% | L44Q (CTG→CAG) | PP_3645 → | L‑proline amide hydrolase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,141,212 | 0 | T | A | 22.7% | 37.6 / 9.9 | 22 | L44Q (CTG→CAG) | PP_3645 | L‑proline amide hydrolase |
Reads supporting (aligned to +/- strand): ref base T (7/10); new base A (2/3); total (9/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01 |
GCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACC > NC_002947/4141073‑4141351 | gcgcCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCAGGCCAACGaa < 2:113167/150‑2 (MQ=255) cTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTACCACGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACcg > 1:573901/1‑150 (MQ=255) ccaaTGATTCAGACTCTGGATGACCCCACCCCTGCCGGCAGCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACcgc < 2:287523/147‑1 (MQ=255) tttCAGCATGTCGCTGAACCTACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATc < 2:573901/149‑1 (MQ=255) cTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGc < 1:56452/150‑1 (MQ=255) aCGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCgg < 2:484079/150‑1 (MQ=255) acacCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCAc > 1:295719/1‑150 (MQ=255) ccTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCAGGCCAACGa < 1:342044/113‑1 (MQ=255) ccTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCAGGCCAACGa > 2:342044/1‑113 (MQ=255) tCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGgc < 2:295719/150‑1 (MQ=255) gCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTAGCGCGTCATCGCCTACGACATGGTAGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGa > 1:177611/1‑150 (MQ=255) cTGATCCACGGCGTGGGACGGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGAc < 2:156086/150‑1 (MQ=255) gTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgc < 2:511957/99‑1 (MQ=255) gTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCcgc > 1:511957/1‑99 (MQ=255) cTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCta < 1:120096/150‑2 (MQ=255) aGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGaccacc > 2:229466/1‑150 (MQ=255) gAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCAGCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTg > 2:276754/1‑150 (MQ=255) ggCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGca < 1:166904/150‑1 (MQ=255) gTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCAc < 1:229466/150‑1 (MQ=255) ttCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCAcc > 2:80181/1‑150 (MQ=255) tCGTTGGCCAGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTAt > 2:431696/1‑98 (MQ=255) tCGTTGGCCAGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTAt < 1:431696/98‑1 (MQ=255) | GCGCCCTCCATGATTCAGCCTGTCGCTGAACGCACACCTGCCGGCACCAGCTATCTGGATGTCGGCCAGGGCCAACCCGTGGTGCTGATCCACGGCGTGGGCCTGAACAAGGAAATGTGGGGCGGTCAGTTCGTTGGCCTGGCCAACGACTACCGCGTCATCGCCTACGACATGCTCGGCCACGGTCAAAGCCGCGTACCGGCCGCCGACACACCGCTGGAAGGCTATGCCGACCAGCTCGCCGAACTGCTTGACCACCTGCAAATTGCACAAGCCACC > NC_002947/4141073‑4141351 |
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: ‑ |