SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Questions on BWA zeam Bioinformatics 7 04-11-2011 10:03 AM
BWA questions totalnew Bioinformatics 16 06-25-2010 07:04 AM
MAQ easyrun problem and questions. answersseq Bioinformatics 0 01-04-2010 02:04 PM
bowtie and maq questions dukevn Bioinformatics 5 05-19-2009 08:47 AM
bwa seeding questions yxi Bioinformatics 0 05-18-2009 01:38 PM

Reply
 
Thread Tools
Old 08-04-2009, 08:14 PM   #1
baohua100
Senior Member
 
Location: Canada

Join Date: Jun 2008
Posts: 103
Default BWA-MAQ questions?

We have RNA-Seq paired-end data, the insert size is about 200bp , and the read length is different.

4_1.fq 50bp
4_2.fq 36bp

5_1.fq 75bp
5_2.fq 50bp

Can i run

bwa aln database.fasta 4_1.fq 4_2.fq 5_1.fq 5_2.fq > aln_sa.sai


How BWA do paired-end alignment for different read length?

And there are many situations like : one read is in one exon ,the other read is in other exon, the mapping position of these two reads on genome is 500, 3600. The distance is 3100 (there is a intron there) and BWA will discard these paired-end reads alignment because their distance is signicicant larger than 200bp?

Or i should map these reads as single reads?

Last edited by baohua100; 08-04-2009 at 08:44 PM. Reason: Questions about aligning paired-end RNA-seq reads using BWA‏
baohua100 is offline   Reply With Quote
Old 08-06-2009, 08:05 AM   #2
totalnew
Member
 
Location: Canada

Join Date: Apr 2009
Posts: 46
Default

Firstly, you will have to align fastq file individually to get correspond suffix array files, as

bwa aln database.fasta 4_1.fq > 4_1.fq.sai
bwa aln database.fasta 4_2.fq > 4_2.fq.sai
bwa aln database.fasta 5_1.fq > 5_1.fq.sai
bwa aln database.fasta 5_2.fq > 5_2.fq.sai

Then use sampe to generate results in sam format, I never align PE reads with different length, but I guess you may be able to get result sam files, just displaying different length at SEQ and QUAL fields. (Try and let us know what you get)

I have no idea about the insert size issue, as bwa is a new tool, everyone is doing experiements on that.
totalnew is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 10:53 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO