Reducing filesize for homozygosity mapper
0
0
Entering edit mode
10.7 years ago

Dear all,

I wanted the homozygosity mapper tool for finding homozygous regions in two affected individuals in a family of four. The data can only be entered in a combined vcf for which the original bam files need to be used. the command mentioned for the purpose is

samtools mpileup -D -gf /path/to/genome.fa *.bam | bcftools view -c -g - > filename.vcf 

I have four bam files (~3-4GB) which generate a combined vcf of 8 GB which cannot be uploaded. Any suggestions?

Thanks :)

homozygosity-mapper samtools • 1.5k views
ADD COMMENT

Login before adding your answer.

Traffic: 3569 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6