I'm trying to use cutadapt to trim my data. The adaptor sequences from Illumina only has a universal adapter, then one sequence for each index.
Example here:
TruSeq Universal Adapter
5’ AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT
TruSeq Adapter, Index 1 5
5’ GATCGGAAGAGCACACGTCTGAACTCCAGTCACATCACGATCTCGTATGCCGTCTTCTGCTTG
TruSeq Adapter, Index 2
5’ GATCGGAAGAGCACACGTCTGAACTCCAGTCACCGATGTATCTCGTATGCCGTCTTCTGCTTG
TruSeq Adapter, Index 3
5’ GATCGGAAGAGCACACGTCTGAACTCCAGTCACTTAGGCATCTCGTATGCCGTCTTCTGCTTG
When I read the guide from cutadapt for pair-end data, they suggest a file named ADAPTER_FWD (should be trimmed from the forward reads (filename.1.fastq.gz) and a file named ADAPTER_REV from the reverse reads (filename.2.fastq.gz).
So my question is: What sequence do I put in ADAPTER_FWD, and what sequence in ADAPTER_REV?
The command is:
cutadapt -a ADAPTER_FWD -A ADAPTER_REV -o out.1.fastq -p out.2.fastq reads.1.fastq reads.2.fastq
Example here:
TruSeq Universal Adapter
5’ AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT
TruSeq Adapter, Index 1 5
5’ GATCGGAAGAGCACACGTCTGAACTCCAGTCACATCACGATCTCGTATGCCGTCTTCTGCTTG
TruSeq Adapter, Index 2
5’ GATCGGAAGAGCACACGTCTGAACTCCAGTCACCGATGTATCTCGTATGCCGTCTTCTGCTTG
TruSeq Adapter, Index 3
5’ GATCGGAAGAGCACACGTCTGAACTCCAGTCACTTAGGCATCTCGTATGCCGTCTTCTGCTTG
When I read the guide from cutadapt for pair-end data, they suggest a file named ADAPTER_FWD (should be trimmed from the forward reads (filename.1.fastq.gz) and a file named ADAPTER_REV from the reverse reads (filename.2.fastq.gz).
So my question is: What sequence do I put in ADAPTER_FWD, and what sequence in ADAPTER_REV?
The command is:
cutadapt -a ADAPTER_FWD -A ADAPTER_REV -o out.1.fastq -p out.2.fastq reads.1.fastq reads.2.fastq
Comment