Difference between revisions of "Program to find probes"

From Ucsbgalaxy
Jump to: navigation, search
(Created page with "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 2…")
 
Line 2: Line 2:
  
 
Program flow:
 
Program flow:
1. Creates a BLAST+ database using sequence.fasta
+
*1. Creates a BLAST+ database using sequence.fasta
2. Reads in sequence.fasta sequence by sequence  
+
*2. Reads in sequence.fasta sequence by sequence  
2b. creates a hash of all sequence ids to false
+
**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
+
*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
+
**3b. For each hit, set the hash with the display_id as key to true

Revision as of 13:10, 20 February 2014

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