Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 6,049,349 | T→C | 10.9% | V52A (GTG→GCG) ‡ | spoT → | bifunctional (p)ppGpp synthase/hydrolase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 6,049,349 | 0 | T | C | 10.9% | 148.2 / 8.0 | 58 | V52A (GTG→GCG) ‡ | spoT | bifunctional (p)ppGpp synthase/hydrolase |
Reads supporting (aligned to +/- strand): ref base T (22/29); new base C (3/3); total (26/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01 |
CTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCA > NC_002947/6049212‑6049498 | cTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgggtg < 1:258538/150‑6 (MQ=255) ttGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTc > 2:1117951/1‑150 (MQ=255) ttGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTc < 2:549851/150‑1 (MQ=255) gAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGAc > 1:1242011/1‑150 (MQ=255) aaCGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACa < 2:1242011/150‑1 (MQ=255) aCGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACAt > 1:321856/1‑150 (MQ=255) ggCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGc > 2:1406604/1‑150 (MQ=255) tttCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGcaca > 1:616826/1‑150 (MQ=255) tttCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGCCGTGCCTATTTCTACCCCGAACAGGCCCACAACGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGcaca > 1:320384/1‑150 (MQ=255) aCCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCCCCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGAc > 1:683392/1‑150 (MQ=255) gTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGcc > 2:807088/1‑150 (MQ=255) aaCCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCCGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCAtg > 2:711340/1‑150 (MQ=255) aaCCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgggtg > 2:101855/1‑102 (MQ=255) aaCCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgggtg < 1:101855/107‑6 (MQ=255) ttCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACg > 2:1019531/1‑150 (MQ=255) gCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATc < 1:1019531/150‑1 (MQ=255) ccTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCAGATGGCGGCCATGCTGCACGATGTGATCg > 1:767061/1‑150 (MQ=255) tttCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATc < 2:737759/145‑1 (MQ=255) tttCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATc > 1:737759/1‑145 (MQ=255) ggCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCaa < 2:683392/150‑1 (MQ=255) gCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTg < 1:850791/60‑1 (MQ=255) gCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTg > 2:850791/1‑60 (MQ=255) ccACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGAtg > 1:1161017/1‑122 (MQ=255) ccACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGAtg < 2:1161017/122‑1 (MQ=255) ccACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCg < 1:1166054/137‑1 (MQ=255) ccACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCg > 2:1166054/1‑137 (MQ=255) ccACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGa > 2:183641/1‑150 (MQ=255) ccACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGAGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGa > 2:920427/1‑150 (MQ=255) aCGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAg < 2:1174290/150‑1 (MQ=255) cGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGc < 1:477735/150‑1 (MQ=255) gATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGcc < 1:371568/150‑1 (MQ=255) ggCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTca < 2:616826/150‑1 (MQ=255) ggCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTca < 1:881737/150‑1 (MQ=255) gCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTcag > 2:1060105/1‑150 (MQ=255) aCGCCGCCGGAGTGGCGAGCCCTACGGGACGCAACCGCTGGCCGTGGGCAGCATCGTGGCCGACATGCACATGGACAGTCAGAGGCGGATGGCGGCCATGCTGCACGATGTGATCGAGGTCACCGGCATCGCGAAGGAAGCCCTcagcca < 2:320384/150‑1 (MQ=255) aCGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTcagcca < 2:899288/150‑1 (MQ=255) gccgccgcAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTcagccagc > 1:250688/1‑150 (MQ=255) ccgccgcAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCcagca < 1:1060105/150‑1 (MQ=255) ccgcAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATggcg < 1:869740/88‑1 (MQ=255) ccgcAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATggcg > 2:869740/1‑88 (MQ=255) gCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCa < 2:434278/63‑1 (MQ=255) gCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCa > 1:434278/1‑63 (MQ=255) gCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGt < 2:946738/150‑1 (MQ=255) cccTACGTGACCCACCCGCTGGCCGTGGCCAGCATCACCGCACCAATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCgcg < 2:944719/150‑3 (MQ=255) cccccccCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGTTATTCTGCACGATGTGATCGAAGACACCGGCATCGCCAATGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGt < 1:734604/147‑1 (MQ=255) ccacccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGt < 2:250688/150‑1 (MQ=255) gatgCTGGCCACGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGa > 1:306777/4‑150 (MQ=255) cccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAATCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGa < 1:524049/150‑1 (MQ=255) cccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGa < 1:1406604/150‑1 (MQ=255) ccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGAt < 2:321856/150‑1 (MQ=255) tgCTGGCCACGGCCAGCATCCTCGCCGACATGCACATGGAc > 1:174442/2‑41 (MQ=255) tgCTGGCCACGGCCAGCATCCTCGCCGACATGCACATGGAc < 2:174442/40‑1 (MQ=255) gCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGAtgt < 1:1303107/75‑1 (MQ=255) gCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATgg > 1:546842/1‑150 (MQ=255) gCTGGACGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCTGCCATGCTGCACGAtgt > 2:1303107/1‑75 (MQ=255) cTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATTCTGCACTATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATggg < 1:807088/150‑1 (MQ=255) ggCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGt < 1:616749/150‑1 (MQ=255) ggCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGt < 1:1117951/150‑1 (MQ=255) tGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCa > 1:1132599/1‑150 (MQ=255) | CTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCA > NC_002947/6049212‑6049498 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |