Liston:Computer Scripts - Revision history

Zachary S. L. Foster: /* Python Scripts */

Zachary S. L. Foster — Thu, 23 Sep 2010 00:07:33 GMT

Python Scripts

←Older revision		Revision as of 00:07, 23 September 2010
Line 32:		Line 32:
	<tr><th>Script Name</th><th>Description</th><th>Input File Format</th><th>Output File Format</th></tr>		<tr><th>Script Name</th><th>Description</th><th>Input File Format</th><th>Output File Format</th></tr>
	<tr><td>[[alignreads.py]] v2.22</td>		<tr><td>[[alignreads.py]] v2.22</td>
-	<td>A pipeline for combining the free aligner NUCmer and the free short read assembler YASRA so full alignments can be made with one command. Alignreads uses the following scripts from the Liston lab: runyasra.py, sumqual.py, and qualtofa.py. It also uses yasra and its binaries, including lastz. It must be possible for the python interpreter to find these scripts to run alignreads; this can be done saving copies of scripts in the Python26 folder, your bin folder, or modifying your .chsrc file to search where they are saved on your system (nessesary if lastz binaries cannot be found). alignreads preserves almost all of the functionality and output of the five programs it uses, so there are many options (over 40). Because there are so many modifiers, the help menu reached by not providing arguments only has what are expected to be the options used most often; to get all of the options type 'python alignreads.py -H' or 'python alignreads.py --advanced-help'. Since there are many files associated with the alignment, the output is a folder with numerous files organized into subfolders. alignreads.py can also start the pipeline from a folder produced by a previous run of alignreads.py or runyasra.py. In that case, the output is saved in the input folder and separate subfolders are created for each alignment. A more complete discription of alignreads can be found [[~~HERE~~]]. </td>	+	<td>A pipeline for combining the free aligner NUCmer and the free short read assembler YASRA so full alignments can be made with one command. Alignreads uses the following scripts from the Liston lab: runyasra.py, sumqual.py, and qualtofa.py. It also uses yasra and its binaries, including lastz. It must be possible for the python interpreter to find these scripts to run alignreads; this can be done saving copies of scripts in the Python26 folder, your bin folder, or modifying your .chsrc file to search where they are saved on your system (nessesary if lastz binaries cannot be found). alignreads preserves almost all of the functionality and output of the five programs it uses, so there are many options (over 40). Because there are so many modifiers, the help menu reached by not providing arguments only has what are expected to be the options used most often; to get all of the options type 'python alignreads.py -H' or 'python alignreads.py --advanced-help'. Since there are many files associated with the alignment, the output is a folder with numerous files organized into subfolders. alignreads.py can also start the pipeline from a folder produced by a previous run of alignreads.py or runyasra.py. In that case, the output is saved in the input folder and separate subfolders are created for each alignment. A more complete discription of alignreads can be found here: [[alignreads.py README]]. </td>
	<td>A FASTA file containing reads, and A FASTA file containing a reference OR a folder from a previous run of alignreads or runyasra.</td>		<td>A FASTA file containing reads, and A FASTA file containing a reference OR a folder from a previous run of alignreads or runyasra.</td>
	<td>A folder containing all of the output from every part of the pipeline organized in subfolders based on each programs' output.</td></tr>		<td>A folder containing all of the output from every part of the pipeline organized in subfolders based on each programs' output.</td></tr>

Zachary S. L. Foster: /* Python Scripts */

Zachary S. L. Foster — Wed, 22 Sep 2010 23:57:46 GMT

Python Scripts

←Older revision		Revision as of 23:57, 22 September 2010
Line 32:		Line 32:
	<tr><th>Script Name</th><th>Description</th><th>Input File Format</th><th>Output File Format</th></tr>		<tr><th>Script Name</th><th>Description</th><th>Input File Format</th><th>Output File Format</th></tr>
	<tr><td>[[alignreads.py]] v2.22</td>		<tr><td>[[alignreads.py]] v2.22</td>
-	<td>A pipeline for combining the free aligner NUCmer and the free short read assembler YASRA so full alignments can be made with one command. Alignreads uses the following scripts from the Liston lab: runyasra.py, sumqual.py, and qualtofa.py. It also uses yasra and its binaries, including lastz. It must be possible for the python interpreter to find these scripts to run alignreads; this can be done saving copies of scripts in the Python26 folder, your bin folder, or modifying your .chsrc file to search where they are saved on your system (nessesary if lastz binaries cannot be found). alignreads preserves almost all of the functionality and output of the five programs it uses, so there are many options (over 40). Because there are so many modifiers, the help menu reached by not providing arguments only has what are expected to be the options used most often; to get all of the options type 'python alignreads.py -H' or 'python alignreads.py --advanced-help'. Since there are many files associated with the alignment, the output is a folder with numerous files organized into subfolders. alignreads.py can also start the pipeline from a folder produced by a previous run of alignreads.py or runyasra.py. In that case, the output is saved in the input folder and separate subfolders are created for each alignment. A more complete discription fo alignreads can be found [[HERE]]. </td>	+	<td>A pipeline for combining the free aligner NUCmer and the free short read assembler YASRA so full alignments can be made with one command. Alignreads uses the following scripts from the Liston lab: runyasra.py, sumqual.py, and qualtofa.py. It also uses yasra and its binaries, including lastz. It must be possible for the python interpreter to find these scripts to run alignreads; this can be done saving copies of scripts in the Python26 folder, your bin folder, or modifying your .chsrc file to search where they are saved on your system (nessesary if lastz binaries cannot be found). alignreads preserves almost all of the functionality and output of the five programs it uses, so there are many options (over 40). Because there are so many modifiers, the help menu reached by not providing arguments only has what are expected to be the options used most often; to get all of the options type 'python alignreads.py -H' or 'python alignreads.py --advanced-help'. Since there are many files associated with the alignment, the output is a folder with numerous files organized into subfolders. alignreads.py can also start the pipeline from a folder produced by a previous run of alignreads.py or runyasra.py. In that case, the output is saved in the input folder and separate subfolders are created for each alignment. A more complete discription of alignreads can be found [[HERE]]. </td>
	<td>A FASTA file containing reads, and A FASTA file containing a reference OR a folder from a previous run of alignreads or runyasra.</td>		<td>A FASTA file containing reads, and A FASTA file containing a reference OR a folder from a previous run of alignreads or runyasra.</td>
	<td>A folder containing all of the output from every part of the pipeline organized in subfolders based on each programs' output.</td></tr>		<td>A folder containing all of the output from every part of the pipeline organized in subfolders based on each programs' output.</td></tr>

Zachary S. L. Foster: /* Python Scripts */

Zachary S. L. Foster — Wed, 22 Sep 2010 23:53:36 GMT

Python Scripts

←Older revision		Revision as of 23:53, 22 September 2010
Line 32:		Line 32:
	<tr><th>Script Name</th><th>Description</th><th>Input File Format</th><th>Output File Format</th></tr>		<tr><th>Script Name</th><th>Description</th><th>Input File Format</th><th>Output File Format</th></tr>
	<tr><td>[[alignreads.py]] v2.22</td>		<tr><td>[[alignreads.py]] v2.22</td>
-	<td>A pipeline for combining the free aligner NUCmer and the free short read assembler YASRA so full alignments can be made with one command. Alignreads uses the following scripts from the Liston lab: runyasra.py, sumqual.py, and qualtofa.py. It also uses yasra and its binaries, including lastz. It must be possible for the python interpreter to find these scripts to run alignreads; this can be done saving copies of scripts in the Python26 folder, your bin folder, or modifying your .chsrc file to search where they are saved on your system (nessesary if lastz binaries cannot be found). alignreads preserves almost all of the functionality and output of the five programs it uses, so there are many options (over 40). Because there are so many modifiers, the help menu reached by not providing arguments only has what are expected to be the options used most often; to get all of the options type 'python alignreads.py -H' or 'python alignreads.py --advanced-help'. Since there are many files associated with the alignment, the output is a folder with numerous files organized into subfolders. alignreads.py can also start the pipeline from a folder produced by a previous run of alignreads.py or runyasra.py. In that case, the output is saved in the input folder and separate subfolders are created for each alignment. </td>	+	<td>A pipeline for combining the free aligner NUCmer and the free short read assembler YASRA so full alignments can be made with one command. Alignreads uses the following scripts from the Liston lab: runyasra.py, sumqual.py, and qualtofa.py. It also uses yasra and its binaries, including lastz. It must be possible for the python interpreter to find these scripts to run alignreads; this can be done saving copies of scripts in the Python26 folder, your bin folder, or modifying your .chsrc file to search where they are saved on your system (nessesary if lastz binaries cannot be found). alignreads preserves almost all of the functionality and output of the five programs it uses, so there are many options (over 40). Because there are so many modifiers, the help menu reached by not providing arguments only has what are expected to be the options used most often; to get all of the options type 'python alignreads.py -H' or 'python alignreads.py --advanced-help'. Since there are many files associated with the alignment, the output is a folder with numerous files organized into subfolders. alignreads.py can also start the pipeline from a folder produced by a previous run of alignreads.py or runyasra.py. In that case, the output is saved in the input folder and separate subfolders are created for each alignment. A more complete discription fo alignreads can be found [[HERE]]. </td>
	<td>A FASTA file containing reads, and A FASTA file containing a reference OR a folder from a previous run of alignreads or runyasra.</td>		<td>A FASTA file containing reads, and A FASTA file containing a reference OR a folder from a previous run of alignreads or runyasra.</td>
	<td>A folder containing all of the output from every part of the pipeline organized in subfolders based on each programs' output.</td></tr>		<td>A folder containing all of the output from every part of the pipeline organized in subfolders based on each programs' output.</td></tr>

Zachary S. L. Foster: /* Python Scripts */

Zachary S. L. Foster — Wed, 22 Sep 2010 23:33:36 GMT

Python Scripts

←Older revision		Revision as of 23:33, 22 September 2010
Line 32:		Line 32:
	<tr><th>Script Name</th><th>Description</th><th>Input File Format</th><th>Output File Format</th></tr>		<tr><th>Script Name</th><th>Description</th><th>Input File Format</th><th>Output File Format</th></tr>
	<tr><td>[[alignreads.py]] v2.22</td>		<tr><td>[[alignreads.py]] v2.22</td>
-	<td>A pipeline for combining the free aligner NUCmer and the free short read assembler YASRA so full alignments can be made with one command. Alignreads uses the following scripts from the Liston lab: runyasra.py, sumqual.py, and qualtofa.py. It must be possible for the python interpreter to find these scripts to run alignreads; this can be done saving copies of scripts in the Python26 folder, your bin folder, or modifying your .chsrc file to search where they are saved on your system. alignreads preserves almost all of the functionality and output of the five programs it uses, so there are many options (over 40). Because there are so many modifiers, the help menu reached by not providing arguments only has what are expected to be the options used most often; to get all of the options type 'python alignreads.py -H' or 'python alignreads.py --advanced-help'. Since there are many files associated with the alignment, the output is a folder with numerous files organized into subfolders. alignreads.py can also start the pipeline from a folder produced by a previous run of alignreads.py or runyasra.py. In that case, the output is saved in the input folder and separate subfolders are created for each alignment. </td>	+	<td>A pipeline for combining the free aligner NUCmer and the free short read assembler YASRA so full alignments can be made with one command. Alignreads uses the following scripts from the Liston lab: runyasra.py, sumqual.py, and qualtofa.py. It also uses yasra and its binaries, including lastz. It must be possible for the python interpreter to find these scripts to run alignreads; this can be done saving copies of scripts in the Python26 folder, your bin folder, or modifying your .chsrc file to search where they are saved on your system (nessesary if lastz binaries cannot be found). alignreads preserves almost all of the functionality and output of the five programs it uses, so there are many options (over 40). Because there are so many modifiers, the help menu reached by not providing arguments only has what are expected to be the options used most often; to get all of the options type 'python alignreads.py -H' or 'python alignreads.py --advanced-help'. Since there are many files associated with the alignment, the output is a folder with numerous files organized into subfolders. alignreads.py can also start the pipeline from a folder produced by a previous run of alignreads.py or runyasra.py. In that case, the output is saved in the input folder and separate subfolders are created for each alignment. </td>
	<td>A FASTA file containing reads, and A FASTA file containing a reference OR a folder from a previous run of alignreads or runyasra.</td>		<td>A FASTA file containing reads, and A FASTA file containing a reference OR a folder from a previous run of alignreads or runyasra.</td>
	<td>A folder containing all of the output from every part of the pipeline organized in subfolders based on each programs' output.</td></tr>		<td>A folder containing all of the output from every part of the pipeline organized in subfolders based on each programs' output.</td></tr>

Zachary S. L. Foster: /* Python Scripts */

Zachary S. L. Foster — Tue, 14 Sep 2010 23:07:28 GMT

Python Scripts

←Older revision		Revision as of 23:07, 14 September 2010
Line 58:		Line 58:
	<td>GSS Basepile output. 1000 bases per contig</td>		<td>GSS Basepile output. 1000 bases per contig</td>
	<td>tab-deliminated list of the statistics for each contig</td></tr>		<td>tab-deliminated list of the statistics for each contig</td></tr>
		+	<tr><td>[[coveragestats.py]] v1.12</td>
		+	<td>Calculates statistics for each contig in a Tablet output file and outputs them to a tab-delimitated '.txt' file following the name of their respective contig. The statistics calculated are: number of coverage values, average coverage depth, percent of bases with a coverage greater than zero, and maximum coverage depth. coveragestats.py requires that contig names have at least one non-numerical character.</td>
		+	<td>Tablet coverage output file (.txt) or one of the same format</td>
		+	<td>Tab-delimitated '.txt' file</td></tr>
	<tr><td>[[fatomum.py]] v1.0</td>		<tr><td>[[fatomum.py]] v1.0</td>
	<td>Creates a MUMmer-style output from a FASTA file of aligned sequences. The first line is always considered the reference and is not included in the output. The start index in terms of the reference, and in terms of the contig, as well as the length are recorded for each matching section of the contig.</td>		<td>Creates a MUMmer-style output from a FASTA file of aligned sequences. The first line is always considered the reference and is not included in the output. The start index in terms of the reference, and in terms of the contig, as well as the length are recorded for each matching section of the contig.</td>

Zachary S. L. Foster: /* Python Scripts */

Zachary S. L. Foster — Wed, 08 Sep 2010 19:04:19 GMT

Python Scripts

Zachary S. L. Foster at 21:34, 9 July 2010

Zachary S. L. Foster — Fri, 09 Jul 2010 21:34:35 GMT

←Older revision		Revision as of 21:34, 9 July 2010
Line 1:		Line 1:
-	This page contains the source code for some of the bioinformatics scripts used by the Liston Lab. Most of the scripts are written in Python 2.6.4 and are designed for Unix systems. A few are written as a list of unix commands designed to be executables.	+	This page contains the source code for some of the bioinformatics scripts used by the Liston Lab. Most of the scripts are written in Python 2.6.4 and are designed for Unix systems, although most of them can be run on windows systems. A few are written as a list of unix commands designed to be executables.

Zachary S. L. Foster: /* Python Scripts */

Zachary S. L. Foster — Fri, 09 Jul 2010 21:30:45 GMT

Python Scripts

←Older revision		Revision as of 21:30, 9 July 2010
Line 32:		Line 32:
	<tr><th>Script Name</th><th>Discription</th><th>Input File Format</th><th>Output File Format</th></tr>		<tr><th>Script Name</th><th>Discription</th><th>Input File Format</th><th>Output File Format</th></tr>
	<tr><td>[[alignreads.py]] v2.0</td>		<tr><td>[[alignreads.py]] v2.0</td>
-	<td>A pipeline for combining the free aligner NUCmer and the free short read assembler YASRA so full alignments can be made with one command. Alignreads uses the following scripts from the Liston lab: runyasra.py, sumqual.py, and qualtofa.py. It must be possible for the python interpreter to find these scripts to run alignreads; this can be done saving copies of scripts in the Python26 folder, your bin folder, or modifying ~~you~~ .chsrc file to search where they are saved on your system. alignreads preserves almost all of the functionality and output of the 5 programs it uses, so there are many options (over 40). Because there are so many modifiers, the help menu reached by not providing arguments only has what are expected to be the options used most often; to get all of the options type 'python alignreads.py -H' or 'python alignreads.py --advanced-help'. Since there are many files associated with the alignment the output is a folder with numerous files organized into subfolders. alignreads.py can also start the pipeline from a folder produced by a previous run of alignreads.py or runyasra.py. In that case, the output is saved in the input folder. </td>	+	<td>A pipeline for combining the free aligner NUCmer and the free short read assembler YASRA so full alignments can be made with one command. Alignreads uses the following scripts from the Liston lab: runyasra.py, sumqual.py, and qualtofa.py. It must be possible for the python interpreter to find these scripts to run alignreads; this can be done saving copies of scripts in the Python26 folder, your bin folder, or modifying your .chsrc file to search where they are saved on your system. alignreads preserves almost all of the functionality and output of the 5 programs it uses, so there are many options (over 40). Because there are so many modifiers, the help menu reached by not providing arguments only has what are expected to be the options used most often; to get all of the options type 'python alignreads.py -H' or 'python alignreads.py --advanced-help'. Since there are many files associated with the alignment the output is a folder with numerous files organized into subfolders. alignreads.py can also start the pipeline from a folder produced by a previous run of alignreads.py or runyasra.py. In that case, the output is saved in the input folder. </td>
	<td>A FASTA file containing reads, and A FASTA file containg a reference OR a folder from a previous run of alingreades or runyasra.</td>		<td>A FASTA file containing reads, and A FASTA file containg a reference OR a folder from a previous run of alingreades or runyasra.</td>
	<td>A folder containing all of the ouput from every part of the pipeline organized in subfolders based on each programs' output.</td></tr>		<td>A folder containing all of the ouput from every part of the pipeline organized in subfolders based on each programs' output.</td></tr>
Line 67:		Line 67:
	<td>FASTA file same number of sequences as the input file but, with gaps removed</td></tr>		<td>FASTA file same number of sequences as the input file but, with gaps removed</td></tr>
	<tr><td>[[runyasra.py]] v1.0</td>		<tr><td>[[runyasra.py]] v1.0</td>
-	<td>Runs the free short read aligner YASRA and orgainzes its output into a folder.</td>	+	<td>Runs the free short read aligner YASRA and orgainzes its output into a folder. Links to the input data are also saved in the output folder.</td>
	<td>A FASTA file containing reads, and a FASTA file containg a reference OR a folder from a previous run of alingreades or runyasra.</td>		<td>A FASTA file containing reads, and a FASTA file containg a reference OR a folder from a previous run of alingreades or runyasra.</td>
	<td>A folder containing all of YASRA's output as well as links to the input data.</td></tr>		<td>A folder containing all of YASRA's output as well as links to the input data.</td></tr>

Zachary S. L. Foster: /* Python Scripts */

Zachary S. L. Foster — Fri, 09 Jul 2010 21:25:52 GMT

Python Scripts

←Older revision		Revision as of 21:25, 9 July 2010
Line 34:		Line 34:
	<td>A pipeline for combining the free aligner NUCmer and the free short read assembler YASRA so full alignments can be made with one command. Alignreads uses the following scripts from the Liston lab: runyasra.py, sumqual.py, and qualtofa.py. It must be possible for the python interpreter to find these scripts to run alignreads; this can be done saving copies of scripts in the Python26 folder, your bin folder, or modifying you .chsrc file to search where they are saved on your system. alignreads preserves almost all of the functionality and output of the 5 programs it uses, so there are many options (over 40). Because there are so many modifiers, the help menu reached by not providing arguments only has what are expected to be the options used most often; to get all of the options type 'python alignreads.py -H' or 'python alignreads.py --advanced-help'. Since there are many files associated with the alignment the output is a folder with numerous files organized into subfolders. alignreads.py can also start the pipeline from a folder produced by a previous run of alignreads.py or runyasra.py. In that case, the output is saved in the input folder. </td>		<td>A pipeline for combining the free aligner NUCmer and the free short read assembler YASRA so full alignments can be made with one command. Alignreads uses the following scripts from the Liston lab: runyasra.py, sumqual.py, and qualtofa.py. It must be possible for the python interpreter to find these scripts to run alignreads; this can be done saving copies of scripts in the Python26 folder, your bin folder, or modifying you .chsrc file to search where they are saved on your system. alignreads preserves almost all of the functionality and output of the 5 programs it uses, so there are many options (over 40). Because there are so many modifiers, the help menu reached by not providing arguments only has what are expected to be the options used most often; to get all of the options type 'python alignreads.py -H' or 'python alignreads.py --advanced-help'. Since there are many files associated with the alignment the output is a folder with numerous files organized into subfolders. alignreads.py can also start the pipeline from a folder produced by a previous run of alignreads.py or runyasra.py. In that case, the output is saved in the input folder. </td>
	<td>A FASTA file containing reads, and A FASTA file containg a reference OR a folder from a previous run of alingreades or runyasra.</td>		<td>A FASTA file containing reads, and A FASTA file containg a reference OR a folder from a previous run of alingreades or runyasra.</td>
-	<td>A folder containing all of the ouput from every part of the pipeline</td></tr>	+	<td>A folder containing all of the ouput from every part of the pipeline organized in subfolders based on each programs' output.</td></tr>
	<tr><td>[[allcomb.py]] v1.2</td>		<tr><td>[[allcomb.py]] v1.2</td>
	<td>Finds all possible combinations of elements (or group of elements) for every line of a tab-delimitated text file and output every combination, on its own line, to a .txt file.</td>		<td>Finds all possible combinations of elements (or group of elements) for every line of a tab-delimitated text file and output every combination, on its own line, to a .txt file.</td>
Line 58:		Line 58:
	<td>GSS Basepile output. 1000 bases per contig</td>		<td>GSS Basepile output. 1000 bases per contig</td>
	<td>tab-deliminated list of the statistics for each contig</td></tr>		<td>tab-deliminated list of the statistics for each contig</td></tr>
		+	<tr><td>[[fatomum.py]] v1.0</td>
		+	<td>Creates a MUMmer-style output from a FASTA file of aligned sequences. The first line is always considered the reference and is not included in the output. The start index in terms of the reference, and in terms of the contig, as well as the length are recorded for each matching section of the contig.</td>
		+	<td>FASTA file containing a reference and its aligned contigs</td>
		+	<td>MUMmuer-like output</td></tr>
	<tr><td>[[gapstrip.py]] v1.0</td>		<tr><td>[[gapstrip.py]] v1.0</td>
	<td>Removes all of the gaps (i.e. '-') that are common (i.e. at the same index) to all of the sequences in an aligned FASTA file. By default, the first sequence is considered the reference and is excluded from the analysis, but the number of sequences that are treated as such can be changed.</td>		<td>Removes all of the gaps (i.e. '-') that are common (i.e. at the same index) to all of the sequences in an aligned FASTA file. By default, the first sequence is considered the reference and is excluded from the analysis, but the number of sequences that are treated as such can be changed.</td>
	<td>FASTA file containing two or more aligned sequences.</td>		<td>FASTA file containing two or more aligned sequences.</td>
	<td>FASTA file same number of sequences as the input file but, with gaps removed</td></tr>		<td>FASTA file same number of sequences as the input file but, with gaps removed</td></tr>
		+	<tr><td>[[runyasra.py]] v1.0</td>
		+	<td>Runs the free short read aligner YASRA and orgainzes its output into a folder.</td>
		+	<td>A FASTA file containing reads, and a FASTA file containg a reference OR a folder from a previous run of alingreades or runyasra.</td>
		+	<td>A folder containing all of YASRA's output as well as links to the input data.</td></tr>
	<tr><td>[[sumqual.py]] v1.7</td>		<tr><td>[[sumqual.py]] v1.7</td>
	<td>A script used in conjunction with the free aligners NUCmer and the free short read assembler YASRA. Creates a file with a similar to the format to that of the .qual file YASRA outputs, except each base can have multiple sets of quality values. Uses a .qual file containing the quality scores of multiple contigs and the output file of NUCmer to make a "consensus" .qual file, which contains the quality values in terms of the reference. Given the reference used with NUCmer to get the alignment, sumqual can output a fasta file containing all the sequence for which there is no quality data.</td>		<td>A script used in conjunction with the free aligners NUCmer and the free short read assembler YASRA. Creates a file with a similar to the format to that of the .qual file YASRA outputs, except each base can have multiple sets of quality values. Uses a .qual file containing the quality scores of multiple contigs and the output file of NUCmer to make a "consensus" .qual file, which contains the quality values in terms of the reference. Given the reference used with NUCmer to get the alignment, sumqual can output a fasta file containing all the sequence for which there is no quality data.</td>
Line 72:		Line 80:
	<td>quality file (.qual); accepts sumqual.py's output</td>		<td>quality file (.qual); accepts sumqual.py's output</td>
	<td>FASTA file</td></tr>		<td>FASTA file</td></tr>
-	~~<tr><td>[[fatomum.py]] v1.0</td>~~
-	<td>Creates a MUMmer-style output from a FASTA file of aligned sequences. The first line is always considered the reference and is not included in the output. The start index in terms of the reference, and in terms of the contig, as well as the length are recorded for each matching section of the contig.</td>
-	~~<td>FASTA file containing a reference and its aligned contigs</td>~~
-	~~<td>MUMmuer-like output</td></tr>~~
	</table>		</table>

Zachary S. L. Foster: /* Python Scripts */

Zachary S. L. Foster — Fri, 09 Jul 2010 21:12:32 GMT

Python Scripts

←Older revision		Revision as of 21:12, 9 July 2010
Line 31:		Line 31:
	<table border=2>		<table border=2>
	<tr><th>Script Name</th><th>Discription</th><th>Input File Format</th><th>Output File Format</th></tr>		<tr><th>Script Name</th><th>Discription</th><th>Input File Format</th><th>Output File Format</th></tr>
		+	<tr><td>[[alignreads.py]] v2.0</td>
		+	<td>A pipeline for combining the free aligner NUCmer and the free short read assembler YASRA so full alignments can be made with one command. Alignreads uses the following scripts from the Liston lab: runyasra.py, sumqual.py, and qualtofa.py. It must be possible for the python interpreter to find these scripts to run alignreads; this can be done saving copies of scripts in the Python26 folder, your bin folder, or modifying you .chsrc file to search where they are saved on your system. alignreads preserves almost all of the functionality and output of the 5 programs it uses, so there are many options (over 40). Because there are so many modifiers, the help menu reached by not providing arguments only has what are expected to be the options used most often; to get all of the options type 'python alignreads.py -H' or 'python alignreads.py --advanced-help'. Since there are many files associated with the alignment the output is a folder with numerous files organized into subfolders. alignreads.py can also start the pipeline from a folder produced by a previous run of alignreads.py or runyasra.py. In that case, the output is saved in the input folder. </td>
		+	<td>A FASTA file containing reads, and A FASTA file containg a reference OR a folder from a previous run of alingreades or runyasra.</td>
		+	<td>A folder containing all of the ouput from every part of the pipeline</td></tr>
	<tr><td>[[allcomb.py]] v1.2</td>		<tr><td>[[allcomb.py]] v1.2</td>
	<td>Finds all possible combinations of elements (or group of elements) for every line of a tab-delimitated text file and output every combination, on its own line, to a .txt file.</td>		<td>Finds all possible combinations of elements (or group of elements) for every line of a tab-delimitated text file and output every combination, on its own line, to a .txt file.</td>

←Older revision		Revision as of 19:04, 8 September 2010
Line 30:		Line 30:

	<table border=2>		<table border=2>
-	<tr><th>Script Name</th><th>~~Discription~~</th><th>Input File Format</th><th>Output File Format</th></tr>	+	<tr><th>Script Name</th><th>Description</th><th>Input File Format</th><th>Output File Format</th></tr>
-	<tr><td>[[alignreads.py]] v2.0</td>	+	<tr><td>[[alignreads.py]] v2.22</td>
-	<td>A pipeline for combining the free aligner NUCmer and the free short read assembler YASRA so full alignments can be made with one command. Alignreads uses the following scripts from the Liston lab: runyasra.py, sumqual.py, and qualtofa.py. It must be possible for the python interpreter to find these scripts to run alignreads; this can be done saving copies of scripts in the Python26 folder, your bin folder, or modifying your .chsrc file to search where they are saved on your system. alignreads preserves almost all of the functionality and output of the 5 programs it uses, so there are many options (over 40). Because there are so many modifiers, the help menu reached by not providing arguments only has what are expected to be the options used most often; to get all of the options type 'python alignreads.py -H' or 'python alignreads.py --advanced-help'. Since there are many files associated with the alignment the output is a folder with numerous files organized into subfolders. alignreads.py can also start the pipeline from a folder produced by a previous run of alignreads.py or runyasra.py. In that case, the output is saved in the input folder. </td>	+	<td>A pipeline for combining the free aligner NUCmer and the free short read assembler YASRA so full alignments can be made with one command. Alignreads uses the following scripts from the Liston lab: runyasra.py, sumqual.py, and qualtofa.py. It must be possible for the python interpreter to find these scripts to run alignreads; this can be done saving copies of scripts in the Python26 folder, your bin folder, or modifying your .chsrc file to search where they are saved on your system. alignreads preserves almost all of the functionality and output of the five programs it uses, so there are many options (over 40). Because there are so many modifiers, the help menu reached by not providing arguments only has what are expected to be the options used most often; to get all of the options type 'python alignreads.py -H' or 'python alignreads.py --advanced-help'. Since there are many files associated with the alignment, the output is a folder with numerous files organized into subfolders. alignreads.py can also start the pipeline from a folder produced by a previous run of alignreads.py or runyasra.py. In that case, the output is saved in the input folder and separate subfolders are created for each alignment. </td>
-	<td>A FASTA file containing reads, and A FASTA file ~~containg~~ a reference OR a folder from a previous run of ~~alingreades~~ or runyasra.</td>	+	<td>A FASTA file containing reads, and A FASTA file containing a reference OR a folder from a previous run of alignreads or runyasra.</td>
-	<td>A folder containing all of the ~~ouput~~ from every part of the pipeline organized in subfolders based on each programs' output.</td></tr>	+	<td>A folder containing all of the output from every part of the pipeline organized in subfolders based on each programs' output.</td></tr>
	<tr><td>[[allcomb.py]] v1.2</td>		<tr><td>[[allcomb.py]] v1.2</td>
	<td>Finds all possible combinations of elements (or group of elements) for every line of a tab-delimitated text file and output every combination, on its own line, to a .txt file.</td>		<td>Finds all possible combinations of elements (or group of elements) for every line of a tab-delimitated text file and output every combination, on its own line, to a .txt file.</td>
Line 40:		Line 40:
	<td>tab-delimitated .txt file</td></tr>		<td>tab-delimitated .txt file</td></tr>
	<tr><td>[[baseanno.py]] v1.0</td>		<tr><td>[[baseanno.py]] v1.0</td>
-	<td>Converts a file containing a list of annotations, as well as each of their respective start and stop indices, into a file containing a list of base indices, each followed by any annotations that apply at that specific base. Each line of the input file is expected to be whitespace-~~deliminated~~, however if your annotations have spaces in them, the script can be made to enforce tab-delimitation. The output file is always tab-delaminated.</td>	+	<td>Converts a file containing a list of annotations, as well as each of their respective start and stop indices, into a file containing a list of base indices, each followed by any annotations that apply at that specific base. Each line of the input file is expected to be whitespace-delimitated, however if your annotations have spaces in them, the script can be made to enforce tab-delimitation. The output file is always tab-delaminated.</td>
	<td>[Annotation text] [Start Index] [End Index]</td>		<td>[Annotation text] [Start Index] [End Index]</td>
-	<td>[Base Index] [Annotation1] [Annotation2] ... [~~AnnotationX~~]</td></tr>	+	<td>[Base Index] [Annotation1] [Annotation2] ... [Annotation]</td></tr>
	<tr><td>[[basediff.py]] v1.2</td>		<tr><td>[[basediff.py]] v1.2</td>
-	<td>Finds base differences between multiple aligned sequences in a single FASTA file and output a tab-~~deliminated~~ txt file containing the base values for all the sequences at the index where the difference occurred. The script has many modifiers that change what is considered a difference.</td>	+	<td>Finds base differences between multiple aligned sequences in a single FASTA file and output a tab-delimitated txt file containing the base values for all the sequences at the index where the difference occurred. The script has many modifiers that change what is considered a difference.</td>
	<td>FASTA file containing two or more aligned sequences</td>		<td>FASTA file containing two or more aligned sequences</td>
	<td>.txt file in the following format: [Base index] [Seq 1 value] [Seq 2 value] ... [Seq N value]</td></tr>		<td>.txt file in the following format: [Base index] [Seq 1 value] [Seq 2 value] ... [Seq N value]</td></tr>
Line 66:		Line 66:
	<td>FASTA file containing two or more aligned sequences.</td>		<td>FASTA file containing two or more aligned sequences.</td>
	<td>FASTA file same number of sequences as the input file but, with gaps removed</td></tr>		<td>FASTA file same number of sequences as the input file but, with gaps removed</td></tr>
-	<tr><td>[[runyasra.py]] v1.0</td>	+	<tr><td>[[runyasra.py]] v1.02</td>
-	<td>Runs the free short read aligner YASRA and ~~orgainzes~~ its output into a folder. Links to the input data are also saved in the output folder.</td>	+	<td>Runs the free short read aligner YASRA and organizes its output into a folder. Links to the input data are also saved in the output folder.</td>
	<td>A FASTA file containing reads, and a FASTA file containg a reference OR a folder from a previous run of alingreades or runyasra.</td>		<td>A FASTA file containing reads, and a FASTA file containg a reference OR a folder from a previous run of alingreades or runyasra.</td>
	<td>A folder containing all of YASRA's output as well as links to the input data.</td></tr>		<td>A folder containing all of YASRA's output as well as links to the input data.</td></tr>
-	<tr><td>[[sumqual.py]] v1.7</td>	+	<tr><td>[[sumqual.py]] v2.12</td>
-	<td>A script used in conjunction with the free aligners NUCmer and the free short read assembler YASRA. Creates a file with a similar to the format to that of the .qual file YASRA outputs, except each base can have multiple sets of quality values. Uses a .qual file containing the quality scores of multiple contigs and the output file of NUCmer to make a "consensus" .qual file, which contains the quality values in terms of the reference. ~~Given the reference used with NUCmer~~ to ~~get the alignment, sumqual can output~~ a fasta ~~file containing all the sequence for which there is no quality data~~.</td>	+	<td>A script used in conjunction with the free aligners NUCmer and the free short read assembler YASRA. Creates a file with a similar to the format to that of the .qual file YASRA outputs, except each base can have multiple sets of quality values. The format preserves all of the alignment, indel, quality, and reference data from the input files. Uses a .qual file containing the quality scores of multiple contigs and the output file of NUCmer to make a "consensus" .qual file, which contains the quality values in terms of the reference. Sumqual's output can be converted to a fasta format using qualtofa.py. td>
	<td><ul><li>FASTA file containing the sequence of multiple contigs</li>		<td><ul><li>FASTA file containing the sequence of multiple contigs</li>
	<li>NUCmer output file</li>		<li>NUCmer output file</li>
	<li>Reference file used with Nucmer</li></ul></td>		<li>Reference file used with Nucmer</li></ul></td>
	<td>similar to .qual input file, except each base can hold the information from more than one contig</td></tr>		<td>similar to .qual input file, except each base can hold the information from more than one contig</td></tr>
-	<tr><td>[[qualtofa.py]] v2.3</td>	+	<tr><td>[[qualtofa.py]] v2.52</td>
-	<td>Selectively extracts the sequence from a quality file and outputs it to a FASTA file. It is designed to be used in conjunction with sumqual.py and accepts its consensus-style output. Conflicting bases from a consensus-style output are condensed into IUPAC ambiguity codes or gaps ~~are~~ added to the reference to ~~accomodate~~ repeated or unmatched sequence. It has ~~an option to mask all bases with a quality value of less than a user-specified threshold. Another threshold can be assigned for SNPs if a reference is given~~. </td>	+	<td>Selectively extracts the sequence from a quality file and outputs it to a FASTA file. It is designed to be used in conjunction with sumqual.py and accepts its consensus-style output, which preserves all of the alignment, quality and reference information. Conflicting bases on overlapping matches from a consensus-style output are condensed into IUPAC ambiguity codes, or gaps can be added to the reference to accommodate repeated or unmatched sequence if a option is used. It has several masking and trimming options. </td>
	<td>quality file (.qual); accepts sumqual.py's output</td>		<td>quality file (.qual); accepts sumqual.py's output</td>
	<td>FASTA file</td></tr>		<td>FASTA file</td></tr>
	</table>		</table>