academic/spidey: Added x86_64 source.

Signed-off-by: Willy Sudiarto Raharjo <willysr@slackbuilds.org>

5 files changed, 349 insertions, 25 deletions

diff --git a/academic/spidey/README b/academic/spidey/README
index 0912181c25..6e7c198210 100644
--- a/academic/spidey/README
+++ b/academic/spidey/README
@@ -1,9 +1,15 @@
-Spidey is an mRNA-to-genomic alignment program. For a complete 
-description of how Spidey works, visit 
-http://www.ncbi.nlm.nih.gov/spidey/spideydoc.html.
+Spidey: an mRNA-to-genomic alignment program.
 
-This is just repackaging of the ready binary for x86 and will not run 
-on x86_64. It will probably work just fine on a Slackware multilib 
-box, but we do not support that ;). If you want to build spidey from 
-source, you should download and compile the NCBI toolkit. For more 
-information: http://www.ncbi.nlm.nih.gov/spidey/spideysource.html
+Spidey is a tool for aligning one or more mRNA  sequences to a given
+genomic sequence. It was written with two main goals in mind:
+1) find good alignments regardless of intron size
+2) avoid getting confused by nearby pseudogenes and paralogs.
+
+The following programs provide a GUI to run spidey:
+-ugene
+-perlprimer
+
+This is just repackaging of precompiled binaries:
+- x86 platform: the executable is provided by upstream (NCBI).
+- x86_64 platform: the executable is kindly provided by the UniPro
+Ugene project, where it is part of their External Tools meta-package.
diff --git a/academic/spidey/slack-desc b/academic/spidey/slack-desc
index cd570e4133..3689502f74 100644
--- a/academic/spidey/slack-desc
+++ b/academic/spidey/slack-desc
@@ -8,12 +8,12 @@
       |-----handy-ruler------------------------------------------------------|
 spidey: spidey (mRNA-to-genomic alignment)
 spidey:
-spidey: Spidey is an mRNA-to-genomic alignment program.
-spidey:
+spidey: Spidey is a tool for aligning one or  more  mRNA  sequences
+spidey: to a given genomic sequence.
 spidey:
+spidey: Home: http://www.ncbi.nlm.nih.gov/spidey/index.html
 spidey:
 spidey:
 spidey:
 spidey:
 spidey:
-spidey: Home: http://www.ncbi.nlm.nih.gov/spidey/index.html
diff --git a/academic/spidey/spidey.1 b/academic/spidey/spidey.1
new file mode 100644
index 0000000000..f3e2836f95
--- /dev/null
+++ b/academic/spidey/spidey.1
@@ -0,0 +1,279 @@
+.TH SPIDEY 1 2005-01-25 NCBI "NCBI Tools User's Manual"
+.SH NAME
+spidey \- align mRNA sequences to a genome
+.SH SYNOPSIS
+.B spidey
+[\|\fB\-\fP\|]
+[\|\fB\-F\fP\ \fIN\fP\|]
+[\|\fB\-G\fP\|]
+[\|\fB\-L\fP\ \fIN\fP\|]
+[\|\fB\-M\fP\ \fIfilename\fP\|]
+[\|\fB\-N\fP\ \fIfilename\fP\|]
+[\|\fB\-R\fP\ \fIfilename\fP\|]
+[\|\fB\-S\fP\ \fIp/m\fP\|]
+[\|\fB\-T\fP\ \fIN\fP\|]
+[\|\fB\-X\fP\|]
+[\|\fB\-a\fP\ \fIfilename\fP\|]
+[\|\fB\-c\fP\ \fIN\fP\|]
+[\|\fB\-d\fP\|]
+[\|\fB\-e\fP\ \fIX\fP\|]
+[\|\fB\-f\fP\ \fIX\fP\|]
+[\|\fB\-g\fP\ \fIX\fP\|]
+\fB\-i\fP\ \fIfilename\fP
+[\|\fB\-j\fP\|]
+[\|\fB\-k\fP\ \fIfilename\fP\|]
+[\|\fB\-l\fP\ \fIN\fP\|]
+\fB\-m\fP\ \fIfilename\fP
+[\|\fB\-n\fP\ \fIN\fP\|]
+[\|\fB\-o\fP\ \fIstr\fP\|]
+[\|\fB\-p\fP\ \fIN\fP\|]
+[\|\fB\-r\fP\ \fIc/d/m/p/v\fP\|]
+[\|\fB\-s\fP\|]
+[\|\fB\-t\fP\ \fIfilename\fP\|]
+[\|\fB\-u\fP\|]
+[\|\fB\-w\fP\|]
+.SH DESCRIPTION
+\fBspidey\fP is a tool for aligning one or more mRNA sequences to a
+given genomic sequence.  \fBspidey\fP was written with two main goals
+in mind: find good alignments regardless of intron size; and avoid
+getting confused by nearby pseudogenes and paralogs.  Towards the
+first goal, \fBspidey\fP uses BLAST and Dot View (another local
+alignment tool) to find its alignments; since these are both local
+alignment tools, \fBspidey\fP does not intrinsically favor shorter or
+longer introns and has no maximum intron size.  To avoid mistakenly
+including exons from paralogs and pseudogenes, \fBspidey\fP first
+defines windows on the genomic sequence and then performs the
+mRNA-to-genomic alignment separately within each window.  Because of
+the way the windows are constructed, neighboring paralogs or
+pseudogenes should be in separate windows and should not be included
+in the final spliced alignment.
+.SS Initial alignments and construction of genomic windows
+\fBspidey\fP takes as input a single genomic sequence and a set of
+mRNA accessions or FASTA sequences.  All processing is done one mRNA
+sequence at a time.  The first step for each mRNA sequence is a
+high-stringency BLAST against the genomic sequence.  The resulting
+hits are analyzed to find the genomic windows.
+.PP
+The BLAST alignments are sorted by score and then assigned into
+windows by a recursive function which takes the first alignment and
+then goes down the alignment list to find all alignments that are
+consistent with the first (same strand of mRNA, both the mRNA and
+genomic coordinates are nonoverlapping and linearly consistent).  On
+subsequent passes, the remaining alignments are examined and are put
+into their own nonoverlapping, consistent windows, until no alignments
+are left.  Depending on how many gene models are desired, the
+top \fIn\fP windows are chosen to go on to the next step and the others
+are deleted.
+.SS Aligning in each window
+Once the genomic windows are constructed, the initial BLAST alignments
+are freed and another BLAST search is performed, this time with the
+entire mRNA against the genomic region defined by the window, and at a
+lower stringency than the initial search.  \fBspidey\fP then uses a
+greedy algorithm to generate a high-scoring, nonoverlapping subset of
+the alignments from the second BLAST search.  This consistent set is
+analyzed carefully to make sure that the entire mRNA sequence is
+covered by the alignments.  When gaps are found between the
+alignments, the appropriate region of genomic sequence is searched
+against the missing mRNA, first using a very low-stringency BLAST and,
+if the BLAST fails to find a hit, using DotView functions to locate
+the alignment.  When gaps are found at the ends of the alignments, the
+BLAST and DotView searches are actually allowed to extend past the
+boundaries of the window.  If the 3' end of the mRNA does not align
+completely, it is first examined for the presence of a poly(A) tail.
+No attempt is made to align the portion of the mRNA that seems to be a
+poly(A) tail; sometimes there is a poly(A) tail that does align to the
+genomic sequence, and these are noted because they indicate the
+possibility of a pseudogene.
+.PP
+Now that the mRNA is completely covered by the set of alignments, the
+boundaries of the alignments (there should be one alignment per exon
+now) are adjusted so that the alignments abut each other precisely and
+so that they are adjacent to good splice donor and acceptor sites.
+Most commonly, two adjacent exons' alignments overlap by as much as 20
+or 30 base pairs on the mRNA sequence.  The true exon boundary may lie
+anywhere within this overlap, or (as we have seen empirically) even a
+few base pairs outside the overlap.  To position the exon boundaries,
+the overlap plus a few base pairs on each side is examined for splice
+donor sites, using functions that have different splice matrices
+depending on the organism chosen.  The top few splice donor sites (by
+score) are then evaluated as to how much they affect the original
+alignment boundaries.  The site that affects the boundaries the least
+is chosen, and is evaluated as to the presence of an acceptor site.
+The alignments are truncated or extended as necessary so that they
+terminate at the splice donor site and so that they do not overlap.
+.SS Final result
+The windows are examined carefully to get the percent identity per
+exon, the number of gaps per exon, the overall percent identity, the
+percent coverage of the mRNA, presence of an aligning or non-aligning
+poly(A) tail, number of splice donor sites and the presence or absence
+of splice donor and acceptor sites for each exon, and the occurrence
+of an mRNA that has a 5' or 3' end (or both) that does not align to
+the genomic sequence.  If the overall percent identity and percent
+length coverage are above the user-defined cutoffs, a summary report
+is printed, and, if requested, a text alignment showing identities and
+mismatches is also printed.
+.SS Interspecies alignments
+\fBspidey\fP is capable of performing interspecies alignments.  The
+major difference in interspecies alignments is that the mRNA-genomic
+identity will not be close to 100% as it is in intraspecies
+alignments; also, the alignments have numerous and lengthy gaps.  If
+\fBspidey\fP is used in its normal mode to do interspecies alignments,
+it produces gene models with many, many short exons.  When the
+interspecies flag is set, \fBspidey\fP uses different BLAST parameters
+to encourage longer and more gaps and to not penalize as heavily for
+mismatches.  This way, the alignments for the exons are much longer
+and more closely approximate the actual gene structure.
+.SS Extracting CDS alignments
+When \fBspidey\fP is run in network-aware mode or when ASN.1 files are
+used for the mRNA records, it is capable of extracting a CDS alignment
+from an mRNA alignment and printing the CDS information also.  Since
+the CDS alignment is just a subset of the mRNA alignment, it is
+relatively straightforward to truncate the exon alignments as
+necessary and to generate a CDS alignment.  Furthermore, the
+untranslated regions are now defined, so the percent identity for the
+5' and 3' untranslated regions is also calculated.
+.PP
+.SH OPTIONS
+A summary of options is included below.
+.TP
+\fB\-\fP
+Print usage message.
+.TP
+\fB\-F\fP\ \fIN\fP
+Start of genomic interval desired (from; 0-based).
+.TP
+\fB\-G\fP
+Input file is a GI list.
+.TP
+\fB\-L\fP\ \fIN\fP
+The extra-large intron size to use (default = 220000).
+.TP
+\fB\-M\fP\ \fIfilename\fP
+File with donor splice matrix.
+.TP
+\fB\-N\fP\ \fIfilename\fP
+File with acceptor splice matrix.
+.TP
+\fB\-R\fP\ \fIfilename\fP
+File (including path) to repeat blast database for filtering.
+.TP
+\fB\-S\fP\ \fIp/m\fP
+Restrict to plus (p) or minus (m) strand of genomic sequence.
+.TP
+\fB\-T\fP\ \fIN\fP
+Stop of genomic interval desired (to; 0-based).
+.TP
+\fB\-X\fP
+Use extra-large intron sizes (increases the limit for initial and
+terminal introns from 100kb to 240kb and for all others from 35kb to
+120kb); may result in significantly longer compute times.
+.TP
+\fB\-a\fP\ \fIfilename\fP
+Output file for alignments when directed to a separate file with
+\fB-p\ 3\fP (default = spidey.aln).
+.TP
+\fB\-c\fP\ \fIN\fP
+Identity cutoff, in percent, for quality control purposes.
+.TP
+\fB\-d\fP
+Also try to align coding sequences corresponding to the given mRNA
+records (may require network access).
+.TP
+\fB\-e\fP\ \fIX\fP
+First-pass e-value (default = 1.0e-10).  Higher values increase speed
+at the cost of sensitivity.
+.TP
+\fB\-f\fP\ \fIX\fP
+Second-pass e-value (default = 0.001).
+.TP
+\fB\-g\fP\ \fIX\fP
+Third-pass e-value (default = 10).
+.TP
+\fB\-i\fP\ \fIfilename\fP
+Input file containing the genomic sequence in ASN.1 or FASTA format.
+If your computer is running on a network that can access GenBank, you
+can substitute the desired accession number for the filename.
+.TP
+\fB\-j\fP
+Print ASN.1 alignment?
+.TP
+\fB\-k\fP\ \fIfilename\fP
+File for ASN.1 output with \fB-k\fP (default = spidey.asn).
+.TP
+\fB\-l\fP\ \fIN\fP
+Length coverage cutoff, in percent.
+.TP
+\fB\-m\fP\ \fIfilename\fP
+Input file containing the mRNA sequence(s) in ASN.1 or FASTA format,
+or a list of their accessions (with \fB-G\fP).  If your computer is
+running on a network that can access GenBank, you can substitute a
+single accession number for the filename.
+.TP
+\fB\-n\fP\ \fIN\fP
+Number of gene models to return per input mRNA (default = 1).
+.TP
+\fB\-o\fP\ \fIstr\fP
+Main output file (default = stdout; contents controlled by \fB-p\fP).
+.TP
+\fB\-p\fP\ \fIN\fP
+Print alignment?
+.RS
+.PD 0
+.IP \fB0\fP
+summary and alignments together (default)
+.IP \fB1\fP
+just the summary
+.IP \fB2\fP
+just the alignments
+.IP \fB3\fP
+summary and alignments in different files
+.PD
+.RE
+.TP
+\fB\-r\fP\ \fIc/d/m/p/v\fP
+Organism of genomic sequence, used to determine splice matrices.
+.RS
+.PD 0
+.IP \fBc\fP
+C. elegans
+.IP \fBd\fP
+Drosophila
+.IP \fBm\fP
+Dictyostelium discoideum
+.IP \fBp\fP
+plant
+.IP \fBv\fP
+vertebrate (default)
+.PD
+.RE
+.TP
+\fB\-s\fP
+Tune for interspecies alignments.
+.TP
+\fB\-t\fP\ \fIfilename\fP
+File with feature table, in 4 tab-delimited columns:
+.RS
+.PD 0
+.IP \fIseqid\fP
+(e.g., \fBNM_04377.1\fP)
+.IP \fIname\fP
+(only \fBrepetitive_region\fP is currently supported)
+.IP \fIstart\fP
+(0-based)
+.IP \fIstop\fP
+(0-based)
+.PD
+.RE
+.TP
+\fB\-u\fP
+Make a multiple alignment of all input mRNAs (which must overlap on
+the genomic sequence).
+.TP
+\fB\-w\fP
+Consider lowercase characters in input FASTA sequences to be masked.
+.SH AUTHOR
+Sarah Wheelan and others at the National Center for Biotechnology
+Information; Steffen Moeller contributed to this documentation.
+.SH SEE ALSO
+.BR blast (1),
+<http://www.ncbi.nlm.nih.gov/spidey>
diff --git a/academic/spidey/spidey.SlackBuild b/academic/spidey/spidey.SlackBuild
index 443940eb11..5a31c4f26a 100644
--- a/academic/spidey/spidey.SlackBuild
+++ b/academic/spidey/spidey.SlackBuild
@@ -1,25 +1,57 @@
 #!/bin/sh
 
-# Slackware build script for spidey
-# Written by Petar Petrov, <ppetrov@paju.oulu.fi> and
-# hereby submitted to the public domain
-
-# THIS SLACKBUILD IS DISTRIBUTETD IN THE HOPE OF BEING
-# USEFUL BUT WITHOUT ANY WARRANTY. THE AUTHOR IS _NOT_
-# RESPONSIBLE FOR ANY DAMAGE OR DATA LOSS CAUSED BY IT.
+# Copyright 2011-2015 Petar Petrov, petar.petrov@student.oulu.fi
+# All rights reserved.
+#
+# Redistribution and use of this script, with or without modification, is
+# permitted provided that the following conditions are met:
+#
+# 1. Redistributions of this script must retain the above copyright
+#    notice, this list of conditions and the following disclaimer.
+#
+#  THIS SOFTWARE IS PROVIDED BY THE AUTHOR "AS IS" AND ANY EXPRESS OR IMPLIED
+#  WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF
+#  MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO
+#  EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+#  SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
+#  PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS;
+#  OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
+#  WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR
+#  OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF
+#  ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 
 PRGNAM=spidey
-VERSION=${VERSION:-20060601}
-BUILD=${BUILD:-1}
+VERSION=${VERSION:-20060601} # Keep the date of the 32bit binary as version.
+BUILD=${BUILD:-2}
 TAG=${TAG:-_SBo}
 
-ARCH=i386
+if [ -z "$ARCH" ]; then
+  case "$( uname -m )" in
+    i?86) ARCH=i386 ;;
+    arm*) ARCH=arm ;;
+       *) ARCH=$( uname -m ) ;;
+  esac
+fi
 
 CWD=$(pwd)
 TMP=${TMP:-/tmp/SBo}
 PKG=$TMP/package-$PRGNAM
 OUTPUT=${OUTPUT:-/tmp}
 
+
+if [ "$ARCH" != "i386" ] && [ "$ARCH" != "x86_64" ]; then
+  printf "\n\n$ARCH is not supported... \n"
+  exit 1
+fi
+
+# Determine the source arch. Many thanks to the Ugene project for the
+# 64bit executable!
+if [ "$ARCH" = "x86_64" ]; then
+  SRCARCH=".64"
+else
+  SRCARCH=""
+fi
+
 set -e
 
 rm -rf $PKG
@@ -28,13 +60,20 @@ cd $TMP
 rm -rf $PRGNAM-$VERSION
 mkdir $PRGNAM-$VERSION
 cd $PRGNAM-$VERSION
-gunzip -c $CWD/$PRGNAM.linux.gz > spidey
+
+gunzip -c $CWD/$PRGNAM.linux${SRCARCH}.gz > spidey
 
 install -D -m755 spidey $PKG/usr/bin/spidey
 
+mkdir -p $PKG/usr/man/man1
+cp $CWD/$PRGNAM.1 $PKG/usr/man/man1/$PRGNAM.1
+
 find $PKG -print0 | xargs -0 file | grep -e "executable" -e "shared object" | grep ELF \
   | cut -f 1 -d : | xargs strip --strip-unneeded 2> /dev/null || true
 
+find $PKG/usr/man -type f -exec gzip -9 {} \;
+for i in $( find $PKG/usr/man -type l ) ; do ln -s $( readlink $i ).gz $i.gz ; rm $i ; done
+
 mkdir -p $PKG/usr/doc/$PRGNAM-$VERSION
 cat $CWD/$PRGNAM.SlackBuild > $PKG/usr/doc/$PRGNAM-$VERSION/$PRGNAM.SlackBuild
 
diff --git a/academic/spidey/spidey.info b/academic/spidey/spidey.info
index 324a7da4a9..3a21a46c7a 100644
--- a/academic/spidey/spidey.info
+++ b/academic/spidey/spidey.info
@@ -3,8 +3,8 @@ VERSION="20060601"
 HOMEPAGE="http://www.ncbi.nlm.nih.gov/spidey/index.html"
 DOWNLOAD="ftp://ftp.ncbi.nih.gov/pub/wheelan/Spidey/spidey.linux.gz"
 MD5SUM="2e56ef2e4fcf57eca266fb1b3bb56c7e"
-DOWNLOAD_x86_64="UNSUPPORTED"
-MD5SUM_x86_64="UNSUPPORTED"
+DOWNLOAD_x86_64="http://www.student.oulu.fi/~ppetrov/source/spidey.linux.64.gz"
+MD5SUM_x86_64="79f1f95976346e0d0f5c7f717deac176"
 REQUIRES=""
 MAINTAINER="Petar Petrov"
-EMAIL="ppetrov@paju.oulu.fi"
+EMAIL="petar.petrov@student.oulu.fi"