# HG changeset patch # User fubar # Date 1708749156 0 # Node ID 53f601fb8664dc445115ca450e5605e17c6ccbc3 # Parent 27df186d544659da1c65b49793b88eabe1ccd972 planemo upload diff -r 27df186d5446 -r 53f601fb8664 mashmap.xml --- a/mashmap.xml Sat Feb 24 04:11:27 2024 +0000 +++ b/mashmap.xml Sat Feb 24 04:32:36 2024 +0000 @@ -1,11 +1,12 @@ - + - + Fast local alignment boundaries mashmap + samtools - + 1: @@ -17,6 +18,7 @@ #end if #end for #end if +samtools faidx '$query' && mashmap --pi '$perc_identity' -s '$seqLength' -f '$filtermode' $dense \ #if int($sketchSize) > 0: -J '$sketchSize' \ @@ -70,8 +72,7 @@ mashmap -r reference.fna -q query.fa - The output is a paf format file (https://github.com/lh3/miniasm/blob/master/PAF.md). - Thi is space-delimited with each line consisting of query name, length, 0-based start, end, strand, target name, length, start, end and mapping nucleotide identity. + The output is space-delimited with each line consisting of query name, length, 0-based start, end, strand, target name, length, start, end and mapping nucleotide identity. Map set of query seqences against a list of reference genomes: