Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 40,391 | G→T | 11.1% | A248E (GCG→GAG) | mecA ← | penicillin binding protein 2 prime |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 40,391 | 0 | G | T | 11.1% | 118.3 / 5.0 | 36 | A248E (GCG→GAG) | mecA | penicillin binding protein 2 prime |
Reads supporting (aligned to +/- strand): ref base G (20/12); new base T (2/2); total (22/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.34e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01 |
CTTTTTATCGTAAAGTTTTTCGAGTCCCTTTTTACCAATAACTGCATCATCTTTATAGCCTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGTTTTAAGTGGAACGAAGGTATCATCTTGTACCCAATTTTGA > CP000730/40264‑40528 | cTTTTTATCGTAAAGTTTTTCGAGTCCCTTTTTACCAATAACTGCATCATCTTTATAGCCTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAgg > 1:171113/1‑141 (MQ=255) cTTTTTATCGTAAAGTTTTTCGAGTCCCTTTTTACCAATAACTGCATCATCTTTATAGCCTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAgg > 1:38141/1‑141 (MQ=255) cTTTTTATCGTAAAGTTTTTCGAGTCCCTTTTTACCAATAACTGCATCATCTTTATAGCCTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAgg > 1:199754/1‑141 (MQ=255) tCGTAAAGTTTTTCGAGTCCCTTTTTACCAATAACTGCATCATCTTTATAGCCTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTa < 1:30452/141‑1 (MQ=255) ttCGATTCCCTTTTTACCAATAACTGCATCATCTTTATAGCCTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATACCCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGt < 2:229380/141‑1 (MQ=255) tCGAGTCCCTTTTTACCAATAACTGCATCATCTTTATAGCCTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGtt < 2:1262/141‑1 (MQ=255) aaTAACTGCATCATCTTTATAGCCTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATg > 1:175433/1‑141 (MQ=255) aTAACTGCATCATCTTTATAGCCTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGa > 2:235946/1‑141 (MQ=255) atcatcTTTATAGCCTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAAtttttt > 2:14456/1‑141 (MQ=255) cTTTATAGCCTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGa > 1:205558/1‑141 (MQ=255) cTTTATAGCCTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGa > 1:138583/1‑141 (MQ=255) taGCCTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGGCCAACATAACCTAATAGATGTGAAGTCTCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCa > 2:160561/1‑141 (MQ=255) taGCCTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCTCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCa > 2:160571/1‑141 (MQ=255) cTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTa < 1:48262/141‑1 (MQ=255) cTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTa < 1:27796/141‑1 (MQ=255) ttATATTCTGTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTaaa < 2:103066/141‑1 (MQ=255) cTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGCGGATAGCTCAGACCATACGTTTCACCAGTTGTAAGATGAAATTTTTTTGCGACATCACTTAAATATTCAt > 1:145383/1‑141 (MQ=255) ttGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCAt > 1:220287/1‑141 (MQ=255) ttGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCAt > 1:67030/1‑141 (MQ=255) cttcAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACgg < 1:235946/141‑1 (MQ=255) cAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGttt > 2:89962/1‑141 (MQ=255) cAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCAATTTTTTTATCGGttt > 2:185115/1‑141 (MQ=255) gTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGTTTTAAg > 1:63520/1‑141 (MQ=255) aaCATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGTTTTAAGTGGAACGAAGGt > 1:27778/1‑141 (MQ=255) tAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGTTTTAAGTGGAACGAAGGTatca < 1:6407/141‑1 (MQ=255) aaTAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGTTTTAAGTGGAACGAAGGTATCATCTTGt > 1:86410/1‑141 (MQ=255) aaTAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTGGCGAAATCACTGAAACATGCATCCATTTTTTTAACGGTTTTAAGTGGAACGAAGGTATCATCTTGt < 2:145383/141‑1 (MQ=255) tAGATGTGAAGTCTCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGTTTTAAGTGGAACGAAGGTATCATCTTGTAc < 1:160571/141‑1 (MQ=255) tAGATGTGAAGTCTCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGTTTTAAGTGGAACGAAGGTATCATCTTGTAc < 1:160561/141‑1 (MQ=255) tgAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGTTATAAGTGGAACGAAGGTATCATCTTGTACCCAAtt > 1:245149/1‑141 (MQ=255) aaGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGTTTTAAGTGGAACGAAGGTATCATCTTGTACCCCAtttt > 2:137354/1‑141 (MQ=255) gTCGCTTTTCCTATAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGTTTTAAGTGGAACGAAGGTATCATCTTGTACCCAATTTTGa < 1:91326/141‑1 (MQ=255) gTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGTTTTAAGTGGAACGAAGGTATCATCTTGTACCCAATTTTGa > 1:61817/1‑141 (MQ=255) gTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGTTTTAAGTGGAACGAAGGTATCATCTTGTACCCAATTTTGa > 1:4684/1‑141 (MQ=255) gTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGTTTTAAGTGGAACGAAGGTATCATCTTGTACCCAATTTTGa < 1:141036/141‑1 (MQ=255) gTCGCTTTTCCTAGAGGAGAGTTACGACTTTCTGTTTCATTAGTTGGAAAAGGAAATTTTGTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGTTTTAAGTGGAACGAAGGTATCATCTTGTACCCAATTTTGa < 1:157656/141‑1 (MQ=255) | CTTTTTATCGTAAAGTTTTTCGAGTCCCTTTTTACCAATAACTGCATCATCTTTATAGCCTTTATATTCTTTTTGTTTTAATTCTTCAGAGTTAATGGGACCAACATAACCTAATAGATGTGAAGTCGCTTTTCCTAGAGGATAGTTACGACTTTCTGTTTCATTAGTTGTAAGATGAAATTTTTTTGCGAAATCACTTAAATATTCATCCATTTTTTTAACGGTTTTAAGTGGAACGAAGGTATCATCTTGTACCCAATTTTGA > CP000730/40264‑40528 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |