Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,276,134 | C→T | 6.9% | R222C (CGC→TGC) | yjcC → | putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,276,134 | 0 | C | T | 6.9% | 156.1 / 2.6 | 58 | R222C (CGC→TGC) | yjcC | putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase |
Reads supporting (aligned to +/- strand): ref base C (39/15); new base T (2/2); total (41/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.73e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01 |
AATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCT > NC_000913/4275988‑4276280 | aaTGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGCCAAAATGc > 1:799304/1‑149 (MQ=255) aaTGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAATGc > 1:358264/1‑149 (MQ=255) gtgtgtATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAAtt > 1:116193/1‑149 (MQ=255) aaggTGACGAAAACCTTTTTCACGTAAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCAATTGAAGGATTTAACCGTACAAACAAATGGCTATTAATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCAatgt < 2:739008/147‑1 (MQ=255) ggTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGTATTAAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCAttgttg < 1:627049/149‑1 (MQ=255) gTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGc < 2:570263/149‑1 (MQ=255) aaaaCCTTTTTCTCGTTAATCAAATAGGCCTCGGAAGCAACGTTTTCTCCGCTGATTCATTTGATGGAGTTAACCGTACAATGAAATGTCTTTTTATATGCGACAGTTTATTCGACAAAACGCCCAATTTCAGCCATTGTTGCGTCTTc > 2:510456/1‑149 (MQ=255) tctcGTTAAGCAAAGAGGCCTCGGCAGCAACTTCTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGt < 2:604428/149‑1 (MQ=255) aaaGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCa > 2:483256/1‑149 (MQ=255) tagGCCTCGGCAGCAACGTTTTCTCCGCTTATTCATTTGAAGGATTTAACCGTACAAAGAAATCGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATGGTTGCGACTTCATATCAACGTCTTATAACCCAttt < 2:714495/148‑1 (MQ=255) gCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTta > 2:529730/1‑149 (MQ=255) ccTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTtat > 1:499683/1‑149 (MQ=255) ccTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTtaa > 2:444270/1‑148 (MQ=255) aaCGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTa > 2:421646/1‑149 (MQ=255) tctcCGCTGATTCATTTTAAGTATTTAAGCGTACAAAGAAATGGCTTTTTATATGAGAAAGTTTATTCGACAAAACGCCCTATTGCAGCAATTGTTGCGACTTCATATCAACGTCTTATAACCCATTATTATAATCATCTTATTTTTGc < 1:853312/149‑1 (MQ=255) ctcCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCg < 2:992442/149‑1 (MQ=255) ctcCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGATGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCg > 2:425288/1‑149 (MQ=255) tcCGCTGATTCATTTGAAGGATTTAACCGTACTAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGt > 1:949538/1‑149 (MQ=255) aTTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTTTATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAAACCATTTTTATAATCTACTTATTTTTGCGTTGCccga > 2:1002070/1‑148 (MQ=255) aTTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAATCGCCCTATTGCAGCCATTGTTGCGAGTTCTTATCAACGGCATATACCCCATATTTAGAACAATCATATTTTTGCGTTCCccgc > 1:448835/1‑149 (MQ=255) ttCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCccgcc > 1:650243/1‑149 (MQ=255) tCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCccgccg > 1:448264/1‑149 (MQ=255) gAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTAttt > 2:382227/1‑149 (MQ=255) aaGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAGCCCAGTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTAtttt > 2:505449/1‑149 (MQ=255) gtattAACCGTACAAAGAAATGGCTATTTATAAGCGACATTTTATTCGACAAAACGCCCAATTTCACCCATTGTTGCGTCTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTggg < 1:444270/146‑1 (MQ=255) aTTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTgggg > 2:753091/1‑149 (MQ=255) tttAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATCACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGa > 1:775533/1‑149 (MQ=255) ttAAACGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAg > 2:539490/1‑149 (MQ=255) aaCCGTACAAAGAAATTTCTATTTATATGCGACAGTTTATTCGACAAAATGCCGAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATATTTGCGTTGCCCGCCGGTATTTTGGGGAggc < 2:799304/149‑3 (MQ=255) aaCCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAATGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTc < 2:358264/149‑1 (MQ=255) aaCCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATGTTATTTTTGCGTTGCCCGCGGGTATTTTGGGGAGtg > 2:164604/1‑148 (MQ=255) aaCCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACACATTTTTATAATCATCTTATTTTTGCTTTGCCCGCCGGTATTTTGGGGAGTc > 1:34306/1‑149 (MQ=255) aaCCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCACAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTCATAACCCATGTTTATAATCATCTTTTTTTTGCGTTGCACGCCGGGAGTTTGGGGAggc > 2:656433/1‑147 (MQ=255) aaCCCTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTc < 1:51024/149‑1 (MQ=255) cGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTg > 1:174851/1‑149 (MQ=255) gTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGt > 1:725058/1‑149 (MQ=255) aaaGAAATGGCTATTTATATGCGCCAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTAGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTAGTTTAGCGTTGCACGCCGGTATTTTGGGGAGTCTTGTTCTg > 2:37114/1‑149 (MQ=255) aaaGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTg > 1:278298/1‑149 (MQ=255) gAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTCTATCCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCtt < 2:266677/149‑2 (MQ=255) aaaTGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTAt > 1:289465/1‑149 (MQ=255) gCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCGATTTTTATAATGATCTTATTTTTGCGTTGCCCGCCGTTATTTTGGGGAGTCTTGTTATGCTATTaata > 2:586240/1‑146 (MQ=255) tATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTTCGTTGCCAGCCGTTATTTTGGGGAGTCTTGTTCTTCTATTACTCTg > 1:578068/1‑149 (MQ=255) atatGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGGAATGTGGGGAGTCTTGTTCTGCTATTACTCTGGCTAc < 1:347525/149‑1 (MQ=255) tatGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACg > 2:625901/1‑149 (MQ=255) tatGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGAGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTAATTTTGCGTTGCCAGCCGGGATTTTGGGGAGTCTTGGTCTGCTATTACTCTGGCTACg > 2:734540/1‑149 (MQ=255) gCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTAt > 1:293093/1‑149 (MQ=255) gACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTc > 1:1026079/1‑149 (MQ=255) aCAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCg > 2:288247/1‑149 (MQ=255) tttATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACaa > 1:466219/1‑149 (MQ=255) tATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACaaaa < 1:529730/149‑1 (MQ=255) aatCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGGCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAAc < 2:116193/147‑1 (MQ=255) ttCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTAGTTCTCTGGCTACTGATTCGACACAACt > 2:958577/1‑149 (MQ=255) gACAAACCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGGTATTACTCTGGCTACGTATTCGACAAAACTAtt < 2:949538/149‑1 (MQ=255) gACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACTAAACTAtt < 2:650243/149‑1 (MQ=255) aCAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTAttt > 2:854228/1‑149 (MQ=255) aaaaCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTAt > 2:858446/1‑149 (MQ=255) aaaaCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTAGTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACGATTTAt > 2:743157/1‑149 (MQ=255) aaCGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTAtct > 1:316461/1‑149 (MQ=255) | AATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCT > NC_000913/4275988‑4276280 |
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: ‑ |