Hi, I have a set of predicted orfs from metagenomics bins after assembly. I have extracted a set of orfs that are single copy genes , reflecting the diverstity of my sample.
How can i get the copy number for each of this genes (actually 139 genes) in my assembly.
My first idea is to map my single copy genes to my assembly and get coverage for each. Is that a suitable approach or is there a better way to get copy number for each my single copy genes ?