Hi everyone,
I have a very basic question, I have metagenomic samples and I want to extract reads from a particular gene (amoA gene) to see the diversity of this gene between samples but I am not sure what would be the best approach. Should I just blast the metagenomic reads against amoA sequences? or maybe it is better to assemble first and then blast and extract the sequences?
Thank you for your time!