View Single Post
Old 03-11-2014, 10:07 AM   #1
NGS_New_User
Member
 
Location: USA

Join Date: Sep 2012
Posts: 41
Exclamation Problem running bedtools (fastFromBed)

Hi Everyone!
I have a gff file of the drosophila r5.48 genome, and the corresponding all chromosome fasta files.
I want to extract the DNA sequences into a fasta file based on the coordinates from the gff file. I have tried using the command 'fastaFromBed' using bedtools, and it starts running well, but then suddenly stops saying that it is "unable to find FASTA index entry for '2L'"
Can anyone please help me out on what I might be doing wrong?

The command I am using is:
fastaFromBed -fi dmel_all_chr.fa -bed dmel_all.gff -fo out.fa

I will appreciate any ideas and/or suggestions
NGS_New_User is offline   Reply With Quote