# HG changeset patch # User fubar # Date 1708762807 0 # Node ID dc53eb4354a60ad01ed09cd9e1bcfa4aa855a532 # Parent 9ba0184870ef9b8114347c51fd63205546feb5cf planemo upload diff -r 9ba0184870ef -r dc53eb4354a6 mashmap.bash.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/mashmap.bash.txt Sat Feb 24 08:20:07 2024 +0000 @@ -0,0 +1,21 @@ +ln -s '$query' 'query' && +#if len($reflist) > 1: +#for i, mash in enumerate($reflist): + #if i == 0: + echo '$mash' > 'reflist' && + #else: + echo '$mash' >> 'reflist' && + #end if +#end for +#end if +samtools faidx 'query' && +mashmap --pi '$perc_identity' -s '$seqLength' -f '$filtermode' $dense \ +#if int($sketchSize) > 0: +-J '$sketchSize' \ +#end if +#if len($reflist) == 1: + -r '$reflist' -q 'query' && +#else +--rl 'reflist' -q 'query' && +#end if +cp 'mashmap.out' '$mashout' diff -r 9ba0184870ef -r dc53eb4354a6 mashmap.xml --- a/mashmap.xml Sat Feb 24 06:58:16 2024 +0000 +++ b/mashmap.xml Sat Feb 24 08:20:07 2024 +0000 @@ -1,6 +1,6 @@ - + Fast local alignment boundaries mashmap @@ -9,7 +9,8 @@ - 1: + 1: #for i, mash in enumerate($reflist): #if i == 0: echo '$mash' > 'reflist' && @@ -18,15 +19,15 @@ #end if #end for #end if -samtools faidx '$query' && +samtools faidx 'query' && mashmap --pi '$perc_identity' -s '$seqLength' -f '$filtermode' $dense \ #if int($sketchSize) > 0: -J '$sketchSize' \ #end if #if len($reflist) == 1: - -r '$reflist' -q '$query' && + -r '$reflist' -q 'query' && #else ---rl 'reflist' -q '$query' && +--rl 'reflist' -q 'query' && #end if cp 'mashmap.out' '$mashout']]>