Skip to content


Switch branches/tags

Name already in use

A tag already exists with the provided branch name. Many Git commands accept both tag and branch names, so creating this branch may cause unexpected behavior. Are you sure you want to create this branch?

Latest commit


Git stats


Failed to load latest commit information.
Latest commit message
Commit time

Kallisto (v2.0)

This is an implementation of Kallisto, intended to be used for RNA quantitation against a human index, Gencode release 37 (GRCh38.p13), or a mouse index, Gencode release M26 (GRCm39).

Module Author: Edwin Juarez, Thorin Tabor, Edwin Huang


Algorithm Version: Kallisto 0.46.1


Used verion 0.46.1, which is the most recent as of 2021-06-08

Technical notes:



  • Default: [blank]
  • Required: Yes
  • Description: The fastq or fastq.gz files to be quantitated and translated into RNA seq transcript and gene counts


  • Default: Human
  • Required: Yes
  • Description: Which transcriptome to use for the pseudoalignment (Human or Mouse).


  • Default: YES
  • Required: YES
  • Description: Perform sequence based bias correction


  • Default: NO
  • Required: Yes
  • Description: Include reads where unobserved rest of fragment is predicted to lie outside a transcript


  • Default: not set
  • Required: Yes
  • Description: For strand specific reads, first read forward, first read reverse or unspecified


  • Default: expression_data
  • Required: Yes
  • Description: The base name of the files to be created


  • Default: 2
  • Required: Yes
  • Description: Number of bootstrap samples


  • Default: 42
  • Required: Yes
  • Description: Seed for the bootstrap sampling


  • Default: NO
  • Required: Yes
  • Description: Quantify single-end reads


  • Default:
  • Required: No
  • Description: Estimated average fragment length (float value)

  • Default:
  • Required: No
  • Description: Estimated standard deviation of fragment length; values are estimated with paired end data, but are required when using single-end reads (float value)


  • Default: NO
  • Required: Yes
  • Description: Save pseudoalignments to transcriptome to a BAM file


  • Default: NO
  • Required: Yes
  • Description: Project pseudoalignments to genome sorted BAM file


  • Default:
  • Required: No
  • Description: GTF file for transcriptome information, required for projecting pseudoalignments to genome sorted BAM file


  • Default:
  • Required: No
  • Description: Tab separated file with chromosome names and lengths, optional for projecting pseudoalignments to genome sorted BAM file, but recommended

Output Files

  • module_log.txt: A list of the commands ran by the module.
  • RNASeq_quant/<output_filename>_abundance.h5: transcript-level abundance. Use this file if you inted to use tximport.
  • RNASeq_quant/<output_filename>_abundance.tsv: transcript-level abundance (same as the h5 file).
  • RNASeq_quant/<output_filename>_normalized_gene_level.csv: normalized gene-level abundance.
  • RNASeq_quant/<output_filename>_raw_estimated_counts.csv: raw gene-level estimated counts.
  • RNASeq_quant/<output_filename>_transcript_level.csv: transcript-level abundance.
  • RNASeq_quant/run_info.json: kallisto's parameters used.
  • stdout.txt: A list of non-essential messages printed by Kallisto, this may be helpful to debug should any errors occur.
  • stderr.txt: If there were errors in completing the job, they'd appear here. Note that the mere presence of a file with this name does not imply that there were errors. Here is an example of the contents of the stderr.txt file from our test run:
[quant] fragment length distribution is truncated gaussian with mean = 200, sd = 30
[index] k-mer length: 31
[index] number of targets: 234,485
[index] number of k-mers: 143,310,200
[index] number of equivalence classes: 1,002,767
[quant] running in single-end mode
[quant] will process file 1: /opt/gpbeta_2/gp_home/users/.cache/uploads/cache/
[quant] finding pseudoalignments for the reads ... done
[quant] processed 50,000 reads, 192 reads pseudoaligned
[   em] quantifying the abundances ... done
[   em] the Expectation-Maximization algorithm ran for 501 rounds
[bstrp] running EM for the bootstrap: 1
[bstrp] running EM for the bootstrap: 2

'gene_mode' is TRUE. Sleuth will do counts aggregation at the gene level for downstream normalization, transformation, and modeling steps, as well as for plotting and results.
reading in kallisto results
dropping unused factor levels
normalizing est_counts
6 targets passed the filter
normalizing tpm
merging in metadata
aggregating by column: gene_name
6 genes passed the filter
summarizing bootstraps
Warning messages:
1: In sleuth_prep(s2c, ~condition, target_mapping = t2g, aggregation_column = "gene_name",  :
  There is only one sample present, but you also provided a model. The model will be set to NULL to prevent downstream errors.
The sample can be viewed using sleuth_live after preparation, but you need more than one sample to run the other aspects of Sleuth.
2: `select_()` is deprecated as of dplyr 0.7.0.
Please use `select()` instead.
This warning is displayed once every 8 hours.
Call `lifecycle::last_warnings()` to see where this warning was generated. 
3: In data.table::melt(bs_df, id.vars = "bootstrap_num", = "target_id",  :
  The melt generic in data.table has been passed a data.frame and will attempt to redirect to the relevant reshape2 method; please note that reshape2 is deprecated, and this redirection is now deprecated as well. To continue using melt methods from reshape2 while both libraries are attached, e.g. melt.list, you can prepend the namespace like reshape2::melt(bs_df). In the next version, this warning will become an error.
4: In data.table::dcast(scaled_bs, sample ~ target_id, value.var = "scaled_reads_per_base") :
  The dcast generic in data.table has been passed a data.frame and will attempt to redirect to the reshape2::dcast; please note that reshape2 is deprecated, and this redirection is now deprecated as well. Please do this redirection yourself like reshape2::dcast(scaled_bs). In the next version, this warning will become an error.
Warning message:
In check_quant_mode(obj, which_units) :
  your sleuth object is in gene mode, but you selected 'est_counts'. Selecting 'scaled_reads_per_base'...


Kallisto itself and this GenePattern module are distributed under a modified BSD license. This module's license is available at

Version Comments Version Release Date Description 2 2021-06-08 Adding Mouse index and upgrading Human index 1 2019-04-22 Initial release of Kallisto