Converting PE to SE
0
0
Entering edit mode
7.2 years ago
Picasa ▴ 650

Hello,

I have RNAseq paired end in fr-firststrand (as described here http://chipster.csc.fi/manual/library-type-summary.html).

I have also RNA reads in single end (reverse orientation).

I need to perform a mapping but since it cannot handle a mix of PE and SE, I have to transform my PE into SE (and make it to the same orientation).

I have a doubt, is it correct that I just have to reverse complement my R2 (my PE) ?

PE SE • 1.5k views
ADD COMMENT
1
Entering edit mode

Why don't you map the PE and SE reads separately, then merge the resulting two bam files ?

ADD REPLY
1
Entering edit mode

BBMap can take in PE and SE data in at the same time. bbwrap.sh in1=read1.fq.gz,singletons.fq.gz in2=read2.fq.gz,null out=mapped.sam append

ADD REPLY
0
Entering edit mode

Thanks for your answer, but I need to be sure that I understand how R2 is :

the reverse complement or only reverse ?

ADD REPLY
0
Entering edit mode

No change. read2.fq.gz is used as is.

ADD REPLY

Login before adding your answer.

Traffic: 2435 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