annotate commandline_sample_STR-FM_reference_profiling @ 12:a3113043abb0 draft

Uploaded
author arkarachai-fungtammasan
date Sun, 24 Jul 2016 17:55:50 -0400
parents d5ed5c2e25c3
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
2
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
1 ## This is a sample PBS script for profiling STR from reference genome using STR-FM version 1.0.0 (April 20, 2014)
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
2 ##
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
3 ##requirement
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
4 ##1 reference genome in FASTA format --> ${INPUT}.fa
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
5 ##
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
6 echo " "
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
7 echo " "
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
8 echo "Job started on `hostname` at `date`"
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
9 cd /working/directory/
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
10 echo " "
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
11 echo " detect STR in reference genome" ## See detail in microsatellite.xml on https://github.com/Arkarachai/STR-FM
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
12 python microsatellite.py ${INPUT}.fa --fasta --period=1 --partialmotifs --minlength=4 --prefix=0 --suffix=0 --hamming=0 --multipleruns --flankdisplay=0 --splitbyvalidity >${INPUT}.mono.out
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
13 python microsatellite.py ${INPUT}.fa --fasta --period=2 --partialmotifs --minlength=6 --prefix=0 --suffix=0 --hamming=0 --multipleruns --flankdisplay=0 --splitbyvalidity >${INPUT}.di.out
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
14 python microsatellite.py ${INPUT}.fa --fasta --period=3 --partialmotifs --minlength=6 --prefix=0 --suffix=0 --hamming=0 --multipleruns --flankdisplay=0 --splitbyvalidity >${INPUT}.tri.out
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
15 python microsatellite.py ${INPUT}.fa --fasta --period=4 --partialmotifs --minlength=8 --prefix=0 --suffix=0 --hamming=0 --multipleruns --flankdisplay=0 --splitbyvalidity >${INPUT}.tetra.out
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
16
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
17 echo "formatting"
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
18 cat ${INPUT}.mono.out | awk 'BEGIN{FS="\t";OFS="\t"};{print $6,$2,$2+$1,$4,$1,length($4) }' > ${INPUT}.mono.TR
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
19 cat ${INPUT}.di.out | awk 'BEGIN{FS="\t";OFS="\t"};{print $6,$2,$2+$1,$4,$1,length($4) }' > ${INPUT}.di.TR
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
20 cat ${INPUT}.tri.out | awk 'BEGIN{FS="\t";OFS="\t"};{print $6,$2,$2+$1,$4,$1,length($4) }' > ${INPUT}.tri.TR
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
21 cat ${INPUT}.tetra.out | awk 'BEGIN{FS="\t";OFS="\t"};{print $6,$2,$2+$1,$4,$1,length($4) }' > ${INPUT}.tetra.TR
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
22
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
23
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
24
d5ed5c2e25c3 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
25 echo "Job end on `hostname` at `date`"