How to draw SNP density diagram
This article mainly shows you "how to draw SNP density map", the content is easy to understand, clear, hope to help you solve your doubts, the following let the editor lead you to study and learn "how to draw SNP density map" this article.
Script running example: Rscript snp.density.map.R-I snp.vcf-n out-s 1000000-c "darkgreen,yellow,red"
-I: follow the input file in vcf file format
-n: sets the prefix of the output file name, and outputs the file to the current directory
-s: sets the window size
-c: set the color gradient. You can set one or more colors. It is recommended to set two or three colors, separated by ",".
The code is as follows: # Beijing formation Biotechnology Co., Ltd. # email: wangq@biomics.com.cnsource ("https://raw.githubusercontent.com/YinLiLin/CMplot/master/R/CMplot.r")library(getopt)#+--------------------# get options#+-spec