This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to visualize a set of sequence alignments that all have the same reference sequence?

I have a handful of similar sequences with a reference that I have aligned each one to with Needleman. As oppose to looking at each alignment one at a time, is it possible to view them all in relation to the reference since they all share the same reference? I'm hoping for it to look like a multiple sequence alignment essentially (I've tried doing MSA on this set but sequences with large indels do not align correctly).

alignment
• 2,426 views • link
updated by bharata1803 bharata1803 • written by nchuang nchuang

1 answer

Maybe you can use IGV application. In that application, you can load the reference and load all of your alignment fils in bam format. This is the example : https://www.dropbox.com/s/65podoee1sdvyzr/ENST00000309035.png?dl=0

• 1 views • link
written by bharata1803 bharata1803

oh why didn't I think of that.thank you. What aligner would I use to generate global alignments in bam?

• 1 views • link
written by nchuang nchuang

There are several tools for aligning and almost every aligner will produce bam. As for me, I use bowtie2 aligner. Bowtie website

• 1 views • link
written by bharata1803 bharata1803

thank you

i ended up scripting something in python but using Needleman

• 1 views • link
written by nchuang nchuang

Log in to answer this question.

More posts like this