SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
samtools sort stdout pederworning Bioinformatics 4 04-09-2011 09:02 AM
@HD SO field after samtools sort telos Bioinformatics 3 11-10-2010 07:12 AM
Cufflinks SAM file sort problem AdamB RNA Sequencing 2 09-14-2010 02:44 AM
SAMtools merge problem cliff Bioinformatics 1 06-12-2010 10:07 PM
samtools sort running extremely slow tsucheta Bioinformatics 2 06-11-2010 07:30 AM

Reply
 
Thread Tools
Old 03-12-2012, 11:41 PM   #1
pandafengye
Junior Member
 
Location: Hangzhou, China

Join Date: Jun 2011
Posts: 5
Default Postfix problem in samtools merge/sort

I am novice to NGS data analysis. I use SHRiMP2 to align SOLiD mate pair reads. Because I don't need the mate-pair information but just need the mapping information, I did the alignment of F3 / R3 reads separately. Each read has only the best hit recorded, with the parameter '--strata -o 1'. The read name in the obtained sam files didnot have the postfix, i.e., '_F3' and '_R3'. Then I used 'samtools view -bS' to convert SAM to BAM format and used 'samtools merge' and 'samtools sort' to do the merge and sort process.

The question is, the size of the merged file is nearly two fold of the sorted file. I guess that, because the reads of the same mate pair have identical names, only one of them is left in the sorted BAM file, am I right? How can I keep both reads in the sorted file?
pandafengye 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 12:23 PM.


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