It's not clear to me what you mean by filtering; do you mean trimming off low-quality bases and dropping short/N-containing reads, or something else?
You might want to know what those 16 % non-mapped reads are - if they contain adapter sequences it is lost data that could possibly be recovered.
It is not really necessary to trim off low quality 3' ends with RSEM, since it allows for many mismatches outside Bowtie's seed region (unless they've changed that recently). However, I'd clip off adapter sequences (if any suspected when examining the data with e.g. FastQC) and remove rRNA reads to avoid total read count biases (rRNA tend to be abundant and sometimes at quite different amounts in different samples).
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Is Data filtering step mandatory?
Hi All,
I have been trying to use Tophat and also RSEM on illumina data. I always filter the data first when I use Tophat and I have always wondered whether it does its own filteration step or it would be better to filter the data yourself and then trim it using FASTX toolkit etc. Using this approach and then using cuffdiff the overlap with qrtpcr has not been good. So I started using RSEM and I thought I made a mistake by not filtering the reads but then I found that RSEM does neglects the reads by itself as is evident from its output.
# reads processed: 54476831
# reads with at least one reported alignment: 45346199 (83.24%)
# reads that failed to align: 9054059 (16.62%)
# reads with alignments suppressed due to -m: 76573 (0.14%)
Thus I was wondering whether my analysis with RSEM would be correct or wrong. As ideally such tools should automatically filter the reads if they are not bad.Tags: None
Latest Articles
Collapse
-
by seqadmin
The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...-
Channel: Articles
07-08-2024, 03:19 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 07-25-2024, 06:46 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
07-25-2024, 06:46 AM
|
||
Started by seqadmin, 07-24-2024, 11:09 AM
|
0 responses
26 views
0 likes
|
Last Post
by seqadmin
07-24-2024, 11:09 AM
|
||
Started by seqadmin, 07-19-2024, 07:20 AM
|
0 responses
160 views
0 likes
|
Last Post
by seqadmin
07-19-2024, 07:20 AM
|
||
Started by seqadmin, 07-16-2024, 05:49 AM
|
0 responses
127 views
0 likes
|
Last Post
by seqadmin
07-16-2024, 05:49 AM
|
Leave a comment: