<p>try with https instead of ftp ?</p>
<pre><code>$ tabix -h https://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20100804/ALL.2of4intersection.20100804.genotypes.vcf.gz 2:39967768-39967768  | cut -f1-8
##fileformat=VCFv4.0
##filedat=20101112
##datarelease=20100804
##samples=629
##description=&quot;Where BI calls are present, genotypes and alleles are from BI.  In there absence, UM genotypes are used.  If neither are available, no genotype information is present and the alleles are from the NCBI calls.&quot;
##FORMAT=&lt;ID=AD,Number=.,Type=Integer,Description=&quot;Allelic depths for the ref and alt alleles in the order listed&quot;&gt;
##FORMAT=&lt;ID=DP,Number=1,Type=Integer,Description=&quot;Read Depth (only filtered reads used for calling)&quot;&gt;
##FORMAT=&lt;ID=GL,Number=3,Type=Float,Description=&quot;Log-scaled likelihoods for AA,AB,BB genotypes where A=ref and B=alt; not applicable if site is not biallelic&quot;&gt;
##FORMAT=&lt;ID=GQ,Number=1,Type=Float,Description=&quot;Genotype Quality&quot;&gt;
##FORMAT=&lt;ID=GT,Number=1,Type=String,Description=&quot;Genotype&quot;&gt;
##FORMAT=&lt;ID=GD,Number=1,Type=Float,Description=&quot;Genotype dosage.  Expected count of non-ref alleles [0,2]&quot;&gt;
##FORMAT=&lt;ID=OG,Number=1,Type=String,Description=&quot;Original Genotype input to Beagle&quot;&gt;
##INFO=&lt;ID=AF,Number=.,Type=Float,Description=&quot;Allele Frequency, for each ALT allele, in the same order as listed&quot;&gt;
##INFO=&lt;ID=DP,Number=1,Type=Integer,Description=&quot;Total Depth&quot;&gt;
##INFO=&lt;ID=CB,Number=.,Type=String,Description=&quot;List of centres that called, UM (University of Michigan), BI (Broad Institute), BC (Boston College), NCBI&quot;&gt;
##INFO=&lt;ID=EUR_R2,Number=1,Type=Float,Description=&quot;R2 From Beagle based on European Samples&quot;&gt;
##INFO=&lt;ID=AFR_R2,Number=1,Type=Float,Description=&quot;R2 From Beagle based on AFRICAN Samples&quot;&gt;
##INFO=&lt;ID=ASN_R2,Number=1,Type=Float,Description=&quot;R2 From Beagle based on Asian Samples&quot;&gt;
#CHROM  POS ID  REF ALT QUAL    FILTER  INFO
2   39967768    rs11124691  T   A   .   PASS    DP=1991;AF=0.435;CB=UM,BI,BC,NCBI;EUR_R2=0.651;AFR_R2=0.742
</code></pre>

CC BY-SA 4.0
2024-07-02T16:18:09.957Z
Weird, not working for me same or similar error.
N Brouwer
CC BY-SA 4.0
2024-07-02T16:36:01.693Z
may be your proxy settings have changed ?
Pierre
CC BY-SA 4.0
2024-07-02T15:24:25.170Z
2024-07-02T15:24:25.170Z
Pierre
<p>I am trying to access <strong>1000 Genomes</strong> (1KG) data using tabix as per the 1KG <code>tabix</code> documentation <a href="https://www.internationalgenome.org/category/tabix/" rel="nofollow noreferrer">&quot;How do I get a genomic region sub-section of your files?&quot;</a> but am getting errors.</p>
<p><strong>Command I'm running</strong>:</p>
<p>I am getting errors when I run the sample command provided in the docs:</p>
<pre><code>tabix -h ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20100804/ALL.2of4intersection.20100804.genotypes.vcf.gz 2:39967768-39967768
</code></pre>
<p><strong>Error message</strong>:</p>
<p>The <strong>error</strong> I'm getting is:</p>
<pre><code>[E::bgzf_read] bgzf_read_block error -1 after 50219 of 52392 bytes
Could not load .tbi/.csi index of [...file name...]
</code></pre>
<p><strong>Manual download</strong></p>
<p>I tried to download a <code>.tbi</code> file and then run <code>tabix</code> as per this <a href="https://github.com/samtools/htslib/issues/1033#issuecomment-597723270" rel="nofollow noreferrer">comment on GitHub</a> but got an error</p>
<p>I ran:</p>
<pre><code>wget ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606_b151_GRCh37p13/VCF/00-All.vcf.gz.tbi
</code></pre>
<p>which seemed to work, but this:</p>
<pre><code>tabix ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606_b151_GRCh37p13/VCF/00-All.vcf.gz 1:17375-17375
</code></pre>
<p>returned an error:</p>
<pre><code>hts_close returned non-zero status:
</code></pre>
<p><strong>Other information</strong></p>
<p>This worked for me a few days ago with the command above and my own queries, but I am now having issues.</p>
<p>I'm not very savvy with command line tools so it could very well be me, though the GUI <a href="https://grch37.ensembl.org/Homo_sapiens/Tools/DataSlicer" rel="nofollow noreferrer"><strong>Data Slicer</strong> tool</a>  for accessing 1KG files doesn't seem to be working right now, so perhaps it's a server-side thing.</p>
<p><strong>Software</strong></p>
<ul>
<li>I'm on a PC using the Windows Linux Subsystem</li>
<li>I installed <code>tabix</code> via <code>bioconda</code> after having issues compiling it directly as specified in <a href="https://www.htslib.org/download/" rel="nofollow noreferrer">docs</a></li>
<li>Since having this issue I have re-installed <code>anaconda</code>, <code>bioconda</code>, <code>tabix</code>, <code>samtools</code>, and <code>htslib</code>.</li>
</ul>
<p><strong>Tabix info</strong></p>
<pre><code>which tabix
/home/brouwern/anaconda3/bin/tabix
</code></pre>
<pre><code>tabix --version

tabix (htslib) 1.3.1
Copyright (C) 2016 Genome Research Ltd.
</code></pre>

<vcf><samtools><1000genomes><tabix>
tabix errors when accessing 1000 Genomes data: "[E::bgzf_read] bgzf_read_block error -1 after 50219 of 52392 bytes" and Could not load .tbi/.csi index
# tabix errors when accessing 1000 Genomes data: "[E::bgzf_read] bgzf_read_block error -1 after 50219 of 52392 bytes" and Could not load .tbi/.csi index

(Asked by: username_0 on 2024-07-02)

I am trying to access 1000 Genomes (1KG) data using tabix as per the 1KG tabix documentation "How do I get a genomic region sub-section of your files?" but am getting errors.
Command I'm running:
I am getting errors when I run the sample command provided in the docs:
tabix -h ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20100804/ALL.2of4intersection.20100804.genotypes.vcf.gz 2:39967768-39967768

Error message:
The error I'm getting is:
[E::bgzf_read] bgzf_read_block error -1 after 50219 of 52392 bytes
Could not load .tbi/.csi index of [...file name...]

Manual download
I tried to download a .tbi file and then run tabix as per this comment on GitHub but got an error
I ran:
wget ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606_b151_GRCh37p13/VCF/00-All.vcf.gz.tbi

which seemed to work, but this:
tabix ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606_b151_GRCh37p13/VCF/00-All.vcf.gz 1:17375-17375

returned an error:
hts_close returned non-zero status:

Other information
This worked for me a few days ago with the command above and my own queries, but I am now having issues.
I'm not very savvy with command line tools so it could very well be me, though the GUI Data Slicer tool  for accessing 1KG files doesn't seem to be working right now, so perhaps it's a server-side thing.
Software

I'm on a PC using the Windows Linux Subsystem
I installed tabix via bioconda after having issues compiling it directly as specified in docs
Since having this issue I have re-installed anaconda, bioconda, tabix, samtools, and htslib.

Tabix info
which tabix
/home/brouwern/anaconda3/bin/tabix

tabix --version

tabix (htslib) 1.3.1
Copyright (C) 2016 Genome Research Ltd.

**Tags:** vcf, samtools, 1000genomes, tabix

---

## 1 answers

### Answer by username_1 on 2024-07-02. Score: 2

try with https instead of ftp ?
$ tabix -h https://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20100804/ALL.2of4intersection.20100804.genotypes.vcf.gz 2:39967768-39967768  | cut -f1-8
##fileformat=VCFv4.0
##filedat=20101112
##datarelease=20100804
##samples=629
##description="Where BI calls are present, genotypes and alleles are from BI.  In there absence, UM genotypes are used.  If neither are available, no genotype information is present and the alleles are from the NCBI calls."
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##FORMAT=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
##INFO=
#CHROM  POS ID  REF ALT QUAL    FILTER  INFO
2   39967768    rs11124691  T   A   .   PASS    DP=1991;AF=0.435;CB=UM,BI,BC,NCBI;EUR_R2=0.651;AFR_R2=0.742

> **Comments:**
> **username_0**: Weird, not working for me same or similar error.  
> **username_1**: may be your proxy settings have changed ?  

---

