SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   Problem running bedtools (fastFromBed) (http://seqanswers.com/forums/showthread.php?t=41585)

NGS_New_User 03-11-2014 10:07 AM

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

TiborNagy 03-12-2014 05:23 AM

I think there is a record in your GFF that contains a contig name that can not be found in your FASTA file. Remove this record.

NGS_New_User 03-18-2014 11:06 AM

Thanks, I will try to look for it


All times are GMT -8. The time now is 03:48 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.