Program to find probes

From Ucsbgalaxy
Revision as of 13:10, 20 February 2014 by 128.111.227.66 (talk)
Jump to: navigation, search

project in /home/likewise-open/ADS/oakley/labdata/OpsinGenomes

Program flow:

  • 1. Creates a BLAST+ database using sequence.fasta
  • 2. Reads in sequence.fasta sequence by sequence
    • 2b. creates a hash of all sequence ids to false
  • 3. Reads through the sequence.fasta file again, this time running blastn on each sequence
    • 3b. For each hit, set the hash with the display_id as key to true