Hello,
I have run featureCounts on aligned bam files (DNA from bacterial genomes).
I have a table with counts per sample as expected. However i would like to normalize by gene length since genes of different length can have more reads mapped then genes that are shorter.
I was wondering what would be the best way to do that ?? Any tool or script ? Once i have normalized by gene length i would like to feed it to DESEq2
Thanks
1 answer
DESeq2 takes care of all necessary steps internally. There is no need to manipulate the data prior to creating your DESeq2 object. Raw counts are expected. Please read the manual.
Not really. There is no need to account for gene length because gene length is identical for all conditions. See also as a starting point for further online searches: https://support.bioconductor.org/p/67131/ and DESeq2 normalisation: is the size of the gene taken into account?
How do you mean it's not possible? DESeq will do a far more sophisticated normalization than you can think off, and you can get the normalized counts. Then you can rank whatever you want. And now, before we help you any further, you are going to read all the links ATpoint gave to you.
Log in to answer this question.
-
Appropiate normalization method for OTU table built from MAGs abundances estimations?
written by v.berriosfarias •
-
Normalization to Read Depth After Alignment to Database
written by Nat •
-
How to normalize gene expression with CPM
written by concetta •
-
Accounting for differences in FASTQ-file size, when comparing metagenomic gene abundance between sa…
written by Hansen_869 •
-
Normalizing gene counts in a quantitative metagenomic DNAseq analysis
written by allerdrengen55
-
DESeq multiple comparisons post normalization
written by zapsep •
-
RPM data analysis with DESeq2
written by ruhollah •
-
Compare samples from bacterial genome
written by David
-
Alignment-free tools (like Kallisto or Salmon): useful for bacterial transcriptomes?
written by arctan •
-
Insert size on PE 2*250
written by David
