Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,787,945 | T→C | 12.9% | L230P (CTG→CCG) | pyrH → | uridylate kinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,787,945 | 0 | T | C | 12.9% | 79.2 / 3.6 | 35 | L230P (CTG→CCG) | pyrH | uridylate kinase |
Reads supporting (aligned to +/- strand): ref base T (13/16); new base C (2/2); total (17/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01 |
CGAGAAATTCGATCACCTGACCTACGATGAGGTTCTGGATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGATCAACGACATCAAGAAAGACGCGCAGGAGCGCATGGGCAAGTCCATCGAAGCCCTGTCCCGCAACCTGGCTGCAATCCGCACCGGTCGCGCTC > NC_002947/1787808‑1788093 | cGAGAAATTCGATCACCTGACCTACGATGAGGTTCTGGATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTAGCGCCCTGAGGAACTTTgtggt > 2:48559/1‑150 (MQ=255) gaAATTCGATCACCTGACCTACGATGAGGTTCTGGATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTggtggg < 1:355378/150‑1 (MQ=255) aaaTTCGATCACCTGACCTACGATGAGGTTCTGGATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTggtgggt < 1:188608/150‑1 (MQ=255) gATCACCTGACCTACGATGAGGTTCTGGATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTggtggg > 2:402687/1‑143 (MQ=255) gATCACCTGACCTACGATGAGGTTCTGGATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTggtggg < 1:402687/143‑1 (MQ=255) aTCACCTGACCTACGATGAGGTTCTGGATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCCGAACATTgtggt < 2:122443/139‑1 (MQ=255) aTCACCTGACCTACGATGAGGTTCTGGATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCCGAACATTgtggt > 1:122443/1‑139 (MQ=255) tGACCTACGATGAGGTTCTGGATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGca > 1:336531/1‑150 (MQ=255) cTACGATGAGGTTCTGGATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCCGAACATTGTGGTGGGTGGCGCTGAAGGCACACt > 1:211268/1‑150 (MQ=255) gTTCTGGATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGtggtgg < 2:87974/121‑1 (MQ=255) gTTCTGGATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGtggtgg > 1:87974/1‑121 (MQ=255) gTTCTGGATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGaa > 1:329671/1‑150 (MQ=255) ttCTGGATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAg > 1:235156/1‑150 (MQ=255) ggATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCa > 1:491493/1‑150 (MQ=255) ggTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGa < 1:241668/128‑1 (MQ=255) ggTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGa > 2:241668/1‑128 (MQ=255) gggtgcgACCACAAGCCGCCTGTGTCGTCAACAACAGCTGCCATAGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCCGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGATCAACGACATCaag < 2:211268/144‑1 (MQ=255) tCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGATCAACGACATCaagaaa < 1:257698/150‑1 (MQ=255) tGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAACGCACACTGATCGAGGAAGGCCAAGCATGATCAACGACATCAAGAAAGACGCGCAGGAgc < 2:111169/150‑1 (MQ=255) gCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTggtggg > 1:354598/1‑79 (MQ=255) gCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTggtggg < 2:354598/79‑1 (MQ=255) gCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGATCAACGACATCAAGAAAGACGCGCAGGAgcg < 2:144809/150‑1 (MQ=255) tgCCGTGACCACAAGATGCCATTGCGGGTATTCAACATTAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGATCAACGACATCAAGAAAGACGCGCAGGAGCGCATgg < 1:48559/150‑1 (MQ=255) tgCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGATCAACGACATCAAGAAAGACGCGCAGGAGCGCATgg < 2:336531/150‑1 (MQ=255) gATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGATCAACGACATCAAGAAAGACGCGCAGGAGCGCATGGGCAAGTCCATCGaa < 1:490374/150‑1 (MQ=255) gCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGATCAACGACATCAAGAAAGACGCGCAGGAGCGCATGGGCAAGTCCATCGAAGcc < 2:235156/150‑1 (MQ=255) aTGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGgcg > 1:379/1‑43 (MQ=255) aTGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGgcg < 2:379/43‑1 (MQ=255) aTGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGAt > 1:179/1‑79 (MQ=255) aTGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGAt < 2:179/79‑1 (MQ=255) aTGAACAAGCCTGGCGCCCTGCAGAAACTCACTGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATTATAAACGACATCAAGAAAGACGCGCAGGAGCGCATGGGCAAGTCCATCGAAGCCATGTCCCGCAACCTGGCTg > 1:395426/1‑150 (MQ=255) gAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGATCAACGACATCAAGAAAGACGCGCAGGAGCGCATGGGCAAGTCCATCGAAGCCCTGTCCCGCAACCTGGCTGCa > 2:159220/1‑150 (MQ=255) cAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGATCAACGACATCAAGAAAGACGCGCAGGAGCGCATGGGCAAGTCCATCGAAGCCCt < 1:178313/128‑1 (MQ=255) cAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGATCAACGACATCAAGAAAGACGCGCAGGAGCGCATGGGCAAGTCCATCGAAGCCCt > 2:178313/1‑128 (MQ=255) ccTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGATCAACGACATCAAGAAAGACGCGCAGGAGCGCATGGGCAAGTCCATCGAAGCCCTGTCCCGCAACCTGGCTGCAATCCGCa > 1:127996/1‑150 (MQ=255) ctgAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGATCAACGACATCAAGAAAGACGCGCAGGAGCGCATGGGCAAGTCCATCGAAGCCCTGTCCCGCAACCTGGCTGCAATCCGCACCGGTCGCGCTc > 1:89552/1‑150 (MQ=255) | CGAGAAATTCGATCACCTGACCTACGATGAGGTTCTGGATCGCAAGCTGGGTGTGATGGATCTGACCGCAATCTGCCTGTGCCGTGACCACAAGATGCCATTGCGGGTATTCAACATGAACAAGCCTGGCGCCCTGCTGAACATTGTGGTGGGTGGCGCTGAAGGCACACTGATCGAGGAAGGCCAAGCATGATCAACGACATCAAGAAAGACGCGCAGGAGCGCATGGGCAAGTCCATCGAAGCCCTGTCCCGCAACCTGGCTGCAATCCGCACCGGTCGCGCTC > NC_002947/1787808‑1788093 |
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: ‑ |