Skip to content

Latest commit

 

History

History
11 lines (11 loc) · 602 Bytes

StarshipExtraAccessoryGenes.md

File metadata and controls

11 lines (11 loc) · 602 Bytes

Identification of Interesting genes in other starship copies

  1. BLAST fgenesh proteins against the 211 B71 predicted proteins
blastp -query fgenesh.pep -subject 211B71.pep -evalue 1e-5 -outfmt '6 qseqid sseqid qlen length' | awk '$4 > $3*0.8' . myMatchingProteins.pep
  1. Identify proteins that show alignment across 80% of their lengths
  2. Create a list of matching proteins
  3. Create table of annoation data (protein ID, motifs and functions) for the proteins found in other starship copies
  4. Use R to Filter the table to remove those proteins that shiwed alignment to the B71