I am trying to make a genome index for mapping and I keep getting this error:
I already followed the instructions given here:
However, the error still persists.
Here is the script I used:
STAR --runThreadN 11 \
--limitBAMsortRAM 75000000000 \
--runMode genomeGenerate \
--genomeDir "$Output/$Species" \
--genomeFastaFiles "$Reference_genome" \
--sjdbGTFfile "$Annotation_file" \
--sjdbGTFtagExonParentTranscript Parent \
--limitSjdbInsertNsj 4000000 \
--sjdbOverhang 149
The annotation file I used is a .gff3, below is the snippet of the annotation file.
Hsteno_chr1H Liftoff gene 26124 32988 . - . ID=Hsteno-30543
Hsteno_chr1H Liftoff mRNA 26124 32988 . - . ID=Hsteno-30543.1;Parent=Hsteno-30543;Dbxref=InterPro:IPR039199,InterPro:IPR028097,InterPro:IPR028091,MobiDBLite:mobidb-lite,PANTHER:PTHR28441,Pfam:PF14648,Pfam:PF14647
Hsteno_chr1H Liftoff exon 26124 26596 . - . ID=Hsteno-30543.1-exon1;Parent=Hsteno-30543.1
Hsteno_chr1H Liftoff exon 26717 28301 . - . ID=Hsteno-30543.1-exon2;Parent=Hsteno-30543.1
Hsteno_chr1H Liftoff exon 28510 28874 . - . ID=Hsteno-30543.1-exon3;Parent=Hsteno-30543.1
Hsteno_chr1H Liftoff exon 29021 29081 . - . ID=Hsteno-30543.1-exon4;Parent=Hsteno-30543.1
Hsteno_chr1H Liftoff exon 29280 29433 . - . ID=Hsteno-30543.1-exon5;Parent=Hsteno-30543.1
Any help would be appreciated. Thank you very much.
Best,
Sandy