Apologies if this has been asked a million times, a link to a previous post would be much appreciated if so.

I am looking at orthologues of human genes from other species. One part of the gene is well-conserved and easy to search for but the rest is more variable. I would like to use the conserved section to find relevant RNA reads and then extend the read by overlapping with other reads from an RNASeq library. I was thinking about how best to do this with Biopython but then came across Galaxy. Does anybody know if there's a standard workflow for this process or any Tools which might be handy?

Thanks a lot,

Joe