maybe samtools mpileup or GATK modify the associated R script
maybe samtools mpileup or GATK
modify the associated R script