Skip to content
Snippets Groups Projects
Commit 09678a53 authored by Helene Rimbert's avatar Helene Rimbert
Browse files

add log file to gmap step

parent 09c1a530
No related branches found
No related tags found
No related merge requests found
...@@ -48,10 +48,10 @@ do ...@@ -48,10 +48,10 @@ do
targetindexdir=$gmapdbdir targetindexdir=$gmapdbdir
targetgff=$rootdir'/gmap.wholeGenome.gff3' targetgff=$rootdir'/gmap.wholeGenome.gff3'
gmapexe='gmapl' gmapexe='gmapl'
eval $gmapexe --min-identity 0.90 --min-trimmed-coverage 0.80 --ordered -n 1 -d $targetindex -D $targetindexdir -f 2 $rootdir'/query.mrna.fasta' 1> $targetgff eval $gmapexe --min-identity 0.90 --min-trimmed-coverage 0.80 --ordered -n 1 -d $targetindex -D $targetindexdir -f 2 $rootdir'/query.mrna.fasta' 1> $targetgff 2> $rootdir'/gmap.wg.log'
else else
echo "running gmap on target" echo "running gmap on target"
eval gmap --ordered -n 1 -g $rootdir'/target.fasta' -f 2 $rootdir'/query.mrna.fasta' 1> $targetgff eval gmap --ordered -n 1 -g $rootdir'/target.fasta' -f 2 $rootdir'/query.mrna.fasta' 1> $targetgff 2> $rootdir'/gmap.target.log'
fi fi
eval gmap --ordered -n 1 -g $rootdir'/query.fasta' -f 2 $rootdir'/query.mrna.fasta' 1> $rootdir'/gmap.query.gff3' eval gmap --ordered -n 1 -g $rootdir'/query.fasta' -f 2 $rootdir'/query.mrna.fasta' 1> $rootdir'/gmap.query.gff3'
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment