'\" t
.\" Title: gt-csa
.\" Author: [FIXME: author] [see http://docbook.sf.net/el/author]
.\" Generator: DocBook XSL Stylesheets v1.78.1
.\" Date: 07/29/2016
.\" Manual: GenomeTools Manual
.\" Source: GenomeTools 1.5.9
.\" Language: English
.\"
.TH "GT\-CSA" "1" "07/29/2016" "GenomeTools 1\&.5\&.9" "GenomeTools Manual"
.\" -----------------------------------------------------------------
.\" * Define some portability stuff
.\" -----------------------------------------------------------------
.\" ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.\" http://bugs.debian.org/507673
.\" http://lists.gnu.org/archive/html/groff/2009-02/msg00013.html
.\" ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
.ie \n(.g .ds Aq \(aq
.el .ds Aq '
.\" -----------------------------------------------------------------
.\" * set default formatting
.\" -----------------------------------------------------------------
.\" disable hyphenation
.nh
.\" disable justification (adjust text to left margin only)
.ad l
.\" -----------------------------------------------------------------
.\" * MAIN CONTENT STARTS HERE *
.\" -----------------------------------------------------------------
.SH "NAME"
gt-csa \- Transform spliced alignments from GFF3 file into consensus spliced alignments\&.
.SH "SYNOPSIS"
.sp
\fBgt csa\fR [option \&...] [GFF3_file]
.SH "DESCRIPTION"
.PP
\fB\-join\-length\fR [\fIvalue\fR]
.RS 4
set join length for the spliced alignment clustering (default: 300)
.RE
.PP
\fB\-v\fR [\fIyes|no\fR]
.RS 4
be verbose (default: no)
.RE
.PP
\fB\-o\fR [\fIfilename\fR]
.RS 4
redirect output to specified file (default: undefined)
.RE
.PP
\fB\-gzip\fR [\fIyes|no\fR]
.RS 4
write gzip compressed output file (default: no)
.RE
.PP
\fB\-bzip2\fR [\fIyes|no\fR]
.RS 4
write bzip2 compressed output file (default: no)
.RE
.PP
\fB\-force\fR [\fIyes|no\fR]
.RS 4
force writing to output file (default: no)
.RE
.PP
\fB\-help\fR
.RS 4
display help and exit
.RE
.PP
\fB\-version\fR
.RS 4
display version information and exit
.RE
.SH "EXAMPLE:"
.sp
Let\(cqs assume we have a GFF3 file \fIcsa_example_spliced_alignments\&.gff3\fR containing the following four overlapping spliced alignments (represented as genes with exons as children):
.sp
.if n \{\
.RS 4
.\}
.nf
##gff\-version 3
##sequence\-region seq 1 290
seq \&. gene 1 209 \&. + \&. ID=gene1
seq \&. exon 1 90 \&. + \&. Parent=gene1
seq \&. exon 110 190 \&. + \&. Parent=gene1
seq \&. exon 201 209 \&. + \&. Parent=gene1
###
seq \&. gene 1 290 \&. + \&. ID=gene2
seq \&. exon 1 90 \&. + \&. Parent=gene2
seq \&. exon 101 190 \&. + \&. Parent=gene2
seq \&. exon 201 290 \&. + \&. Parent=gene2
###
seq \&. gene 10 290 \&. + \&. ID=gene3
seq \&. exon 10 90 \&. + \&. Parent=gene3
seq \&. exon 110 190 \&. + \&. Parent=gene3
seq \&. exon 201 290 \&. + \&. Parent=gene3
###
seq \&. gene 181 290 \&. + \&. ID=gene4
seq \&. exon 181 190 \&. + \&. Parent=gene4
seq \&. exon 201 290 \&. + \&. Parent=gene4
###
.fi
.if n \{\
.RE
.\}
.sp
To compute the consensus spliced alignments we call:
.sp
.if n \{\
.RS 4
.\}
.nf
$ gt csa csa_example_spliced_alignments\&.gff3
.fi
.if n \{\
.RE
.\}
.sp
Which returns:
.sp
.if n \{\
.RS 4
.\}
.nf
##gff\-version 3
##sequence\-region seq 1 290
seq gt csa gene 1 290 \&. + \&. ID=gene1
seq gt csa mRNA 1 290 \&. + \&. ID=mRNA1;Parent=gene1
seq gt csa exon 1 90 \&. + \&. Parent=mRNA1
seq gt csa exon 110 190 \&. + \&. Parent=mRNA1
seq gt csa exon 201 290 \&. + \&. Parent=mRNA1
seq gt csa mRNA 1 290 \&. + \&. ID=mRNA2;Parent=gene1
seq gt csa exon 1 90 \&. + \&. Parent=mRNA2
seq gt csa exon 101 190 \&. + \&. Parent=mRNA2
seq gt csa exon 201 290 \&. + \&. Parent=mRNA2
###
.fi
.if n \{\
.RE
.\}
.sp
As one can see, they have been combined into a consensus spliced alignment (represented as genes with mRNAs as children which in turn have exons as children) with two alternative splice forms\&. The first and the third spliced alignment have been combined into the first alternative splice form (mRNA1) and the the second and the fourth spliced alignment into the second alternative splice form (mRNA2)\&.
.sp
As one can see, the second exon from the first alternative splice form is shorter than the corresponding exon from the second alternative splice form\&.
.SH "REPORTING BUGS"
.sp
Report bugs to \&.