How should I use skewer to trim paired end reads?
1
0
Entering edit mode
3.0 years ago
jrleary ▴ 190

I'm attempting to use skewer to trim paired end ChIPseq fastq files. I'm not solid on how to actually format the command used to run skewer as there's little in the way of user guides or documentation available. Here's what I'm using right now:

skewer -f auto -m pe -z -o ${sample}_trimmed.fastq.gz $read1 $read2

where $sample is the sample name, and $read1 and $read2 are the first and second paired end reads, respectively. This command only outputs one trimmed fastq file, and I believe it should output two, one for the first read and one for the second. Any tips would be greatly appreciated.

ChIP-Seq skewer • 1.5k views
ADD COMMENT
0
Entering edit mode

Are you open to using other trimming programs that have plenty of documentation? BBduk from BBMap suite, Trimmomatic or cutadapt.

ADD REPLY
0
Entering edit mode

I like skewer because of fast it is, but I'm open to using other tools if the documentation / speed tradeoff is worth it

ADD REPLY
1
Entering edit mode

bbduk is plenty fast and is multi-threaded to boot. You can use multiple cores.

ADD REPLY
1
Entering edit mode

I'll second this suggestion and throw in an additional plug for fastp, which is also multi-threaded.

ADD REPLY
0
Entering edit mode

I'll give both of these a look and compare their performance with that of skewer now that I've got it working. Thanks for the options.

ADD REPLY
1
Entering edit mode
3.0 years ago
GenoMax 123k

As I see it, you need to provide a base name option (e.g. -o trimmed_files) to have trimmed files saved with that basename for skewer.

ADD COMMENT
0
Entering edit mode

This rectified the problem, thank you!

ADD REPLY

Login before adding your answer.

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