Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,807,986 | G→A | 15.1% | intergenic (+121/+34) | 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,807,986 | 0 | G | A | 15.1% | 72.1 / 6.0 | 33 | intergenic (+121/+34) | mrcA/maeB | penicillin‑insensitive transglycosylase/penicillin‑sensitive transpeptidase/malic enzyme B |
Reads supporting (aligned to +/- strand): ref base G (16/12); new base A (2/3); total (18/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.39e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01 |
GACGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACG > NC_002947/5807844‑5808128 | gACGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGCGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGcc < 8:325084/149‑1 (MQ=255) gCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAATCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑ggg > 2:321944/1‑149 (MQ=255) gCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCGCATGAGCC‑ggg < 1:204237/149‑1 (MQ=255) gCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAACCCCCGGCTCGCATGAGCC‑ggg < 2:352108/149‑1 (MQ=255) gTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑g > 1:346781/1‑135 (MQ=255) gTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑g < 2:346781/135‑1 (MQ=255) ttCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCtttttt < 6:82835/144‑1 (MQ=255) tAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCGCATGAGCC‑gggg > 6:221783/1‑134 (MQ=255) tAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCGCATGAGCC‑gggg < 5:221783/134‑1 (MQ=255) aGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTg < 2:303312/147‑1 (MQ=255) aTCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTGATGCGACCC‑GGGGCTTTTTTATGTCGCATGCTCCAATCAGCCGTTGCACACTTCATCCAAGCTTGt > 7:179416/1‑147 (MQ=255) cTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGc > 6:28106/1‑149 (MQ=255) tGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCg > 4:7888/1‑149 (MQ=255) tAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGt < 1:321944/149‑1 (MQ=255) gcgcCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAacac > 7:103880/1‑123 (MQ=255) gcgcCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAacac < 8:103880/123‑1 (MQ=255) gcgcCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAAt > 5:41227/1‑149 (MQ=255) cgcCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATg < 6:41227/149‑1 (MQ=255) cccTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAg > 1:377589/1‑149 (MQ=255) gACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGtt < 7:15090/101‑1 (MQ=255) gACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCAAGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGtt > 8:15090/1‑101 (MQ=255) acacATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGa < 5:28106/148‑1 (MQ=255) acaTTCCCTAGCACAATATAAAAGCCCCGGCTCATGCGATCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATGATCCGTTGAAAACTTCATCCACGCTGTTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATg < 8:179416/148‑1 (MQ=255) acaTTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGCGCCGGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCg > 1:282070/1‑117 (MQ=255) aGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTgg < 2:377589/149‑1 (MQ=255) aTAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGGTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAg > 7:165559/1‑149 (MQ=255) aaGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGAt > 4:141475/1‑126 (MQ=255) ccccGGCTCGCATGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGaaa > 2:258087/1‑149 (MQ=255) ccccGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGaaa > 6:90753/1‑149 (MQ=255) ccGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAac > 8:327379/1‑148 (MQ=255) gCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAacgacg > 4:204622/1‑149 (MQ=255) cTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGc > 4:297213/1‑129 (MQ=255) cTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGc < 3:297213/129‑1 (MQ=255) | GACGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCC‑GGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACG > NC_002947/5807844‑5808128 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |