anchor.diff:               allocate(seq(nv),stat=err)
                allocate(seqs(nv,ns),stat=err)
             end if
  -          call fasta_string2seq(string,seq,error_code,error_string)
  +          call fasta_string2seq(string,symbols_set,seq,error_code)
             if (error_code > 0) return
             seqs(:,k) = seq
             string = ""
  
positive.code:        ['A','C','G','T','U','-']
    character(len=kprot) :: protein_set='ACDEFGHIKLMNPQRSTVWY-'
    character(len=knuc)  :: nuc_acid_set='ACGTU-'
  
  contains
  
  
negative.code:              allocate(seq(nv),stat=err)
               allocate(seqs(nv,ns),stat=err)
            end if
            call fasta_string2seq(string,symbols_set,seq,error_code)
            if (error_code > 0) return
            seqs(:,k) = seq
            exit 
  