SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
How to compare two WIG files?? tfcheng Bioinformatics 1 08-28-2013 01:23 AM
Compare two bam files. fabrice Bioinformatics 4 10-03-2011 03:20 PM
compare a few sra files rafi bondi Bioinformatics 0 05-16-2011 03:49 AM
Compare mapped reads from different aligners epigen Bioinformatics 3 06-25-2010 08:06 AM
Compare different files of SNVs Lien Bioinformatics 1 06-01-2010 06:10 PM

Reply
 
Thread Tools
Old 07-06-2010, 11:53 AM   #1
guykol
Junior Member
 
Location: Israel

Join Date: May 2010
Posts: 1
Default Compare results from different aligners (SAM files)

I would like to take 2 SAM files that are the results of alignment of the same input by different aligners and compare them.
Does anyone know of a code that already does that?

Thanks,

Guy.
guykol is offline   Reply With Quote
Old 07-07-2010, 05:08 AM   #2
epigen
Senior Member
 
Location: Germany

Join Date: May 2010
Posts: 101
Default

I already asked that and noone could provide me with any useful code, see
http://seqanswers.com/forums/showthread.php?t=5666
So I hope someone might give a positive answer here.
epigen is offline   Reply With Quote
Old 07-07-2010, 05:24 AM   #3
Dethecor
Member
 
Location: Germany

Join Date: May 2010
Posts: 24
Default Define "Compare"

Maybe if you could specify what you want to compare there would be a bigger chance that someone will provide more useful information.

To me it sounds like something you would do on your own with a little script in whatever programming language you like best
You could for example use samtools to sort both files and then use awk or whatever you like to check which reads are aligned in the same spot and which are not o.O
Dethecor is offline   Reply With Quote
Old 07-08-2010, 07:09 AM   #4
Bruins
Member
 
Location: Groningen

Join Date: Feb 2010
Posts: 78
Default

I'm interested too.

I will be looking at this sometime in the near future (before August), I'll post whatever methods I use/come up with then.
Bruins is offline   Reply With Quote
Old 04-11-2012, 12:42 PM   #5
mistro
Junior Member
 
Location: Cold Spring Harbor, New York

Join Date: Nov 2011
Posts: 1
Default

I came across this utility
http://genome.sph.umich.edu/wiki/BamUtil:_diff

$bam diff --in1 file1.sam/bam --in2 file2.sam/bam --out file1.fil2.diff.sam/bam

This will create three sam or bam files. proceeding from there...
mistro 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 02:40 AM.


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