Trimmomatic trimming 3' bases
2
0
Entering edit mode
3.5 years ago
rob.costa1234 ▴ 310

I want to trim 3' - 20 bases of FASTq in PE file but could not figure out the option @ http://www.usadellab.org/cms/index.php?page=trimmomatic

Any pointer will be highly appreciated.

Thanks

sequence • 1.4k views
ADD COMMENT
1
Entering edit mode
3.5 years ago

Your question is rather terse, but if you want the last 20 bases gone, it looks like the option you want is CROP.

ADD COMMENT
1
Entering edit mode
3.5 years ago
shelkmike ★ 1.2k

Trimmomatic has no option to trim a fixed number of bases from the 3' end. The option CROP trims the read to a specific length. For example, CROP:20 will make all reads 20 bp long if they were longer than that.

To do what you want, you can use the program Cutadapt with the following command:

cutadapt -u -20 -o reads_trimmed_R1.fastq -p reads_trimmed_R2.fastq reads_initial_R1.fastq reads_initial_R2.fastq
ADD COMMENT

Login before adding your answer.

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