Hi,
Lately, I was playing a little bit with STAR and cufflinks to perform transcriptome assembly. I have a gtf file but it is not the best one (I study a non-well annotated species ). Here is the parameters I used for STAR. FYI, I have 2x50bp strand specific reads (~15 samples).
STAR --genomeDir $stargenomeDir --outFilterIntronMotifs RemoveNoncanonicalUnannotated --alignSJDBoverhangMin 1 --outFilterType BySJout --outFilterMismatchNmax 5 --readFilesIn $r1 $r2 --readFilesCommand zcat
and for cufflinks :
cufflinks -g $gtf -b $genomeFile -u --library-type fr-firststrand aligned.bam
After that, I want to merge the assembly using cuffmerge.
Anyone has advices to better tune STAR and cufflinks paramteres for transcriptome assembly ? or maybe other idea/tools to perform that on RNA-Seq data
Thanks
N.