# HG changeset patch # User boris # Date 1370455516 14400 # Node ID 663fc7a8549b0e91115f4811727300c77cab08df # Parent e0638858594fbdb099fbe697cf77224d55a2246e uploaded tool xml diff -r e0638858594f -r 663fc7a8549b phylorelatives.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/phylorelatives.xml Wed Jun 05 14:05:16 2013 -0400 @@ -0,0 +1,113 @@ + + Relatedness of minor alelle sequences in NJ tree + + R + rpy2 + dendropy + ape + + phylorelatives.py + #for $sequence in $sequences: + -i "${sequence.sequence}" + #end for + -m "${multifasta}" + -b "$iterations" + $pairwise + #if str($root.root_selector)=="upload_root": + -r "${root.root_fa}" + #end if + --relatives-out "$siblings" + --newick-out "$newick" + --trees-out "$plot" + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +(boris-at-bx.psu.edu) + + +----- + +.. class:: infomark + +**What it does** + +Constructs relatedness of a set of sequences based on the pairwise proportion of different sites. +One or more test sequences are accepted as long as their name include the strict suffix "_minor" or "_test" (i.e. >seq1_minor). +It returns the FASTA multiple alignment used, and reports a table with the relatives of the minor allele sequence in a NJ tree, the tree plot and newick string. + +----- + +.. class:: warningmark + +**Note** + +This tools DOES NOT align the sequences. +The FASTA multiple alignment file is fabricated by concatenating same length homologous sequences. +An independently aligned FASTA multiple alignment file can be used as well. + +----- + +.. class:: infomark + +**About formats** + +**FASTA multiple alignment** -- + +describe format shortly + +----- + +**Example** + +- For the following dataset:: + >sample1_major + >sample1_minor + >sample2_major + >sample3_major + >sample4_major + + +- running this tool with **b = 1000** **m = yes** **root = RSRS**, will return:: + newick + relatives + plot + + + \ No newline at end of file