Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,508,089 | A→C | 12.5% | intergenic (+19/+43) | sspB → / ← PP_1322 | ClpXP protease specificity‑enhancing factor/lipoprotein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,508,089 | 0 | A | C | 12.5% | 93.3 / 6.2 | 40 | intergenic (+19/+43) | sspB/PP_1322 | ClpXP protease specificity‑enhancing factor/lipoprotein |
Reads supporting (aligned to +/- strand): ref base A (19/16); new base C (3/2); 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.38e-01 |
CCGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCCAGCAG > NC_002947/1507945‑1508236 | ccGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCggatggat < 1:661049/150‑1 (MQ=255) cTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGcc < 1:746733/150‑1 (MQ=255) tggatggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATTACGATGGTCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑tttttt < 2:81469/150‑1 (MQ=255) tggatggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCATCCGGATCGCC‑tttttt < 1:735407/150‑1 (MQ=255) atggCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTtttgtt < 1:230828/150‑1 (MQ=255) gCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTgg < 1:948201/150‑1 (MQ=255) gCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTgg > 2:102185/1‑150 (MQ=255) acgaGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAACGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCt > 2:57430/1‑150 (MQ=255) tGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAg > 1:1015575/1‑150 (MQ=255) tGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAAC‑AAAAAAGGCGATCCATCCGGATCGCC‑tttttttgtta > 1:821766/1‑134 (MQ=37) tGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAAC‑AAAAAAGGCGATCCATCCGGATCGCC‑tttttttgtta < 2:821766/139‑6 (MQ=37) gAAGACGAGGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAGGTGGTCAGGTAACAAAAAAAGGCGATCCGGATGGATCGAC‑TATTTTCTTTGGCCTGTCCAGGc > 2:565518/1‑150 (MQ=255) aCGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCAAAGGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATc > 1:900222/1‑150 (MQ=255) gatgatGTGCAGCCGGATGACGATGGCACGCCAGACGATGGAGGGCAGCCGCCGCGTACGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTGGTTTGGCCTGTCCAGGCACTCAc > 2:119803/1‑150 (MQ=255) tgatgTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTg > 2:393996/1‑150 (MQ=255) gCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCaa < 1:186853/150‑1 (MQ=255) aGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATc < 1:152578/150‑1 (MQ=255) gATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTTCAGGCAATCACTGAGGCAATCGATGTACTCGa > 1:692665/1‑150 (MQ=255) aTGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGaa < 2:23168/150‑1 (MQ=255) gACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCggatgg > 1:37845/1‑82 (MQ=255) gACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCggatgg < 2:37845/82‑1 (MQ=255) tGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGAtcttc < 1:565518/150‑1 (MQ=255) gccggccAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGcaccac > 1:794087/1‑150 (MQ=255) cggccAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACgg < 2:1015575/150‑1 (MQ=255) cAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGt > 2:836859/1‑150 (MQ=255) gAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAg > 2:34226/1‑150 (MQ=255) aaGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGtt < 1:102185/150‑1 (MQ=255) gTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTgg < 1:34226/150‑1 (MQ=255) tGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGc > 2:751844/1‑150 (MQ=255) aaGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGcttctt > 1:915470/1‑150 (MQ=255) aaacaaaaaaGGCGATCCATCCGGATCGCCTTTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGGACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGAc > 1:531159/5‑150 (MQ=38) aaacaaaaaaGGCGATCCATCCGGATCGCCTTTTTTTGATTGGGCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGAc > 1:985680/5‑150 (MQ=38) aCAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACgg < 1:836859/150‑1 (MQ=255) aCAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACgg > 1:349126/1‑150 (MQ=255) aCAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTATTTGAGCTGTCTAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCAGGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACgg < 1:119803/150‑1 (MQ=255) cAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACgg > 2:264689/1‑149 (MQ=255) cAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACgg < 1:264689/149‑1 (MQ=255) aaaGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGAcc > 2:370000/1‑150 (MQ=255) aaaGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGAcc < 1:393996/150‑1 (MQ=255) aaGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGa < 2:380561/81‑1 (MQ=255) aaGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGa > 1:380561/1‑81 (MQ=255) gCGATCCGGATGGATCGCC‑TTTTTTGTTTGGTCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGGTTGACAATCTTCTGCACTCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTACCGGGTGACCAGGc > 2:74001/1‑150 (MQ=255) aTCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACAAGGCCca > 1:382279/1‑150 (MQ=255) ggatggatCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCcagcag < 1:370000/150‑1 (MQ=255) ggatggatCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCcagcag > 2:1078058/1‑150 (MQ=255) ggatggatCGCC‑TTTTTTGTTTGGCCTGTCCAGGACCCGACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCcagcag < 1:57430/150‑1 (MQ=255) | CCGTCTTCGATGGATGGCGACGAGCTGGAAGACGATGATGTGCAGCCGGATGACGATGGCCCGCCAGACGATGGAGGGCAGCCGCCGCGTCCGACTGGCCGGCCAAGCCTGAAAGTGGTCAAGTAACAAAAAAAGGCGATCCGGATGGATCGCC‑TTTTTTGTTTGGCCTGTCCAGGCAATCACTGAGTCAATCGATGTACTCGAACAGCTTGACGATCTTCTGCACGCCCGAAACACCCTGCACCACGGCAGTGGCCGAGTTGGCTTCTTGACGGGTGACCAGGCCCAGCAG > NC_002947/1507945‑1508236 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |