# HG changeset patch # User davidvanzessen # Date 1481896900 18000 # Node ID 54f6756bacb1b0bcb558e86f3b71668838cffe4e # Parent d001d0c05dbe8d918de49e3435d82aa3f50e3329 Uploaded diff -r d001d0c05dbe -r 54f6756bacb1 complete.sh --- a/complete.sh Fri Dec 16 07:30:32 2016 -0500 +++ b/complete.sh Fri Dec 16 09:01:40 2016 -0500 @@ -46,7 +46,7 @@ bash ${dir}/imgt_loader/imgt_loader.sh $current $parsedFileName "${fileName}" else echo "Sample $count of patient $id is not a zip file so assuming fasta/fastq, using igBLASTn" >> $html - bash ${dir}/igblast/igblast.sh $current $species $locus $parsedFileName + bash ${dir}/igblast/igblast.sh $current "$species" $locus $parsedFileName fi mergerInput+=($parsedFileName) count=$((count+1)) diff -r d001d0c05dbe -r 54f6756bacb1 igblast/igblast.sh --- a/igblast/igblast.sh Fri Dec 16 07:30:32 2016 -0500 +++ b/igblast/igblast.sh Fri Dec 16 09:01:40 2016 -0500 @@ -9,10 +9,16 @@ declare -A speciesdict -speciesdict=(["Bos taurus functional"]="BosTaurus" ["Bos taurus non-functional"]="BosTaurus" ["Camelus dromedarius functional"]="CamelusDromedarius" ["Camelus dromedarius non-functional"]="CamelusDromedarius" ["Canis lupus familiaris functional"]="CanisLupusFamiliaris" ["Canis lupus familiaris non-functional"]="CanisLupusFamiliaris" ["Danio rerio non-functional"]="DanioRerio" ["Danio rerio functional"]="DanioRerio" ["Homo sapiens functional"]="HomoSapiens" ["Homo sapiens non-functional"]="HomoSapiens" ["Macaca mulatta non-functional"]="MacacaMulatta" ["Macaca mulatta functional"]="MacacaMulatta" ["Mus musculus functional"]="MusMusculus" ["Mus musculus non-functional"]="MusMusculus" ["Mus spretus functional"]="MusSpretus" ["Mus spretus non-functional"]="MusSpretus" ["Oncorhynchus mykiss functional"]="OncorhynchusMykiss" ["Oncorhynchus mykiss non-functional"]="OncorhynchusMykiss" ["Oryctolagus cuniculus functional"]="OryctolagusCuniculus" ["Oryctolagus cuniculus non-functional"]="OryctolagusCuniculus" ["Rattus norvegicus functional"]="RattusNorvegicus" ["Rattus norvegicus non-functional"]="RattusNorvegicus" ["Sus scrofa functional"]="SusScrofa" ["Sus scrofa non-functional"]="SusScrofa" ["BosTaurus"]="BosTaurus" ["CamelusDromedarius"]="CamelusDromedarius" ["CanisLupusFamiliaris"]="CanisLupusFamiliaris" ["DanioRerio"]="DanioRerio" ["HomoSapiens"]="HomoSapiens" ["MacacaMulatta"]="MacacaMulatta" ["MusMusculus"]="MusMusculus" ["MusSpretus"]="MusSpretus" ["OncorhynchusMykiss"]="OncorhynchusMykiss" ["OryctolagusCuniculus"]="OryctolagusCuniculus" ["RattusNorvegicus"]="RattusNorvegicus" ["SusScrofa"]="SusScrofa") +speciesdict=(["Rattus norvegicus functional"]="rat" ["Rattus norvegicus non-functional"]="rat" ["Oryctolagus cuniculus functional"]="rabbit" ["Oryctolagus cuniculus non-functional"]="rabbit" ["Mus musculus functional"]="mouse" ["Mus musculus non-functional"]="mouse" ["Homo sapiens functional"]="human" ["Homo sapiens non-functional"]="human" ["Macaca mulatta non-functional"]="rhesus_monkey" ["Macaca mulatta functional"]="rhesus_monkey") + +echo "Species: $species ${speciesdict[$species]}" + +species="${speciesdict[$species]}" echo "$input $species $locus $output" java -Xmx64G -jar $IGBLASTWRP/igblastwrp.jar -p 4 -S $species -R $locus ${input} $PWD/blasted_output 2>&1 +#java -Xmx64G -jar $dir/igblastwrp.jar -p 4 -S $species -R $locus ${input} $PWD/blasted_output 2>&1 +#java -Xmx64G -jar $dir/migmap-1.0.1.jar -p 4 -S "$species" -R $locus ${input} $PWD/blasted_output 2>&1 Rscript --verbose $dir/igblast.r "$PWD/blasted_output.L2.txt" "$output" 2>&1