Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,524,063 | T→A | 18.3% | V398E (GTG→GAG) | mrcA → | peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,524,063 | 0 | T | A | 18.3% | 34.9 / 2.6 | 22 | V398E (GTG→GAG) | mrcA | peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA |
Reads supporting (aligned to +/- strand): ref base T (7/8); new base A (2/2); total (9/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01 |
GCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCG > NC_000913/3523858‑3524277 | accgaaaaaaacaacaacAGTACCGAAACACTGAAGGCGCGGCAAACCTATGGTCCGCTGCGGCCGGCCGCGATCACCAGGCCCAATCCTCAGCAAGCTAAGGCGAGGCTGGCGGACGGGTCGACCGTCGCAGTGAGTATGGAAGGCGTACGCTGGGCGCGTCCCTAACGTTCGGATACTCAGCTAGGACCGACGCCGCGTAAACTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGTTTCGTCAGGTTGGCGATGCATGGTGGCTGGcaca < 2:91786/266‑1 (MQ=255) tgtgATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCGGCCGCAGGCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGGCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAGGGGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAg < 1:210274/271‑1 (MQ=255) aGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGAGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGtggctgg > 1:91786/1‑255 (MQ=255) gacaggCTGGAGGCGCTGCCAACCGTGGGCCCGCGGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCGGCAAGCGACGGCGATGCGGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGgcgc < 2:51885/269‑1 (MQ=255) aGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGt < 1:169909/172‑1 (MQ=255) tGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCg > 1:51885/1‑254 (MQ=255) cAACCTATGGTCGCCTGCCGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGAGGCGGGCGGACGGGTCGACCGTCGCATTGAGTATAGAAGACGTTCGCTGGGCGCGTCATTACCGTCCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAg < 1:8312/188‑1 (MQ=255) tATGTTCCTCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACgg < 1:169628/275‑1 (MQ=255) tGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGAAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAGGTGACATGGGCGCTGGTGTTGATCAATCCGCAAAACGGTg > 1:57171/1‑275 (MQ=255) ctgctgCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCCCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCCAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCCCAAGTGCCGGAAGTGAAATAGCCCCTGGTGTCGACACATCCCGCAAACGGGGCCCTTAt > 1:218329/1‑275 (MQ=255) gctgcCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGAGGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGAc < 2:101189/150‑1 (MQ=255) ggggggCGGGGGGGGCGACCGGAGCGTGGAGTAGGGAAGGGCTGCGCTGGGCGCGTCCTTACCGTTCGGATACGCAGCGAGGACCGACGCCGCGTAAAGGGACCGATGTTCTGCAAACGGGTCAGCAAGTCTGGGTTCGTCAGGTTGTCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCgcgc < 2:218329/270‑1 (MQ=255) cGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCGTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGGGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACgg < 1:99276/207‑1 (MQ=255) tGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCAGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTa < 1:75602/192‑1 (MQ=255) gCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATTCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTAAATCAGAGCAAGTTTAACCGCGCACCCCAGGCCCTGCGTCAGGTGGGTTCCAACATCAAACCg > 1:13995/1‑275 (MQ=255) ttCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGtggct > 2:92285/1‑99 (MQ=255) ttCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTg > 1:98497/1‑118 (MQ=255) gATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGca < 1:92285/99‑1 (MQ=255) gATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACt < 2:98497/118‑1 (MQ=255) cgTAAAGAGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGtggctgg > 1:176705/1‑73 (MQ=255) cgTAAAGAGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGtggctgg > 1:176704/1‑73 (MQ=255) aaGAGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGcaca < 2:176704/73‑1 (MQ=255) aaGAGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGcaca < 2:176705/73‑1 (MQ=255) | GCGTGGGATAACAACAAGATTACCGATACGCTGAAGGCGCTGCCAACCTATGGTCCGCTGCTGCCTGCCGCAGTCACCAGCGCCAATCCTCAGCAAGCGACGGCGATGCTGGCGGACGGGTCGACCGTCGCATTGAGTATGGAAGGCGTTCGCTGGGCGCGTCCTTACCGTTCGGATACTCAGCAAGGACCGACGCCGCGTAAAGTGACCGATGTTCTGCAAACGGGTCAGCAAATCTGGGTTCGTCAGGTTGGCGATGCATGGTGGCTGGCACAAGTGCCGGAAGTGAACTCGGCGCTGGTGTCGATCAATCCGCAAAACGGTGCCGTTATGGCGCTGGTCGGTGGCTTTGATTTCAATCAGAGCAAGTTTAACCGCGCCACCCAGGCACTGCGTCAGGTGGGTTCCAACATCAAACCG > NC_000913/3523858‑3524277 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |