This is a test version of Biostars. For the public version, visit https://www.biostars.org.
bam to vcf to .bed

Is the following the correct process? I do not believe I can "convert" .bam to .bed directly. I am converting .bam to vcf using bcftools, and then using plink to prepare .bed and .fam files.

bed to vcf using bcftools

vcf to bed using plink

bed + fam using plink

Thank you.

bed bam vcf
• 2,961 views • link
updated by finswimmer finswimmer • written by genomics Newbie • genomics Newbie

I think you're using the word "convert" a little incorrectly here. One cannot "convert" a BAM to a VCF because they do not represent the same information. You can obtain a VCF by processing a BAM with a (reference) FASTA, but cannot convert the BAM to a VCF. You can also extract a BED out of a VCF, but you don't convert it there either. I'd only use "convert" if the information content in the formats is comparable tending to equivalence.

• 1 views • link
written by Ram Ram

What do you want your bed file to contain? Coordinates of variants? Coordinates of mapped reads?

• 1 views • link
written by swbarnes2 swbarnes2

genomics Newbie : Don't forget to follow up on your threads.

If an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. You can accept more than one if they all work.
Upvote|Bookmark|Accept

Please do the same for your previous posts as well.

• 1 views • link
written by GenoMax GenoMax

1 answer

Important: there are two types of bed files containing entirety different information. The type you want (plink) is not the same as mentioned in the other comments. What you are doing is fine.

• 1 views • link
written by WouterDeCoster WouterDeCoster

Log in to answer this question.

More posts like this