def vcf_annotator(self):
        """Vcf annotator"""

        tstart = datetime.now()

        # python ../scripts/annotate_vcfs.py -i mm13173_14.ug.target1.vcf -r 1000genomes dbsnp138 clinvar esp6500 -a ../data/1000genomes/ALL.wgs.integrated_phase1_v3.20101123.snps_indels_sv.sites.vcf.gz ../data/dbsnp138/00-All.vcf.gz ../data/dbsnp138/clinvar_00-latest.vcf.gz ../data/ESP6500/ESP6500.vcf.gz

        # command = 'python %s/vcf_annotator_parallel.py -n %s -i sanity_check/checked.vcf -r 1000genomes dbsnp clinvar esp6500 -a %s %s %s %s 2>log/pynnotator.log' % (scripts_dir, pynnotator_cores, genomes1k, dbsnp, clinvar, esp)
        # self.shell(command)

        resources = "genomes1k dbsnp clinvar esp6500 ensembl_phen ensembl_clin hgmd"  # 
        resources = resources.split(' ')
        annfiles = [
            "%s/1000genomes/%s" % (settings.data_dir, settings.genomes1k_file),
            "%s/dbsnp/%s" % (settings.data_dir, settings.dbsnp_file),
            "%s/dbsnp/%s" % (settings.data_dir, settings.clinvar_file),
            "%s/esp6500/%s" % (settings.data_dir, settings.esp_final_file),
            "%s/ensembl/%s" % (settings.data_dir, settings.ensembl_phenotype_file),
            "%s/ensembl/%s" % (settings.data_dir, settings.ensembl_clinically_file),
        ]

        if os.path.exists(settings.hgmd):
            annfiles.append("%s/hgmd/%s" % (settings.data_dir, settings.hgmd),)
        #



        # annfiles = " ".join(annfiles)
        # annfiles = ["%s/1000genomes/%s" % (settings.data_dir, settings.genomes1k_file)]

        annotator_obj = vcf_annotator.VCF_Annotator(self.vcf_file, annfiles, resources, settings.vcf_annotator_cores)

        annotator_obj.run()

        tend = datetime.now()
        execution_time = tend - tstart