Hi, To build a pairwise assembly alignments, do I need to go through all the steps as if I was building a msa? Or, is there a simpler way to do this? More specifically, can I skip the ref_assembly_self alignment?
BTW, I think there is a minor mistake on the http://code.google.com/p/pygr/wiki/NlmsaFromAxtPairwise wiki page. Shouldn't the seqlist be a seqdict? # a dictionary of genomic sequences seqdict = {'hg18':hg18, 'panTro2':panTro2, 'mm8':mm8, 'rn4':rn4, 'canFam2':canFam2} for k, v in seqdict.iteritems(): genomes[k] = seqdb.SequenceFileDB(v) Best wishes, Jenny --~--~---------~--~----~------------~-------~--~----~ You received this message because you are subscribed to the Google Groups "pygr-dev" group. To post to this group, send email to pygr-dev@googlegroups.com To unsubscribe from this group, send email to pygr-dev+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/pygr-dev?hl=en -~----------~----~----~----~------~----~------~--~---
