Assembly Scavenger HuntΒΆ

This exercise is meant to bring together knowledge from the whole week, and also just be fun. I’ve taken some pain text, embedded it in DNA with a simple algorithm along with some random sequence, and fragmented it to produce reads. Your job will be to assemble the reads and put the results back through the script to retrieve the original message.

We’ll need a few things:

cd /mnt
mkdir scavenger-hunt
cd scavenger-hunt

curl -O
curl -O
curl -O

You should have velvet already, but if not:

cd /root
curl -O
tar xzf velvet_1.2.10.tgz
cd velvet_1.2.10
cp velvet? /usr/local/bin

cd /mnt/scavenger-hunt

Of course, sequencing chemistry is always improving, so you may want to use:

curl -O

You’ll want to use velvet to assemble scavenger_reads.fa. They’re 36-base singled-ended reads, and like an actual metagenome, have variable coverage. This means that you might need to do some parameter exploration to get the contigs you want out of it; I would recommend looking at the exp_cov parameter of velvetg in particular.

To decode your results, make use of the script. It’s usage is:

python contigs.fa > contigs.text

Which you can then look at with less:

less contigs.text

Further, you might want to make use of an ipython notebook which plots k-mer abundance distributions, and could help you with parameters:

curl > /usr/local/notebooks/2013-caltech-scavenger-hunt.ipynb

Which you can then access by going to https://ec2-???????????

Happy hunting!

comments powered by Disqus