Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,808,040 | A→G | 26.7% | intergenic (+119/+36) | mrcA → / ← maeB | penicillin‑insensitive transglycosylase/penicillin‑sensitive transpeptidase/malic enzyme B |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,808,040 | 0 | A | G | 26.7% | 18.3 / 5.3 | 15 | intergenic (+119/+36) | mrcA/maeB | penicillin‑insensitive transglycosylase/penicillin‑sensitive transpeptidase/malic enzyme B |
Reads supporting (aligned to +/- strand): ref base A (6/5); new base G (2/2); total (8/7) | |||||||||||
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.60e-01 |
GACGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCAT‑AAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTG > NC_002947/5807900‑5808187 | gACGCCGATGGATCTGTTCTAACGCCGGGCGTCTGACTCAGGTTGTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTACCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCAT‑AAAAAAGCCCCGGCTCATGCGAGCCg < 2:544759/150‑1 (MQ=255) tGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCAT‑AAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGt > 1:480853/1‑150 (MQ=255) aaGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACGTTCCCAAGCCCCAT‑GAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTtg < 2:137573/150‑1 (MQ=255) cgcTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCAT‑AAAAAAGCCCCGGCTCGCATGAGCCGGGGCTTTTTTATg > 1:530279/1‑110 (MQ=37) cgcTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCAT‑AAAAAAGCCCCGGCTCGCATGAGCCGGGGCTTTTTTATg < 2:530279/110‑1 (MQ=37) tAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTg < 1:379339/140‑1 (MQ=255) tAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTg > 2:379339/1‑140 (MQ=255) aaGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCAT‑AAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGa < 1:62513/150‑1 (MQ=255) aCCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCAT‑AAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAggg < 2:480853/150‑1 (MQ=255) cAACACCTTGCAGCCCTGACACATTCCCAAGCCCCAT‑AAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTa > 2:53356/1‑150 (MQ=255) cAT‑AAAAAAGCCCCGGCTCGCATGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGc < 1:192940/118‑1 (MQ=37) cAT‑AAAAAAGCCCCGGCTCGCATGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGc > 2:192940/1‑118 (MQ=37) cAT‑AAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGCACCCCCGCCGGATTCGATGCCAGCCCTGccccccc > 1:13488/1‑146 (MQ=255) aaGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGa > 2:266573/1‑150 (MQ=255) tCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTg > 1:122653/1‑150 (MQ=255) | GACGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCAT‑AAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTG > NC_002947/5807900‑5808187 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |