|
8
|
1 # automate rnastar index generation
|
|
9
|
2 # from existing fasta genomes using limited ram - cut genomeChrbinNbits down to 13, genomeSAindexNbases 12 genomeSAsparseD 1
|
|
|
3
|
|
8
|
4 INDEXROOT=/mnt/galaxyIndices/genomes
|
|
|
5 # fix that, the list of genomes below and SOURCE below to suit your needs. This is for GVL august 2013
|
|
|
6 for genome in Arabidopsis_thaliana_TAIR9 ce10 danRer7 dm3 hg19 hg_g1k_v37 mm9 phiX rn5 sacCer2 sacCer3 xenTro2
|
|
4
|
7 do
|
|
8
|
8 SOURCE=${INDEXROOT}/${genome}/seq/
|
|
4
|
9 TARGET=${INDEXROOT}/${genome}/rnastar
|
|
|
10 mkdir -p $TARGET
|
|
|
11 cd $TARGET
|
|
9
|
12 /mnt/galaxy/tools/rnastar/2.3.0e/fubar/rgrnastar_203e/9425706f4e97/bin/STAR --runMode genomeGenerate --genomeDir $TARGET --genomeFastaFiles ${SOURCE}/${genome}.fa --runThreadN 2 genomeLoad=NoSharedMemory --genomeChrBinNbits 13 genomeSAindexNbases 12 genomeSAsparseD 1
|
|
4
|
13 done
|
|
|
14
|
|
9
|
15 echo "#add these lines TABS not spaces please to rnastar_indexes.loc in tool-data" > rnastar_indexes.loc
|
|
|
16 for genome in Arabidopsis_thaliana_TAIR9 ce10 danRer7 dm3 hg19 hg_g1k_v37 mm9 phiX rn5 sacCer2 sacCer3 xenTro2
|
|
|
17 do
|
|
|
18 SOURCE=${INDEXROOT}/${genome}/seq/
|
|
|
19 TARGET=${INDEXROOT}/${genome}/rnastar
|
|
|
20 echo ${genome}'\t'${genome}'\t'${genome}'\t'/mnt/galaxyIndices/genomes/${genome}/rnastar >> rnastar_indexes.loc
|
|
|
21 done
|
