Change some blastx parametter to get better matches by taking into

account intron size and the good genetic code

Former-commit-id: 6600123fbdce2070058074e82c791c7fc260c39b
Former-commit-id: ac413cc4a49844d4fa4087107aa84680d36f3df1
This commit is contained in:
2016-10-06 12:36:43 -03:00
parent e4f3081fa8
commit 4992483b80

View File

@ -12,6 +12,9 @@ unsetenv ORG_SOURCED
setenv ORG_HOME `dirname $0`/../../.. setenv ORG_HOME `dirname $0`/../../..
source $ORG_HOME/scripts/csh_init.sh source $ORG_HOME/scripts/csh_init.sh
set ParamsDir = $PROG_DIR/../params
set ModelsDir = $PROG_DIR/../models
NeedArg 2 NeedArg 2
set GenoFile = $Argv[1]; Shift set GenoFile = $Argv[1]; Shift
@ -19,6 +22,14 @@ set ProtFile = $Argv[1]; Shift
NeedFile $GenoFile NeedFile $GenoFile
NeedFile $ProtFile NeedFile $ProtFile
NeedFile $ParamsDir/default
#
# general parameters
#
source $ParamsDir/default
set IDMIN = 70 set IDMIN = 70
set NBMIN = 50 set NBMIN = 50
@ -43,7 +54,9 @@ endif
# #
Notify " blasting $ProtFile" Notify " blasting $ProtFile"
blastx -query $GenoFile -db $ProtFile -outfmt 7 |\ blastx -query $GenoFile -db $ProtFile -outfmt 7 \
-query_gencode $PASS1_BLASTX_FILTER_GCODE \
-max_intron_length $PASS1_MAX_INTRON |\
$AwkCmd -v IDMIN=$IDMIN \ $AwkCmd -v IDMIN=$IDMIN \
-v NBMIN=$NBMIN \ -v NBMIN=$NBMIN \
-v NBMAX=$NBMAX \ -v NBMAX=$NBMAX \