Question: running bwa in pair end reads
0
gravatar for luzglongoria
8 months ago by
luzglongoria40
luzglongoria40 wrote:

Hi friends,

I would like to run bwa in my reads (12 pair end reads = 24 samples) and create a .sam file as an out put.

I want to use the command:

bwa sampe ref.fa aln_sa1.sai aln_sa2.sai read1.fq read2.fq > aln-pe.sam

My question is if I can use an unique command for all my 12 pair reads or I have to run 12 different commands (one per pair end read).

Any help is more than welcome :)

bwa rna-seq aligment sorted • 255 views
ADD COMMENTlink modified 8 months ago by Devon Ryan95k • written 8 months ago by luzglongoria40
1
gravatar for Devon Ryan
8 months ago by
Devon Ryan95k
Freiburg, Germany
Devon Ryan95k wrote:

You need a separate command per pair of reads.

ADD COMMENTlink written 8 months ago by Devon Ryan95k

Thanks !

The thing here is that I want a final .bam file for working with freebayes. Reading the freebayes manual it seems like I need a unique .bam file. Shall I merge the .sam files resulting from bwa, create a .bam file and then work with it on freebayes?

Does it make sense?

ADD REPLYlink written 8 months ago by luzglongoria40

Shall I merge the .sam files resulting from bwa, create a .bam file and then work with it on freebayes?

yes

ADD REPLYlink written 8 months ago by Pierre Lindenbaum129k
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: 1597 users visited in the last hour