Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Comparing sets of contigs LizBent Bioinformatics 1 05-17-2012 02:39 AM
How are you comparing assemblies? Hobbe Bioinformatics 0 02-17-2011 10:38 PM
Comparing Transcriptomes poregide Bioinformatics 1 08-25-2010 11:23 PM
beginner comparing exomes drillbit Bioinformatics 0 06-16-2010 08:09 AM
Comparing bwa with shrimp jmj1091 Bioinformatics 3 09-24-2009 12:25 PM

Thread Tools
Old 10-11-2012, 01:11 AM   #1
Junior Member
Location: Utah

Join Date: Sep 2010
Posts: 1
Default comparing alignments

Hi all,

I have a question regarding what I though would be a simple thing to do, but cannot find an efficient method to implement it...

I have several individual alignments that were each assembled to the same annotated reference genome (large gaps in each as I targeted genes for sequencing), and I am trying to find a way to compare the various alignments to one another to see whether or not the same gene (or other) regions were recovered in each assembly (i.e. repeatability of my sequencing strategy). I don't care about minor differences (e.g. SNPs) between the individuals, just whether or not they 'hit' the same genomic region. It seems that the information is in the SAM/BAM file I have produced, but I can't seem to find an efficient way to extract and compare each region for all of the individuals.

Any help would be much appreciated.
sooper is offline   Reply With Quote
Old 10-11-2012, 04:46 AM   #2
Location: freiburg

Join Date: Apr 2010
Posts: 25

I do not know if this of great help, but if you are a little bit familiar with python (really just a little bit), the HTSEQ ( package could be of help. It delivers (in my opinion) very good and convenient methods to parse bam files and extract for instance the regions, where your reads aligned to. You can also very easily create .bed files that can be loaded into a genome browser. On the webpage you find a very good tutorial where the most common applications are explained in ready to use examples.
moritzhess is offline   Reply With Quote

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 06:11 AM.

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