Saturday, April 19, 2014

Bedtools: Extract FASTA sequences based on a BED file

here is how:

$ bedtools getfasta -fi <input FASTA> -bed <BED/GFF/VCF> -fo <output FASTA>
 
For example:
 
$ bedtools getfasta -fi hg19.fa -bed wzy.bed -fo out.fasta

No comments:

Post a Comment