Scroll to navigation

SEGEMEHL(1) User Commands SEGEMEHL(1)

NAME

segemehl - Heuristic mapping of short sequences

SYNOPSIS

segemehl [-besVOc] -d <file> [<file>] [-q <file>] [-p <file>] [-i <file>] [-j <file>] [-x <file>] [-y <file>] [-G <file>] [-g <string>] [-t <n>] [-o <string>] [-u <file>] [-B <string>] [-F <n>] [-S [<basename>]] [-A <n>] [-D <n>] [-E <double>] [-H] [-m <n>] [-Z <n>] [-W <n>] [-U <n>] [-l <f>] [-w <double>] [-X <n>] [-J <n>] [-I <n>] [-M <n>] [-n <n>] [-r <n>] [--skipidcheck] [--showalign] [--nohead]

DESCRIPTION

Segemehl is a software to map short sequencer reads to reference genomes. Segemehl implements a matching strategy based on enhanced suffix arrays (ESA). Segemehl accepts fasta and fastq queries (gzip’ed and bgzip'ed). In addition to the alignment of reads from standard DNA- and RNA-seq protocols, it also allows the mapping of bisulfite converted reads (Lister and Cokus) and implements a split read mapping strategy. The output of segemehl is a SAM or BAM formatted alignment file. In the case of split-read mapping, additional BED files are written to the disc. These BED files may be summarized with the postprocessing tool haarz. In the case of the alignment of bisulfite converted reads, raw methylation rates may also be called with haarz.

In brief, for each suffix of a read, segemehl aims to find the best-scoring seed. Seeds might contain insertions, deletions, and mismatches (differences). The number of differences allowed within a single seed is user-controlled and is crucial for the runtime of the program. Subsequently, seeds that undercut the user-defined E-value are passed on to an exact semi-global alignment procedure. Finally, reads with a minimum accuracy of percent are reported to the user.

OPTIONS

INPUT

list of path/filename(s) of fasta database sequence(s)
path/filename of query sequences (default:none)
path/filename of mate pair sequences (default:none)
path/filename of db index (default:none)
path/filename of second db index (default:none)
generate db index and store to disk (default:none)
generate second db index and store to disk (default:none)
filename to read @RG header (default:none)
read group id (default:none)
start <n> threads (default:1)

OUTPUT

outputfile (default:none)
generate a bam output (-o <filename> required)
filename for unmatched reads (default:none)
brief cigar string (M vs X and =)
show a progress bar
file bins with basename <string> for easier data handling (default:none)
output MEOP field for easier variance calling in SAM (XE:Z:)

ALIGNMENT

bisulfite aln with methylC-seq/Lister et al. (=1) or bs-seq/Cokus et al. protocol (=2) (default:0)
detect split/spliced reads. (default:none)
min percentage of matches per read in semi-global alignment (default:90)
search seeds initially with <n> differences (default:1)
max evalue (default:5.000000)
report only best scoring hits (=1) or all (=0) (default:1)
minimum length of queries (default:12)
min length of a spliced fragment (default:20)
min coverage for spliced transcripts (default:80)
min score of a spliced fragment (default:18)
report spliced alignment with score s only if <f>*s is larger than next best spliced alignment (default:0.900000)
max evalue for splits (default:50.000000)

SPECIAL

dropoff parameter for extension (default:8)
search seeds with jump size <n> (0=automatic) (default:0)
sorts the output by chromsome and position (might take a while!)
maximum size of the inserts (paired end) in case of multiple hits (default:200000)
maximum width of a suffix array interval, i.e. a query seed will be omitted if it matches more than <n> times (default:100)
check index
penalty for a mismatch during extension (default:4)
maximum number of alignments that will be reported. If set to zero, all alignments will be reported (default:0)
do not check whether the fastq ids of mates / paired ends match. Instead, the first mate (-q) will be used for output only.
show alignments
do not output header

BUGS

Please report bugs to steve@bioinf.uni-leipzig.de

SEE ALSO

http://www.bioinf.uni-leipzig.de/Software/segemehl/

REFERENCES

2008 Bioinformatik Leipzig
2018 Leibniz Institute on Aging (FLI)

AUTHOR

This manpage was written by Andreas Tille for the Debian distribution and can be used for any other usage of the program.

October 2018 segemehl 0.3