Indexing with STAR
0
0
Entering edit mode
5 months ago
Fizzah ▴ 20

Hello, I am working with RNA seq data and creating an index of reference genome Gossypium hirsutum by using STAR. STAR asks GTF annotation format while my file is GFF3. According to literature, in order to run GFF file I need to remove --sjdbOverhang 50 and also need to replace -sjdbGTFfile with sjdbGTFfeatureExon.

My question is:

  1. Removing "sjdbOverhang 50" from command will affect the index efficacy? As sjdbOverhang option needs to be specified for detecting possible splicing sites.

  2. Is replacing sjdbGTFfile with sjdbGTFfeatureExon mean indexing will be of exons only not of whole genome? Please explain what "sjdbGTFfeatureExon" meant?

This is the command I used:

STAR --runMode genomeGenerate --genomeDir indexes/Gh                                                                                                 --genomeFastaFiles /mnt/e/fizza\ data/S017679/trimmed/Ghirsutumv1.1_genome.fasta                                                  
--sjdbGTFfeatureExon exonfile /mnt/e/fizza\ data/S017679/trimmed/Ghirsutumv1.1_genome_repeat.gff3                      
--outFileNamePrefix Gh_align

Regards Fizzah

STAR RNAseq indexing • 258 views
ADD COMMENT

Login before adding your answer.

Traffic: 1315 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6