SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Can Cuffdiff treat paired-end and single-end reads at the same time? zun RNA Sequencing 3 06-12-2012 05:37 PM
Can paired-end mapping produce more reads than single-end ? warrenemmett Bioinformatics 13 03-20-2012 11:10 PM
RNA-seq: Replicates, single-end, paired-end story pasta Bioinformatics 2 07-04-2011 11:51 PM
Sam flags for bwa-aligned paired end reads with identical + / - strand coordinates spark Bioinformatics 0 03-09-2011 04:00 AM
Does Cufflinks support single-end and paired end data together ? ersenkavak Bioinformatics 1 10-22-2010 07:26 AM

Reply
 
Thread Tools
Old 07-20-2010, 03:09 AM   #1
ramouz87
Member
 
Location: Doha, Qatar

Join Date: Oct 2009
Posts: 35
Default Paired-end Bam from single-end aligned sam

Hi,
I have 2 sam(reverse and farward) files generated by bowtie.
I would like to generate a paired-end bam file from the 2 sam files.
using samtools merge -nr generate only a file with the concatenation of the singles reads from bothb file but do not do the pairing
tried also samtools fixmate but it's not working.
I've looked at the spec for SAM format, in 1.4.4. Storing paired-end reads the following is stated : "A tool is also provided to reconstruct mating information from BAM, although this is done at the cost of intensive computation and large disk space."
Is there any tool that allow to do the pairing of single reads based on the read ID ?
Thanks in advance for your help,
Regards,
Ramzi

Last edited by ramouz87; 07-20-2010 at 05:20 AM.
ramouz87 is offline   Reply With Quote
Old 08-17-2011, 06:44 AM   #2
oghabian
Junior Member
 
Location: Helsinki

Join Date: Jun 2011
Posts: 2
Default

Hi everyone!
Did you manage to get an answer for your question dear ramouz87?
I have a similar questions except that I have two .bam files that include the pairs that are mapped, separately. I want to associate the pairs and have the results as a single .bam file.

Does anyone know how that is possible?
oghabian is offline   Reply With Quote
Old 08-17-2011, 10:17 AM   #3
efoss
Member
 
Location: Seattle

Join Date: Jul 2011
Posts: 98
Default

I just came to this site to ask the same question - I have two sam files from Bowtie that come from paired end reads and I want to create one sam or bam file that incorporates the paired end information.

Eric
efoss is offline   Reply With Quote
Old 08-17-2011, 12:22 PM   #4
lh3
Senior Member
 
Location: Boston

Join Date: Feb 2008
Posts: 693
Default

Please redo the alignment with a mapper supporting proper paired-end mapping. Merging single-end alignment afterwards always gives you worse results.
lh3 is offline   Reply With Quote
Old 08-17-2011, 12:55 PM   #5
dreesbl
Registered Vendor
 
Location: Seattle, WA

Join Date: Nov 2009
Posts: 6
Default

@ efoss, oghabian, ramouz87

I don't think there's any reasonable way to get from one to the other. The alignments you got by treating the reads as singles are not the same as what you'll get if you align them as paired. You can either redo the alignment as paired (the best solution) or forge ahead with what you've got treating the reads as single end.

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


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