Aligning All Samples Exercise


15 Minutes


We just learned about how to use RSEM & STAR, but now we need to align all of the rest of our samples to the reference genome. In this breakout exercise, we’ll build upon some concepts we’ve learned previously.


Instructions:


  • Work independently in the main room, posting any questions that arise to slack.
  • Recommendations for writing your own code:
    • Read function documentation
    • Test out ideas - it’s okay to make mistakes and generate errors
    • Use a search engine to look up errors or recommended solutions using keywords
  • We’ll review possible solutions after time is up as a group.


  • Review what we’ve learned about running RSEM + STAR, to determine an appropriate command for aligning our samples.
  • Using what we’ve learned previously, create a script using this command to quickly and easily align the rest of our samples.
  • Run the script, view the output, and verify that we have the files we need.


LS0tCnRpdGxlOiAiRXhlcmNpc2UgMDIiCmF1dGhvcjogIlVNIEJpb2luZm9ybWF0aWNzIENvcmUiCm91dHB1dDoKICAgICAgICBodG1sX2RvY3VtZW50OgogICAgICAgICAgICBpbmNsdWRlczoKICAgICAgICAgICAgICAgIGluX2hlYWRlcjogaGVhZGVyLmh0bWwKICAgICAgICAgICAgdGhlbWU6IHBhcGVyCiAgICAgICAgICAgIGZpZ19jYXB0aW9uOiB0cnVlCiAgICAgICAgICAgIG1hcmtkb3duOiBHRk0KICAgICAgICAgICAgY29kZV9kb3dubG9hZDogdHJ1ZQotLS0KPHN0eWxlIHR5cGU9InRleHQvY3NzIj4KYm9keXsgLyogTm9ybWFsICAqLwogICAgICBmb250LXNpemU6IDE0cHQ7CiAgfQpwcmUgewogIGZvbnQtc2l6ZTogMTJwdAp9Cjwvc3R5bGU+Cgo8YnI+CgojIyBBbGlnbmluZyBBbGwgU2FtcGxlcyBFeGVyY2lzZQoKPGJyPgoKKioxNSBNaW51dGVzKioKCjxicj4KCldlIGp1c3QgbGVhcm5lZCBhYm91dCBob3cgdG8gdXNlIFJTRU0gJiBTVEFSLCBidXQgbm93IHdlIG5lZWQgdG8gYWxpZ24gYWxsIG9mIHRoZSByZXN0IG9mIG91ciBzYW1wbGVzIHRvIHRoZSByZWZlcmVuY2UgZ2Vub21lLiBJbiB0aGlzIGJyZWFrb3V0IGV4ZXJjaXNlLCB3ZSdsbCBidWlsZCB1cG9uIHNvbWUgY29uY2VwdHMgd2UndmUgbGVhcm5lZCBwcmV2aW91c2x5LgoKPGJyPgoKIyMjIEluc3RydWN0aW9uczoKCjxicj4KCi0gV29yayBpbmRlcGVuZGVudGx5IGluIHRoZSBtYWluIHJvb20sIHBvc3RpbmcgYW55IHF1ZXN0aW9ucyB0aGF0IGFyaXNlIHRvIHNsYWNrLgotIFJlY29tbWVuZGF0aW9ucyBmb3Igd3JpdGluZyB5b3VyIG93biBjb2RlOgogIC0gUmVhZCBmdW5jdGlvbiBkb2N1bWVudGF0aW9uCiAgLSBUZXN0IG91dCBpZGVhcyAtIGl0J3Mgb2theSB0byBtYWtlIG1pc3Rha2VzIGFuZCBnZW5lcmF0ZSBlcnJvcnMKICAtIFVzZSBhIHNlYXJjaCBlbmdpbmUgdG8gbG9vayB1cCBlcnJvcnMgb3IgcmVjb21tZW5kZWQgc29sdXRpb25zIHVzaW5nIGtleXdvcmRzCi0gV2UnbGwgcmV2aWV3IHBvc3NpYmxlIHNvbHV0aW9ucyBhZnRlciB0aW1lIGlzIHVwIGFzIGEgZ3JvdXAuCgo8YnI+CgotIFJldmlldyB3aGF0IHdlJ3ZlIGxlYXJuZWQgYWJvdXQgcnVubmluZyBSU0VNICsgU1RBUiwgdG8gZGV0ZXJtaW5lIGFuIGFwcHJvcHJpYXRlIGNvbW1hbmQgZm9yIGFsaWduaW5nIG91ciBzYW1wbGVzLgotIFVzaW5nIHdoYXQgd2UndmUgbGVhcm5lZCBwcmV2aW91c2x5LCBjcmVhdGUgYSBzY3JpcHQgdXNpbmcgdGhpcyBjb21tYW5kIHRvIHF1aWNrbHkgYW5kIGVhc2lseSBhbGlnbiB0aGUgcmVzdCBvZiBvdXIgc2FtcGxlcy4KLSBSdW4gdGhlIHNjcmlwdCwgdmlldyB0aGUgb3V0cHV0LCBhbmQgdmVyaWZ5IHRoYXQgd2UgaGF2ZSB0aGUgZmlsZXMgd2UgbmVlZC4KCjxicj4K