Open Euphrasiologist opened 3 months ago
You don't need to reverse complement the FASTA file yourself. But how you do it will depend on the tool that you use to fetch the subsequence(s) from the FASTA file. Different tools use different conventions for coordinates, especially for the reverse complement strand.
If you're using our esl-sfetch
tool, if you give the coordinates as 151479..151422
(i.e. with start > end), it fetches the reverse complemented subsequence from the FASTA file.
Okay thank you so much, very helpful! I didn't know about the easel tools (I see you wrote more about this here: http://cryptogenomicon.org/extracting-hmmer-results-to-sequence-files-easel-miniapplications.html).
Hi, thanks for this great tool!
I have a question about extracting sequences from the input fasta. I've got an example output below:
For the first line, the extraction is easy because the strand is "+". I just search the input fasta for the position 32213024 to 32213361. Great. For the second line the strand is "-". This might be a silly question, but do I need to reverse complement the input fasta before I do the search (151422 to 151479)?
Cheers, M