view CSP2/CSP2_env/env-d9b9114564458d9d-741b3de822f2aaca6c6caa4325c4afce/opt/mummer-3.23/run-mummer1 @ 69:33d812a61356

planemo upload commit 2e9511a184a1ca667c7be0c6321a36dc4e3d116d
author jpayne
date Tue, 18 Mar 2025 17:55:14 -0400
parents
children
line wrap: on
line source
#! -f
#
# **SEVERELY** antiquated script for running the mummer 1 suite
# -r option reverse complements the query sequence, coordinates of the reverse
# matches will be relative to the reversed sequence
#

set ref = $1
set qry = $2
set pfx = $3
set rev = $4

set bindir = /mnt/c/Users/crash/Documents/BobLiterman/CSP2_Galaxy/CSP2/CSP2_env/env-d9b9114564458d9d-741b3de822f2aaca6c6caa4325c4afce/opt/mummer-3.23

if($ref == '' || $qry == '' || $pfx == '') then
    echo "USAGE: $0 <fasta reference> <fasta query> <prefix> [-r]"
    exit(-1)
endif

echo "Find MUMs"
$bindir/mummer -mum -l 20 $rev $ref $qry | tail +2 > $pfx.out
echo "Determine gaps"
$bindir/gaps $ref $rev < $pfx.out > $pfx.gaps
echo "Align gaps"
$bindir/annotate $pfx.gaps $qry > $pfx.align
mv witherrors.gaps $pfx.errorsgaps