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
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
Comment