Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 851,009 | C→G | 6.7% | intergenic (‑589/‑216) | pcs ← / → hemA | phosphatidylcholine synthase/glutamyl‑tRNA reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 851,009 | 0 | C | G | 6.7% | 182.2 / 5.0 | 60 | intergenic (‑589/‑216) | pcs/hemA | phosphatidylcholine synthase/glutamyl‑tRNA reductase |
Reads supporting (aligned to +/- strand): ref base C (28/28); new base G (2/2); total (30/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01 |
TACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAA > NC_002947/850865‑851157 | tACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTgc < 1:780856/150‑1 (MQ=255) tACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTgc < 1:417767/150‑1 (MQ=255) tGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCAATTCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCa > 1:73470/1‑150 (MQ=255) gTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCaa > 2:689139/1‑150 (MQ=255) gCTTCTAAGCGCGCGGTTGATCCGGAGATGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGc < 2:316494/150‑1 (MQ=255) ttCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATTCCCGTCACCTCCCTGCGCAAGGCCt > 1:854830/1‑150 (MQ=255) tCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGCCCTg < 1:481075/150‑1 (MQ=17) aGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGtt > 1:852750/1‑150 (MQ=255) gcgcgcGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTc < 1:689139/150‑1 (MQ=255) gcgcgGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCgg < 1:272195/150‑1 (MQ=255) tGATCAGGCGCTGCAGATTGCGACTTCAGCAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCTTGCGGCCTCCTg > 1:549606/1‑150 (MQ=255) aTCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGaa > 1:732526/1‑150 (MQ=255) tGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCAAAGCGCTTGCCGGGACACCAGGGCCGTTACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaag > 1:142007/1‑150 (MQ=21) aTTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTTGGCCTCCTGAAGTCGCGaagatcaa < 1:695749/150‑1 (MQ=255) gACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAagatcaagatca > 1:579098/1‑150 (MQ=255) cctcgaggcggCCGAGCAGAGGGCGGGCTCTCTGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAg < 2:560400/145‑1 (MQ=18) ttCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGa > 2:64906/1‑150 (MQ=255) ttCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCACGa < 1:1207617/150‑1 (MQ=255) tCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGAt < 2:672567/150‑1 (MQ=255) cgagattAGAGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCa < 1:473389/144‑1 (MQ=255) aGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGa < 2:535963/150‑1 (MQ=255) gCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGaa < 2:925227/115‑1 (MQ=255) gCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGaa > 1:925227/1‑115 (MQ=255) ggCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAg > 1:322537/1‑150 (MQ=255) cTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCCAGCTCAGGATGAAAAGCCAGATTaaaattc < 2:620706/150‑4 (MQ=18) tGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCgg < 2:579098/150‑1 (MQ=255) tGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGACACAGCGGTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCgg < 2:73470/150‑1 (MQ=255) gaCCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGaaa > 2:638736/1‑150 (MQ=255) gaCCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGaaa > 1:438179/1‑150 (MQ=255) aCCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGaaaa > 1:885270/1‑150 (MQ=255) aCCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGaaaa < 2:885270/150‑1 (MQ=255) gcgcTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGCAAacac < 2:322537/150‑1 (MQ=255) gctgaggccccaggatggggccttcgGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaa < 1:114088‑M1/70‑1 (MQ=255) gctgaggccccaggatggggccttcgGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaa > 2:114088‑M1/27‑96 (MQ=255) ctgaggccccaggatggggccttcgGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAgg > 1:894477‑M1/26‑67 (MQ=255) ctgaggccccaggatggggccttcgGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAgg < 2:894477‑M1/42‑1 (MQ=255) tGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatc > 2:1010607/1‑97 (MQ=255) tGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatc < 1:1010607/97‑1 (MQ=255) tGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGAACAGAAAACTCCGcg > 2:688711/1‑149 (MQ=255) gCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTa < 1:377656/150‑1 (MQ=255) aGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAgg > 1:562900/1‑150 (MQ=255) ccccATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCa > 1:130029/1‑150 (MQ=255) ccATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAgg > 2:574212/1‑150 (MQ=255) gggCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGaaa > 2:929455/1‑150 (MQ=255) cacaGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAgg < 2:562900/150‑1 (MQ=255) aGCGTTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGa < 2:732526/150‑1 (MQ=255) aGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGAc > 1:123297/1‑114 (MQ=255) aGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGAc < 2:123297/114‑1 (MQ=255) GCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACTGCTAGGGCAGGCTACGAAAAGGCGAt < 2:502271‑M1/150‑1 (MQ=255) gggCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgc > 1:494205/1‑150 (MQ=255) gCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgcgc > 2:966020/1‑150 (MQ=255) cATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCAAGATCAACAGCCGGCAAGTACTGTAATGACCGCTATCGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCa < 1:929455/150‑1 (MQ=25) aTGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAg < 2:852750/150‑1 (MQ=255) ccGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTAc > 1:48009/1‑118 (MQ=255) ccGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTAc < 2:48009/118‑1 (MQ=255) gTCACCTCCCTGCGCAAGGCCTGCGTTTGGCCTCCTAAAGTCGCGAAGATCAAGATAAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTAGGAAAAGACGATAGCCGTTGCGCTATCAGCACAAt > 2:520927/1‑150 (MQ=38) cTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAacac > 1:344761/1‑95 (MQ=255) cTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAacac < 2:344761/95‑1 (MQ=255) tCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGaa > 1:336071/1‑150 (MQ=255) tCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGaa > 1:60460/1‑150 (MQ=255) | TACTGTCGCTTCTAAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAA > NC_002947/850865‑851157 |
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: ‑ |