Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,807,985 | T→C | 21.7% | intergenic (+120/+35) | 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,985 | 0 | T | C | 21.7% | 54.6 / 12.3 | 33 | intergenic (+120/+35) | mrcA/maeB | penicillin‑insensitive transglycosylase/penicillin‑sensitive transpeptidase/malic enzyme B |
Reads supporting (aligned to +/- strand): ref base T (12/13); new base C (3/4); total (16/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 = 5.97e-01 |
GACGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAAC > NC_002947/5807844‑5808124 | 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) | GACGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGC‑CCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAAC > NC_002947/5807844‑5808124 |
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: ‑ |