SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem running vcfutils.pl Balat Bioinformatics 11 11-11-2014 10:11 AM
Which cloumn to take from SV Detect to make bed format to make it running on BEdtools vishal.rossi Bioinformatics 1 05-16-2013 03:59 AM
Problem to use Bedtools after filtering uniquely mapped reads with samtools eilosei Bioinformatics 2 12-21-2011 03:51 PM
problem with running MAQ seq_lover Bioinformatics 0 10-30-2011 09:31 PM
tophat running problem yasu Bioinformatics 4 02-08-2010 10:54 PM

Reply
 
Thread Tools
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
Old 03-12-2014, 05:23 AM   #2
TiborNagy
Senior Member
 
Location: Budapest

Join Date: Mar 2010
Posts: 329
Default

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.
TiborNagy is offline   Reply With Quote
Old 03-18-2014, 11:06 AM   #3
NGS_New_User
Member
 
Location: USA

Join Date: Sep 2012
Posts: 41
Default

Thanks, I will try to look for it
NGS_New_User is offline   Reply With Quote
Reply

Tags
bedtools, fastafrombed

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 03:13 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO