I'm running a Perl script from the EMBL (found here https://github.com/EMBL-EBI-GCA/reseqtrack/blob/master/scripts/variation_data/calculate_allele_frq_from_vcf.pl) Under Ubuntu 16.10 I have installed Vcftools and Tabix as is required and both have been tested to work accordingly. I execute the script with the following command:
perl /home/[user]/calculate_allele_frq_from_vcf.pl \ -vcf /home/[user]/share/1000genomesPhase3VCFFiles/ALL.chr1.phase3_shap +eit2_mvncall_integrated_v5a.20130502.genotypes.vcf.gz \ -sample_panel /home/[user]/share/1000genomesPhase3VCFFiles/integrated_ +call_samples_v3.20130502.ALL.panel \ -out_dir /home/[user]/Desktop/AlleleFrequency \ -tabix /home/[user]/tabix-0.2.6 \ -vcftools_dir /home/[user]/vcftools_0.1.13 \ -region 1:1-10000 \ -pop CEU,FIN \
This returns the following error
Smartmatch is experimental at /home/[user]/calculate_allele_frq_from_v +cf.pl line 133. Smartmatch is experimental at /home/[user]/calculate_allele_frq_from_v +cf.pl line 144. sh: 1: /home/[user]/tabix-0.2.6: Permission denied Broken VCF header, no column names? at /home/[user]/vcftools_0.1.13/perl/Vcf.pm line 172. Vcf::throw(Vcf4_2=HASH(0x55761601d320), "Broken VCF header, no col +umn names?") called at /home/[user]/vcftools_0.1.13/perl/Vcf.pm line +866 VcfReader::_read_column_names(Vcf4_2=HASH(0x55761601d320)) called +at /home/[user]/vcftools_0.1.13/perl/Vcf.pm line 601 VcfReader::parse_header(Vcf4_2=HASH(0x55761601d320)) called at /ho +me/[user]/vcftools_0.1.13/perl/vcf-subset line 121 main::vcf_subset(HASH(0x557615b24558)) called at /home/[user]/vcft +ools_0.1.13/perl/vcf-subset line 12 Broken VCF header, no column names? at /home/[user]/vcftools_0.1.13/perl/Vcf.pm line 172. Vcf::throw(Vcf4_2=HASH(0x55fa53436f60), "Broken VCF header, no col +umn names?") called at /home/[user]/vcftools_0.1.13/perl/Vcf.pm line +866 VcfReader::_read_column_names(Vcf4_2=HASH(0x55fa53436f60)) called +at /home/[user]/vcftools_0.1.13/perl/Vcf.pm line 601 VcfReader::parse_header(Vcf4_2=HASH(0x55fa53436f60)) called at /ho +me/[user]/vcftools_0.1.13/perl/fill-an-ac line 45 main::fill_an_ac(undef) called at /home/[user]/vcftools_0.1.13/per +l/fill-an-ac line 9 sh: 1: /home/[user]/tabix-0.2.6: Permission denied Broken VCF header, no column names? at /home/[user]/vcftools_0.1.13/perl/Vcf.pm line 172. Vcf::throw(Vcf4_2=HASH(0x5607d059a190), "Broken VCF header, no col +umn names?") called at /home/[user]/vcftools_0.1.13/perl/Vcf.pm line +866 VcfReader::_read_column_names(Vcf4_2=HASH(0x5607d059a190)) called +at /home/[user]/vcftools_0.1.13/perl/Vcf.pm line 601 VcfReader::parse_header(Vcf4_2=HASH(0x5607d059a190)) called at /ho +me/[user]/vcftools_0.1.13/perl/vcf-subset line 121 main::vcf_subset(HASH(0x5607d00a1558)) called at /home/[user]/vcft +ools_0.1.13/perl/vcf-subset line 12 Broken VCF header, no column names? at /home/[user]/vcftools_0.1.13/perl/Vcf.pm line 172. Vcf::throw(Vcf4_2=HASH(0x556b8e952110), "Broken VCF header, no col +umn names?") called at /home/[user]/vcftools_0.1.13/perl/Vcf.pm line +866 VcfReader::_read_column_names(Vcf4_2=HASH(0x556b8e952110)) called +at /home/[user]/vcftools_0.1.13/perl/Vcf.pm line 601 VcfReader::parse_header(Vcf4_2=HASH(0x556b8e952110)) called at /ho +me/[user]/vcftools_0.1.13/perl/fill-an-ac line 45 main::fill_an_ac(undef) called at /home/[user]/vcftools_0.1.13/per +l/fill-an-ac line 9 /home/[user]/Desktop/AlleleFrequency/calculated_fra.process3085.1.1-10 +000.CEU_FIN
It appears that the script doesn't have permission to access Tabix for some reason. I've given the folder full read/write permission for the user executing the script (me). Any ideas?

In reply to Perl script cannot access folder by Billebaars

Title:
Use:  <p> text here (a paragraph) </p>
and:  <code> code here </code>
to format your post, it's "PerlMonks-approved HTML":



  • Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!
  • Titles consisting of a single word are discouraged, and in most cases are disallowed outright.
  • Read Where should I post X? if you're not absolutely sure you're posting in the right place.
  • Please read these before you post! —
  • Posts may use any of the Perl Monks Approved HTML tags:
    a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr
  • You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)
            For:     Use:
    & &amp;
    < &lt;
    > &gt;
    [ &#91;
    ] &#93;
  • Link using PerlMonks shortcuts! What shortcuts can I use for linking?
  • See Writeup Formatting Tips and other pages linked from there for more info.