<<

NAME

frag_lens - Calculate paired end fragment lengths from read alignment

SYNOPSIS

frag_lens [--sam <in.sam>] OR [--ref <cons.fa> --forward <R1.fq> --reverse <R2.fq>] [options] > frag_lens.txt

DESCRIPTION

Calculates library fragment lengths based on paired-end read alignment. Takes as input either a preprepared SAM alignment or a reference and read files from which it produces an alignment. Outputs calculated fragment lengths, one per line.

PREREQUISITES

Requires the following binaries:

OPTIONS

Input option one

--sam filename

Path to input SAM alignment.

Input option two

--ref filename

Path to reference sequence (e.g. assembly)

--forward filename

Forward reads in FASTQ format

--reverse filename

Reverse reads in FASTQ format

Configuration

--max_align

Maximum number of alignment records to read as input. Used to limit run times.

--threads

Number of threads to use for alignment (ignored if --sam is given)

CAVEATS AND BUGS

Please submit bug reports to the issue tracker in the distribution repository.

AUTHOR

Jeremy Volkening (jeremy.volkening@base2bio.com)

LICENSE AND COPYRIGHT

Copyright 2014-23 Jeremy Volkening

This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program. If not, see <http://www.gnu.org/licenses/>.

<<