Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,314,956 | A→G | 11.8% | S36P (TCA→CCA) | msrP ← | methionine sulfoxide reductase catalytic subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,314,956 | 0 | A | G | 11.8% | 81.5 / 3.5 | 34 | S36P (TCA→CCA) | msrP | methionine sulfoxide reductase catalytic subunit |
Reads supporting (aligned to +/- strand): ref base A (14/16); new base G (2/2); total (16/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 = 3.23e-01 |
GTCCTTCACCGTCACCGCCTGCCAGCGCGTGGCCGCGAGCTTGTCAGTGAACCAGGCCGGCGCGGCCCCAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGACCTGGGAAGCTTGATAAGCATGAACGGCTCCATAGCACTGAGTAACGGATGTACCAGTAAGACTATGGAGCCGTGAG > NC_002947/5314807‑5315105 | gTCCTTCACCGTCACCGCCTGCCAGCGCGTGGCCGCGAGCTTGTCAGTGAACCAGGCCGGCGCGGCCCCAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGa > 2:289051/1‑150 (MQ=255) gtcaccgCCTGCCAGCGCGTGGCCGCGAGCTTGTCAGTGAACCAGGCCGGCGCGGCCCCAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCcagc > 1:336816/1‑150 (MQ=255) ccgCCTGCCAGCGCGTGGCCGCGAGCTTGTCAGTGAATCAGGCCGGCGCGGCCCCAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAggg > 1:80598/1‑150 (MQ=255) gcgcGTGGCCGCGAGCTTGTCAGTGAACCAGGCCGGCGCGGCCCCAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGa > 1:1857/1‑150 (MQ=255) gcgTGGCCGCGAGCTTGTCAGTGAACCAGGCCGGCGCGGCCCCAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTATGACGGGaaa < 1:42694/150‑1 (MQ=255) tGGCCGCGAGCTTGTCAGTGAACCAGGCCGGCGCGGCCCCAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGt < 1:306518/150‑1 (MQ=255) gCTTGTCAGTGAACCAGGCCGGCGCGGCCCCAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCtt > 2:473484/1‑150 (MQ=255) tGTCAGTGAACCAGGCCGGCGCGGCCCCAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCCCCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCgg < 1:412473/150‑1 (MQ=255) tCAGTGAACCAGGCCGGCGCGGCCCCAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAggg > 1:349223/1‑121 (MQ=255) tCAGTGAACCAGGCCGGCGCGGCCCCAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAggc < 2:349223/121‑2 (MQ=255) tCAGTGAACCAGGCCGGCGCGGCCCCAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCgtgg < 1:507339/150‑1 (MQ=255) gTGAACCAGGCCGGCGCGGCCCCAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTg > 2:48961/1‑150 (MQ=255) ccAGGCCGGCGCGGCCCCAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCTGGGGTGAtctc > 1:459791/1‑150 (MQ=255) cccAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCcagc > 2:407590/1‑94 (MQ=255) cccAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCcagc < 1:407590/94‑1 (MQ=255) cccAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTc > 2:361913/1‑150 (MQ=255) caTCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGACCTgg < 1:473484/150‑1 (MQ=255) cgACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGGGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGAc > 1:345224/1‑133 (MQ=255) cgACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGGGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGAc > 1:345175/1‑133 (MQ=255) cgACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGGGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGAc < 2:345175/133‑1 (MQ=255) cgACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGGGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGAc < 2:345224/133‑1 (MQ=255) cgACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGACCTGGGAAGCTTGATAAg < 2:80598/150‑1 (MQ=255) gACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGa > 1:431996/1‑83 (MQ=255) gACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGa < 2:431996/83‑1 (MQ=255) gACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGACCTGGGAAGCTTGATAAGc < 2:336816/150‑1 (MQ=255) ttGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGACCTGGGAAGCTTGATAAGCATGaa < 2:1857/150‑1 (MQ=255) cccgcccgGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGACCTGGGAAGCTTGATAAGCATGAACGGCTCCATAGCACt > 1:8161/1‑150 (MQ=255) aaTGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGACCTGGGAAGCTTGATAAGCATGAACGGCTCCATAGCACTGAGTAAc < 1:309813/147‑1 (MQ=255) aaTGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGACCTGGGAAGCTTGATAAGCATGAACGGCTCCATAGCACTGAGTAAc > 2:309813/1‑147 (MQ=255) gcCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGACCTGGGAAGCTTGATAAGCATGAACGGCTCCATAGCACTGAGTAACGGATGTAc < 2:8161/150‑1 (MQ=255) cgcCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGCAATCTATATGTCCTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGACCTGGGAAGCTTGATAAGCATGAACGGCTCCATAGCACTGAGTAACGGATGTACCAGTaa < 1:289051/150‑1 (MQ=255) gcCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGACCTGGGAAGCTTGATAAGCATGAACGGCTCCATAGCACTGAGTAACGGa > 2:418282/1‑138 (MQ=255) gcCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCAATCGGCAGACCTGGGAAGCTTGATAAGCATGAACGGCTCCATAGCACTGAGTAACGGa < 1:418282/138‑1 (MQ=255) agccaATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGACCTGGGAAGCTTGAt < 2:257197/96‑1 (MQ=255) agccaATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGACCTGGGAAGCTTGAt > 1:257197/1‑96 (MQ=255) aGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCGGCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGACCTGGGAAGCTTGATAAGCATGAACGGCTCCATAGCACTGAGTAACGGATGTACCAGTAAGACTATGGAGCCGTGAg < 2:120302/150‑1 (MQ=255) | GTCCTTCACCGTCACCGCCTGCCAGCGCGTGGCCGCGAGCTTGTCAGTGAACCAGGCCGGCGCGGCCCCAGCCTGCACATCGGCATAACGCGACATTTGCGCCGCGCCCACCCCGCCCGGCAATGCGCCTAGCGCCAAGCCAGCCAATGAGCCACCCAGCAGGGTACGACGGGAAAGGTAAATGCCTTCGGGGGTGATCTCCGACGCCTTGCACTCGGAAGACCTGGGAAGCTTGATAAGCATGAACGGCTCCATAGCACTGAGTAACGGATGTACCAGTAAGACTATGGAGCCGTGAG > NC_002947/5314807‑5315105 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |