Skip to content

Ruby script for calling bacterial probes. Bacterial genome files cut up to 40mers, screened for suitability, and blasted with a local copy of NCBI Blast, results narrow down to lists of good candidates for each genome in the queue that are unique at the species or genus level.

Notifications You must be signed in to change notification settings

edallen/blaster

Repository files navigation

Blaster README
The Blaster project is copyright 2008, 2009, Stanford University, released as open source under the MIT license. 

Will make a proper README later. For now, the blaster_poster.ppt file
 (I know, ugh, Powerpoint...) has the best documentation on the project. The
file blaster_plan.txt has some older working notes to myself that may be 
useful as well.

As will be evident in the code, this was a "my first big ruby script outside
of Rails" and I would have refactored the code quite a bit, finished more tests,
etc, had my time at Stanford continued long enough for me to go on with it. 
But it does run and called a lot of good probes that Richard Hyman had synthesized
 and used successfully. Releasing as open source at Richard's request.

When the Other Genomic database splits like the human one did, possibly already, since this code was written in early 2008, this project will take some modification, mostly in bacterium.rb and an extensions list for it in config.yaml,  like the one for the human database there to iterate over the parts as it does for the Human database.

A more sophisticated script running on a tighter schedule would parallelize to threads feeding multiple blast instances, possibly on multiple machines. This simple implementation, running blast on a dedicated Linux PC locally, generally got through two or three genomes a day, adequate to our purposes.

Ed Allen
[email protected] 

About

Ruby script for calling bacterial probes. Bacterial genome files cut up to 40mers, screened for suitability, and blasted with a local copy of NCBI Blast, results narrow down to lists of good candidates for each genome in the queue that are unique at the species or genus level.

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages