SEQanswers

Go Back   SEQanswers > Applications Forums > Epigenetics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Bowtie2 output with paired end reads sampie Bioinformatics 5 08-15-2014 05:35 AM
Using Bowtie2 for to align paired end reads sewellh Bioinformatics 5 06-05-2014 10:50 AM
Paired End Tags in Bowtie2 sameeranc Bioinformatics 2 05-09-2014 03:00 AM
Bowtie2 mapping and paired-end reads westeros Bioinformatics 3 03-22-2014 01:02 AM
Bowtie2- --al-conc: Writting to screen and file?? WiscoADL Bioinformatics 1 10-23-2012 10:39 AM

Reply
 
Thread Tools
Old 08-18-2016, 12:38 AM   #1
tesa
Junior Member
 
Location: The Netherlands

Join Date: Apr 2008
Posts: 6
Default Bowtie2 paired end --un-conc

Hi all,
As part of a normalization strategy for a ChIP-seq experiment, I am using Bowtie2 to align paired end ChiP-seq reads first to the fly genome and then take the unaligned reads and align to the human genome.

In order to output the reads that fail to align to fly concordantly I'm using --un-conc like this:

bowtie2 -p 8 --no-unal --un-conc unaligned.fq -x DMBDGP6.dna.primary_assembly -1 sample1_R1.fastq.gz -2 sample2_R2.fastq.gz -S output.sam

This produces two fastq files called unaligned.1.fq and unaligned.2.fq (.1 and .2 added to make per-mate filenames). I want to use these files downstream to align to the human as paired-end data. But, the output is different in length (different numbers of reads), which causes Bowtie2 to fail when mapping paired-end with -1 and -2.

Is there a way to have Bowtie2 output only the paired reads to these two files that fail to align concordantly? Or do I need to align downstream as independent unpaired reads, which completely defeats the purpose of using PE reads in the first place.

Thanks for the advice!
T
tesa is offline   Reply With Quote
Reply

Tags
alignment, mapping

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 08:53 AM.


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