Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 901,418 | G→A | 10.3% | intergenic (‑137/‑194) | cysZ ← / → trxB | sulfate transport protein CysZ/thioredoxin reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 901,418 | 0 | G | A | 10.3% | 154.7 / 6.2 | 58 | intergenic (‑137/‑194) | cysZ/trxB | sulfate transport protein CysZ/thioredoxin reductase |
Reads supporting (aligned to +/- strand): ref base G (26/26); new base A (4/2); total (30/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.71e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01 |
GCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGC > NC_002947/901274‑901565 | gCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGc > 6:232284/1‑149 (MQ=255) tGCAGCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTa < 4:21746/149‑1 (MQ=255) aTCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGa > 4:278650/1‑149 (MQ=255) cTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAgg > 2:63860/1‑149 (MQ=255) agagACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTg > 8:50086/1‑149 (MQ=255) gagaCTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTgg < 3:278650/149‑1 (MQ=255) gagaCTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTg > 5:47264/1‑148 (MQ=255) agaCTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGa > 1:343978/1‑149 (MQ=255) ccGCAGAGAAGGAACACGCGCTGACCTTACCGAGTTTGCCGCGCCGACCAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTc < 7:66549/149‑1 (MQ=255) ctcAGAGAATGAAATCTCGCTGTCCTTAACGAGTTTGCCGCTCCGTAGAAATGTCACTAGCCGTTGGCCTACAGTAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGt < 1:279795/144‑1 (MQ=255) cGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGATTGTTTATAGATGAAATCGCAGCATAGGCGATGCATATGAGGTGGATTGTc > 8:68847/1‑148 (MQ=255) cAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCt > 1:203907/1‑147 (MQ=255) gagaAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCATCGCCTATgctgc > 3:312269/1‑132 (MQ=255) gagaAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCATCGCCTATgctgc < 4:312269/136‑5 (MQ=255) aaGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACAGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGata > 1:276574/1‑149 (MQ=255) cgcTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAAt > 7:194223/1‑149 (MQ=255) aCCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTTTAAGGATATTTCCTTAATCttt < 7:50086/148‑1 (MQ=255) ccTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTg < 4:32671/148‑1 (MQ=255) ccTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGc < 1:63860/149‑1 (MQ=255) ccGAGTTTGCCGTGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCt > 2:328370/1‑149 (MQ=255) acgaAA‑GCCAC‑‑GCGGTGGGC‑‑‑GGGAAAGGGGATGAAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCaa < 2:343978/149‑1 (MQ=255) gaAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACt > 7:29031/1‑149 (MQ=255) gaAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACt > 6:95648/1‑149 (MQ=255) aa‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGt > 1:167712/1‑137 (MQ=255) aa‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGt < 2:167712/137‑1 (MQ=255) a‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCttt > 1:241470/1‑120 (MQ=255) a‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCttt < 2:241470/120‑1 (MQ=255) gCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCg > 7:285201/1‑149 (MQ=255) aC‑‑GGCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCt < 5:27056/149‑1 (MQ=255) tGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCgg < 7:68847/149‑1 (MQ=255) tGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCgg < 2:217308/149‑1 (MQ=255) gggAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTg > 8:196809/1‑149 (MQ=255) ggATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTAAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGtttt < 2:276574/149‑1 (MQ=255) gATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTAAATCTTTGCGACCTCGCTACTGTGCGCTCACCTTTTCGCTTTGCGGGCCTGCGCGTTTTa < 6:47264/149‑1 (MQ=255) aaCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGc > 8:138002/1‑93 (MQ=255) aaCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGc < 7:138002/93‑1 (MQ=255) aCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTa > 6:135787/1‑143 (MQ=255) aCAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTa < 5:135787/143‑1 (MQ=255) aaaaGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTccc > 6:267729/1‑149 (MQ=255) aGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCaa < 2:203907/149‑1 (MQ=255) cTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTgg < 5:267729/149‑1 (MQ=255) gATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCttttt > 2:58635/1‑148 (MQ=255) gAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTc < 5:78106/120‑1 (MQ=255) gAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTc > 6:78106/1‑120 (MQ=255) aaTCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTa < 4:319152/144‑1 (MQ=255) aaTCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTa > 3:319152/1‑144 (MQ=255) ggagctTACGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTACTTAATCTTTGCGATCTCTCTACAGTGCGCTGAACTTTTCTCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCTAGTGCTGCGTGGGTCCTTTTTATTTTCTtg < 7:196809/146‑1 (MQ=255) gcagcaTAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtgg < 5:232284/149‑1 (MQ=255) ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAAttt > 5:240557/3‑142 (MQ=255) ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAAttt < 6:240557/140‑1 (MQ=255) ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtggct > 7:89424/3‑149 (MQ=255) ctcaTAGGCATCGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtggct > 4:56377/3‑149 (MQ=255) gcaTAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTgg > 3:331082/1‑127 (MQ=255) gcaTAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTgg < 4:331082/127‑1 (MQ=255) caTAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTtggctgg < 1:79070/149‑1 (MQ=255) gCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAg > 6:246212/1‑149 (MQ=255) gCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAg > 8:101177/1‑149 (MQ=255) cGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGc < 8:194223/149‑1 (MQ=255) | GCTTGCATCTAGAGACTCCGCAGAGAAGGAAAACGCGCTGACCTTACCGAGTTTGCCGCGCCGACGAAA‑GGCAC‑‑GCCGTGGGC‑‑‑GGGAAACGGGATGTAACAAAAGCGCCCGGGGCTTGTTTATCGATGAAATCGCAGCATAGGCGATGCCTATGAGGTGGATTGTCTAAGGATATTTCCTTAATCTTTGCGACCTCGCTACAGTGCGCTCAACTTTTCGCTTTGCGGGCCTGCGCGTTTTAGCCTTCCCCAAGTGCTGCGTGGGTCCTTTTTAATTTCTTGGCTGGCGCAGC > NC_002947/901274‑901565 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |