SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Star with cufflinks Nino Bioinformatics 20 12-12-2016 09:06 AM
STAR: Alignment reported with less matches than parameter --outFilterMatchNmin AnneBiton Bioinformatics 0 07-29-2014 11:45 AM
STAR 2-pass ruggedtextile Bioinformatics 1 06-06-2014 01:27 AM
STAR vs Tophat (2.0.5/6) dvanic Bioinformatics 44 05-21-2014 07:08 AM
STAR-alignment/confusion with mismatches & multi-mapping parameter babi2305 Bioinformatics 9 05-07-2013 12:45 PM

Reply
 
Thread Tools
Old 09-08-2014, 07:10 PM   #1
sbdk82
Member
 
Location: USA

Join Date: Jul 2014
Posts: 25
Default STAR for alignment

I am trying to use STAR instead of Tophat, but I am bit confused. I already have the reference file and GTF file. So should I skip the "generate genome" and "generate genome with annotation" steps (4 and 5 in manual)?

If I directly use the following, should I keep my reference file under "/path/to/GenomeDir" ?

HTML Code:
STAR --genomeDir /path/to/GenomeDir --readFilesIn /path/
to/read1 [/path/to/read2] --runThreadN <n> --<inputParameterName> <input 
parameter value(s)>

I have paired end reads, reference file and GTF file. Please help me which command I need to use.
sbdk82 is offline   Reply With Quote
Old 09-08-2014, 11:21 PM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,478
Default

The "generate genome" and "generate genome with annotation" create an index that STAR uses for mapping. You can't skip those steps.

"/path/to/GenomeDir" is just an example name, you can keep things where ever you want and then just change that to match.
dpryan is offline   Reply With Quote
Old 09-09-2014, 08:22 AM   #3
sbdk82
Member
 
Location: USA

Join Date: Jul 2014
Posts: 25
Default

Thanks Devon !!

So I should use the following to build the index

STAR --runMode genomeGenerate --genomeDir <output directory> --genomeFastaFiles /path/reference.fa --runThreadN 5 --sjdbGTFfile <GTF file name>.
sbdk82 is offline   Reply With Quote
Reply

Tags
rna-seq, star

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 01:53 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO