Jump to content

  • Log in with Facebook Log in with Twitter Log In with Google      Sign In   
  • Create Account

- - - - -

Sequence alignment: one large sequence (e.g. 1200 bp) and many small (e.g. 30 bp


  • Please log in to reply
3 replies to this topic

#1 seanspotatobusiness

seanspotatobusiness

    Veteran

  • Active Members
  • PipPipPipPipPipPipPipPipPipPip
  • 115 posts
1
Neutral

Posted 25 May 2011 - 06:15 AM

I would like to align a large sequence with many small ones (ASCII schematic below). MegAlign software tried to align all sequences to each other but I only want the small ones aligned to the large one. I could do them all separately, but that's a lot messier than it needs to be.

This is the sort of thing I'm trying to achieve:

---------------------------------------------------------------------------------------------------------------------------------------
-------------------
		   ----------------------
																  ----------------------------

Is LaserGene's MegAlign can't do it, can some other software?

#2 Micro

Micro

    Enthusiast

  • Active Members
  • PipPipPipPipPip
  • 56 posts
4
Neutral

Posted 06 June 2011 - 08:30 PM

Sequencher if you happen to know someone that has it, or you could try the freeware program ClustalX, or alternatively interfaces that use ClustalW eg. BioEdit or MEGA. In any of these programs you will need to reduce the "gap opening penalty" and/or the "gap extension penalty" in order to help it along to create the large gaps needed for this alignment.

Otherwise if this doesn't work (and assuming you don't have hundreds of fragments to align) you could select your long sequence and one of the short sequences and then align, repeat for each short sequence.

Cheers
M  :)

#3 Hamish McWilliam

Hamish McWilliam

    member

  • Members
  • Pip
  • 4 posts
0
Neutral

Posted 27 February 2012 - 11:49 PM

This is a typical situation in genome assembly where you have a reference sequence and many smaller sequences to be aligned against it. Multiple Sequence Alignment (MSA) tools do not behave well in a case such as this since they assume that all the sequences have a mutual relationship, which is not the case here. There are many specialist tools to deal with this, to begin with I would try the est2genome program in the EMBOSS suite. You can find online services providing this tool.

#4 bob1

bob1

    Thelymitra pulchella

  • Global Moderators
  • PipPipPipPipPipPipPipPipPipPip
  • 4,427 posts
233
Excellent

Posted 28 February 2012 - 11:51 AM

You can also try making what is known as a contig.  I think lasergene has a tool for that.




Home - About - Terms of Service - Privacy - Contact Us

©1999-2012 Protocol Online, All rights reserved.