Go Back   SEQanswers > Applications Forums > RNA Sequencing

Similar Threads
Thread Thread Starter Forum Replies Last Post
RNA-seq read alignment best practices for low-input Smart-seq2 samples Jon Goldberg RNA Sequencing 1 01-10-2018 12:43 PM
views Bowtie2 : aligned result read count not matching with sam Read count vinumanikandan Bioinformatics 0 05-13-2015 03:00 AM
HiSeq Low Read Count - Barcode in Read1 tmm343 Illumina/Solexa 9 12-04-2014 05:45 AM
DESeq2: dealing with count outliers and interpretation of the results Tatsiana_by Bioinformatics 13 11-10-2014 04:33 PM
dealing with homopolymer/low complexity errors The_Roads Ion Torrent 0 07-12-2012 07:09 AM

Thread Tools
Old 09-10-2018, 12:55 PM   #1
Junior Member
Location: CDC,Atlanta, GA

Join Date: Feb 2012
Posts: 1
Default Dealing with low read count in some samples

Greeting Members

I have an experiment with 4 time points and 5 biological replicates and each biological replicate has 5 technical replicates (100 samples thus)

The first two time points are really early is the organism's life cycle and hence we often expect less reads for some genes:

After summing the reads from the technical replicates, I get this profile (OBVIOUSLY AN ARTEFACT)
0h GeneA 3 reads
0h GeneA 2 reads
0h GeneA 0 reads
0h GeneA 1 reads
0h GeneA 3 reads

24h GeneA 0 reads
24h GeneA 1 reads
24h GeneA 2 reads
24h GeneA 0 reads
24h GeneA 1 reads

But DeSeq2 analysis of DE of these genes gives me a profile like this:
baseMean log2FoldChange lfcSE stat pvalue padj
31.88739159 -4.026822896 0.943450084 -4.268188603 1.97E-05 0.007173227

How do I filter for such artefacts ? Though I give an example here, my results are severely afflicted with these kind of low read count but highly significant statistics.
Any suggestion is highly appreciated...

I cannot set a minreadcount filter as these genes will have higher read count as the time line progresses (and will also skew the stats if I only remove it for one time point)
sganesh02 is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 09:47 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO