comparison rgrnastar/prepStar.sh @ 9:c3ac5a2a2bfb draft

Uploaded
author fubar
date Sat, 10 Aug 2013 08:58:12 -0400
parents 6b232288f519
children 1e40b977aa58
comparison
equal deleted inserted replaced
8:6b232288f519 9:c3ac5a2a2bfb
1 # automate rnastar index generation 1 # automate rnastar index generation
2 # from existing fasta genomes 2 # from existing fasta genomes using limited ram - cut genomeChrbinNbits down to 13, genomeSAindexNbases 12 genomeSAsparseD 1
3
3 INDEXROOT=/mnt/galaxyIndices/genomes 4 INDEXROOT=/mnt/galaxyIndices/genomes
4 # fix that, the list of genomes below and SOURCE below to suit your needs. This is for GVL august 2013 5 # fix that, the list of genomes below and SOURCE below to suit your needs. This is for GVL august 2013
5 for genome in Arabidopsis_thaliana_TAIR9 ce10 danRer7 dm3 hg19 hg_g1k_v37 mm9 phiX rn5 sacCer2 sacCer3 xenTro2 6 for genome in Arabidopsis_thaliana_TAIR9 ce10 danRer7 dm3 hg19 hg_g1k_v37 mm9 phiX rn5 sacCer2 sacCer3 xenTro2
6 do 7 do
7 SOURCE=${INDEXROOT}/${genome}/seq/ 8 SOURCE=${INDEXROOT}/${genome}/seq/
8 TARGET=${INDEXROOT}/${genome}/rnastar 9 TARGET=${INDEXROOT}/${genome}/rnastar
9 mkdir -p $TARGET 10 mkdir -p $TARGET
10 cd $TARGET 11 cd $TARGET
11 /mnt/galaxy/tools/rnastar/2.3.0e/fubar/rgrnastar_203e/33b5c8dc8e60/bin/STAR --runMode genomeGenerate --genomeDir $TARGET --genomeFastaFiles ${SOURCE}/seq/${genome}.fa --runThreadN 2 genomeLoad=NoSharedMemory 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
12 --genomeChrBinNbits 12
13 done 13 done
14 14
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