Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,653,708 | T→A | 27.6% | intergenic (+104/‑21) | PP_2325 → / → PP_2326 | 2‑dehydropantoate 2‑reductase/universal stress protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,653,708 | 0 | T | A | 27.6% | 26.2 / 7.6 | 18 | intergenic (+104/‑21) | PP_2325/PP_2326 | 2‑dehydropantoate 2‑reductase/universal stress protein |
Reads supporting (aligned to +/- strand): ref base T (5/8); new base A (2/3); total (7/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAGCGCTGAGTTTCATCGACTGCAGCGACCGGCCGGCCTGAGGGTGGCCGAGGCTTCGGCGGGGCCCTCCAGCCCACGGCAGCAGCGCGCCGGACGGCAAAGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCTGCGGCAGTCGAGGAGGTCGAATGATCCGCTCCATGCTGTACGCCACCGACCTCGGTGTCTACGCCCCTTTCGTGATGCAACATGCGCTGGCCCTGGCGCGCACGTTCGGTGCCGAGTTGTATGTGATCCACGCGGTCGAGC > NC_002947/2653564‑2653849 | aaGCGCTGAGTTTCATCGACTGCAGCGACCGGCCGGCCTGAGGGTGGCCGAGGCTTCGGCGGGGCCCTCCAGCCCACGGCAGCAGCGCGCCGGACGGCAAAGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCTGCGGc > 2:65492/1‑150 (MQ=255) tGAGTTTCATCGACTGCAGCGACCGGCCGGCCTGAGGGTGGCCGAGGCTTCGGCGGGGCCCTCCAGCCCACGGCAGCAGCGCGCCGGATGGCAAAGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCAGCGGCAgagca < 1:204933/150‑5 (MQ=255) tGAGTTTCATCGACTGCAGCGACCGGCCGGCCTGAGGGTGGCCGAGGCTTCGGCGGGGCCCTCCAGCCCACGGCAGCAGCGCGCCGGACGGCAAAGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCTGCGGCAGTCga < 1:65492/150‑1 (MQ=255) gAGTTCCATCGACTGCAGCGACCGGCCGGCCTGAGGGTGGCCGACGGCCCGCCGGGGCCCTCCAGCCCACGGCAGCAGCGCGCCGGACGGCAAAGCGAGCGTCCGGCGCGCTTCTAAGCTGAAGCTTGCTCTGCCGCTGCGGCAGTCgag < 1:691801/150‑1 (MQ=255) tttCATCGACTGCAGCGACCGGCCGGCCTGAGGGTGGCCGAGGCTTCGGCGGGGCCCTCCAGCCCACGGCAGCAGCGCGCCGGACGGCAAAGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCTGCGGCAGTCgaggag < 2:506981/150‑1 (MQ=255) tCATCGACTGCAGCGACCGGCCGGCCTGAGGGTGGCCGAGGCTTCGGCGGGGCCCTCCAGCCCACGGCTGCAGCGCGCCGGACGGCAAAGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCTGCGGCAGTCGAGGAGGt < 1:248728/150‑1 (MQ=255) tGAGGGTGGCCGAGGCTTCGGCGGGGCCCTCCAGCCCACGGCAGCAGCGCGCCGGACGGCAAAGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCTGCGGCAGTCGAGGAGGTCGAATGATCCGCTCCATGCTGTACGc > 1:230799/1‑150 (MQ=255) ggTGGCCGAGGCTTCGGCGGGGCCCTCCAGCCCACGGCAGCAGCGCGCCGGACGGCAAAGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCTGCGGCAGTCGAGGAGGTCGAATGATCCGCTCCATGCTGTACGCCAcc < 2:470104/150‑1 (MQ=255) cGAGGCTTCGGCGGGGCCCTCCAGCCCACGGCAGCAGCGCGCCGGACGGCAAAGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCTGCGGCAGTCGAGGAGGTCGAATGATACGCTCCATGCTGTACGCCACCGACCTc > 2:101258/1‑150 (MQ=255) tCCAGCCCACGGCAGCAGCGCGCCGGACGGCAAAGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCAGCGGCAga > 1:208844/1‑85 (MQ=255) tCCAGCCCACGGCAGCAGCGCGCCGGACGGCAAAGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCAGCGGCAga < 2:208844/86‑2 (MQ=255) ccAGCCCACGGCAGCAGCGCGCCGGACGGCAAAGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCTGCGGCAGTCGAGGAGGTCGAATGATCCGCTCCATGCTGTACGCCACCGACCTCTGTGTCTACGCCCCTTTCGt < 1:528267/150‑1 (MQ=255) ccAGCCCACGGCAGCAGCGCGCCGGACGGCAAAGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCTGCGGCAGTCGAGGAGGTCGAATGATCCGCTCCATGCTGTACGCCACCGACCTCGGTGTCTACGCCCCTTTCGt < 2:202256/150‑1 (MQ=255) aaaGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCTGCGGCAGTCGAGGAGGTCGAATGATCCGCTCCATGCTGTACGCCACCGACCTCGGTGTCTACGCCCCTTTCGTGATGCAACATGCGCTGGCCCTg > 1:826413/1‑142 (MQ=255) aaaGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCTGCGGCAGTCGAGGAGGTCGAATGATCCGCTCCATGCTGTACGCCACCGACCTCGGTGTCTACGCCCCTTTCGTGATGCAACATGCGCTGGCCCTg < 2:826413/142‑1 (MQ=255) tcgactGCCGCAGCGGCAGTCGAGGAGGTCGAATGATCCGCTCCATGCTGTACGCCACCGACCTCGGTGTCTACGCCCCTTTCGTGATGCAACATGCGCTGGCCCTGGCGCGCACGTTCGGTGCCGAGTTGTATGTGATCCACGCgg < 1:775997/143‑1 (MQ=255) tcgactGCCGCAGCGGCAGTCGAGGAGGTCGAATGATCCGCTCCATGCTGTACGCCACCGACCTCGGTGTCTACGCCCCTTTCGTGATGCAACATGCGCTGGCCCTGGCGCGCACGTTCGGTGCCGAGTTGTATGTGATCCACGCgg > 2:775997/5‑147 (MQ=255) tctGCCGCTGCGGCAGTCGAGGAGGTCGAATGATCCGCTCCATGCTGTACGCCACCGACCTCGGTGTCTACGCCCCTTTCGTGATGCAACATGCGCTGGCCCTGGCGCGCACGTTCGGTGCCGAGTTGTATGTGATCCACGCGGTCGAGc > 2:261223/1‑150 (MQ=255) | AAGCGCTGAGTTTCATCGACTGCAGCGACCGGCCGGCCTGAGGGTGGCCGAGGCTTCGGCGGGGCCCTCCAGCCCACGGCAGCAGCGCGCCGGACGGCAAAGCGAGCGTCCGGCGCGCTGCTAAGCTGAAGCTTGCTCTGCCGCTGCGGCAGTCGAGGAGGTCGAATGATCCGCTCCATGCTGTACGCCACCGACCTCGGTGTCTACGCCCCTTTCGTGATGCAACATGCGCTGGCCCTGGCGCGCACGTTCGGTGCCGAGTTGTATGTGATCCACGCGGTCGAGC > NC_002947/2653564‑2653849 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |