Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,270,127 | T→C | L50P (CTG→CCG) | mepS → | peptidoglycan DD‑endopeptidase/peptidoglycan LD‑carboxypeptidase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,270,127 | 0 | T | C | 97.5% | 136.1 / ‑4.2 | 40 | L50P (CTG→CCG) | mepS | peptidoglycan DD‑endopeptidase/peptidoglycan LD‑carboxypeptidase |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (21/18); total (22/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01 |
AACCGATTTTGAGATATATCTTGCGCGGGATTCCCGCGATTGCAGTAGCGGTTCTGCTTTCTGCATGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACTGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTGGGCGGCAGCACTAAAAAAGGTATCGATTGTTCT > NC_000913/2269992‑2270263 | aaCCGATTTTGAGATATATCTTGCGCGGGATTCCCGCGATTGCAGTAGCGGTTCTGCTTTCTGCATGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCaa > 1:633013/1‑140 (MQ=255) aCCGATTTTGAGATATATCTTGCGCGGGATTCCCGCGATTGCAGTAGCGGTTCTGCTTTCTGCATGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAg > 1:578488/1‑140 (MQ=255) aCCGATTTTGAGATATATCTTGCGCGGGATTCCCGCGATTGCAGTAGCGGTTCTGCTTTCTGCATGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAg < 2:565281/140‑1 (MQ=255) tGCGCGGGATTCCCGCGATTGCAGTAGCGGTTCTGCTTTCTGCATGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGaa > 1:423964/1‑140 (MQ=255) gggATTCCCGCGATTGCAGTAGCGGTTCTGCTTTCTGCATGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACACCATCACCACAAGCTACTCAGGCTGAATTTGAAAACCt > 1:361333/1‑140 (MQ=255) ggATTCCCGCGATTGCAGTAGCGGTTCTGCTTTCTGCATGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTg < 2:422012/140‑1 (MQ=255) aTTGCAGTAGCGGTTCTGCTTTCTGCATGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGt > 2:374147/1‑140 (MQ=255) ttCTGCTTTCTGCATGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATcg < 1:405871/140‑1 (MQ=255) gCTTTCTGCATGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAACGTCGACGTAAACTCGCGaa < 2:423964/140‑1 (MQ=255) ttCTGCATGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTa < 2:578488/140‑1 (MQ=255) cTGCATGTAGTGCAAATAACACCGCTAAGTATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATg < 1:335511/140‑1 (MQ=255) cTGCATGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATg < 1:48153/140‑1 (MQ=255) cTGCATGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATg < 1:643289/140‑1 (MQ=255) tGCATGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATgg > 1:102264/1‑140 (MQ=255) aTGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATc < 1:374147/140‑1 (MQ=255) aTGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATc < 2:84344/140‑1 (MQ=255) gTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAAtt > 1:53336/1‑128 (MQ=255) gTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAAtt < 2:53336/128‑1 (MQ=255) tGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGCCGTGAATTCGCGACTTATGGATCAGTATg > 1:261349/1‑140 (MQ=255) cacCGCAAAGAATATGCATCCTCAGACACGTGCAGTGGGCAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGa < 2:228689/140‑1 (MQ=255) acCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACTGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCATAATGTCGGCGTAAAAACGCGAATTATGGAGCAGTATGCTGACTGGaa > 2:221404/1‑140 (MQ=255) cAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATg > 1:104760/1‑126 (MQ=255) cAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATg < 2:104760/126‑1 (MQ=255) cAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGc < 2:633013/140‑1 (MQ=255) aaaGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCg > 2:581033/1‑140 (MQ=255) cTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTg > 2:23431/1‑140 (MQ=255) cTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTg > 2:728812/1‑140 (MQ=255) cTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTg > 1:701625/1‑140 (MQ=255) cTGAGACACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTg > 2:650611/1‑140 (MQ=255) gagaCACGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTggg > 2:701399/1‑140 (MQ=255) cGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTGGGCGgcag > 1:328139/1‑140 (MQ=255) cGTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTGGGCGgcag < 2:102264/140‑1 (MQ=255) gTGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTGGGCGgcagc > 2:15226/1‑140 (MQ=255) tGCAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTGGGCGgcagca > 1:427970/1‑140 (MQ=255) cAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTGGGCGGCAGCACt > 1:658981/1‑140 (MQ=255) cAGTGGGTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTGGGCGGCAGCACt < 1:650611/140‑1 (MQ=255) gggTAGTGAAACATCATCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTGGGCGGCAGCACTaaaa < 2:40772/140‑1 (MQ=255) catcatCACCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGAAAAGGCGTCCGTGATCGTCTGGGCGGCAGCTCTAAAAAAGGTATCGAt > 2:209051/1‑140 (MQ=255) tcaCCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTGGGCGGCAGCACTAAAAAAGGTATCGATTGTTc < 1:701399/140‑1 (MQ=255) tcaCCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGACGTACGTTTTCGTCAGTCCGGCAGCTCTAAGAAAGGTATAGATTGTTc > 2:360738/1‑140 (MQ=255) caCCGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTGGGCGGCAGCACTAAAAAAGGTATCGATTGTTCt > 1:546465/1‑140 (MQ=255) | AACCGATTTTGAGATATATCTTGCGCGGGATTCCCGCGATTGCAGTAGCGGTTCTGCTTTCTGCATGTAGTGCAAATAACACCGCAAAGAATATGCATCCTGAGACACGTGCAGTGGGTAGTGAAACATCATCACTGCAAGCTTCTCAGGATGAATTTGAAAACCTGGTTCGTAATGTCGACGTAAAATCGCGAATTATGGATCAGTATGCTGACTGGAAAGGCGTACGTTATCGTCTGGGCGGCAGCACTAAAAAAGGTATCGATTGTTCT > NC_000913/2269992‑2270263 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |