This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to get the diversity one SNP loci by one using the sequenceDiversity function of vcflib

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?

sequence diversity nucleotide diversity pi vcflib

0 answers

No answers yet.

Log in to answer this question.