Im using Bio::Index::Fasta in order to retrieve sequences from my indexed fasta file but I keep seeing MSG: Did not provide a valid Bio::PrimarySeqI object when I call fetch followed by write_seq() on a Bio::SeqIO handle. Why?
It’s likely that fetch didn’t retrieve a Bio::Seq object. There are few possible explanations but the most common cause is that the id you’re passing to fetch is not the key to that sequence in the index. For example, if the FASTA header is >gi|12366 and your id is 12366 then fetch won’t find the sequence, it expects to see gi|12366.
Related Questions
- Im using Bio::Index::Fasta in order to retrieve sequences from my indexed fasta file but I keep seeing MSG: Did not provide a valid Bio::PrimarySeqI object when I call fetch followed by write_seq() on a Bio::SeqIO handle. Why?
- If you parse a FASTA sequence format file with Bio::SeqIO the sequences won have the accession number. What to do?
- How to get a fasta file of multiple DNA/RNA sequences?
