AdapterRemoval Tool line error
0
0
Entering edit mode
22 months ago
Sarah • 0

Hi, everyone.

I am very new to bioinformatics and struggling with AdapterRemoval. I am on a MacBook Pro and using a conda environment in the terminal. I have two documents that I am trying to trim and pair. My code is below:

adapterremoval \
  --file1 /Users/*******/Desktop/PRO_CAN_NK15018_20M_R1.fastq.gz \
  --file2 /Users/*******/Desktop/PRO_CAN_NK15018_20M_R2.fastq.gz \
  --basename /Users/*******/Desktop/PRO_CAN_PAIRED_R3.fastq.gz

The code runs, but then stops with the error:

Error reading FASTQ record at line 14821369; aborting:
    partial FASTQ record; cut off after sequence
Aborting thread due to error.

It seems there is something wrong with the fastq.gz files I am inputting, but I'm not sure how to fix the error. The files I am inputting have been checked for quality with FastQC and downsampled using Reformat.sh in bbmap.

Thanks for any help.

bbmap adapterremoval fastqc • 905 views
ADD COMMENT
0
Entering edit mode

Have you looked at the 14821369th line?

ADD REPLY
0
Entering edit mode

Is there a way to view a specific line in the terminal? I'm not sure how to maneuver such a large file.

ADD REPLY
0
Entering edit mode

You can do so with zcat and sed

zcat my.fastq.gz | sed -n 14821360,14821380p
ADD REPLY

Login before adding your answer.

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