anchor.diff:      logical :: res
   
       res = .false.
  -    if (verify(string,nuc_acid_set) == 0) res = .true.
  +    if (verify(trim(string),trim(nuc_acid_set)) == 0) res = .true.
       
     end function is_nuc_acid
     
  
positive.code:   character(len=1) :: nuc_acid_alphabet(6) = &
         ['A','C','G','T','U','-']
    character(len=21) :: protein_set='ACDEFGHIKLMNPQRSTVWY-'
    character(len=5)  :: nuc_acid_set='ACGTU-'
    integer          :: aamap(239) ! 239 is n. of ascii codes (0 is null)
  
  contains
  
negative.code:   implicit none
    private 
    public :: protein_alphabet
    public :: nuc_acid_alphabet
    public :: fasta_read
    character(len=1) :: protein_alphabet(21) = &
         ['A','C','D','E','F','G','H','I','K',&
  