This is a test version of Biostars. For the public version, visit https://www.biostars.org.
the latest human reference genome fasta file

Hi all

I would like to download the latest human reference genome (GRCH38) in fasta and gtf format for my RNA seq analysis. I would like to know which database is the beast,Genbank version 21 or ensemble?

Where can I get the fasta file of whole genome of ensemble version?

Is the below link below contains this file?

ftp://ftp.ensembl.org/pub/release-77/fasta/homo_sapiens/dna/Homo_sapiens.GRCh38.dna.toplevel.fa.gz

Is there any difference between GenBank and Ensembl's alignment and annotation output result ?

Thanks in advance

rna-seq
• 26,467 views • link
updated by Ram Ram • written by hana hana

Hi all

I would like to download the latest human reference genome (GRCH38) in fasta and gtf format for my RNA seq analysis. I would like to know which database is the beast, Genbank version 21 or ensemble?

where can I get the fasta file of whole genome of Ensembl version?

Is the below link below contains this file?

ftp://ftp.ensembl.org/pub/release-77/fasta/homo_sapiens/dna/Homo_sapiens.GRCh38.dna.toplevel.fa.gz

Is there any difference between Genbank and Ensembl's alignment and annotation output result?

Thanks in advance

• 0 views • link
updated by Ram Ram • written by hana hana

1 answer

database is the beast?????

Yes, Its the one from ensembl.

You can download it from here, same way as you previously downloaded hg19 from UCSC

whole genome fasta

http://hgdownload.soe.ucsc.edu/goldenPath/hg38/bigZips/

GTFs from

http://genome.ucsc.edu/cgi-bin/hgTables?

• 0 views • link
updated by Ram Ram • written by Manvendra Singh Manvendra Singh

Which file is contained the whole genome file? hg38.2bit?

Thank you

• 0 views • link
updated by Ram Ram • written by hana hana

Yes, you need to convert it to fasta

You can get the utility program TwoBitToFa from here.

Once you downloaded it, you must change permissions first to allow it to be executed as a program.

Then you execute it from a terminal without arguments to see the options:

$ /path/to/twoBitToFa

twoBitToFa - Convert all or part of .2bit file to fasta
usage:
twoBitToFa input.2bit output.fa
• 0 views • link
updated by Ram Ram • written by Manvendra Singh Manvendra Singh

thank you very much

• 0 views • link
written by hana hana

Just directly download the fasta file. There's no need to deal with 2bit.

• 0 views • link
written by Devon Ryan Devon Ryan

Exactly, hana was asking about .2bit so I wrote that he can convert them as well.

• 0 views • link
written by Manvendra Singh Manvendra Singh

Sorry @Devon, from where I could find fasta file for each individual chromosome for hs37d5.fa ?

• 0 views • link
written by zizigolu zizigolu

Hi

I have already download the fasta and gtf files of hg38 from USCS database and run tophat.

But I have a problem with running cuffllinks . I got the below error

cufflinks  --GTF   genome.gtf   -o   /home/ra/cufflinks    /home/ra/accepted_hits.bam

[20:55:32] Loading reference annotation.
Error parsing strand (1) from GFF line:
uc001aaa.3    chr1    +    11873    14409    11873    11873    3    11873,12612,13220,    12227,12721,14409,        uc001aaa.3

Can you please tell me what dose it mean and how can I solve it

Thank you

• 0 views • link
updated by Ram Ram • written by hana hana

That's not a GTF file. You have to explicitly set the output format to GTF, otherwise you'll get all of the table columns as is.

• 0 views • link
updated by Ram Ram • written by Devon Ryan Devon Ryan

Log in to answer this question.

More posts like this
  • is there any difference between size of RefGenome in UCSC and Ensembl?
    written by star star
    I have downloaded the Human Refrence Genome from Ensembl and UCSC like below: UCSC: http://hgdownload.soe.ucsc.edu/goldenPath/hg38/bigZips/hg38.fa.gz Ensembl: ftp://ftp.ensembl.org/pub/release-99/fasta/homo_sapiens/dna/Homo_sapiens.GRCh38.dna.toplevel.fa.gz and when I make the uncompressed using gunzip, …
  • Where can I download the gff3 file for a specific human genome build? (hg38 release 77)
    written by O.rka O.rka
    I found this: ftp://ftp.ensembl.org/pub/release-98/gff3/homo_sapiens but there are only GTF files for release 77: ftp://ftp.ensembl.org/pub/release-77 Does anyone know where I can find this information?
  • convert GRCh38 FASTA .fa to .vcf
    written by LTDavid • LTDavid
    How do you convert the GRCh38.fa FASTA file to a VCF file? The file is from Ensembl primary_assembly [ ftp://ftp.ensembl.org/pub/release-97/fasta/homo_sapiens/dna/][1] [1]: ftp://ftp.ensembl.org/pub/release-97/fasta/homo_sapiens/dna/
  • Which Ensembl chromosome sequence file is used for the gff3 files?
    written by O.rka O.rka
    **I can't which of the following files were used to create the GFF3 files?** ftp://ftp.ensembl.org/pub/grch37/current/fasta/homo_sapiens/dna/Homo_sapiens.GRCh37.dna.toplevel.fa.gz ftp://ftp.ensembl.org/pub/grch37/current/fasta/homo_sapiens/dna/Homo_sapiens.GRCh37.dna.primary_assembly.fa.gz or an entirely different one from the ftp: ftp://ftp.ensembl.org/pub/grch37/current/fasta/homo_sapiens/dna/ …
  • Fasta length does not match Ensembl info
    written by firestar firestar
    I am calculating the total genome length from a fasta file using the following code `zcat genome.fa.gz | grep -v ">" | wc | awk …
  • Fasta file and GTF file for STAR alignment
    written by snp87 • snp87
    Hello, this is a very basic question but I was wondering if someone could help me understand if I've used the correct GTF file and …
  • choice of human genome before RNASeq mapping
    written by Omics data mining Omics data mining
    Dear all, I have RNASEQ data and want to map against the reference genome. I download human genome from ftp://ftp.ensembl.org/pub/release-91/fasta/homo_sapiens/dna/ Toplevel sequences unmasked: Homo_sapiens.GRCh38.dna.toplevel.fa.gz Here, …
  • hg38 Ensembl GTF file
    written by Gregor Rot Gregor Rot
    Dear all, what is the difference between the hg38 Ensembl GTF file: ftp://ftp.ensembl.org/pub/release-90/gtf/homo_sapiens/Homo_sapiens.GRCh38.90.gtf.gz and ftp://ftp.ensembl.org/pub/release-90/gtf/homo_sapiens/Homo_sapiens.GRCh38.90.chr.gtf.gz (note the .chr. addition in the name) Thank you, Gregor
  • No XIST in Emsembl cdna.fa
    written by WouterDeCoster WouterDeCoster
    I noticed that Salmon didn't count XIST in my RNA-seq data, and tracked it down to the transcriptome fasta from Ensembl which was used to …
  • which genome sequence fasta file should be used for builindg bowtie2 index ?
    written by Pei Pei
    <p>Hi all:</p> <p>this is a simple quesion:</p> <p>when building bowtie2 index for tophat2 analysis, which genome sequence file should I used:</p> <p>ftp://ftp.ensembl.org/pub/release-83/fasta/homo_sapiens/dna/<a href="ftp://ftp.ensembl.org/pub/release-83/fasta/homo_sapiens/dna/Homo_sapiens.GRCh38.dna.toplevel.fa.gz">Homo_sapiens.GRCh38.dna.toplevel.fa.gz</a></p> <p>or</p> <p>ftp://ftp.ensembl.org/pub/release-83/fasta/homo_sapiens/dna/<a …