I have a metagenomic dataset across multiple samples that has been grouped and binned (we did this due to poor sequence quallity).
It is several samples taken from termite guts of differing castes, the samples have been grouped to make binning easier due to poor sequence suitability. They have been binned via ASVs, assembling reads from all samples to generate the putative bacterial species. So I have fasta files of putative species that have been annotated via prodigal. I want to now calculate the abundance of each bin in a given sample. I think the best way to do this would be at the gene level by looking at the number of reads that map to a gene in any given sample. I am just unsure how to do this.
I have annotated the bins and now want to calculate the abundance of each bin in a sample but looking at the number of reads from a sample that map to each gene in a bin.
Is this a viable method? If so how would I do it and if not what else could I do?