Hi, I am struggling with the equation provided in this publication (https://www.nature.com/articles/ismej201559#Sec2). The authors have presented a formula for calculating the abundance of ARGs per copy of 16S rRNA genes :
I annotated my metagenomic reads (fastq) against the CARD database using DIAMOND Blastx, and determined the number of 16S rRNA genes identified using Metaxa2. My question regarding the formula is whether I should use the average length of the total reads or the average length of the reads carrying only one specific ARG type since the read length is not unique like in the case of Illumina sequencing? Another question is whether I should use the same number of 16S rRNA genes or the total number of 16S rRNA genes identified in my data?
