Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,883,770 | C→A | 18.1% | R167S (CGC→AGC) | pucL → | 2‑oxo‑4‑hydroxy‑4‑carboxy‑5‑ureidoimidazoline decarboxylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,883,770 | 0 | C | A | 18.1% | 60.6 / 3.2 | 22 | R167S (CGC→AGC) | pucL | 2‑oxo‑4‑hydroxy‑4‑carboxy‑5‑ureidoimidazoline decarboxylase |
Reads supporting (aligned to +/- strand): ref base C (6/12); new base A (2/2); total (8/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.02e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01 |
CAGTTCCCGTTCATCATGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATGCGCACCCTTATGATCGAGCCCTTGACCAAAGAAGCCTTCGCCCAGTTCGGAGACGTGA > NC_002947/4883635‑4883897 | cAGTTCCCGTTCATCATGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTcc < 1:13947/136‑1 (MQ=255) cccGTTCATCATGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCACCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTtg < 1:95831/136‑1 (MQ=255) ttCATCATGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTtgctgc > 2:304694/1‑136 (MQ=255) tcatcaTGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCa > 1:151188/1‑136 (MQ=255) tcatcaTGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCa > 1:457548/1‑136 (MQ=255) atcatGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGc < 1:263051/136‑1 (MQ=255) ggCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTaa < 1:33645/136‑1 (MQ=255) cACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAg < 1:106905/136‑1 (MQ=255) cAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGaaa < 2:50482/136‑1 (MQ=255) tctcGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGaa < 1:304694/136‑1 (MQ=255) gccgccTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAAt < 2:357438/136‑1 (MQ=255) cATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCAt < 1:215105/136‑1 (MQ=255) aTCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATg < 1:427970/136‑1 (MQ=255) aTCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATg < 1:238488/136‑1 (MQ=255) aTCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATg < 2:4804/136‑1 (MQ=255) tgccgatgccgaGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCAGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATg < 1:449990/125‑1 (MQ=255) tgccgatgccgaGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCAGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATg < 1:452102/125‑1 (MQ=255) tgccgatgccgaGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCAGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATg > 2:449990/1‑125 (MQ=255) tgccgatgccgaGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCAGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATg > 2:452102/1‑125 (MQ=255) cgatgccgaGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATGCGCACCCTTATGAt > 1:215036/1‑136 (MQ=255) cGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGATCAATAGAAGAGATAACCGCATGCGCACCCTTATGATCGAGCCCTTGACCAAAGAAGCCTTCGCCCAGTTCGGAGACGt < 2:476368/136‑1 (MQ=255) cccTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATGCGCACCCTTATGATCGAGCCCTTGACCAAAGAAGCCTTCGCCCAGTTCGGAGACGTGa > 1:481714/1‑136 (MQ=255) cccTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATGCGCACCCTTATGATCGAGCCCTTGACCAAAGAAGCCTTCGCCCAGTTCGGAGACGTGa > 2:321999/1‑136 (MQ=255) | CAGTTCCCGTTCATCATGGCGGTGAAGGGCAGCAACCGGCACCAGATTCTCGCCGCCTTTGAAAAGCGCATCCACAACGATGCCGATGCCGAGTTCAAGGAAGCCCTGGCACAGATCAACCTGATCGCCCTGTTCCGCTTGCTGCAGCTGTAAGCAACAGCCGAACCTTATTCAGAAAACGAACAATAGAAGAGATAACCGCATGCGCACCCTTATGATCGAGCCCTTGACCAAAGAAGCCTTCGCCCAGTTCGGAGACGTGA > NC_002947/4883635‑4883897 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |