How to ensure that the read names are identical except for /1 or /2
1
0
Entering edit mode
3.7 years ago

I am trying to do a de novo transcriptome with paired-end reads using Trinity. I have received the following error:

 pairs.K25.stats is empty.  Be sure to check your fastq reads and ensure that the read names are identical except for th
e /1 or /2 designation. at /usr/local/trinity/trinityrnaseq-Trinity-v2.4.0/util/insilico_read_normalization.pl line 904.
CMD finished (38 seconds)

My files are quite large and I'm using quite a few, I was wondering if anyone had any advice on how to identify where this error is originating from without manually scanning all of my files

next-gen assembly software error RNA-Seq • 664 views
ADD COMMENT
1
Entering edit mode
3.7 years ago

I manually looked through each read using:

head read1.fq
head read2.fq

I did this for each pair & noticed one pair were both coded as read 1. I observed the contents of the sample file throughout the preprocessing steps, saw where the error occurred and applied all the necessary changes from that point. I ran trinity again with no problems.

ADD COMMENT
0
Entering edit mode

You should have used a utility like repair.sh from BBTools. It is not practical to manually inspect things like large sequence data files.

ADD REPLY

Login before adding your answer.

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