RNA seq library size
I apologize for the silly question, but how do I find the library size of rna seq data if I have the count table? The count table was created using htseq.
Is it just the sum of the counts of genes?
• 4,597 views
•
link
written
by
felipead66
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Normalization of ChIP-Seq count matrix ??
written by Ankit
-
strand-specific RNA-seq analysis
written by 17318598206 •
-
STAR Aligner --quantMode GeneCounts
written by Seigfried
-
RNA-Seq, HTSeq-Counts, counts are ~1M less compared to the reference research paper
written by gdm01 •
-
Cannot find a gene in the counts table
written by felipead66
-
How to normalize gene expression with CPM
written by concetta •
-
Calculating FPKM and TPM by hand from htseq-count output?
written by Joel Wallenius
-
Normalization of read counts from HTSeq by library size
written by Susmita Mandal
-
How to separate htseq-count table to coding RNAs and non-coding RNAs
written by Naresh D J •
-
Normalization Procedure For The Rna-Seq Data Based On Counts
written by mssestak •
the problem is not the question, but the term vague "library size". It depends on what you want to do with the estimate. For diff expression normalisation one normally uses the sum of counts in genes. If you mean the count of all the reads that came out of the sequencer for your library? Then depending on the preprocessing (e.g. adapter trimming, removal of small fragments, removal of contaminants) you might not have all the needed info. Then the best approximation you have is sum of count in genes + the sum of the last 4? rows which is the unassigned or ambigous reads.
The actual "usable" (and therefore meaningful) library size that goes into normalization and DEG analysis is the column sum for each sample considering the genes you quantified against excluding rRNAs imho.