Topic: Linux Basics

  1. Download Halobacterium proteome and inspect it
     wget # downloads code from this slide
     module load ncbi-blast/2.2.26
     gunzip GCA_000069025.1_ASM6902v1_protein.faa.gz
     mv GCA_000069025.1_ASM6902v1_protein.faa halobacterium.faa
     less halobacterium.faa # press q to quit
  2. How many protein sequences are stored in the downloaded file?
     grep '>' halobacterium.faa | wc
     grep '^>' halobacterium.faa --count
  3. How many proteins contain the pattern WxHxxH or WxHxxHH?
     egrep 'W.H..H{1,2}' halobacterium.faa --count
  4. Use less to find IDs for pattern matches or use awk
     awk --posix -v RS='>' '/W.H..(H){1,2}/ { print ">" $0;}' halobacterium.faa | less
     awk --posix -v RS='>' '/W.H..(H){1,2}/ { print ">" $0;}' halobacterium.faa | grep '^>' | cut -c 2- | cut -f 1 -d\ > myIDs
  5. Create a BLASTable database with formatdb
     formatdb -i halobacterium.faa -p T -o
  6. Query BLASTable database by IDs stored in a file (e.g. myIDs)
     fastacmd -d halobacterium.faa -i myIDs > myseq.fasta
  7. Run BLAST search for sequences stored in myseq.fasta
     blastall -p blastp -i myseq.fasta -d halobacterium.faa -o blastp.out -e 1e-6 -v 10 -b 10
     blastall -p blastp -i myseq.fasta -d halobacterium.faa -m 8 -e 1e-6 >

Additional exercise material in Linux Manual

Homework assignment

Perform above analysis on the following Escherichia coli strain: Record result from final BLAST command (with -m 8) in text file.

Homework submission

Upload result file to your private course GitHub repository under Homework/HW2/HW2.txt.

Due date

Most homeworks will be due one week after they are assigned. This one is due on Thu, April 13th at 6:00 PM.