running bwa in pair end reads
1
0
Entering edit mode
17 months ago
luzglongoria ▴ 40

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 :)

RNA-Seq bwa sorted aligment • 456 views
ADD COMMENT
1
Entering edit mode
17 months ago

You need a separate command per pair of reads.

ADD COMMENT
0
Entering edit mode

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 REPLY
0
Entering edit mode

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

yes

ADD REPLY

Login before adding your answer.

Traffic: 1735 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6