Question: merge .fastq reads of different lengths and align
gravatar for User6891
4.2 years ago by
User6891250 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 • 2.3k views
ADD COMMENTlink modified 4.2 years ago by Devon Ryan88k • written 4.2 years ago by User6891250
gravatar for Devon Ryan
4.2 years ago by
Devon Ryan88k
Freiburg, Germany
Devon Ryan88k 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 4.2 years ago by Devon Ryan88k
gravatar for David Langenberger
4.2 years ago by
David Langenberger8.7k 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 4.2 years ago by David Langenberger8.7k
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: 921 users visited in the last hour