SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
samtools merge bair Bioinformatics 5 10-10-2012 12:51 PM
samtools merge nans_bn Bioinformatics 2 07-22-2011 03:45 AM
samtools merge error underestimating the sequence counts Seq1234 Bioinformatics 0 11-22-2010 02:42 PM
samtools view error: CIGAR and sequence length are inconsistent (tophat/bowtie) glacierbird Bioinformatics 2 06-29-2010 01:58 AM
Samtools merge wangzkai Bioinformatics 1 05-01-2010 12:35 PM

Reply
 
Thread Tools
Old 01-17-2011, 09:35 PM   #1
Fumumumu
Junior Member
 
Location: Japan

Join Date: Jan 2011
Posts: 3
Default Error in samtools merge; different target sequence name

I am now making a script with bowtie and samtools, and
troubled in merging bam files from bowtie.
1) bowtie with every chromosome ewbt; output is sam file
2) sam to bam (every chromosome sam file)
3)* merge all bams
4) bam sort
5) bam index

#!/bin/sh
date=`date +"%Y%m%d_%H%M"`
samtools merge ${date}.bam ${date}_bam_chr1.bam, ${date}_bam_chr2.bam,... ,${date}_bam_chrY.bam

error massage was
[bam_merge_core] different target sequence name : 'chr1' != 'chr2' in file 20110115_0050_bam_chr2.bam

In this case, date was 20110115_0050.
What was wrong??
Fumumumu is offline   Reply With Quote
Old 01-26-2011, 05:22 AM   #2
jmrosa
Junior Member
 
Location: Madrid

Join Date: Sep 2010
Posts: 3
Default

Hi,

I had the same problem.

Check that chromosomes headers are in the same order in the headers of all the bam files.

Cheers
jmrosa is offline   Reply With Quote
Reply

Tags
samtools bowtie

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


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