Cutadapt All Samples Exercise
15 Minutes
Now that we’ve learned the basics of running Cutadapt, we need to
trim all the rest of our samples. If you remember from the Computational
Foundations course, we learned about using bash variables. Let’s try an
exercise where we use a bash variable to trim each one of our FASTQ
files.
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 Cutadapt’s help page and choose the proper arguments for our
Cutadapt command(s).
- Use a bash variable along with Cutadapt to trim all remaining FASTQ
files.
- Confirm that we have all of our expected output files.
Hint: Using a bash variable allows us to quickly change some
arguments in a repeated command, e.g. :
noun="World"
echo "Hello, $noun!"
noun="Class"
echo "Hello, $noun!"
LS0tCnRpdGxlOiAiRXhlcmNpc2UgMDEiCmF1dGhvcjogIlVNIEJpb2luZm9ybWF0aWNzIENvcmUiCm91dHB1dDoKICAgICAgICBodG1sX2RvY3VtZW50OgogICAgICAgICAgICBpbmNsdWRlczoKICAgICAgICAgICAgICAgIGluX2hlYWRlcjogaGVhZGVyLmh0bWwKICAgICAgICAgICAgdGhlbWU6IHBhcGVyCiAgICAgICAgICAgIGZpZ19jYXB0aW9uOiB0cnVlCiAgICAgICAgICAgIG1hcmtkb3duOiBHRk0KICAgICAgICAgICAgY29kZV9kb3dubG9hZDogdHJ1ZQotLS0KPHN0eWxlIHR5cGU9InRleHQvY3NzIj4KYm9keXsgLyogTm9ybWFsICAqLwogICAgICBmb250LXNpemU6IDE0cHQ7CiAgfQpwcmUgewogIGZvbnQtc2l6ZTogMTJwdAp9Cjwvc3R5bGU+Cgo8YnI+CgojIyBDdXRhZGFwdCBBbGwgU2FtcGxlcyBFeGVyY2lzZQoKPGJyPgoKKioxNSBNaW51dGVzKioKCjxicj4KCk5vdyB0aGF0IHdlJ3ZlIGxlYXJuZWQgdGhlIGJhc2ljcyBvZiBydW5uaW5nIEN1dGFkYXB0LCB3ZSBuZWVkIHRvIHRyaW0gYWxsIHRoZSByZXN0IG9mIG91ciBzYW1wbGVzLiBJZiB5b3UgcmVtZW1iZXIgZnJvbSB0aGUgQ29tcHV0YXRpb25hbCBGb3VuZGF0aW9ucyBjb3Vyc2UsIHdlIGxlYXJuZWQgYWJvdXQgdXNpbmcgYmFzaCB2YXJpYWJsZXMuIExldCdzIHRyeSBhbiBleGVyY2lzZSB3aGVyZSB3ZSB1c2UgYSBiYXNoIHZhcmlhYmxlIHRvIHRyaW0gZWFjaCBvbmUgb2Ygb3VyIEZBU1RRIGZpbGVzLgoKPGJyPgoKIyMjIEluc3RydWN0aW9uczoKCjxicj4KCi0gV29yayBpbmRlcGVuZGVudGx5IGluIHRoZSBtYWluIHJvb20sIHBvc3RpbmcgYW55IHF1ZXN0aW9ucyB0aGF0IGFyaXNlIHRvIHNsYWNrLgotIFJlY29tbWVuZGF0aW9ucyBmb3Igd3JpdGluZyB5b3VyIG93biBjb2RlOgogIC0gUmVhZCBmdW5jdGlvbiBkb2N1bWVudGF0aW9uCiAgLSBUZXN0IG91dCBpZGVhcyAtIGl0J3Mgb2theSB0byBtYWtlIG1pc3Rha2VzIGFuZCBnZW5lcmF0ZSBlcnJvcnMKICAtIFVzZSBhIHNlYXJjaCBlbmdpbmUgdG8gbG9vayB1cCBlcnJvcnMgb3IgcmVjb21tZW5kZWQgc29sdXRpb25zIHVzaW5nIGtleXdvcmRzCi0gV2UnbGwgcmV2aWV3IHBvc3NpYmxlIHNvbHV0aW9ucyBhZnRlciB0aW1lIGlzIHVwIGFzIGEgZ3JvdXAuCgo8YnI+CgotIFJldmlldyBDdXRhZGFwdCdzIGhlbHAgcGFnZSBhbmQgY2hvb3NlIHRoZSBwcm9wZXIgYXJndW1lbnRzIGZvciBvdXIgQ3V0YWRhcHQgY29tbWFuZChzKS4KLSBVc2UgYSBiYXNoIHZhcmlhYmxlIGFsb25nIHdpdGggQ3V0YWRhcHQgdG8gdHJpbSBhbGwgcmVtYWluaW5nIEZBU1RRIGZpbGVzLgotIENvbmZpcm0gdGhhdCB3ZSBoYXZlIGFsbCBvZiBvdXIgZXhwZWN0ZWQgb3V0cHV0IGZpbGVzLgoKPGJyPgoKPiBIaW50OiBVc2luZyBhIGJhc2ggdmFyaWFibGUgYWxsb3dzIHVzIHRvIHF1aWNrbHkgY2hhbmdlIHNvbWUgYXJndW1lbnRzIGluIGEgcmVwZWF0ZWQgY29tbWFuZCwgZS5nLiA6Cj4KPiB+fn4KPiBub3VuPSJXb3JsZCIKPiBlY2hvICJIZWxsbywgJG5vdW4hIgo+IG5vdW49IkNsYXNzIgo+IGVjaG8gIkhlbGxvLCAkbm91biEiCj4gfn5+Cg==