anchor.diff:      logical :: res
   
       res = .false.
  -    if (verify(string,protein_set) == 0) res = .true.
  +    if (verify(trim(string),trim(protein_set)) == 0) res = .true.
       
     end function is_protein
   
  
positive.code:     logical :: res
  
      res = .false.
      if (verify(string,nuc_acid_set) == 0) res = .true.
      
    end function is_nuc_acid
    
  
negative.code:   subroutine fasta_init(data_type)
      character(len=*), intent(in) :: data_type
      ! set the aa -> class mapping
      integer :: k
      
      aamap = 0
      select case(data_type)
      case('protein')
         do k = 1,21
            aamap(iachar(protein_alphabet(k))) = k
         end do
      case('nuc_acid')
         do k = 1,6
            aamap(iachar(nuc_acid_alphabet(k))) = k
         end do
      end select
      
  