Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,850,810 | (C)5→4 | intergenic (+117/‑50) | sppA → / → ansA | protease IV (signal peptide peptidase)/cytoplasmic L‑asparaginase 1 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,850,806 | 0 | C | . | 98.4% | 252.5 / ‑2.5 | 63 | intergenic (+113/‑54) | sppA/ansA | protease IV (signal peptide peptidase)/cytoplasmic L‑asparaginase 1 |
Reads supporting (aligned to +/- strand): ref base C (0/1); new base . (30/32); total (30/33) | |||||||||||
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.15e-01 |
TGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGA > NC_000913/1850676‑1850899 | tGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCt < 2:1667657/139‑1 (MQ=255) tGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCt < 2:1193625/139‑1 (MQ=255) gcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCAc < 1:1323913/139‑1 (MQ=255) aCATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGtat > 2:1631700/1‑139 (MQ=255) gTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGct < 1:1589393/139‑1 (MQ=255) tACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCttt < 1:376764/137‑1 (MQ=255) tACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCttt > 2:376764/1‑137 (MQ=255) tGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGata < 1:1435475/139‑1 (MQ=255) aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑cccctgc < 2:765073/108‑4 (MQ=255) aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑cccctgc > 1:765073/1‑105 (MQ=255) aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCAc < 2:1064292/113‑1 (MQ=255) aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCAc > 1:1064292/1‑113 (MQ=255) aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtc > 1:947611/1‑139 (MQ=255) aGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtc > 2:50041/1‑139 (MQ=255) ggCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt < 2:947611/139‑1 (MQ=255) ccGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt < 1:360681/139‑1 (MQ=255) aCAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatat < 1:50041/139‑1 (MQ=255) aTCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGctc > 2:1093258/1‑117 (MQ=255) aTCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGctc < 1:1093258/117‑1 (MQ=255) ggCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATg < 1:1631700/139‑1 (MQ=255) acaTTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCaaa > 1:228486/1‑139 (MQ=255) caTTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCg > 2:752345/1‑111 (MQ=255) caTTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCg < 1:752345/111‑1 (MQ=255) tACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCt > 2:1675212/1‑77 (MQ=255) tACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCt < 1:1675212/77‑1 (MQ=255) aCGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCt < 2:929279/76‑1 (MQ=255) aCGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCt > 1:929279/1‑76 (MQ=255) gTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt > 2:376870/1‑106 (MQ=255) gTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt < 1:376870/106‑1 (MQ=255) cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCt < 1:198265/72‑1 (MQ=255) cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCt > 2:198265/1‑72 (MQ=255) cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGta > 1:1473393/1‑78 (MQ=255) cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGta < 2:1473393/78‑1 (MQ=255) cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt > 2:426471/1‑104 (MQ=255) cgacgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt < 1:426471/104‑1 (MQ=255) acgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTg > 1:158734/1‑139 (MQ=255) cgaATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGc < 2:158734/139‑1 (MQ=255) tctTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt > 2:323389/1‑94 (MQ=255) tctTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt < 1:323389/94‑1 (MQ=255) ttCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACAcggg > 1:1001691/1‑139 (MQ=255) gCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCt > 2:555781/1‑55 (MQ=255) gCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCt < 1:555781/55‑1 (MQ=255) gCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGa > 1:243902/1‑139 (MQ=255) gCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTAACGTTGCCTACAcgggcg > 2:1055239/1‑136 (MQ=255) gCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAAGTTACGTTGCCTACAcgggcg < 1:1055239/136‑1 (MQ=255) tGTTTCGCCCCATATTTCCTTATCCTCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGct < 2:1725026/71‑1 (MQ=255) tGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGct > 1:1725026/1‑71 (MQ=255) gTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCt < 1:1002654/52‑1 (MQ=255) gTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCt > 2:1002654/1‑52 (MQ=255) aTCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg > 2:971588/1‑98 (MQ=255) aTCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg < 1:971588/98‑1 (MQ=255) tCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGtata > 2:332056/1‑39 (MQ=255) tCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGtata < 1:332056/39‑1 (MQ=255) tCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg > 2:1768656/1‑97 (MQ=255) tCCGCATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg < 1:1768656/97‑1 (MQ=255) cATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcattca < 2:153958/61‑1 (MQ=255) cATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcattca > 1:153958/1‑61 (MQ=255) cATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt > 2:1177220/1‑59 (MQ=255) cATAGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt < 1:1177220/59‑1 (MQ=255) aGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt > 1:1594418/1‑56 (MQ=255) aGTATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcatt < 2:1594418/56‑1 (MQ=255) tATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTaca < 2:1053073/97‑1 (MQ=255) tATCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTaca > 1:1053073/1‑97 (MQ=255) tCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg > 1:778651/1‑86 (MQ=255) tCAGGTGCGCT‑‑CCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACg < 2:778651/86‑1 (MQ=255) | TGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCT‑CCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAAAAGAAATCAATTTACGTTGCCTACACGGGCGGGA > NC_000913/1850676‑1850899 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |