|
I am trying to fetch promoter sequences from a bacterial genome, retrieving only those sequences that don't fall into cds. I think the problem is that some sequences fall out of the FASTA sequences (2 chromosomes, 2 plasmids) that are not recognized as circular giving me an error from scanFa. Here there's my code
this is the error I get from scanFa:
this is the seqinfo of phy.fasta
I have tried also with isCircular but this method doesn't work with FAfile I hope that I have been clear, and thanks in advance for your time and your precious advises! |