| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,807,987 | C→T | 21.1% | intergenic (+122/+33) | 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,987 | 0 | C | T | 21.1% | 62.1 / 11.0 | 34 | intergenic (+122/+33) | mrcA/maeB | penicillin‑insensitive transglycosylase/penicillin‑sensitive transpeptidase/malic enzyme B |
| Reads supporting (aligned to +/- strand): ref base C (13/13); new base T (3/4); total (17/17) | |||||||||||
| 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.96e-01 | |||||||||||
GACGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTGAT > NC_002947/5807844‑5808133 | gACACCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGc < 5:17924/148‑1 (MQ=255) cGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCgg < 8:27527/149‑1 (MQ=255) gCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCCGCCCTGACACATTCCCAAGCCCCATCAAAAAGCCCCGGCTCATGCGAGCCggg < 5:80743/149‑1 (MQ=255) gCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCGACATAAAAAAGCCCCGGCTCGCATGAGCCggg < 1:383678/149‑1 (MQ=39) gCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCGCATGAGCCggg < 6:143759/149‑1 (MQ=255) gCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCGCATGAGCCggg < 3:338609/149‑1 (MQ=255) gCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCggg < 6:374140/149‑1 (MQ=255) ccGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCgggg < 1:427212/149‑1 (MQ=255) cGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCGCATGAGCCGGGGc < 3:318347/149‑1 (MQ=255) aTGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CATGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCtt < 7:276630/149‑1 (MQ=255) ttCTAAGGCAGGGCGCCTGACTCAGTTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTTACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCtttttt < 2:269833/144‑1 (MQ=255) tAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCt > 2:363609/1‑149 (MQ=255) ggCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTtgc > 7:264562/1‑149 (MQ=255) gCAGGGCGCCTGACTCAGGTTTTGCGGGGGCGCTCGATCTTATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGTCTCATGCGAGCCGGGGCTTTTTTATGTCGCTtgct < 4:272976/149‑1 (MQ=255) tCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCAGCTGACCAACACATTGCAGCCATGACACATACCAAAGCCCCATAAAAAACCCCCGCAACATGCGAGCCGGGCCTTTTTAATTACCATTGATGCAAACAGCCttt > 1:308137/1‑149 (MQ=255) ttGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCAGGCGACCCGGGGCGTTTTTATGTCGCGTGCTGCAATCAGCCGTTGAacacgt > 1:386813/1‑147 (MQ=255) tGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCACTACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATcc < 1:363609/149‑1 (MQ=255) cTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTc > 7:403159/1‑149 (MQ=255) acacATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTaa < 4:348052/111‑1 (MQ=255) acacATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTaa > 3:348052/1‑111 (MQ=255) acacATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGAt < 8:403159/149‑1 (MQ=255) cacaTTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATg < 8:264562/149‑1 (MQ=255) acaTTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTTATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATgg > 5:372529/1‑149 (MQ=255) ccATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCcac > 3:205525/1‑149 (MQ=255) aaaaaGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGGGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCAt > 5:257246/1‑149 (MQ=255) aaaaGCCCCGGCTCGCATGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATc > 5:38063/1‑149 (MQ=255) aaaGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCg > 1:192684/1‑149 (MQ=255) aGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGAtt > 7:25996/1‑121 (MQ=255) aGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGAtt < 8:25996/121‑1 (MQ=255) gCCCCGGCTCGCATGAGCCGGGGCTTTTTTATGGGGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGa > 7:51971/1‑148 (MQ=38) ccccGGCTCGCATGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGaaa > 6:421791/1‑149 (MQ=255) cccGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAac > 4:353005/1‑149 (MQ=255) ccGGCTCATGATAGCCGGGGCTTTTTAATGTCGCTTGCTGCAATCAGCCGTTGAACACTTAATCCACGCTCGTCACCGGGTAAAGCTTCGCATAAGGCAGCGTAGCAACACCGTATTCGAAGGCATACTTGCCCACAGCATCGAAAac > 8:90013/1‑148 (MQ=255) gCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTGAt > 5:66451/1‑148 (MQ=255) | GACGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTGAT > NC_002947/5807844‑5808133 |
| 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: ‑ |