Answers.Body: <p>Regarding a human genome reference, have a look at CHM13/HG002 from the Telomere-to-Telomere consortium:</p>
  <p><a href="https://github.com/marbl/CHM13#assembly-releases" rel="nofollow noreferrer">https://github.com/marbl/CHM13#assembly-releases</a></p>
  <p>The benefit of the T2T assembly is that it is fully-assembled across all repetitive regions, which improves the mapping accuracy. The pangenome consortium is working on additional T2T assemblies from other people, trying to capture as much human [structural] variation as possible.</p>
  <p><a href="https://humanpangenome.org/" rel="nofollow noreferrer">https://humanpangenome.org/</a></p>
  <p>The v2.0 assembly, adding chrY from HG002, has only just been released. From the looks of it, it is not yet annotated, but there's gene annotation for all other chromosomes via the v1.1 release:</p>
  <p><a href="https://github.com/marbl/CHM13#downloads" rel="nofollow noreferrer">https://github.com/marbl/CHM13#downloads</a></p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2022-02-28T22:05:06.423Z
Answers.LastActivityDate: 2022-02-28T22:05:06.423Z
Answers.UserDisplayName: gringer
Body: <p>When I trying to use bcftools mpileup to convert bam to vcf files,
  I try to use:</p>
  <pre><code>bcftools mpileup -Ob -o resources-broad-hg38-v0-Homo_sapiens_assembly38.fasta.fai input.bam
  </code></pre>
  <p>However, I got this error message:</p>
  <pre><code>[E::fai_build3_core] Failed to open the file resources-broad-hg38-v0-Homo_sapiens_assembly38.fasta.fai
  </code></pre>
  <p>Or, I would get another error message when I use another ref file:</p>
  <pre><code>[E::fai_build_core] Format error, unexpected &quot;c&quot; at line 1
  </code></pre>
  <p>Does anyone know how to solve these problems? Thanks in advance</p>
  <p>Also, could anyone recommend some human gene Reference? The one I am using now is from GATK pipeline (but I don't know if that is good enough).</p>
  <p>The version I am using is:</p>
  <pre><code>Program: bcftools (Tools for variant calling and manipulating VCFs and BCFs)
  License: GNU GPLv3+, due to use of the GNU Scientific Library
  Version: 1.10.2 (using htslib 1.10.2)
  </code></pre>
  
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2022-02-28T21:29:40.183Z
Comments.Text: What are you trying to do exactly? The `-o` option is to set an output file, and you seem to be trying to overwrite an existing fasta index (`.fai`) file. Did you mean to use `-f` to give the genome (which would be `resources-broad-hg38-v0-Homo_sapiens_assembly38.fasta` and not `resources-broad-hg38-v0-Homo_sapiens_assembly38.fasta.fai`)?
Comments.UserDisplayName: terdon
Tags: <bcftools>
Title: bcftools mpileup "Failed to open the file" questions
ThreadText: # bcftools mpileup "Failed to open the file" questions
  
  (Asked by: username_0 on 2022-02-28)
  
  When I trying to use bcftools mpileup to convert bam to vcf files,
  I try to use:
  
  ```
  bcftools mpileup -Ob -o resources-broad-hg38-v0-Homo_sapiens_assembly38.fasta.fai input.bam
  
  ```
  
  However, I got this error message:
  
  ```
  [E::fai_build3_core] Failed to open the file resources-broad-hg38-v0-Homo_sapiens_assembly38.fasta.fai
  
  ```
  
  Or, I would get another error message when I use another ref file:
  
  ```
  [E::fai_build_core] Format error, unexpected "c" at line 1
  
  ```
  
  Does anyone know how to solve these problems? Thanks in advance
  
  Also, could anyone recommend some human gene Reference? The one I am using now is from GATK pipeline (but I don't know if that is good enough).
  
  The version I am using is:
  
  ```
  Program: bcftools (Tools for variant calling and manipulating VCFs and BCFs)
  License: GNU GPLv3+, due to use of the GNU Scientific Library
  Version: 1.10.2 (using htslib 1.10.2)
  
  ```
  
  **Tags:** bcftools
  
  > **Comments:**
  > **username_1**: What are you trying to do exactly? The `-o` option is to set an output file, and you seem to be trying to overwrite an existing fasta index (`.fai`) file. Did you mean to use `-f` to give the genome (which would be `resources-broad-hg38-v0-Homo_sapiens_assembly38.fasta` and not `resources-broad-hg38-v0-Homo_sapiens_assembly38.fasta.fai`)?  
  
  ---
  
  ## 1 answers
  
  ### Answer by username_2 on 2022-02-28. Score: 2
  
  Regarding a human genome reference, have a look at CHM13/HG002 from the Telomere-to-Telomere consortium:
  
  <https://github.com/marbl/CHM13#assembly-releases>
  
  The benefit of the T2T assembly is that it is fully-assembled across all repetitive regions, which improves the mapping accuracy. The pangenome consortium is working on additional T2T assemblies from other people, trying to capture as much human [structural] variation as possible.
  
  <https://humanpangenome.org/>
  
  The v2.0 assembly, adding chrY from HG002, has only just been released. From the looks of it, it is not yet annotated, but there's gene annotation for all other chromosomes via the v1.1 release:
  
  <https://github.com/marbl/CHM13#downloads>
  
  
  ---
  
  