GFF and Genome Fasta file

Hello,

Please let me know if the gff and genome fasta file have to be the same when running a cuffdiff analysis. for example,

can I use these:

For the genome version: ftp://ftp.ensemblgenomes.org/pub/release-28/plants/fasta/arabidopsis_thaliana/dna/Arabidopsis_thaliana.TAIR10.28.dna.toplevel.fa.gz

For the gff files:

Please note that chromosomes are labeled as 'ChrX', while the fasta is labeled as 'X' ftp://ftp.arabidopsis.org/Maps/gbrowse_data/TAIR10/TAIR10_GFF3_genes_transposons.gff

rna-seq gff fasta cuffdiff genome

• 3,067 views • link

updated 5.3 years ago by Ram

• written 10.8 years ago by aj123

1 answer

Yes, the chromosome names must be exact matches (unless things have changed recently), so having "ChrX" vs "X" in different files will be an issue.

The *nix command-line utility sed would be the best option, and it would be simpler (fewer changes) to change the fasta file chromosome names to match the gff names.

• 1 views • link

updated 5.3 years ago by Ram

• written 10.8 years ago by Joseph Pearson

yes, thanks sed and awk seem to work best in such cases!

• 0 views • link

written 10.8 years ago by aj123

More posts like this

indexing the drosophila genome for kallisto
written 2.1 years ago by gogeni5529 •

I want to index the drosophila genome to run kallisto (v. 0.50). I was wondering which fasta file I should use, when downloading it from …
Salmon Index
written 3.2 years ago by ExtentHonest56 •

Hello, I had a question about indexing with Salmon. I saw on the Salmon github pipeline that you can use the cDNA sequence with no …
Index genome for Triticum aestivum
written 4.0 years ago by mostafarafiepour

Hi all, I tried many times to create the wheat genome index file, but due to the size of the wheat genome, I faced a …
Identify transcripts code for longest protein from gene annotation file
written 6.9 years ago by waqaskhokhar999

I have reference annotation file of Arabidopsis thaliana and I am interested to identify extract transcipts that code for longest protein isoform and then extract …
Reference sequence for indica Rice Paired End reads assembly using HISAT2
written 7.5 years ago by sharadlko1 •

I am working on RNA-Seq based analysis of an indica rice variety. For reference based assembly of Paired End reads using HISAT2, I found indicia …
Transcriptome analysis for differential gene expression
written 7.6 years ago by umeshtanwar2 •

Hi. I am new for performing RNA-Seq. I have got the sequencing raw reads (PE) for my samples from Arabidopsis thaliana plants. I have analysed …
Fasta length does not match Ensembl info
written 7.8 years ago by firestar

I am calculating the total genome length from a fasta file using the following code `zcat genome.fa.gz | grep -v ">" | wc | awk …
Drosophila genome indexing using STAR for mapping reads
written 8.2 years ago by ishento •

Could anyone let me know which FASTA files that can be used to build genome index for analysis using STAR. 1- In Ensemble there is …
Fasta file and GTF file for STAR alignment
written 8.2 years ago by 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 …
which genome sequence fasta file should be used for builindg bowtie2 index ?
written 10.5 years ago by Pei

Hi all: this is a simple quesion: when building bowtie2 index for tophat2 analysis, which genome sequence file should I used: 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> or ftp://ftp.ensembl.org/pub/release-83/fasta/homo_sapiens/dna/<a …