Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: Using Bwa Sampe Using Illumna Pair End Reads
0
gravatar for Juliofdiaz
9.0 years ago by
Juliofdiaz • 130
Toronto, Ontario, Canada
Juliofdiaz • 130 wrote:

I am trying to map illumina paired reads using BWA and I am having some trouble understanding the workflow. I have two fastq files from the illumina run (read1.fq & read2.fq) and my reference genome in a fasta file(ref.fa). I have first indexed the reference file using:

bwa index ref.fa

and i got the following files

ref.fa.amb
ref.fa.ann
ref.fa.bwt
ref.fa.pac
ref.fa.rbwt
ref.fa.rpac
ref.fa.rsa
ref.fa.sa

Now I want to map my illumina reads,and I think the tool I'm looking for is bwa sampe. The help for sampe gives this format.

 bwa sampe [options] <prefix> <in1.sai> <in2.sai> <in1.fq> <in2.fq>

I dont quite understand where I'm supposed to get the .sai files from or what those files are. Thanks for any pointers.
illumina paired bwa • 6.0k views
ADD COMMENT • link •
written 9.0 years ago by Juliofdiaz • 130
4
gravatar for Zev.Kronenberg
9.0 years ago by
Zev.Kronenberg ♦ 11k
United States
Zev.Kronenberg ♦ 11k wrote:

The *.sai files are generated by first using:

bwa aln
ADD COMMENT • link written 9.0 years ago by Zev.Kronenberg ♦ 11k
1

Yes you first need to aln both reads independently. That will generate the sai files you need to create the SAM file.

ADD REPLY • link written 9.0 years ago by Zev.Kronenberg ♦ 11k

So first I need to align each read to the ref separately

ADD REPLY • link written 9.0 years ago by Juliofdiaz • 130
Please log in to add an answer.

Similar posts • Search »

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