Question: Convert fastq to plink .bed
0
gravatar for genomics Newbie
2.6 years ago by
genomics Newbie20 wrote:

I’m a genomic newbie, so I believe my question is a basic one.

I need to transform raw fastq files into plink .bed files. From my biostar searchers it appears I cannot directly transform fastq sequence files to a .bed format, as I would first need to first transform the fastq to .bam files. It appears after I generate .bam files I can then convert these .bam to plink .bed files?

Questions:

• Which tools (i.e. galaxy, etc) are best to convert fastq to .bam

• Which tools (i.e. bedTools) can be used to convert .bam to plink .bed

Any additional information regarding these format "conversions" would be helpful. I had difficulty installing Bowtie.

Thank you!

plink bam fastq bed • 2.8k views
ADD COMMENTlink modified 2.6 years ago by WouterDeCoster44k • written 2.6 years ago by genomics Newbie20

Any information about necessary alignments would also be helpful!

ADD REPLYlink written 2.6 years ago by genomics Newbie20

good luck with that.

ADD REPLYlink written 2.6 years ago by Pierre Lindenbaum129k

I'm sure you can come up with a join/comm/awk oneliner for this conversion.

ADD REPLYlink written 2.6 years ago by WouterDeCoster44k
1
gravatar for WouterDeCoster
2.6 years ago by
Belgium
WouterDeCoster44k wrote:

"Converting" fastq to bam is called alignment. You didn't tell us about the data format you have, but if it's just regular DNA sequencing you should be fine with bwa mem.

After alignment, you should do variant calling on the bam to get a vcf file of variants. This could be done using samtools/bcftools or GATK.

The vcf you obtained can then be converted to plink.bed using, well, plink.

ADD COMMENTlink written 2.6 years ago by WouterDeCoster44k

Thank you so very much for your quick response. Yes, you are correct that my original format is just regular DNA sequencing (fastq format).

By bwa mem do you mean Burrow-Wheeler Aligner for short-read alignment (see minimap2 for long-read alignment)? Can i do the alignment with any of the packages in galaxy? I can use the samtools in Galaxy, correct? I will also try the Genome analysis toolkit.

ADD REPLYlink written 2.6 years ago by genomics Newbie20

I don't use galaxy, but yeah, those tools would be fine.

ADD REPLYlink written 2.6 years ago by WouterDeCoster44k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1169 users visited in the last hour