Tutorial: Orchestrating a small, parallel, RNA-seq pre-processing workflow using R
gravatar for Sean Davis
4.3 years ago by
Sean Davis21k
United States
Sean Davis21k wrote:


In this little workflow, we will be using a relatively new technology, pseudoalignment and quantification to deal with RNA-seq data from eight samples. The technical steps are:

  1. Use the SRA SDK to download FASTQ files for each sample
  2. Build a transcriptome index for Kallisto
  3. Pseudoalignment and quantification with Kallisto
  4. Read Kallisto output into a SummarizedExperiment object

Technical skills being showcased include:

  1. Accessing data from SRA using R/SRA SDK
  2. Use of R system() functionality to orchestrate workflows involving command-line programs
  3. Parallel processing with BiocParallel
  4. RNA-seq processing with Kallisto
ADD COMMENTlink modified 4.3 years ago • written 4.3 years ago by Sean Davis21k

Great tutorial. I found that the "abundance.txt" file referred to before the call to runKalllisto() is actually called "abundance.tsv". After making that change everything runs smoothly.

ADD REPLYlink written 4.3 years ago by Diego Diez730
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 337 users visited in the last hour