This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Retrieve specific fasta sequences from a group of assemblies

Hi all,

Sorry if this question has been addressed before but I haven't been able to find a solution to this. I have a lot of assemblies (around 800) and I would like to retrieve the fasta sequence for a specific housekeeping gene which should (in theory) be present in all of them. Is there any tool that can take the fasta assembly as input and retrieve a specific gene with certain % variation to retrieve the gene even if it has mutations? Alternatively it could take a gbk or gff3 as input and use the gene annotation as retrieval criterion.

Thanks a lot!

cds gene sequence
• 1,733 views • link
written by SushiRoll • SushiRoll

Don't know if there is a ready made tool. You will need to align the gene to your assemblies and then it is a matter of parsing the results and retrieving the sequence you need using samtools faidx and similar options.

• 0 views • link
written by GenoMax GenoMax

Great, I'll give it a shot.

Thanks!

• 0 views • link
written by SushiRoll • SushiRoll

Personally, I never worked on similar tasks and thus unfortunately can't provide you with a polished solution, but what you are trying here is to find orthologous genes. Using this keyword, you should find tools suitable for this task, e.g. OrthoFinder showed up in a quick search.

• 0 views • link
written by Matthias Zepper Matthias Zepper

Thanks Matthias, that's a great starting point, I'll check what's out there.

• 0 views • link
written by SushiRoll • SushiRoll

0 answers

No answers yet.

Log in to answer this question.

More posts like this