dividing fastq file into 2 separate files
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.
• 1,902 views
•
link
updated
by
kk.mahsa
•
written
by
abel.arensberg •
•
written
by
abel.arensberg •
1 answer
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
Log in to answer this question.
More posts like this
-
Trimming pair end fastq files merged not interleaved
written by ricfoz •
-
scRNA-seq: Kallisto processing of bioproject data (fastq-dump)
written by bsmith030465
-
Illumina Overlapping R1/R2 reads, error-correction in python
written by Mick •
-
Combine paired-end fastq files
written by 112498262 •
-
strand identification from sam file
written by alireza346 •
-
bamtools stats output
written by GK1610
-
Trouble preparing paired-end (MiSeq 2x250), double-tagged and dual-indexed reads from metabarcoding…
written by jack1120 •
-
should I merge fastq files for different lanes before do QC?
written by Lila M
-
Demultiplex fastq miseq reads
written by bishnuadhikari323 •
-
Rna-Seq Paired End Analysis
written by KS
you can use SRAtoolkit to do that
No you can't use SRAtoolkit for the purpose OP is asking about.