Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_002947 | 5,378,332 | +A | intergenic (‑61/‑132) | grpE ← / → recN | heat shock protein GrpE/DNA repair/recombination protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,378,331 | 1 | . | A | 97.6% | 137.3 / ‑4.7 | 42 | intergenic (‑60/‑133) | grpE/recN | heat shock protein GrpE/DNA repair/recombination protein |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base A (19/22); total (20/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.76e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑T‑ATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGAGCCGACCCCATGCTGG > NC_002947/5378182‑5378470 | cTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCta‑t > 2:1037978/1‑150 (MQ=255) ccAGAACACGGGTATCAGCTGCATTGCACGCACTGGCTTCCTATCTGTCCAGAATGTTCTCGTTCAGCTGCTACTTAGCTATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGc < 1:762766/150‑1 (MQ=255) gTACCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCa > 2:561509/1‑150 (MQ=255) tATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAg < 1:929766/150‑1 (MQ=255) gCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGTGAGCGTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTC‑TAATGGCGCCGGAAAAACCAGGTTCa > 1:333589/1‑150 (MQ=255) gCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCta‑t‑aatg > 1:613284/1‑126 (MQ=255) gCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCta‑t‑aatg < 2:613284/129‑4 (MQ=255) tGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAg < 1:1037978/150‑1 (MQ=255) gCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTa < 1:326990/150‑1 (MQ=255) gTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTc < 1:624752/150‑1 (MQ=255) ttcCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGcc < 2:580523/111‑1 (MQ=255) ttcCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGcc > 1:580523/1‑111 (MQ=255) ttcCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGAc > 1:550891/1‑150 (MQ=255) cACGTTGAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGa < 2:567044/146‑1 (MQ=255) cACGTTGAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGa > 1:567044/1‑146 (MQ=255) cACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGA‑GCGGTGTAGTTTTCAAGACTGc > 1:322831/1‑150 (MQ=255) aCGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAg > 1:492245/1‑118 (MQ=255) aCGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAg < 2:492245/118‑1 (MQ=255) cGTTAAGGTCTTTCTCTTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCACAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCaa < 1:305633/150‑1 (MQ=255) gTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAg > 2:872877/1‑117 (MQ=255) gTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAg < 1:872877/117‑1 (MQ=255) ttCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAgggg > 2:664735/1‑150 (MQ=255) cAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCa < 2:888225/150‑1 (MQ=255) cgatctatcaGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGGTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTCGTCAAGACTGCAACATAATCTGCGAAGGGGCATTg < 1:612836/144‑1 (MQ=255) ctCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCa < 1:169250/150‑1 (MQ=255) tcaGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGGGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGACTATGAGACCTGTACATAATCTGCGAAGGGGCATTggagatcg > 2:612836/1‑144 (MQ=255) aGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTg < 1:664735/150‑1 (MQ=255) ccATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaa > 2:506678/1‑150 (MQ=255) aTGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaaa > 2:406333/1‑150 (MQ=255) tGGGGTCTCTCCTGCGCAATTTTGTGGACTACAGGAGAGGCTCGCTATTGATTTCGTCTACA‑ATGGGGCTGG‑AAAACCAGCTGTAGTGGTAGAAGTGAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTAGCTGGCCTGaaaaa > 1:504757/1‑150 (MQ=255) gggTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAg < 1:561509/150‑1 (MQ=255) tctctcCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGGTA‑TAATGGGGCCGGGAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTAc > 1:889812/1‑150 (MQ=255) cgctcCTGCGCAATTTGGTGGACTGCGAGCCAGGCACGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACt < 2:328299/148‑1 (MQ=255) cAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTataaataa < 1:506678/150‑1 (MQ=255) ttttGTGGACTGTAAGCGAGGCTCGCCTTGGATGTCCGCTA‑TAATATGGACTGCAACAACAAGATCGAGAG‑GCAAGGGAGTTGTGAAGACTGGTACATAATCTGGTAAGGGGCATTGGCAGGCCTGAAAAAAATAAAATAAAAATAACCa > 2:419669/1‑150 (MQ=255) tgtgGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGAc < 2:550891/150‑1 (MQ=255) tGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGAttt > 1:282730/1‑150 (MQ=255) gAGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCg < 2:504757/150‑1 (MQ=255) aGCCAGGCTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCgg > 1:116243/1‑150 (MQ=255) cTCGCCTTGGATGTCGGCTA‑TAATGGGGCCTGAAAAACCAGGTTCAAGGG‑GCAAGGTAGCTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGAGCCGa > 2:322094/1‑150 (MQ=255) cTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAggg < 1:639283/90‑1 (MQ=255) cTCGCCTTGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAggg > 2:639283/1‑90 (MQ=255) ttGGATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGAGCCGACCCCAt < 2:322831/150‑1 (MQ=255) ggATGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGAGCCGACCCCATGc < 1:406333/150‑1 (MQ=255) tGTCGGCTA‑TAATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGAGCCGACCCCATGCTgg > 1:645934/1‑150 (MQ=255) | CTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTA‑T‑ATGGGGCCGGAAAAACCAGGTTCAAGGG‑GCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGAGCCGACCCCATGCTGG > NC_002947/5378182‑5378470 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |