Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,935,901 | C→G | 15.4% | intergenic (+153/‑109) | PP_2568 → / → PP_2569 | ring‑cleaving dioxygenase/MHS family metabolite MFS transporter |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,935,901 | 0 | C | G | 15.4% | 61.6 / 6.5 | 26 | intergenic (+153/‑109) | PP_2568/PP_2569 | ring‑cleaving dioxygenase/MHS family metabolite MFS transporter |
Reads supporting (aligned to +/- strand): ref base C (10/12); new base G (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 = 4.09e-01 |
AAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCC > NC_002947/2935762‑2936047 | aaaaTCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGGCCGCCtcgga < 2:124811/150‑3 (MQ=255) ccAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGTCGGCCCGCCGTTTTGGGCTTCgcgc > 2:286397/1‑150 (MQ=255) cccGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAAcc > 1:128290/1‑150 (MQ=255) cccGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGACTCGTCGTGCTTCgcgcgaaa > 2:284901/1‑146 (MQ=255) gCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCaaa < 2:500044/150‑1 (MQ=255) tATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAaagc < 2:128290/150‑1 (MQ=255) aCCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAggg > 1:506549/1‑150 (MQ=255) cccGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGcc > 2:420023/1‑150 (MQ=255) cGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCTCGGCGGGGGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGGGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTg > 2:1466/1‑150 (MQ=255) cGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTg < 1:145766/150‑1 (MQ=255) gtcCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGaa < 1:258643/150‑1 (MQ=255) tcCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGGCCGCCtcgga < 2:121796/100‑3 (MQ=255) tcCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGGCCGCCtcgga > 1:121796/1‑98 (MQ=255) ccAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCCAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGt > 2:281716/1‑150 (MQ=255) cgccCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTtgt < 2:506549/150‑1 (MQ=255) gccCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTtgtg < 1:422948/150‑1 (MQ=255) ctctCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGCGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAgg < 1:411967/150‑1 (MQ=255) tcGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGc > 1:369459/1‑150 (MQ=255) aGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGc > 2:182822/1‑150 (MQ=255) cccGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGt < 1:284901/150‑1 (MQ=255) gTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGacaac < 2:185257/150‑1 (MQ=255) cAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAg < 1:182822/150‑1 (MQ=255) tCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAAccc > 2:487898/1‑150 (MQ=255) gcGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCAct < 2:407602/150‑2 (MQ=255) cGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAgc < 1:422881/150‑1 (MQ=255) gAGGCGGGCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTACTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAgcg > 1:564214/1‑150 (MQ=255) cggcCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGcc < 2:369459/150‑1 (MQ=255) | AAAATCGTCAGCCAACCCGCTATCGCCCGGCACCCCGGTGCTGCGTGTCGTCCAGGCCGCCCTCTCGTGCAGGGCCCGGCGGGTGGATATCCAGCCAATCCCGATTATGCTTCATGAACCTTGCTGCGATCCGAGGCGGCCCGCCTCGTCGTGCTTCGCGCAACCAAAGCAAGCGCCAGGGCCTGGGTGCGACGGTGAAGTGGTTGTGAGGTGCTCCGGCCTGTATGACAACATGACGGAGGCAACCCATGGCGGTTCTCGACAGCACATCCACGGGCAGCGCGCC > NC_002947/2935762‑2936047 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |