.\" Automatically generated by Pod::Man 2.22 (Pod::Simple 3.07) .\" .\" Standard preamble: .\" ======================================================================== .de Sp \" Vertical space (when we can't use .PP) .if t .sp .5v .if n .sp .. .de Vb \" Begin verbatim text .ft CW .nf .ne \\$1 .. .de Ve \" End verbatim text .ft R .fi .. .\" Set up some character translations and predefined strings. \*(-- will .\" give an unbreakable dash, \*(PI will give pi, \*(L" will give a left .\" double quote, and \*(R" will give a right double quote. \*(C+ will .\" give a nicer C++. Capital omega is used to do unbreakable dashes and .\" therefore won't be available. \*(C` and \*(C' expand to `' in nroff, .\" nothing in troff, for use with C<>. .tr \(*W- .ds C+ C\v'-.1v'\h'-1p'\s-2+\h'-1p'+\s0\v'.1v'\h'-1p' .ie n \{\ . ds -- \(*W- . ds PI pi . if (\n(.H=4u)&(1m=24u) .ds -- \(*W\h'-12u'\(*W\h'-12u'-\" diablo 10 pitch . if (\n(.H=4u)&(1m=20u) .ds -- \(*W\h'-12u'\(*W\h'-8u'-\" diablo 12 pitch . ds L" "" . ds R" "" . ds C` "" . ds C' "" 'br\} .el\{\ . ds -- \|\(em\| . ds PI \(*p . ds L" `` . ds R" '' 'br\} .\" .\" Escape single quotes in literal strings from groff's Unicode transform. .ie \n(.g .ds Aq \(aq .el .ds Aq ' .\" .\" If the F register is turned on, we'll generate index entries on stderr for .\" titles (.TH), headers (.SH), subsections (.SS), items (.Ip), and index .\" entries marked with X<> in POD. Of course, you'll have to process the .\" output yourself in some meaningful fashion. .ie \nF \{\ . de IX . tm Index:\\$1\t\\n%\t"\\$2" .. . nr % 0 . rr F .\} .el \{\ . de IX .. .\} .\" .\" Accent mark definitions (@(#)ms.acc 1.5 88/02/08 SMI; from UCB 4.2). .\" Fear. Run. Save yourself. No user-serviceable parts. . \" fudge factors for nroff and troff .if n \{\ . ds #H 0 . ds #V .8m . ds #F .3m . ds #[ \f1 . ds #] \fP .\} .if t \{\ . ds #H ((1u-(\\\\n(.fu%2u))*.13m) . ds #V .6m . ds #F 0 . ds #[ \& . ds #] \& .\} . \" simple accents for nroff and troff .if n \{\ . ds ' \& . ds ` \& . ds ^ \& . ds , \& . ds ~ ~ . ds / .\} .if t \{\ . ds ' \\k:\h'-(\\n(.wu*8/10-\*(#H)'\'\h"|\\n:u" . ds ` \\k:\h'-(\\n(.wu*8/10-\*(#H)'\`\h'|\\n:u' . ds ^ \\k:\h'-(\\n(.wu*10/11-\*(#H)'^\h'|\\n:u' . ds , \\k:\h'-(\\n(.wu*8/10)',\h'|\\n:u' . ds ~ \\k:\h'-(\\n(.wu-\*(#H-.1m)'~\h'|\\n:u' . ds / \\k:\h'-(\\n(.wu*8/10-\*(#H)'\z\(sl\h'|\\n:u' .\} . \" troff and (daisy-wheel) nroff accents .ds : \\k:\h'-(\\n(.wu*8/10-\*(#H+.1m+\*(#F)'\v'-\*(#V'\z.\h'.2m+\*(#F'.\h'|\\n:u'\v'\*(#V' .ds 8 \h'\*(#H'\(*b\h'-\*(#H' .ds o \\k:\h'-(\\n(.wu+\w'\(de'u-\*(#H)/2u'\v'-.3n'\*(#[\z\(de\v'.3n'\h'|\\n:u'\*(#] .ds d- \h'\*(#H'\(pd\h'-\w'~'u'\v'-.25m'\f2\(hy\fP\v'.25m'\h'-\*(#H' .ds D- D\\k:\h'-\w'D'u'\v'-.11m'\z\(hy\v'.11m'\h'|\\n:u' .ds th \*(#[\v'.3m'\s+1I\s-1\v'-.3m'\h'-(\w'I'u*2/3)'\s-1o\s+1\*(#] .ds Th \*(#[\s+2I\s-2\h'-\w'I'u*3/5'\v'-.3m'o\v'.3m'\*(#] .ds ae a\h'-(\w'a'u*4/10)'e .ds Ae A\h'-(\w'A'u*4/10)'E . \" corrections for vroff .if v .ds ~ \\k:\h'-(\\n(.wu*9/10-\*(#H)'\s-2\u~\d\s+2\h'|\\n:u' .if v .ds ^ \\k:\h'-(\\n(.wu*10/11-\*(#H)'\v'-.4m'^\v'.4m'\h'|\\n:u' . \" for low resolution devices (crt and lpr) .if \n(.H>23 .if \n(.V>19 \ \{\ . ds : e . ds 8 ss . ds o a . ds d- d\h'-1'\(ga . ds D- D\h'-1'\(hy . ds th \o'bp' . ds Th \o'LP' . ds ae ae . ds Ae AE .\} .rm #[ #] #H #V #F C .\" ======================================================================== .\" .IX Title "2NDSCORE 1" .TH 2NDSCORE 1 "2011-02-19" "" "User Contributed Documentation" .\" For nroff, turn off justification. Always turn off hyphenation; it makes .\" way too many mistakes in technical documents. .if n .ad l .nh .SH "NAME" .Vb 1 \& 2ndscore \- find the best hairpin anchored at each position. .Ve .SH "SYNOPSIS" .IX Header "SYNOPSIS" .Vb 1 \& 2ndscore in.fasta > out.hairpins .Ve .SH "DESCRIPTION" .IX Header "DESCRIPTION" For every position in the sequence this will output a line: .PP .Vb 2 \& \-0.6 52 .. 62 TTCCTAAAGGTTCCA GCG CAAAA TGC CATAAGCACCACATT \& (score) (start .. end) (left context) (hairpin) (right contenxt) .Ve .PP For positions near the ends of the sequences, the context may be padded with \&'x' characters. If no hairpin can be found, the score will be 'None'. .PP Multiple fasta files can be given and multiple sequences can be in each fasta file. The output for each sequence will be separated by a line starting with \&'>' and containing the \s-1FASTA\s0 description of the sequence. .PP Because the hairpin scores of the plus-strand and minus-strand may differ (due to \s-1GU\s0 binding in \s-1RNA\s0), by default 2ndscore outputs two sets of hairpins for every sequence: the \s-1FORWARD\s0 hairpins and the \s-1REVERSE\s0 hairpins. All the forward hairpins are output first, and are identified by having the word '\s-1FORWARD\s0' at the end of the '>' line preceding them. Similarly, the \s-1REVERSE\s0 hairpins are listed after a '>' line ending with '\s-1REVERSE\s0'. If you want to search only one or the other strand, you can use: .PP .Vb 2 \& \-\-no\-fwd Don\*(Aqt print the FORWARD hairpins \& \-\-no\-rvs Don\*(Aqt print the REVERSE hairpins .Ve .PP You can set the energy function used, just as with transterm with the \-\-gc, \&\-\-au, \-\-gu, \-\-mm, \-\-gap options. The \-\-min\-loop, \-\-max\-loop, and \-\-max\-len options are also supported. .SS "\s-1FORMAT\s0 \s-1OF\s0 \s-1THE\s0 .BAG \s-1FILES\s0" .IX Subsection "FORMAT OF THE .BAG FILES" The columns for the .bag files are, in order: .PP .Vb 6 \& 1. gene_name \& 2. terminator_start \& 3. terminator_end \& 4. hairpin_score \& 5. tail_score \& 6. terminator_sequence \& \& 7. terminator_confidence: a combination of the hairpin and tail score that \& takes into account how likely such scores are in a random sequence. This \& is the main "score" for the terminator and is computed as described in \& the paper. \& \& 8. APPROXIMATE_distance_from_end_of_gene: The *approximate* number of base \& pairs between the end of the gene and the start of the terminator. This \& is approximate in several ways: First, (and most important) TransTermHP \& doesn\*(Aqt always use the real gene ends. Depending on the options you give \& it may trim some off the ends of genes to handle terminators that \& partially overlap with genes. Second, where the terminator "begins" \& isn\*(Aqt that well defined. This field is intended only for a sanity check \& (terminators reported to be the best near the ends of genes shouldn\*(Aqt be \& _too far_ from the end of the gene). .Ve .SS "\s-1USING\s0 \s-1TRANSTERM\s0 \s-1WITHOUT\s0 \s-1GENOME\s0 \s-1ANNOTATIONS\s0" .IX Subsection "USING TRANSTERM WITHOUT GENOME ANNOTATIONS" TransTermHP uses known gene information for only 3 things: (1) tagging the putative terminators as either \*(L"inside genes\*(R" or \*(L"intergenic,\*(R" (2) choosing the background GC-content percentage to compute the scores, because genes often have different \s-1GC\s0 content than the intergenic regions, and (3) producing slightly more readable output. Items (1) and (3) are not really necessary, and (2) has no effect if your genes have about the same GC-content as your intergenic regions. .PP Unfortunately, TransTermHP doesn't yet have a simple option to run without an annotation file (either .ptt or .coords), and requires at least 2 genes to be present. The solution is to create fake, small genes that flank each chromosome. To do this, make a fake.coords file that contains only these two lines: .PP .Vb 2 \& fakegene1 1 2 chome_id \& fakegene2 L\-1 L chrom_id .Ve .PP where L is the length of the input sequence and L\-1 is 1 less than the length of the input sequence. \*(L"chrom_id\*(R" should be the word directly following the \*(L">\*(R" in the .fasta file containing your sequence. (If, for example, your .fasta file began with \*(L">seq1\*(R", then chrom_id = seq1). .PP This creates a \*(L"fake\*(R" annotation with two 1\-base\-long genes flanking the sequence in a tail-to-tail arrangement: \-\-> <\-\-. TransTermHP can then be run with: .PP .Vb 1 \& transterm \-p expterm.dat sequence.fasta fake.coords .Ve .PP If the G/C content of your intergenic regions is about the same as your genes, then this won't have too much of an effect on the scores terminators receive. On the other hand, this use of TransTermHP hasn't been tested much at all, so it's hard to vouch for its accuracy. .SH "SEE ALSO" .IX Header "SEE ALSO" \&\fItransterm\fR\|(1) .SH AUTHOR Alex Mestiashvili