How should I use skewer to trim paired end reads?
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.
1 answer
Log in to answer this question.
More posts like this
-
Creating fastq subsets from existing files
written by Zack •
-
Issues with adapter trimming when there are no adapters?
written by DareDevil
-
Galaxy RNASeq Paired-End Reads cutadapt adapter sequence setting
written by schen •
-
Explanation for kmer enrichment after skewer trimming of the RNA-seq reads
written by tobytaogf •
-
Coverage Discrepancy for samtools mpileup
written by kspata
-
How to avoid outputting untrimmed second reads when demultiplexing paired-end reads with cutadapt? …
written by ajfaure •
-
Use Bowtie to trim reads in fastq files
written by Jackie
-
beginner question BWA program
written by reza
-
Extracting unmapped paired end from an alignment
written by Picasa
-
bwa sampe error [fread] Unexpected end of file
written by gsuryawanshi •
Are you open to using other trimming programs that have plenty of documentation? BBduk from BBMap suite, Trimmomatic or cutadapt.
I like skewer because of fast it is, but I'm open to using other tools if the documentation / speed tradeoff is worth it
bbdukis plenty fast and is multi-threaded to boot. You can use multiple cores.I'll second this suggestion and throw in an additional plug for fastp, which is also multi-threaded.
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.