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

keep only one hit with gmap (-n1)

parent 7ded6fa6
No related branches found
No related tags found
No related merge requests found
......@@ -61,9 +61,9 @@ do
echo "3. run gmap on target and query"
if [ ! -s "$rootdir/gmap.target.gff3" ]
then
gmap -d target -D $rootdir -f 2 $rootdir'/query.mrna.fasta' > $rootdir'/gmap.target.gff3' &
gmap -n 1 -d target -D $rootdir -f 2 $rootdir'/query.mrna.fasta' > $rootdir'/gmap.target.gff3' &
pid=$!
gmap -d query -D $rootdir -f 2 $rootdir'/query.mrna.fasta' > $rootdir'/gmap.query.gff3' &
gmap -n 1 -d query -D $rootdir -f 2 $rootdir'/query.mrna.fasta' > $rootdir'/gmap.query.gff3' &
pid2=$!
blastn -query $rootdir'/query.fasta' -subject $rootdir'/target.fasta' -outfmt 6 -out $rootdir'/blast_m8.tab' &
pid3=$!
......
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