PhyloSkims/annotate
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

correction of the IR detection

Browse Source
This commit is contained in:
Eric Coissac
2025-05-25 19:38:01 +02:00
parent 478a6bdca7
commit 28df0c35c1
2 changed files with 2 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
detectors/cds/bin/go_cds.sh
View File

@ -111,7 +111,7 @@ if [[ "$cdsdetection_pass3" == "yes" ]] ; then
loginfo $fams loginfo $fams
loginfo "running pass3:$dir exonerate of $Genome on $DbRoot" loginfo "running pass3:$dir exonerate of $Genome on $DbRoot"
for f in $fams ; do for f in $fams ; do
tcsh -f $PROG_DIR/do_exonerate.csh Pass3 $Fasta $f $AnnotFile $DbRoot/models $temp echo tcsh -f $PROG_DIR/do_exonerate.csh Pass3 $Fasta $f $AnnotFile $DbRoot/models $temp
done done
fi fi
done) | parallel -j $Threads done) | parallel -j $Threads

3
detectors/normalize/lib/selectIR.py
View File

@ -120,8 +120,7 @@ for line in repeats:
# pvalue= # pvalue=
# print("c.lsc = %f c.ssc = %f o.lsc = %f o.ssc = %f score = %6.4f (len=%d)" % (c_lsc,c_ssc,o_lsc,o_ssc,score,len1), # print("c.lsc = %f c.ssc = %f o.lsc = %f o.ssc = %f score = %6.4f (len=%d)" % (c_lsc,c_ssc,o_lsc,o_ssc,score,len1),
# file=sys.stderr) # file=sys.stderr)
if (int(score*100) >= int(scoreMax*100)) and ((len1 > len1Max) or (len2 > len2Max)):
if (score >= scoreMax) and ((len1 > len1Max) or (len2 > len2Max)):
scoreMax = score scoreMax = score
pos1Max = pos1 pos1Max = pos1
pos2Max = pos2 pos2Max = pos2