This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Dealing with gff files in STAR

Hi all,

I am really having trouble doing the alignment in STAR using gff file. I tried downloading gtf files but I think finding and getting the files into cluster from Ensembl and USCS was very confusing. I was able to easily download gff files from NCBI but then using it in STAR is more complicated. I know I have to specify extra parameter for gff files and here is the code I am using:

STAR --runThreadN 12 --runMode genomeGenerate \
  --genomeDir /home/n-z/valizad2/Sepsis_RNAseq/HumanData/GenomeDirectory \
  --limitGenomeGenerateRAM 206609344554 \
  --genomeFastaFiles GCF_000001405.38_GRCh38.p12_genomic.fna \
  sjdbGTFtagExonParentTranscript Parent \
  --sjdbGTFfile GCF_000001405.38_GRCh38.p12_genomic.gff \
  --sjdbOverhang 100

Adding sjdbGTFtagExonParentTranscript Parent gives me the error below:

EXITING: FATAL INPUT ERROR: empty value for parameter "sjdbGTFtagExonParentTranscript Parent" in input "Command-Line-Initial"
SOLUTION: use non-empty value for this parameter

I don't know why it's telling me the parameter is empty. I would appreciate your help.

Thanks!

rna-seq
• 6,816 views • link
updated by Biostar Biostar • written by valizad2 • valizad2

Did you really run the command with sjdbGTFtagExonParentTranscript Parent, instead of --sjdbGTFtagExonParentTranscript Parent?

• 1 views • link
written by h.mon h.mon

Yes, that was my mistake! Sorry. I ended up having to download the gtf files because I see in several places that using gff for gene counting is not a good idea.

Do you happen to know why I might have gotten this error when generating the index file using the gtf file?

Segmentation fault?

   Oct 22 18:52:35 ..... started STAR run
Oct 22 18:52:35 ... starting to generate Genome files
Oct 22 18:53:44 ... starting to sort Suffix Array. This may take a long time...
Oct 22 18:53:54 ... sorting Suffix Array chunks and saving them to disk...
Oct 22 21:15:05 ... loading chunks from disk, packing SA...
Oct 22 21:18:26 ... finished generating suffix array
Oct 22 21:18:26 ... generating Suffix Array index
Oct 22 21:21:13 ... completed Suffix Array index
Oct 22 21:21:13 ..... processing annotations GTF
/tmp/slurmd/job1091443/slurm_script: line 17: 3763135 Segmentation fault      STAR --runThreadN 12 --runMode genomeGenerate --genomeDir /home/n-z/valizad2/Sepsis_RNAseq/HumanData/GenomeDirectory --limitGenomeGenerateRAM 206609344554 --genomeFastaFiles /home/n-z/valizad2/Sepsis_RNAseq/HumanData/Human_GeneGTF/Homo_sapiens.GRCh38.dna.primary_assembly.fa --sjdbGTFfile /home/n-z/valizad2/Sepsis_RNAseq/HumanData/Human_GeneGTF/Homo_sapiens.GRCh38.94.gtf --sjdbOverhang 100
• 1 views • link
updated by h.mon h.mon • written by valizad2 • valizad2

Hi valizad2, please:

1) use the code button to format your post

2) do not post questions as answers, either open a new question, or post as a comment to my comment above, or edit your question.

I fixed for you this time.

• 1 views • link
written by h.mon h.mon

Sure, Thanks! Sorry about that.

• 1 views • link
written by valizad2 • valizad2

Try not to mention RAM allocation or check your space. I hope in this way. Thanks

• 1 views • link
written by k.kathirvel93 k.kathirvel93

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • Indexing the reference genome
    written by Yasemin • Yasemin
    ``` nohup STAR \ --runMode genomeGenerate \ --genomeDir /Users/yasi/mockexp/genome/genome_index/ \ --genomeFastaFiles /Users/yasi/mockexp/genome/GCF_000001405.39_GRCh38.p13_genomic.fna \ --sjdbGTFfile /Users/yasi/mockexp/genome/genomic.gtf --sjdbOverhang 80 > star_genome_generate.log 2>&1 ``` What I am missing …
  • STAR Genome index Error
    written by Prasanna • Prasanna
    I tried to run STAR command for RNAseq but I got the following error. /home/pshekar/RNAseq/STAR-2.7.11a/source/STAR --runMode genomeGenerate \ > > --genomeDir GRCh38.79.chrom1 \ > > …
  • STAR index issues to the step 4
    written by aka • aka
    I am trying to index a reference genome on a pipeline with snakemake and I made this rule: ```python rule reference_faidx_star: input: "../resources/reference/Qrob_PM1N.fa" output: "../resources/reference/ref/" …
  • What should be the value of Sjdboverhang in STAR?
    written by DareDevil DareDevil
    I run fastqc on the adapter removed fastq file and gets `Sequence length : 20-150`. Next I want to create index for reference genome using …
  • Creating a non-coding RNA index for STAR
    written by plberry • plberry
    I am trying to generate an index of ncRNAs for use with STAR aligner. I downloaded the long nc RNA GTF and FASTA files from …
  • STAR index file for GRCH37
    written by bioinformatics.queries bioinformatics.queries
    Hello everyone I am trying to generate the index file for STAR alignment using hg19 genome. I used the following commad STAR --runThreadN 30 --runMode …
  • RNA-seq STAR genome mapping
    written by waaaaaa • waaaaaa
    I'm working on the genome index file using STAR tool. Following is my code: `STAR --runMode genomeGenerate --genomeDir /GenomeDirectory --genomeFastaFiles D_genome.fasta --sjdbGTFfile D_assembl.gtf ` However, …
  • Segmentation Fault error with STAR
    written by valizad2 • valizad2
    Hi, I am trying make the index file using STAR. I am getting segmentation fault error. Can anyone help me with that? Thanks! Oct 22 …
  • 2-pass mapping with STAR
    written by Lidia • Lidia
    Hello everyone! I am new in the use of this tool and therefore I have some questions. I am working with RNA paired end read …
  • STAR indexing doubt
    written by Marshall • Marshall
    noobie question: just started an indexing of a genome: with the command STAR --runMode genomeGenerate --runThreadN 4 --genomeDir STAR_index/ --genomeFastaFiles genome.fa --sjdbGTFfile genome.gff3 --sjdbGTFtagExonParentTranscript Parent …