Wednesday, July 2, 2014

Plans for today:

[1] Find missing data for high school student
[2] Figure out how to teach RaxML to a high school student (JW slides) COMPLETE
[3] Get Gazey-Staley algorithm set up in R and running. COMPLETE
      -Okay, so we don't actually have to do this step, as its just a measure of cloning success (and since we didn't clone. We don't really need to measure this. What we do need to measure (but not right away, is the # of sequences (reads) v. # genes (unique contigs). This will tell us if we have sampled enough of the genome for the OR genes given our degenerate primers.

Anyways, to download the R code for the Gazey-Staley algorithm along with a heterogeneity test (with very good documentation), visit:

[4] Figure out most recent problem with babbler paper.
[5] Set up ORA pipeline COMPLETE
To install ORA (after having install BioPerl, HMMER, and FASTA), download the .zip file from:

Unzip and navigate to the ORA folder.

#run the perl script to make the "makefile"

perl Makefile.PL
sudo make install
perl Build.PL

WOOOOHOOOO got it working.
There are two main issues (assuming you have everything installed).
1) You must run the perl script in the /ora-*/scripts folder because it must access the .hmm files
2) When you install ORA, it includes the older version of HMMER. The error message you get when you run "", is: 
binary auxfiles are in an outdated HMMER format (3/b); please hmmpress your HMM file again

To fix this, do the following (assuming you have HMMER installed).
-In the /ora-*/scripts folder, there should be several files related to HMMER
Delete the following: or.hmm.h3f     or.hmm.h3i     or.hmm.h3m     or.hmm.h3p
On the command line, type:
hmmpress or.hmm

That should fix things.

To run the ORA:
perl --sequence ~/Documents/Analyses/Carrollia_OR/GPC-trim.deg.fasta
#still playing with all of the fun parameters

Now I need to trim the data of a certain length.
By the end of the weekend, it would be great to make a tree from the sequences.

No comments:

Post a Comment