Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Mapping SOLiD colorspace paired end reads NestorNotabilis SOLiD 10 12-12-2012 06:14 PM
Can paired-end mapping produce more reads than single-end ? warrenemmett Bioinformatics 13 03-20-2012 11:10 PM
should I reverse paired end reads before mapping? supermario Bioinformatics 2 03-06-2012 08:59 PM
Fastq: Paired end reads and mapping cedance Bioinformatics 7 06-18-2011 12:33 PM
BFAST mapping paired end reads. tanghz Bioinformatics 10 04-29-2011 06:29 AM

Thread Tools
Old 06-02-2011, 01:00 PM   #1
Junior Member
Location: Vancouver

Join Date: Jun 2011
Posts: 6
Default Paired-end reads mapping to different chromosomes...

Hi there,

I've been stumbling through some Illumina HiSeq data for a while, and while I've managed to answer some of the questions we wanted to ask, I'm having trouble with some others. One such question is whether it is possible to identify and print out only those paired-end reads that map to different chromosomes or have a larger insert size than expected. I was hoping to use this for identifying new retrotransposition events and rare fusions.

I have aligned the reads using bwa sampe, and aligned the forward and reverse separately with bwa samse.

My questions are as follows:

1. Can I use my paired-end aligned sam/bam file to pull out the reads where each mate maps to a different chromosome? I'm not sure if bwa biases against these and forces poorer alignments that are spatially close (I know that ambiguously aligned reads have the XA:Z: column that identifies all other possible mapping locations, but not sure if I can use this). If the sam file I have can't be used to answer this question, could anyone suggest a good program (fusionhunter?) that can do this?

2. If I can use the sam file I currently have (and not have to re-align with a specialized program), could anyone suggest basic code to pull out these sequences? I'm thinking of using either awk, Perl or R, but I'm still learning how to write in these languages.

I would want to match lines that have the same first column of the sam file (ignoring the #0 and #0/3 at the end which denotes forward or reverse read), then only pull out the lines where the thrid column (chromosome) is different (or the forth column (location in bp) is greater than 50kb).

Any help or advice would be greatly appreciated.

Thanks fo your time.
rareaquaticbadger is offline   Reply With Quote
Old 06-27-2011, 04:24 AM   #2
Junior Member
Location: china

Join Date: Apr 2011
Posts: 6

you can try deFuse
kissthefuture is offline   Reply With Quote

chromosome, fusion, hiseq, paired end, transposition

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:07 AM.

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