Read1 and Read2 are not consistent

priya · 07-16-2013, 02:09 AM

Hi, I got the sequencing data from Illumina HiSeq2000. Its pairend sequencing.
The two fastq files(read1 and read2) does not contain equal number of lines.
19166068 read1.fastq
19295372 read2.fastq

Is there any way that I can make them equal by removing the lines which are not there in read1 file?

dpryan · 07-16-2013, 02:48 AM

Yeah, have a look at this thread on Biostars for some example commands. Allegedly, you can also use cmpfastq, which is a perl script, to do this too (I've never used this so I can't say that it works properly).

In the future, use a trimmer that handles both reads of a pair at once so you can avoid these issues.

priya · 07-16-2013, 03:42 AM

Quote:

Originally Posted by dpryan

Yeah, have a look at this thread on Biostars for some example commands. Allegedly, you can also use cmpfastq, which is a perl script, to do this too (I've never used this so I can't say that it works properly).

In the future, use a trimmer that handles both reads of a pair at once so you can avoid these issues.

Thanks for your reply.I will look into it..

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
BWA aligned pairs read1=XT:A:U, read2=XT:A:R	jflowers	Bioinformatics	1	07-29-2014 11:06 PM
Paired end- search for complement of Read2 in Read1	flacchy	Introductions	9	05-23-2013 07:55 AM
refGene.txt sizes and frameshifts of exons not consistent	jdenvir	Bioinformatics	0	04-09-2013 12:36 PM
Cufflinks FPKM output not consistent with different thread numbers	Rachel Hillmer	RNA Sequencing	1	07-23-2012 04:27 PM
Current read1, index, and read2 primers for HiSeq2000	SeqVicious	Illumina/Solexa	1	09-26-2011 03:12 PM

07-16-2013, 02:09 AM	#1
priya Member Location: sweden Join Date: Apr 2013 Posts: 57	Read1 and Read2 are not consistent Hi, I got the sequencing data from Illumina HiSeq2000. Its pairend sequencing. The two fastq files(read1 and read2) does not contain equal number of lines. 19166068 read1.fastq 19295372 read2.fastq Is there any way that I can make them equal by removing the lines which are not there in read1 file?

07-16-2013, 02:48 AM	#2
dpryan Devon Ryan Location: Freiburg, Germany Join Date: Jul 2011 Posts: 3,480	Yeah, have a look at this thread on Biostars for some example commands. Allegedly, you can also use cmpfastq, which is a perl script, to do this too (I've never used this so I can't say that it works properly). In the future, use a trimmer that handles both reads of a pair at once so you can avoid these issues.