I want to extract some base pairs from a FASTA file. There is a given sequence, say marker sequence. Every time this marker sequence occurs in the FASTA file, I want to extract n base pairs to the left of it (before it). Is it possible to do so using biopython? If so, please tell.
Thank you.
PS - I know it seems quite simple. But I am very new to python. And I have to do this using biopython only. So finding it very difficult to understand from the cookbook.
•
written
by
priyankakukreja267 •
Why not just use regex capture groups? I believe this would be much more efficient.
•
written
by
Matt Shirley
Agreed.
•
written
by
Matt Shirley
Biopython is the only package that I am acquainted with. That`s why the special emphasis on using that one.
But now, as it turns out I'll have to learn regex too.. Thank you for your help.. :)
•
written
by
priyankakukreja267 •
As mentioned by RamRS, you could use Bio.SeqIO for parsing the FASTA file. Or if you prefer plain strings, from Bio.SeqIO.FastaIO import SimpleFastaParser might be useful?
If you are looking for an exact substring match, both the Python string object and the Biopython Seq (sequence) object both offer a .find(...) method.
If you are looking for a more complicated pattern, then as RamRS suggested the Python regular expression library might be a good choice.
Log in to answer this question.
