Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_002947 | 1,499,508 | +C | intergenic (+85/‑79) | trpS → / → zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,499,506 | 1 | . | C | 98.1% | 158.2 / ‑3.1 | 52 | intergenic (+83/‑81) | trpS/zapE | tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (22/29); total (23/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.42e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01 |
GCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑TCGGG‑CACG‑CCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCA > NC_002947/1499361‑1499639 | gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CAcg‑cc > 6:179437/1‑149 (MQ=37) gttaTGCGGCCATTCAGAGCCTTACCGAAGCCTTGCCACCACAGCAGGAGTGAGGCTGCCAAACGGTGGGAGGGGGCGTGCCCGCGAAACTGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCc‑aca < 8:105391/146‑1 (MQ=18) gTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCc‑a > 2:66192/1‑147 (MQ=25) agagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCAtt < 8:304972/148‑1 (MQ=21) cTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCgg < 3:218303/149‑1 (MQ=21) gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCa > 5:181071/1‑135 (MQ=21) gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCa < 6:181071/135‑1 (MQ=21) aGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCg < 1:66192/149‑1 (MQ=21) ccTTGGCACCACAGCATGTCTGAGGCTTCCAAACGGTGGGAGCTGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGATTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCt < 5:155262/149‑1 (MQ=21) ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTa > 5:350225/1‑148 (MQ=21) caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTg < 1:232526/149‑1 (MQ=21) cacaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCCTTCCCGCGAAACAAGCGCCGCGGTGCATGGCACCCGCTTTGCCGGTGT‑TAGCGGG‑GACGCCCGCTCCC‑CCAGACAGATTTCCAATAAGCCGGCCTATCGCGACAGTGACg > 6:352125/1‑149 (MQ=255) gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCCTGCCCGCCAAACACGCCCCCCGGTGCAAGGCACCCGCATTGCCGGTGT‑TAACGGGCCCCC‑CCCCTCCC‑ACAGACAGAAAGGCAATAAACGGGCCGAAAGCAACAATGAccggcc > 2:107054/1‑144 (MQ=255) agGACTGAGGCTGCCTATCGGGGGGAGCGGGCGTGCCGGCGAAACTGGCGCCGCGGTGCATGGCACCGGCTTGGCCTGTGT‑TCGGGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGcgcgg < 4:329850/149‑2 (MQ=21) aCTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCCCCGTGCATTGCACCCGCTTTGCCGGTGT‑TAGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCAATAAGCGGGCCGATCGCAAAAGTGTCGGCCCGgttt > 8:293932/1‑149 (MQ=18) cTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt > 5:322882/1‑149 (MQ=21) aGGCTGCCAACGGTTGGGAGCGGGCGTGCCCGCGAAACCGGCGCCGCGGGGCATGGCACCGGCTTTGCCGGTTT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt < 3:160093/147‑1 (MQ=21) gCTGCCTAACGGTGGGAGGGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGGGGG‑CACGCCCGCTCCCTACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt < 1:107054/149‑1 (MQ=21) gTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCAGGGCACCGGCTTTGCCGGTGTGC‑GCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACg < 4:58937/149‑1 (MQ=21) gTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc < 6:322882/149‑1 (MQ=37) ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt > 7:92093/1‑147 (MQ=37) ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt > 2:223838/1‑147 (MQ=37) cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt > 4:305147/1‑146 (MQ=37) gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc > 5:288427/1‑143 (MQ=37) gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc < 6:288427/143‑1 (MQ=37) aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt < 6:350225/149‑1 (MQ=37) cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTaa > 8:192370/1‑123 (MQ=25) cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTaa < 7:192370/123‑1 (MQ=25) ggCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc < 7:293932/149‑1 (MQ=37) gcgcCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAGTCATGACTCCCTTAGAACGCTATc > 4:124735/1‑148 (MQ=37) gccgcGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa > 5:309455/1‑136 (MQ=37) gccgcGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa < 6:309455/136‑1 (MQ=37) gccgcGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg > 8:314336/1‑149 (MQ=37) ccgcGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAGTCATGACTCCCTTAGAACGCTATCAAGc < 3:124735/149‑1 (MQ=37) cgcgGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAACATGACTCCCTTAGAACGCTATCAAGCa < 1:223838/149‑1 (MQ=37) aTGGCACCGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGa > 6:9755/1‑148 (MQ=38) tggcaccGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGaaa > 5:168638/7‑149 (MQ=38) ggcaccGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAAc > 7:6290/6‑149 (MQ=38) accGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTc < 7:314336/147‑1 (MQ=38) ccGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc < 8:6290/148‑1 (MQ=39) ccGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc > 7:244503/2‑149 (MQ=39) ccGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc > 6:101688/2‑149 (MQ=39) cGGCTTTGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTccc < 3:226507/149‑1 (MQ=39) caaagCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc < 5:281134/145‑1 (MQ=255) caaagCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGAc < 7:214302/145‑1 (MQ=255) tttGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAAc < 5:204639/140‑1 (MQ=39) tttGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAAc > 6:204639/1‑140 (MQ=39) tGCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActt < 3:48602/148‑1 (MQ=255) gCCGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtc < 4:49889/148‑1 (MQ=255) ccGGTGT‑TCGCGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtctt < 6:168638/149‑1 (MQ=255) gtgttcgcGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTACCTTAGAACGCTATAAAGCAGATCTGAAACGTCCCGACTTCTTCCa > 2:310096/8‑149 (MQ=255) gcGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt < 7:305076/73‑1 (MQ=255) gcGGG‑CACGCCCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt > 8:305076/2‑74 (MQ=255) | GCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑TCGGG‑CACG‑CCGCTCCC‑ACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCA > NC_002947/1499361‑1499639 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |