SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Aligning paired end Illumina data with Bowtie kopardev Bioinformatics 5 03-29-2012 09:46 AM
Simulating paired-end reads and bowtie alignments droog_22 Bioinformatics 0 02-09-2012 08:53 AM
The best output format for paired-end alignment using bowtie chedonat Bioinformatics 4 02-01-2012 02:03 PM
Bowtie output from paired end reads godzilla07 Bioinformatics 0 01-06-2011 12:36 PM
BWA alignment for paired end reads AvinashP Genomic Resequencing 2 06-08-2010 04:11 AM

Reply
 
Thread Tools
Old 09-19-2011, 07:49 AM   #1
empyrean
Member
 
Location: EU

Join Date: Sep 2010
Posts: 52
Default Bowtie Illumina paired end reads alignment

Hello all,

I am new to NGS and trying to align illumina paired end data to a region which is of approximately 100kb where as the sequenced genome is about 1gb. Here i wanted to identify single pairs mapped to end of this region where as its other pair doesnot. How could i achieve this easily? Which is better bowtie or BWA? i used following command in bowtie but not working properly.

Code:
 bowtie -t -a -m 5 -S -p 10 -X 600 -1 001.fastq -2 002.fastq --al aligned_out
empyrean is offline   Reply With Quote
Old 09-20-2011, 07:06 AM   #2
empyrean
Member
 
Location: EU

Join Date: Sep 2010
Posts: 52
Default

can someone help?
empyrean is offline   Reply With Quote
Old 09-20-2011, 08:38 AM   #3
westerman
Rick Westerman
 
Location: Purdue University, Indiana, USA

Join Date: Jun 2008
Posts: 1,104
Default

When you say that something is "not working" it is useful to the rest of us if you can tell us the error message or other reason the program is not working.
westerman is offline   Reply With Quote
Old 09-20-2011, 10:51 AM   #4
empyrean
Member
 
Location: EU

Join Date: Sep 2010
Posts: 52
Default

I should have done that, but now bowtie ran succesfully on my data and gave sam output.

Code:
HWI-ST5 1:Y:0:CGAAGT    77      *       0       0       *       *       0       0       ANAANNNNNGGCTCC   ###########   XM:i:0
I got something like the top output in sam format. But i need the information of a read in a pair matches to the reference where as its another pair doesnot. The third column which says * was supposed to be a query_identifier and when i got the reads which are not equal to *, i got all pairs but not singles which map to reference. Can i get this info from bowtie?
empyrean 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 03:21 AM.


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