SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
How to count number of mapped paired-end and single-end rna-seq reads repinementer Bioinformatics 8 01-06-2013 05:06 AM
50 bp Single end Vs 100 bp Single end Vs 50 bp Paired end dhanapala RNA Sequencing 4 06-08-2012 05:09 PM
Paired-end Bam from single-end aligned sam ramouz87 Bioinformatics 4 08-17-2011 12:55 PM
Does Cufflinks support single-end and paired end data together ? ersenkavak Bioinformatics 1 10-22-2010 07:26 AM

Reply
 
Thread Tools
Old 03-15-2013, 02:41 PM   #1
puggie
Member
 
Location: Europe

Join Date: Nov 2011
Posts: 52
Default Forcing paired-end data mapped as single-end in SAM

Im working with paired-end RNA-seq data, and for my purpose I map the forward and reverse reads individually. What is the easiest way to assemle reads into a single alignment file while forcing the original direction to be preserved? So e.g. I want forward and reverse information of reads in final alignments.
puggie is offline   Reply With Quote
Old 03-16-2013, 10:50 AM   #2
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

You know that aligners that make .sams will natively figure out how to combined your paired end data so that you still know which read cam from which fastq, right?

If you really have to work with these bams, set the flags yourself to indicate which come from read1, and which come from read 2, then merge the .bams.
swbarnes2 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 09:23 AM.


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