Question: STAR GenomeIndex output (which is the BAM file?)
0
gravatar for daniloceschin
23 months ago by
daniloceschin0 wrote:

Hi all, I'm using STAR to generate a Genome Index to then use it for a Genome-guided Trinity analisys. The parameter I have used are:

STAR --runThreadN 20 \
     --runMode genomeGenerate \
     --genomeDir tibGenIdex/ \
     --genomeFastaFiles nanPar1.fa \
     --genomeChrBinNbits 10 \
     --outSAMtype BAM SortedByCoordinate \
     --sjdbGTFfile tibetan.gtf \

The output files are:

115K  - chrLength.txt
457K - chrNameLength.txt
343K - chrName.txt
271K - chrStart.txt
9.1M - exonGeTrInfo.tab
3.8M - exonInfo.tab
635K - geneInfo.tab
2.0G - Genome
567 - genomeParameters.txt
16G - SA
1.5G - SAindex
6.4M - sjdbInfo.txt
5.9M - sjdbList.fromGTF.out.tab
5.9M - sjdbList.out.tab
2.6M - transcriptInfo.tab

Which is supouse to be the BAM file?

There are huge files: Genome: 2.0G; SA: 16G and SAindex: 1.5G. Should I mantain those files?

Thank you all for any help! Best regards,

ADD COMMENTlink modified 23 months ago • written 23 months ago by daniloceschin0
1

This produces a genome index which is used for alignment of regular fastq data independent of Trinity. You can't use this index for Trinity directly.

From Trinity support page.

Users must provide read alignments to Trinity as a coordinate-sorted bam file. Use GSNAP, TopHat, STAR or other favorite RNA-Seq read alignment tool to generate the bam file, and be sure it's coordinate sorted by running 'samtools sort' on it.

ADD REPLYlink modified 23 months ago • written 23 months ago by genomax72k

Thank you! Now I understand how it works. You have made me think! :)

Thank you all for your fast replay.

ADD REPLYlink written 23 months ago by daniloceschin0
2
gravatar for yhoogstrate
23 months ago by
yhoogstrate50
Netherlands
yhoogstrate50 wrote:

genomeGenerate doesn't prdoce a SAM/BAM file, STAR only produces one if you run it in normal mode upon some FASTQ/FASTA data and the generated index.

ADD COMMENTlink modified 23 months ago • written 23 months ago by yhoogstrate50
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1968 users visited in the last hour