edwardslab

  • Increase font size
  • Default font size
  • Decrease font size
Home Kate
Kate's lab blog

Perl Smith Waterman script

The following is a small script to perform a Smith Waterman alignment, and to calculate the percent identity between two sequences.

 

Read more...
 

Calculating the average and standard deviation in perl

The following two subroutines can be used as a drop in replacement for Math::NumberCruncher. This is useful when trying to limit the number of dependencies of your perl script.

Read more...
 

Perl subroutine to extract fasta sequences from a file

The following perl subroutine will read in a fasta formatted file, parse the file, and return all the sequences in a reference to a hash table.

Read more...
 

PHACTS flowchart

I have been working on getting an understandable flowchart constructed for PHACTS.  Because of PHACTS odd complexity it has been difficult to get a meaningful and understandable chart.  Here is the most current version.
 

Varying Variant Counts



I looked at different strains of the same bacterium, to compare how variant counts might vary. This is of concern to me since I have Biolog plates for a species but not a specific taxon. Below is a table of variant counts returned for the compound C00022, pyruvate.(* an integer denotes that many positive variant counts returns for that reaction; while a 0.01 denotes a zero or blank variant count; a 0.0001 denotes that a NULL was returned for that reaction and taxon)
Actinobacillus pleuropneumoniae strains
L20
serovar1 S4074
7.01 0.0001 23.08 0 2 0.0001 2 2 4.01 4.01 0.0001 2 1 0.0001 13.11 0.0001 0.0001 2.01 0.0001 21.05 1 21.05 19.05 0.0001 3.03 12 2 0 11.04 7.01 0.0001 2.01 10.05 0.0001 0.0001 3 0.0001 2 0.0001 1 4.03 0.0001 0.0001 0.01 4.01 0.0001 0.0001 2 23.06 0.0001 0 0.0001 4.01 0.0001 0 4.01 0.0001 0.03 0.0001 4.01 3.02 4.01 0.0001 0.0001 6.05 0.0001 1 2 4 4 3 0.0001 0.0001 0 0 4.01 0.0001 1.01 0 4 0.0001 0.02 0.02 1 0.01 1.01 0.0001 0.0001 0.0001 0.0001 0.01 0.01 0.0001 5.03 0.0001 0.0001 0.0001 6.02 4.01 0 0.0001 0.01 0.0001 1 0 4.01 6.05 1.03 0.0001 1.01 0.0001 0.0001 0.0001 3.02 7.01 0.0001 0.0001 0.0001 1 1.01 0.0001 0.0001 0.0001 0.0001 0.0001 0.0001 0 1 7.01 0.0001 0.02 2.09 2.01 0.0001 0.0001 4 2.02 3 0.01 0.0001 1.01 2.03 0.0001 2.09 2.09 2.09 1 0.0001 0.0001 0.0001 0.0001 0.0001 1 0 0.0001 0.0001 0.0001 0.01
7.01 0.0001 20.09 0 2 0.0001 2 2 5 5 0.0001 2 1 0.0001 20.02 0.0001 0.0001 2.01 0.0001 20.06 1 20.06 18.06 0.0001 3.03 12 2 0 7.08 7.01 0.0001 2.01 12.03 0.0001 0.0001 3 0.0001 2 0.0001 1 4.01 0.0001 0.0001 0 5 0.0001 0.0001 2 22.07 0.0001 0 0.0001 5 0.0001 0 5 0.0001 0.03 0.0001 4.01 4.01 5 0.0001 0.0001 7.02 0.0001 1 2 4 4 3 0.0001 0.0001 0 0 5 0.0001 1.01 0 4 0.0001 2 2 1 0.01 1.01 0.0001 0.0001 0.0001 0.0001 0.01 0.01 0.0001 6 0.0001 0.0001 0.0001 5.02 5 0 0.0001 0.01 0.0001 1 0 5 7.02 0.03 0.0001 2 0.0001 0.0001 0.0001 4.01 7.01 0.0001 0.0001 0.0001 1 2 0.0001 0.0001 0.0001 0.0001 0.0001 0.0001 0 1 7.01 0.0001 2 4.06 1.02 0.0001 0.0001 4 2.02 1.02 0.01 0.0001 2 2.03 0.0001 4.06 4.06 4.06 1 0.0001 0.0001 0.0001 0.0001 0.0001 1 0 0.0001 0.0001 0.0001 0.01
 
  • «
  •  Start 
  •  Prev 
  •  1 
  •  2 
  •  Next 
  •  End 
  • »


Page 1 of 2