Hi,
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?
2 answers
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.
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.
Log in to answer this question.
-
MACS for a paired-end sequencing
written by KT •
-
Help Understanding PairedEnd read merging
written by m.radz •
-
Difference in R1 and R2 FASTQ files from RNA seq paired end
written by jaynishshah •
-
add UMI sequences to fastq read name
written by User6891
-
Merge fastq reads for several samples
written by rim.klabi •
-
hisat2 parameters for paired end and single end reads from same samples
written by bioinfo17 •
-
Merging illumina paired-end data Lxxx and R1/R2
written by euduca •
-
BWA MEM alignment output of splited fastq files differ from the original(unsplit) fastq file
written by manojkumar_bhosale
-
bwa mapping after fastx-trimmer
written by User6891
-
How to deal with demultiplexed Miseq pair-end (2*250bp) 16S data using QIIME?
written by nkuyfq
