Abundance of ARGs per copy of 16S rRNA genes · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Abundance of ARGs per copy of 16S rRNA genes

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 : enter image description here

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?

abundnace args

• 581 views • link

written 3.1 years ago by lea •

0 answers

No answers yet.

Log in to answer this question.

More posts like this

Very high number of unique sequences in PacBio full-length 16S CCS data using DADA2
written 5 months ago by muhammad.ali •

I am analysing PacBio Sequel II full-length 16S rRNA CCS reads (~1450 bp) using the DADA2 long-read workflow and observing an unusually high number of …
STAR aligner can't map too short reads
written 3.6 years ago by Assa Yeroslaviz

For our Ribo-seq data set I tried the star aligner but was able to map only a very small fraction of the reads (<1% in …
question about the equation of counts in STAR ReadsPerGene.out.tab
written 6.7 years ago by boymin2020 •

Hi, Recently I have been evaluating gene expression levels using STAR software. Regarding the counts in the ReadsPerGene.out.tab table, I am confused about the equation …
STAR output - high multi-mapping, low alignment (human)
written 8.2 years ago by anikb •

Hello, I am new to RNA seq analysis and was hoping if the community would help me understand a few things about the analysis. I …
MuTect2 still discard many reads, how to fix?
written 8.5 years ago by Sharon

Hello everyone MuTect2 from GATK3 discarded `87.33%` of the reads. This is during preparing panel of normals. I will use sp1.vcf in the PON creation …
Extracting an information about conservation of 16S rRNA genes of bacteria from multiple sequence a…
written 11.1 years ago by Peter Leontev

Hi. I am trying to get any useful information from 16S rRNA bacterial genes multiple sequence alignment (MSA), i.e. for sake of simplicity, average length …
Problems with STAR mapping output
written 11.1 years ago by manekineko

Hi, I have problems seeing the output from STAR mapper. I have done the mapping, then made an sorted bam for IGV but I see …
Normalization of certain gene abundance in a metagenome
written 11.7 years ago by bioinfo

Hi I have analysed two metagenomes (whole-genome shotgun; 5 and 8 million reads respectively, 101 bp reads) to look for antibiotic resistance genes and will …
Confusion in HTSeq paired-end output count
written 11.9 years ago by Whoknows

Hi friends I used HTSeq for creating count table (outputs) in paired-end file, actually I generated RPKM values from these count tables with below formula …
Rpkm Calculation For Genes
written 13.8 years ago by siddharth.sethi5

<p>Hi everyone</p> <p>I understand the RPKM fomula is as follows:</p> <p>C = Number of reads mapped to a gene</p> <p>N = Total mapped reads in …