1 Introduction

Carrie and Dr. Taneyhill have some new data. I know nothing about it! So I am excited to read some reviews and learn the experimental design. In the mean time, I am going to submit some jobs on the cluster to get a feeling for the data.

cd preprocessing
module add cyoa/202302
start=$(pwd)
for sample in A* B* C*; do
    echo "Starting $sample."
    cd $sample
    cyoa --method prnaseq --species mm39_109 \
         --input $(/bin/ls | tr '\n' ':' | sed 's/:$//g') \
         --gff_tag gene
    cd $start
done
LS0tCnRpdGxlOiAiUHJlcHJvY2Vzc2luZyBzb21lIG5ldyBzYW1wbGVzISIKYXV0aG9yOiAiYXRiIGFiZWxld0BnbWFpbC5jb20iCmRhdGU6ICJgciBTeXMuRGF0ZSgpYCIKb3V0cHV0OgogIGh0bWxfZG9jdW1lbnQ6CiAgICBjb2RlX2Rvd25sb2FkOiB0cnVlCiAgICBjb2RlX2ZvbGRpbmc6IHNob3cKICAgIGZpZ19jYXB0aW9uOiB0cnVlCiAgICBmaWdfaGVpZ2h0OiA3CiAgICBmaWdfd2lkdGg6IDcKICAgIGhpZ2hsaWdodDogemVuYnVybgogICAga2VlcF9tZDogZmFsc2UKICAgIG1vZGU6IHNlbGZjb250YWluZWQKICAgIG51bWJlcl9zZWN0aW9uczogdHJ1ZQogICAgc2VsZl9jb250YWluZWQ6IHRydWUKICAgIHRoZW1lOiByZWFkYWJsZQogICAgdG9jOiB0cnVlCiAgICB0b2NfZmxvYXQ6CiAgICAgIGNvbGxhcHNlZDogZmFsc2UKICAgICAgc21vb3RoX3Njcm9sbDogZmFsc2UKICBybWRmb3JtYXRzOjpyZWFkdGhlZG93bjoKICAgIGNvZGVfZG93bmxvYWQ6IHRydWUKICAgIGNvZGVfZm9sZGluZzogc2hvdwogICAgZGZfcHJpbnQ6IHBhZ2VkCiAgICBmaWdfY2FwdGlvbjogdHJ1ZQogICAgZmlnX2hlaWdodDogNwogICAgZmlnX3dpZHRoOiA3CiAgICBoaWdobGlnaHQ6IHplbmJ1cm4KICAgIHdpZHRoOiAzMDAKICAgIGtlZXBfbWQ6IGZhbHNlCiAgICBtb2RlOiBzZWxmY29udGFpbmVkCiAgICB0b2NfZmxvYXQ6IHRydWUKICBCaW9jU3R5bGU6Omh0bWxfZG9jdW1lbnQ6CiAgICBjb2RlX2Rvd25sb2FkOiB0cnVlCiAgICBjb2RlX2ZvbGRpbmc6IHNob3cKICAgIGZpZ19jYXB0aW9uOiB0cnVlCiAgICBmaWdfaGVpZ2h0OiA3CiAgICBmaWdfd2lkdGg6IDcKICAgIGhpZ2hsaWdodDogemVuYnVybgogICAga2VlcF9tZDogZmFsc2UKICAgIG1vZGU6IHNlbGZjb250YWluZWQKICAgIHRvY19mbG9hdDogdHJ1ZQotLS0KCjxzdHlsZSB0eXBlPSJ0ZXh0L2NzcyI+CmJvZHksIHRkIHsKICBmb250LXNpemU6IDE2cHg7Cn0KY29kZS5yewogIGZvbnQtc2l6ZTogMTZweDsKfQpwcmUgewogZm9udC1zaXplOiAxNnB4Cn0KYm9keSAubWFpbi1jb250YWluZXIgewogIG1heC13aWR0aDogMTYwMHB4Owp9Cjwvc3R5bGU+CgpgYGB7ciBvcHRpb25zLCBpbmNsdWRlPUZBTFNFfQpsaWJyYXJ5KGhwZ2x0b29scykKbGlicmFyeShyZXRpY3VsYXRlKQp0dCA8LSB0cnkoZGV2dG9vbHM6OmxvYWRfYWxsKCJ+L2hwZ2x0b29scyIpKQprbml0cjo6b3B0c19rbml0JHNldCgKICBwcm9ncmVzcyA9IFRSVUUsIHZlcmJvc2UgPSBUUlVFLCB3aWR0aCA9IDkwLCBlY2hvID0gVFJVRSkKa25pdHI6Om9wdHNfY2h1bmskc2V0KAogIGVycm9yID0gVFJVRSwgZmlnLndpZHRoID0gOCwgZmlnLmhlaWdodCA9IDgsIGZpZy5yZXRpbmEgPSAyLAogIGZpZy5wb3MgPSAidCIsIGZpZy5hbGlnbiA9ICJjZW50ZXIiLCBkcGkgPSBpZiAoa25pdHI6OmlzX2xhdGV4X291dHB1dCgpKSA3MiBlbHNlIDMwMCwKICBvdXQud2lkdGggPSAiMTAwJSIsIGRldiA9ICJwbmciLAogIGRldi5hcmdzID0gbGlzdChwbmcgPSBsaXN0KHR5cGUgPSAiY2Fpcm8tcG5nIikpKQpvbGRfb3B0aW9ucyA8LSBvcHRpb25zKGRpZ2l0cyA9IDQsCiAgICAgICAgICAgICAgICAgICAgICAgc3RyaW5nc0FzRmFjdG9ycyA9IEZBTFNFLAogICAgICAgICAgICAgICAgICAgICAgIGtuaXRyLmR1cGxpY2F0ZS5sYWJlbCA9ICJhbGxvdyIpCmdncGxvdDI6OnRoZW1lX3NldChnZ3Bsb3QyOjp0aGVtZV9idyhiYXNlX3NpemUgPSAxMikpCnZlciA8LSAiMjAyMzA1IgpwcmV2aW91c19maWxlIDwtICIiCnZlciA8LSBmb3JtYXQoU3lzLkRhdGUoKSwgIiVZJW0lZCIpCgojI3RtcCA8LSBzbShsb2FkbWUoZmlsZW5hbWU9cGFzdGUwKGdzdWIocGF0dGVybj0iXFwuUm1kIiwgcmVwbGFjZT0iIiwgeD1wcmV2aW91c19maWxlKSwgIi12IiwgdmVyLCAiLnJkYS54eiIpKSkKcm1kX2ZpbGUgPC0gInByZXByb2Nlc3NpbmcuUm1kIgpgYGAKCiMgSW50cm9kdWN0aW9uCgpDYXJyaWUgYW5kIERyLiBUYW5leWhpbGwgaGF2ZSBzb21lIG5ldyBkYXRhLiAgSSBrbm93IF9ub3RoaW5nXyBhYm91dAppdCEgIFNvIEkgYW0gZXhjaXRlZCB0byByZWFkIHNvbWUgcmV2aWV3cyBhbmQgbGVhcm4gdGhlIGV4cGVyaW1lbnRhbApkZXNpZ24uICBJbiB0aGUgbWVhbiB0aW1lLCBJIGFtIGdvaW5nIHRvIHN1Ym1pdCBzb21lIGpvYnMgb24gdGhlCmNsdXN0ZXIgdG8gZ2V0IGEgZmVlbGluZyBmb3IgdGhlIGRhdGEuCgpgYGB7YmFzaCwgZXZhbD1GQUxTRX0KY2QgcHJlcHJvY2Vzc2luZwptb2R1bGUgYWRkIGN5b2EvMjAyMzAyCnN0YXJ0PSQocHdkKQpmb3Igc2FtcGxlIGluIEEqIEIqIEMqOyBkbwogICAgZWNobyAiU3RhcnRpbmcgJHNhbXBsZS4iCiAgICBjZCAkc2FtcGxlCiAgICBjeW9hIC0tbWV0aG9kIHBybmFzZXEgLS1zcGVjaWVzIG1tMzlfMTA5IFwKICAgICAgICAgLS1pbnB1dCAkKC9iaW4vbHMgfCB0ciAnXG4nICc6JyB8IHNlZCAncy86JC8vZycpIFwKICAgICAgICAgLS1nZmZfdGFnIGdlbmUKICAgIGNkICRzdGFydApkb25lCmBgYAo=