cutadapt --cut option for paired-end reads?
2
0
Entering edit mode
8.0 years ago
blur ▴ 280

Hi, I have paired-end read data, and I need to remove 10 bp at the end of R1 ans the first 10bp in R2 - how can I do this without ruining the pairs? I tried cutting each read on its own, and ended up with unmatched sized reads that I couldn't align in bowtie...

cutadapt • 3.2k views
ADD COMMENT
1
Entering edit mode

Did you read the manual?

ADD REPLY
0
Entering edit mode

actually he is repeating the question with different title :)

C: uneven paired end reads in bowtie

ADD REPLY
1
Entering edit mode
7.9 years ago
Marcel M ▴ 100
cutadapt -u 10 -U -10 -o out.1.fastq.gz -p out.2.fastq.gz in.1.fastq.gz in.2.fastq.gz

The -u 10 removes 10 bp from the 5' end of read 1; the -U -10 removes 10 bp from the 3' end of read 2.

ADD COMMENT
0
Entering edit mode
8.0 years ago
Benn 8.3k

I would advice to use Trim Galore instead.

http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/

ADD COMMENT

Login before adding your answer.

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