cutadapt: forward and reverse adapter sequence options
0
0
Entering edit mode
6.3 years ago
umn_bist ▴ 390

So I have a fasta file of all Illumina sequences that I am feeding cutadapt to trim adapters for my pair-end RNA-seq

cutadapt -q 10,10 -a "${adapter}" -o "${file1%_1.fastq}_1_trimmed.fastq" -p "${file2%_2.fastq}_2_trimmed.fastq" "${file1}" "${file2}"

First, if my files are from CGHub/TCGA, does the '_1' and '_2" at the end denote a particular orientation that I should keep in mind when doing QC?

In cutadapt -a option is for forward adapter sequence and -A is for the reverse adapter sequence (for the paired file). I am wondering if cutadapt will know to reverse the adapters sequences I feed automatically or if I have to reverse the adapter sequences manually.

So like thisĀ 

cutadapt -q 10,10 -a "${adapter}" -A "${adapter}" -o "${file1%_1.fastq}_1_trimmed.fastq" -p "${file2%_2.fastq}_2_trimmed.fastq" "${file1}" "${file2}"
RNA-Seq cutadapt • 3.4k views
ADD COMMENT
0
Entering edit mode

_1 = R1 = Read 1

_2 = R2 = Read 2

ADD REPLY
0
Entering edit mode

Right. But does it provide additional information other than read 1 and read 2? Ie, which one was sequenced in forward direction or reverse direction?

ADD REPLY
0
Entering edit mode

DNA is always sequenced in one direction 5'-->3'. You could consider R1 as forward and R2 as reverse, if that is how you want to think about the watson and crick strand.

ADD REPLY

Login before adding your answer.

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