Entering edit mode
8.1 years ago
robertorun
•
0
I want to to get the diversity one SNP loci by one using the sequenceDiversity function of vcflib, but it's failure.
As we all known, there is w argument for windows size (INFO: optional: w,window -- argument: the number of SNPs per window; default is 20) . But if you set "-w 1", it will give the sequence diversity from this SNP loci to the front next SNP loci (aslo the region result). The following result is setting w to 1.
A01 9712 15700 3.40818e-05 0.58816
A01 15700 17403 2.5044e-05 0.913932
A01 17403 24170 3.07121e-05 0.580598
How can i get it?
