wd="/Volumes/web/scaphapoda/Grace/Transcriptomes/mercenaria"
dircode="me"

cd {wd}

!blastx \
-query query.fa \
-db /Volumes/Data/blast_db/uniprot_sprot \
-max_target_seqs 1 \
-max_hsps 1 \
-outfmt 6 \
-num_threads 8 \
-out blast_sprot.tab

!wc -l blast_sprot.tab

!tr '|' "\t" <blast_sprot.tab> blast_sprot.sql.tab

!python /Applications/sqlshare-pythonclient-master/tools/singleupload.py \
-d {dircode}_uniprot \
blast_sprot_sql.tab

!python /Applications/sqlshare-pythonclient-master/tools/singleupload.py \
-s "SELECT Column1, term, GOSlim_bin, aspect, ProteinName FROM [graceac9@washington.edu].[md_uniprot]md left join [samwhite@washington.edu].[UniprotProtNamesReviewed_yes20130610]sp on md.Column3=sp.SPID left join [sr320@washington.edu].[SPID and GO Numbers]go on md.Column3=go.SPID left join [sr320@washington.edu].[Go_to_Goslim]slim on go.GOID=GO_id where aspect like 'P'" \
-f tsv \
-o {dircode}_descriptions.txt

!head {dircode}_descriptions.txt

!egrep --color "male|female|genitalia|gonad|ovarian|reproduction|estrogen|testosterone|gametogenesis|germination|ovulation|penile|prostate|vulval" {dircode}_descriptions.txt / {dircode}_reprot.txt

!head {dircode}_reprot.txt

!wc -l {dircode}_repro.txt