Question: can I map multiple (>2) fastq reads against and assembly using bwa?
0
gravatar for sa91
12 days ago by
sa910
sa910 wrote:

Hi

I want to map the output of 4 separate illumina smrt cells against a pacbio reference genome using BWA. I've indexed the full genome.

I can try run the the equivalent of this line, but using additional reads:

bwa mem ref.fa read1.fq read2.fq > aln-pe.sam

However, if I try

bwa mem ref.fa read1.fq read2.fq read3.fq read4.fq > aln-pe.sam

include read3.fq or more nothing runs.. Instead I get the usage instructions for bwa.

Is there a way to map 4 illumina runs against the reference genome at once, or am I misunderstanding this process?

There are simply 4 runs to increase robustness.

ADD COMMENTlink modified 8 days ago by swbarnes27.8k • written 12 days ago by sa910

It is unclear what you do. Do you have different paired-end samples, all with R1 and R2? Should all the four samples go to the same SAM file or to different ones?

ADD REPLYlink written 12 days ago by ATpoint34k

These are 4 separate runs of illumina output that I am supposed to analyse and compare to my pacbio genome that I sequenced. I want all 4 in the same SAM file.

ADD REPLYlink written 8 days ago by sa910

4 separate illumina smrt cells against a pacbio reference genome using BWA

You are mixing terminologies here. PacBio sequencer uses SMRTCells where as Illumina sequencers use flowcells. Do you have 4 separate samples/lanes/flowcells worth Illumina data? Depending on the answer there would be a different solution.

ADD REPLYlink modified 8 days ago • written 8 days ago by genomax83k

Thanks for correcting me. I was getting confused and using the wrong terminology. It is from flow cells. I've been asked to map illumina short reads against my pacbio genome assembly to examine the coverage.

ADD REPLYlink written 2 days ago by sa910
1
gravatar for swbarnes2
8 days ago by
swbarnes27.8k
United States
swbarnes27.8k wrote:

Your first command lined is for paired end sequencing. Is that what you have?

If you do not have paired end sequencing, just two fastqs of the same single-end library, then your first command line is wrong.

If you have 4 fastqs that are all the same single end library, you can cat them all together before giving them to bwa

cat read*.fq | bwa mem ref.fa - > out.sam
ADD COMMENTlink written 8 days ago by swbarnes27.8k

Thank you! I didn't realise that command was specifically for paired end sequencing. That was where i went wrong. i'll cat the other fastqs together

ADD REPLYlink modified 2 days ago • written 2 days ago by sa910
1
gravatar for gayachit
12 days ago by
gayachit190
India
gayachit190 wrote:

Hi There are some posts similar to yours. Have a look.

bwa for multiple fastq files , BWA mem on multiple samples, perl script for BWA-mem on multiple different files

ADD COMMENTlink written 12 days ago by gayachit190
Please log in to add an answer.

Help
Access

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