dividing fastq file into 2 separate files
1
0
Entering edit mode
6.6 years ago

I have RNAseq data for both strands. how can I separate 2 strands from each other and make 2 fastq files? one for forward and one for reverse.

RNA-Seq • 1.4k views
ADD COMMENT
0
Entering edit mode

you can use SRAtoolkit to do that

ADD REPLY
0
Entering edit mode

No you can't use SRAtoolkit for the purpose OP is asking about.

ADD REPLY
0
Entering edit mode
6.6 years ago
GenoMax 143k

You first have to align the data first and then you can use splitsam.sh tool from BBMap suite to separate reads mapping to two strands. You will need to be mindful of multi-mappers, secondary alignments etc (if you have paired-end data). Finally use reformat.sh tool to recover the fastq files.

splitsam.sh mapped.sam plus.sam minus.sam unmapped.sam
reformat.sh in=plus.sam out=plus.fq
reformat.sh in=minus.sam out=minus.fq rcomp
ADD COMMENT

Login before adding your answer.

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