Liston:Computer Scripts

From OpenWetWare
Revision as of 21:43, 9 January 2010 by Zachary S. L. Foster (talk | contribs) (Python Scripts)
Jump to: navigation, search

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.

Python Script Conventions

The scripts must be compiled using a Python compiler in the following format:

   python [modifiers] <Arguments>. 

For example, in order to run the script one could enter the following into an Unix shell:

   python -c -v ../myQualFile.qual ../myMumFile

This would compile and run the script with the modifiers -c and -v, using myQualFile.qual and myMumFile as arguments. All of the scripts save their output in a file in the current working directory, with a name usually composed of some combination of the arguments and the name of the script. However, one can save the output anywhere, under any name, using the following technique:

   python -c -v ../myQualFile.qual ../myMumFile > ../myOutput.ext

The order in which the modifiers are given is not important, however, the order of the required arguments is important. For Example the above modifiers could be entered in the opposite order (-v -c), but the two file paths need to be in a predetermined order. Some scripts have modifiers that require arguments of their own. These modifier arguments should be written directly after their respective modifier. For example, if the above modifier, -c, had a argument, one would type,

   python -c theArgument -v ../myQualFile.qual ../myMumFile

Every Script has a description of what it does and how/when to use it in its source code. The list all the modifiers that the script supports and what they do is also included. A similar help menu can be viewed by calling the script with no arguments. For example, typing the following,


would cause a help menu to be printed to the screen.

Python Scripts

Script NameDiscriptionInput File FormatOutput File Format 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-delaminated, however if your annotations have spaces in them, the script can be made to enforce tab-delimitation. The output file is always tab-delaminated. [Annotation text] [Start Index] [End Index] [Base Index] [Annotation1] [Annotation2] ... [AnnotationX] Finds base differences between multiple aligned sequences in a single FASTA file and output a tab-delaminated 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. FASTA file containing two or more aligned sequences .txt file in the following format: [Base index] [Seq 1 value] [Seq 2 value] ... [Seq N value] Performs and outputs various statistical tests for each contig in a GSS basepile output. it is assumed that the information for each contig is 1000 bases long. The following statistics are outputted in a tab-delimitated list for each contig by the script:
  • Reference Match Length: The index of the last known base (i.e. not 'N')
  • Target Match Sum: The number of known bases (i.e. not 'N')
  • Coverage Proportion: proportion of 'Target match sum' to ' Reference match length'
  • Average Density: The average density value for all bases within ' Reference match length'
  • Median Density: median density of the entire range of density values
GSS Basepile output. 1000 bases per contig tab-deliminated list of the statistics for each contig 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. FASTA file containing two or more aligned sequences. FASTA file same number of sequences as the input file but, with gaps removed