I am very new to this forum and also new to do variant calling. I am using samtools mpileup|bcftools|vcfutils.pl commands for this.
I need to find the mean position of variant (or positions of the variant) within supporting reads. Does anyone know how to obtain this?
Thanks a lot in advance for any of your comments and suggestions.