Hi,
The VCF output currently has ADr and ADf, a combined AD field would be very useful to potentially have CRISP VCF be compatible with the poolfstat package in R.
Edit: I've written a python script to add an AD field and related counts to the VCF and it is now being accepted by poolfstat.
If the results I obtain are the same as when I use a mpileup then this could make CRISP a really useful tool for these kinds of pool-seq studies as the mpileup file is 2.5T in size whereas a crisp VCF is 14G in size.
Hi, The VCF output currently has ADr and ADf, a combined AD field would be very useful to potentially have CRISP VCF be compatible with the poolfstat package in R.
Edit: I've written a python script to add an AD field and related counts to the VCF and it is now being accepted by poolfstat. If the results I obtain are the same as when I use a mpileup then this could make CRISP a really useful tool for these kinds of pool-seq studies as the mpileup file is 2.5T in size whereas a crisp VCF is 14G in size.