Answers.Body: <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>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2024-07-02T16:18:09.957Z
Answers.Text: Weird, not working for me same or similar error.
Answers.UserDisplayName: N Brouwer
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2024-07-02T16:36:01.693Z
Answers.Text: may be your proxy settings have changed ?
Answers.UserDisplayName: Pierre
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2024-07-02T15:24:25.170Z
Answers.LastActivityDate: 2024-07-02T15:24:25.170Z
Answers.UserDisplayName: Pierre
Body: <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>
  
Tags: <vcf><samtools><1000genomes><tabix>
Title: 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
ThreadText: # 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?"](https://www.internationalgenome.org/category/tabix/) 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](https://github.com/samtools/htslib/issues/1033#issuecomment-597723270) 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](https://grch37.ensembl.org/Homo_sapiens/Tools/DataSlicer) 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](https://www.htslib.org/download/)
  * 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=<ID=AD,Number=.,Type=Integer,Description="Allelic depths for the ref and alt alleles in the order listed">
  ##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Read Depth (only filtered reads used for calling)">
  ##FORMAT=<ID=GL,Number=3,Type=Float,Description="Log-scaled likelihoods for AA,AB,BB genotypes where A=ref and B=alt; not applicable if site is not biallelic">
  ##FORMAT=<ID=GQ,Number=1,Type=Float,Description="Genotype Quality">
  ##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
  ##FORMAT=<ID=GD,Number=1,Type=Float,Description="Genotype dosage.  Expected count of non-ref alleles [0,2]">
  ##FORMAT=<ID=OG,Number=1,Type=String,Description="Original Genotype input to Beagle">
  ##INFO=<ID=AF,Number=.,Type=Float,Description="Allele Frequency, for each ALT allele, in the same order as listed">
  ##INFO=<ID=DP,Number=1,Type=Integer,Description="Total Depth">
  ##INFO=<ID=CB,Number=.,Type=String,Description="List of centres that called, UM (University of Michigan), BI (Broad Institute), BC (Boston College), NCBI">
  ##INFO=<ID=EUR_R2,Number=1,Type=Float,Description="R2 From Beagle based on European Samples">
  ##INFO=<ID=AFR_R2,Number=1,Type=Float,Description="R2 From Beagle based on AFRICAN Samples">
  ##INFO=<ID=ASN_R2,Number=1,Type=Float,Description="R2 From Beagle based on Asian Samples">
  #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 ?  
  
  ---
  
  