Question: merge .fastq reads of different lengths and align
gravatar for User6891
5.6 years ago by
User6891270 wrote:


I have .fastq files from the same sample sequenced in two sequencing runs.

First run was on Illumina Hiseq and generated paired-end 2x100bp reads. Second run was on Miseq and generated paired-end 2x250bp reads. 

Can I just merge the .fastq files of both runs (using cat), keeping off course R1 and R2 separate and then align with BWA? Or are the different read lengths going to pose a problem?

bwa fastq • 3.0k views
ADD COMMENTlink modified 5.6 years ago by Devon Ryan96k • written 5.6 years ago by User6891270
gravatar for Devon Ryan
5.6 years ago by
Devon Ryan96k
Freiburg, Germany
Devon Ryan96k wrote:

Different lengths aren't a problem. Note that it's unclear if these were sequenced from independent libraries or not. If they were, then you might want to keep them separate if you want to do variant calling, since there may be different biases between them.

ADD COMMENTlink written 5.6 years ago by Devon Ryan96k
gravatar for David Langenberger
5.6 years ago by
David Langenberger9.5k wrote:

Of course, you can do that. Different read length do not influence the mapping algorithms. Actually, when clipping the adapter sequences, you end up with a lot reads having different read length.

ADD COMMENTlink written 5.6 years ago by David Langenberger9.5k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1358 users visited in the last hour