[Tutor] Need advise on protein sequence alignment code in python 2.7 giving wrong output

Laura Creighton lac at openend.se
Tue Aug 11 21:14:18 CEST 2015


I wrote:	 
  >>> nwalign.global_align("PLEASANTLY", "MEANLY", gap_open=-5, gap_extend=-5, matrix='./BLOSUM62.txt')
  >>> ('PLEASANTLY', '-MEA--N-LY')

I forgot to mention that I got my BLOSUM62.txt from
  http://www.ncbi.nlm.nih.gov/Class/FieldGuide/BLOSUM62.txt

in case that matters.

Laura


More information about the Tutor mailing list