Scroll to navigation

spaln(1) General Commands Manual spaln(1)

NAME

spaln - mapping and alignment of cDNA/protein sequences onto genomic sequence or rapid homology search against protein sequence database and (semi)global alignment

SYNOPSIS

spaln -WGenome.bk[n|p] -K[D|P] [W_Options] Genome.mfa(.gz)
spaln -WProsDb.bka -KA [W_Options] ProteinSeqenceDB.faa(.gz)
spaln [-Q[0|1|2|3]] [R_options] Genome_segment [cDNA|protein]_queries
spaln -Q[4|5|6|7] -dGenome [R_options] [cDNA|protein]_queries
spaln -Q[4|5|6|7] -aProsDb [R_options] Genomic_segment
spaln -Q[4|5|6|7] -aProsDb [R_options] protein_queries
spaln -Q[4|5|6|7] [R_options] Genome.mfa [cDNA|protein]_queries
spaln -Q[4|5|6|7] [R_options] ProsDb.faa protein_queries

DESCRIPTION

spaln is a stand-alone program that maps and aligns a set of cDNA/protein sequences onto a genomic sequence. From version 1.4, spaln supports a combination of protein sequence database and a given genomic segment. From version 2.2, it also supports rapid similarity search of protein sequences against a protein sequence database followd by ordinary alignment.

Spaln runs with relatively small internal memory, making it feasible to search a whole mammalian genome in a single job on a conventional computer.

EXAMPLES

Make the block index table 'ddignm.bkn' of the genomic sequence 'ddignm.mfa' in multi-fasta format with the word size of 10 nucleotides and the block size of 8192.
Align the cDNA sequence 'cdna.fa' onto the genomic segment of 'chr1.fa' within the range from the 10001-th to the 40000-th nucleotide of the complementary strand in the semi-global mode.
Report the gene structures corresponding to individual cDNAs in 'ddicdna.mfa'. The local alignment mode is used and the results are written in 'ddi.exon' in the exon-oriented format. The computation proceeds with 10 threads in parallel.
Cross-species comparison between the genome 'mmugnm' and the cDNAs in 'ratcdna.mfa' from the first to the 20-th entries. The maximal expected gene size is reset to 1 Mbp. The outputs go to 'mmu.intron' in the intron-oriented format. The tetrapod-specific parameter set is used.
A set of ORFs in the genomic segment are translated and rapidly searched against the protein database, and then the best-hit sequence is used as the template of the spliced alignment to infer the organization of the gene located in the relevant genomic region. The output is shown in Gff3 gene format. The dicot-specific parameter set is used.
Find up to four SwissProt sequences most similar to each query sequence, and report global alignment statistics. To show alignment themselves, use -O1 option in stead of -O0. The computation proceeds with 10 threads in parallel.

OPTIONS

Conventions: #: number; $: string; default values in ()
(default for DNA, with -yX option, default for protein)

FORMATING OPTIONS

Format the genomic sequence for DNA ($=D) or protein ($=P) queries, or format the protein database sequences ($=A) for rapid search of template.

Word size (11 for DNA, 5 for protein)
Block size (4096)
Maximum gene size (262144)
Abundance factor (10)
Shifts between adjacent seeds (k)

RUN TIME OPTIONS

NCBI transl_table number indicating the genetic code (1)
Minimum alignment score for report (35)
Smith-Waterman-type local alignment
Multiple loci (0)
#=empty: Multiple loci maximally up to 4
#=0: Single locus
#=1: Re-search unaligned parts
#>1: Multiple loci maximally up to #
Output format (4)
#=0: GFF3 format (gene) in genome vs [cDNA|protein] mode
#=0: Alignment statistics in protein vs protein mode
#=1: Alignment
#=2: GFF3 format (match)
#=3: Bed format
#=4: Exon-oriented format similar to output of megablast -D 3
#=5: Intron-oriented output
#=6: Concatenated exon sequence
#=7: Translated amino-acid sequence
#=8: Mapping (block) information only. Use with -Q4
#=12: Output the same information as -O4 in binary formats
Select algorithm (3)
0<=#<=3: Genomic segment in the fasta format given by the first argument vs. cdna/protein given by the second argument
4<=#<=7: Genome mapping and alignment
#=0,4: DP procedure without HSP search
#=1-3,5-7: Recursive HSP searches up to the level of (# % 4)
Read block index table from the file $
Specify the orientation of query (0)
0: depend on the query annotation
1: forward direction only
2: reverse direction only
3: both directions;
Specify the species (in combination with the -yS option for cDNA query)
Map/align without splicing
Minimum space to induce the Hirschberg's algorithm (16M)
Write block index table to file $
Paired-end reads as the queries from a single or two files
-ia: 5' and 3' matching pairs must appear alternatively in a file
-ip: 5' and 3' matching pairs must appear in the same order in the two files
Destination of output file (stdout)
Suppress trimming of terminal polyA or polyT sequence
Suppress some outputs to stderr
Report the result irrespective of the alignment score
Gap-extension penalty (3, 2, 2)
Gap-opening penalty (8, 6, 9)
Bit pattern of the seeds used for HSP search at level 1
Bit pattern of the seeds used for HSP search at level 3
Dinucleotide pairs at the ends of an intron (0)
0: canonical only (GT..AG, GC..AG, AT..AC)
1: relaxed to (GT..AG, GC..AG, AT..AN)
2: #=1 + allow 1 mismatch from GT..AG
3: any;
Intron penalty (11, 8, 11)
Incline of long gap penalty (0.6)
Flex point where the incline of gap penalty changes (7)
Double affine gap penalty if #=3; affine penalty otherwise
Score for a nucleotide match (2, 2)
Penalty for a nucleotide mismatch (6, 2)
Penalty for in-frame termination codon (100)
PAM level used in the alignment (third) phase (150)
PAM level used in the second phase (50)
Penalty for a frame shift (100)
Relative contribution of splicing signal (8)
Relative contribution of coding potential (2)
Relative contribution of the translational initiation or termination signal (8)
Relative contribution of branch point signal (0)
Minimum exon length (2)
Intron distribution parameters
Relative contribution of the bonus given to a conserved intron position
The minimum intron length (20)
Percent contribution of the species-specific splice signal. if #=0, only the ubiquitous signal given to the dinucleotide pair at the ends of an intron is used. By default #=0 for DNA and #=100 for protein queries. -yX option automatically sets #=100 for DNA and #=30 for protein. -yS For cDNA queries, use species-specific exon-intron boundary signals. For protein queries, invoke 'salvage' procedure in phase 1
For a DNA query, this option sets parameter values for cross-species comparison. Conversely, this option specifies an intra-species mode for a protein query.
Relative contribution of length-dependent part of intron penalty (8)
Relative contribution of oligomer composition within an intron (0)

REFERENCES

(1) "A Space-Efficient and Accurate Method for Mapping and Aligning cDNA Sequences onto Genomic Sequence", Osamu Gotoh, Nucleic Acid Res., 36 (8), 2630-2638 (2008).
(2) "Direct Mapping and Alignment of Protein Sequences onto Genomic Sequence", Osamu Gotoh, Bioinformatics, 24 (21) 2438-2444 (2008).
(3) "Benchmarking spliced alignment programs including Spaln2, an extended version of Spaln that incorporates additional species-specific features", Iwata, H. and Gotoh, O.", Nucleic Acids Res., 40 (20) e161 (2012).

AUTHOR

Osamu Gotoh <o.gotoh@aist.go.jp>

2018-09-06