SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Tophat/Bowtie2 inconsistency in number of paired reads BenLerch RNA Sequencing 0 07-14-2012 11:52 AM
Reads not mapped in proper pair: Bowtie output RGP Bioinformatics 4 12-14-2011 06:21 AM
Base pair sequence order in paired reads seqmonkey Illumina/Solexa 2 10-07-2011 03:56 AM
Are the mate-pair reads always shadowed by their friends the paired-end (picture) seb567 Illumina/Solexa 1 07-15-2011 09:49 AM
How to check if reads are properly paired in mate-pair data? genepool_bee Bioinformatics 2 02-22-2011 01:07 AM

Reply
 
Thread Tools
Old 09-06-2013, 11:24 AM   #1
shawn.mek
Member
 
Location: Colorado

Join Date: Feb 2013
Posts: 12
Default Bowtie2 output fastq of paired reads where at least one in the pair aligns

Bowtie2 has the option on paired reads to spit out fastq files of things that fail to align concordantly (--un-conc <path>), and it has the option to spit out fastq files of things than align concordantly (--al-conc <path>).

What I want though, is it to spit out paired reads (a 1 and 2 fastq file) for those paired reads where one of the pair aligns. --I will then be aligning those output paired fastq files to a difference reference.

Is there a way to do this, or maybe to run the alignment and use samtools view on the bam file to extract reads that hit one of the pair and turn that back into a fastq or align what's in there to a different reference?

Thanks
shawn.mek is offline   Reply With Quote
Reply

Tags
bowtie 2, fastq, paired end reads, samtools

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 05:30 PM.


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