Sickle is giving empty trimmed file
0
0
Entering edit mode
2.2 years ago
salmon ▴ 10

Hello,

I am using RNA sequence data (forward and reverse fastq sequence) for my pure culture bacteria. I used the following sickle code and I am getting empty trimmed file. I am not sure where I am going wrong

sickle pe -f Data/K1_R1.fastq -r Data/ K1_R2.fastq -t sanger -o K1_trimmed.fastq -p K2_trimmed.fastq -s qtrim.unparied.fastq -q 30 -l 130
Warning: PE file 2 is shorter than PE file 1. Disregarding rest of PE file 1.

Result
FastQ paired records kept: 0 (0 pairs)
FastQ single records kept: 0 (from PE1: 0, from PE2: 0)
FastQ paired records discarded: 0 (0 pairs)
FastQ single records discarded: 0 (from PE1: 0, from PE2: 0)

I checked the Fastqc report of K1_R1 (Forward read) and K1_R2 (reverse read)as attached below for quality check.

enter image description here

................................................................

enter image description here

Sickle • 905 views
ADD COMMENT
0
Entering edit mode

I just want to let you know that sickle seems quite outdated. I’d go for trimmomatic or some other tools.

ADD REPLY
0
Entering edit mode

Agreed. I'd recommend fastp or bbduk.

ADD REPLY
0
Entering edit mode

-t sanger

You are not likely to have sanger formatted data if this is recent data.

Warning: PE file 2 is shorter than PE file 1. Disregarding rest of PE file 1.

That message is odd. FastQC indicates that there are an equal number of reads in your files. I second jkim suggestion of trying a different program.

ADD REPLY
0
Entering edit mode

Thank you so much for your reply.

ADD REPLY

Login before adding your answer.

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