Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 6,170,990 | G→A | 7.5% | T456T (ACC→ACT) | glmS ← | L‑glutamine/D‑fructose‑6‑phosphate aminotransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 6,170,990 | 0 | G | A | 7.5% | 136.1 / 2.1 | 53 | T456T (ACC→ACT) | glmS | L‑glutamine/D‑fructose‑6‑phosphate aminotransferase |
Reads supporting (aligned to +/- strand): ref base G (25/24); new base A (2/2); total (27/26) | |||||||||||
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.94e-01 |
TCGCCTGCTGGGTATGCTTCGGCGTGGATGTAGGAGATTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGACCAAGGGCCAGGGTCAGCAGCATCAGCGAAACCAGCTGAGTGGT > NC_002947/6170847‑6171136 | tCGCCTGCTGGGTATGCTTCGGCGTGGATGTAGGAGATTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACAGTGGCa > 2:195798/1‑150 (MQ=255) tCGCCTGCTGGGTATGCTTCGGCGTGGATGTAGGAGATTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACAGTGGCa > 2:195871/1‑150 (MQ=255) tCGCCTGCTGGGTATGCTTCGGCGTGGATGTAGGAGATTTCTTTAAGCTTGAGAGCACATTCCATCGCCACCGGGAACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCa > 1:410959/1‑150 (MQ=255) ctgctgGGTATGCTTCGGCGTGGATGTAGGAGATTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCa > 1:216984/1‑150 (MQ=255) tgctgGGTATGCTTCGGCGTGGATGTAGGAGATTTCTTTAAGCTTGAGAGCACCTTCCGTCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCaa < 2:17135/150‑2 (MQ=255) ctgGGTATGCTTCGGCGTGGATGTAGGAGATTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATgg > 2:833211/1‑150 (MQ=255) ctgGGTATGCTTCGGCGTGGATGTAGGAGATTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACGGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATgg < 1:1072326/150‑1 (MQ=255) tgGGTATGCTTCGGCGTGGATGTAGGAGATTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCGCGGCAATTTTCTGGACGGTGGCATCCATGGc > 2:217939/1‑150 (MQ=255) tATGCTTCGGCGTGGATGTAGGAGATTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAgg > 2:356445/1‑150 (MQ=255) tGCTTCGGCGTGGATGTAGGAGATTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGc < 2:988680/150‑1 (MQ=255) tGCTTCGGCGTGGATGTAGGAGATTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGc < 2:410959/150‑1 (MQ=255) tGCTTCGGCGTGGATGTAGGAGATTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGCAACCGCGTATGGTGCTTGTCGCCGAACAGCTCGGCAATTGTCTCGACGGTGGCATCCATGGCCAGGGc > 2:386246/1‑150 (MQ=255) ggATGTAGGAGATTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGc > 1:701611/1‑150 (MQ=255) tAGGAGATTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTa < 2:132063/150‑1 (MQ=255) aGGAGATTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAg < 2:505749/150‑1 (MQ=255) gagaTTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGc > 2:186790/1‑150 (MQ=255) gaTTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTc < 1:67431/132‑1 (MQ=255) gaTTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTc > 2:67431/1‑132 (MQ=255) gaTTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACCGCGCGCCGTGTCGGGGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAg > 2:431340/1‑150 (MQ=255) tttAAGCTTGAGAGCACCTTCCATCGCCNCCGGGTACTGCGCGCNGNGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCGATTTTCTCGACGGTGGCATTTTTAGTCAGGGCTTCGCCCAGGCGGGTAGGCAGTCNTCg > 1:608337/1‑150 (MQ=255) tttAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACg < 1:452128/150‑1 (MQ=255) tttAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACg > 1:127411/1‑150 (MQ=255) tttAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGGAGGCAGTCGACg > 1:1075528/1‑150 (MQ=255) cTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGctc > 1:832697/1‑150 (MQ=255) cTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGctc > 2:945948/1‑150 (MQ=255) tGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCtt > 2:396693/1‑150 (MQ=255) gagaGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGc < 2:404013/119‑1 (MQ=255) gagaGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGc > 1:404013/1‑119 (MQ=255) tCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCtt < 2:995359/150‑1 (MQ=255) tCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCtt > 1:1071051/1‑150 (MQ=255) cGCTACTGGGTACTCCGCGCCGCGGCCGAGGAACAGGGTACAATACATATCGTCCGATTTTTAATCAATTTGCTCGAAGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTc < 1:406970/150‑1 (MQ=255) cGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTc < 2:473442/150‑1 (MQ=255) cACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGAc > 1:889232/1‑122 (MQ=255) cACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGAc < 2:889232/122‑1 (MQ=255) gTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGc < 1:1015471/150‑1 (MQ=255) ggCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTAccac > 2:872664/1‑148 (MQ=255) aaCAGGGTATGGTGCCTGTAGGCGAACCGCATGTAAGATTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTCCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGa < 2:675316/150‑1 (MQ=255) ggTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGACCAAGGGCCa > 1:538922/1‑150 (MQ=255) gTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGACCAAGGGCCAg > 1:113432/1‑150 (MQ=255) ggCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGACCAAGGGCCAgg < 1:431340/149‑1 (MQ=255) cGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCAAAATTTACCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGACCAAGGGCCAGGGTcagca < 2:820242/150‑1 (MQ=255) gCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGACCAAGGGCCAGGGTCagcagc > 2:47771/1‑150 (MQ=255) gCGAACAGCTCGGCAATTTTCTCGACAGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGACCAAGGGCCAGGGTCagcagc < 1:195871/150‑1 (MQ=255) gCGAACAGCTCGGCAATTTTCTCGACAGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGACCAAGGGCCAGGGTCagcagc < 1:195798/150‑1 (MQ=255) aaCAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCACACCTGACCAAGGGCCAGGGTCAGCAGCATc > 2:983301/1‑150 (MQ=255) aGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTGGGCAGTGGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGACCAAGGGCCAGGGTCAGCAGCATCAGc > 2:799322/1‑150 (MQ=255) aGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGACCAAGGGCCAGGGTCAGCAGCATCAGc < 1:396693/150‑1 (MQ=255) aGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGACCAAGGGCCAGGGTCAGCAGCATCAGc < 2:1075528/150‑1 (MQ=255) cGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCACGCACCTGACCAAGGGCCAGGGTCAGCAGCATCAGCGaaa < 1:872664/150‑1 (MQ=255) ttttCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGACCAAGGGCCAGGGTCAGCAGCATCAGCGAAACCAGCt < 2:1071051/150‑1 (MQ=255) tcGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGACCAAGGGCCAGTGTCAGCAGCATCAGCGAAACCAGCTGAgtg < 2:127411/150‑1 (MQ=255) gACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGACCAAGGGCCAGGGTCAGCAGCATCAGCGAAACCAGCTGAgnggt < 2:701611/150‑1 (MQ=255) gACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCACACCTGACCAAGGGCCAGGGTCAGCAGCATCAGCGAAACCAGCTGAgtggt < 1:983301/150‑1 (MQ=255) | TCGCCTGCTGGGTATGCTTCGGCGTGGATGTAGGAGATTTCTTTAAGCTTGAGAGCACCTTCCATCGCCACCGGGTACTGCGCGCCGCGGCCGAGGAACAGGGTATGGTGCTTGTCGGCGAACAGCTCGGCAATTTTCTCGACGGTGGCATCCATGGCCAGGGCTTCGCCCAGGCGGGTAGGCAGTCGACGCAGCTCTTCAACCAGCTCAGCTTCGACACCGGCTTCCAGCGTACCGCGCACCTGACCAAGGGCCAGGGTCAGCAGCATCAGCGAAACCAGCTGAGTGGT > NC_002947/6170847‑6171136 |
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: ‑ |