How to deal with single reads from PE reads after trimming?
1
1
Entering edit mode
9.9 years ago
Cacau ▴ 520

I have some single reads from paired-ended reads data after trimming the reads. Should I map these to the genome? For example, to use the same parameters as paired reads?

RNA-Seq alignment ChIP-Seq • 1.9k views
ADD COMMENT
0
Entering edit mode
9.9 years ago

You certainly can map them like the paired-end reads. The question becomes if it's worth your time to bother. If you don't have a lot of orphans like this, then you can just ignore them. If, however, you have a LOT of orphans (e.g., there was a problem with read #2) then you pretty much have to align the orphaned reads. Use the same settings as for the paired-end datasets (within reason).

Edit: I should add that the SE reads might map slightly differently than the PE reads. In many cases this isn't a problem. However if whatever you're doing is very sensitive to any bias then either compare the results with/without including the SE reads or at least ensure that the proportion of SE reads is approximately constant across samples (i.e., try to ensure that all of the samples have a similar bias).

ADD COMMENT
0
Entering edit mode

Thanks a lot!

ADD REPLY

Login before adding your answer.

Traffic: 1003 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6