This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Samtools Mpileup Defined Regions

Hi! I need to get genotypes in a list of positions (for example BED file) having BAM file. Trying to use tool samtools mpileup. I know I can do it by creating consensus fastq using

samtools mpileup -uf hg19.fa file.bam | \
bcftools view -cg - | \
perl vcfutils.pl vcf2fq > consensus.fastq

and then getting genotypes from fastq file, but it's ineffective. Do you know how to get differences from reference, genotype calls at interesting positions simultaneously?

samtools
• 5,165 views • link
updated by swbarnes2 swbarnes2 • written by Vova Naumov Vova Naumov

Please take a closer look at mpileup in samtools. The documentation clearly states how to call variants and use a bed file.

• 139 views • link
written by Zev.Kronenberg Zev.Kronenberg

Thank's a lot but an example would be very helpful

• 0 views • link
written by Vova Naumov Vova Naumov

1 answer

I'd use BEDTools intersectBed to get the intersection between your all-points vcf and your bed file coordinates. That's way simpler than making a consensus fastq.

• 256 views • link
written by swbarnes2 swbarnes2

Log in to answer this question.

More posts like this