#!/bin/bash
printf "Enter ID : "; read id
printf "What panel: "; read panel
cd 'C:\Users\cmccabe\Desktop\annovar'
[ -z "$id" ] && break
[ "$id" = "end" ] && break
OMR=Output_Mutation_Report
perl -aF/\\t/ -lne 'BEGIN{%m=map{chomp;s/\cM|\cJ//g;$p=join("\t",(split/\t/)[4,5]);($p,$_)} <>;$m{"#CHROM\tINFO"}=$m{"Chr\tSegment Position"}};/SEGPOS=(\d+)/ || /\t(INFO)\t/ or next;$p=$F[0]."\t".$1;exists $m{$p} and print join("\t",$_,$m{$p})' ${id}_${panel}_${OMR}.txt < ${id}_${panel}_${OMR}_Filtered.vcf > ${id}_matched.vcf
I am trying to combine shells: bash ~/newbatch.sh
Code:
#!/bin/bash
printf "What is the id of the patient to be matched : "; read id
printf "What panel: "; read panel
case "$id" in
cd 'C:\Users\cmccabe\Desktop\annovar'
OMR=Output_Mutation_Report
$( perl -aF/\\t/ -lne 'BEGIN{%m=map{chomp;s/\cM|\cJ//g;$p=join("\t",(split/\t/)[4,5]);($p,$_)} <>;$m{"#CHROM\tINFO"}=$m{"Chr\tSegment Position"}};/SEGPOS=(\d+)/ || /\t(INFO)\t/ or next;$p=$F[0]."\t".$1;exists $m{$p} and print join("\t",$_,$m{$p})' ${id}_${panel}_${OMR}.txt < ${id}_${panel}_${OMR}_Filtered.vcf > ${id}_matched.vcf )
;;
show_menu=true
menu=( [yY] [nN] )
func_show() { # "MESSAGE STRING"
# Prints passed argument
# Returns nothing
printf '\n\t%s\n\n' "$1"
}
#
# Display & Action
#
clear
while $show_menu
do printf '#----------------------#\n%s\n' "Are there additonal patients to be matched"
select entry in "${menu[@]}" Back
do case "$entry" in
Back) show_menu=false
;;
[yY]) $( perl -aF/\\t/ -lne 'BEGIN{%m=map{chomp;s/\cM|\cJ//g;$p=join("\t",(split/\t/)[4,5]);($p,$_)} <>;$m{"#CHROM\tINFO"}=$m{"Chr\tSegment Position"}};/SEGPOS=(\d+)/ || /\t(INFO)\t/ or next;$p=$F[0]."\t".$1;exists $m{$p} and print join("\t",$_,$m{$p})' ${id}_${panel}_${OMR}.txt < ${id}_${panel}_${OMR}_Filtered.vcf > ${id}_matched.vcf ) (index id 0) :: $entry
;;
[nN]) func_show "printf "Does the file need to be converted : " ; read id :: ${menu[${#menu[@]}-1]}"
while true
do
printf "Enter ID : " ; read id
cd 'C:\Users\cmccabe\Desktop\annovar'
[ -z "$id" ] && break
[ "$id" = "end" ] && break
$( perl convert2annovar.pl -includeinfo -format vcf4old ${id}_matched.vcf > ${id}_matched.avinput )
done
;;
*)
esac
break
done
done
I am sure the code needs work, but why does it throw an error in the underlined part when it was fine before? Thanks .
Code:
bash ~/newbatch.sh
What is the id of the patient to be matched : H62947
What panel: Epilepsy70
/home/cmccabe/newbatch.sh: line 5: syntax error near unexpected token `'C:\Users\cmccabe\Desktop\annovar''
/home/cmccabe/newbatch.sh: line 5: ` cd 'C:\Users\cmccabe\Desktop\annovar''
Hi all,
I'm trying to wirte a small shell script in Linux. My script has the flow like,
cmd1
cmd2
cd testdata
cmd3
After exiting the program, the CWD remains the same as where I execute the program. I need it to be changed to the latest updated directory in the program. How can I do... (1 Reply)
hi,
Iam in directory A. I run a script from there. inside the script i have a command cd B.
When i come out of the script directory is A only.
Even when i come out scrip i want the directory to be B
How to achieve (2 Replies)
I would like to have a script that would change my current working directory. However, any time I execute a 'cd' command in a script, it holds only for the life of that script -- the working directory on exit is the same as when the script was initiated. Is it possible to have the script return... (3 Replies)
I have a directory that is existing under my root dir of the FTP server. The DIR name is 'Software Patch'. I want to move in to that DIR to download some patches. But, when I issued a command 'cd SOftware Patch', the system said that it cannot find the dir 'Software'. I tried all possible ways like... (2 Replies)
I have a directory of files and each file has a random 5 digit string at the beginning that needs to be removed. Plus, there are some files that will be identically named after the 5 digit string is removed and I want those eliminated or moved.
any ideas? (17 Replies)
I have a simple shell script that prompts the user to enter a directory to navigate to.
What i want it to do and i don't know how to do this is if the directory is invalid automatically navigate to the home directory.
echo "enter a directory to navigate to:"
read directory
cd $directory... (6 Replies)
Hi,
I Installed mysql on my CentOS 6.2 Server. But when I tried to change the location of /var/lib/mysql to another directory. I can't start the mysql. Below is what I've done
yum install mysql mysql-server mysql-devel
mkdir /path/to/new/
cp -R /var/lib/mysql /path/to/new
chown -R... (1 Reply)
I am trying to do the following task :
export ENV=aaa
export ENV_PATH=$(cd /apps | ls | grep $ENV)
However, it's not working. What's the way to change to directory and search some file in that directory in single command
Please help. (2 Replies)
Hi All,
There is a code like below in my script
###############################################
###Create Directories and Sub-Directories
###############################################
dpdir=DP_FROM_${from}_TO_${to}
mkdir $dpdir
cd $dpdir
mkdir AWQM WFCONTROLLER PROVCO PRISM
... (1 Reply)