This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Metagenomics normalization for long reads

Hello everyone,

I am currently working on a metagenomics analysis using ONT and I have used epi2me wimp for classification of the reads. What I want is to plot the relative abundance of each species present in the sample.

Epi2me outputs a read count assigned to each taxa so I can produce a plot for the proportion of reads per taxa but since each reads has a different length should I normalize by read length ?

metagenomics normalization nanopore
• 1,543 views • link
updated by colindaven colindaven • written by fonteneaudam • fonteneaudam

1 answer

The more relevant normalization would be

  1. normalize to one million sequenced reads per sample (if you sequenced more reads, you'll get more hits)
  2. normalize to bacterial genome size (if the bacteria has a big genome, you'll get more reads from it compared to one with a small genome)

We cover these in our pipeline Wochenende https://github.com/MHH-RCUG/nf_wochenende (older version https://github.com/MHH-RCUG/Wochenende )

• 0 views • link
written by colindaven colindaven

I get it one read equals one fragment no matter its size but the genome size is of influences on the number of fragment so I have to normalize for that. Thanks

• 0 views • link
written by fonteneaudam • fonteneaudam

Log in to answer this question.

More posts like this