SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Illumina/Solexa (http://seqanswers.com/forums/forumdisplay.php?f=6)
-   -   Read1 and Read2 are not consistent (http://seqanswers.com/forums/showthread.php?t=31973)

priya 07-16-2013 02:09 AM

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?

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 (Post 110586)
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..:)


All times are GMT -8. The time now is 06:46 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.