Versions Compared


  • This line was added.
  • This line was removed.
  • Formatting was changed.

Rationale and background:


Rascaf: A tool for genome assembly improvement with RNA-seq data.



  1. A CyVerse account. (Register for an CyVerse account here -
  2. Input
     a. Mandatory
         -b bamfile(s):  one or multiple BAM file
         -f path_to_assembly: path to the raw assembly FASTA file
     b. Optional
         -ms minreads: minimum support for connecting two contigs (default: 2)
         -ml minlen:  minimum exonic length if no intron (default: 200)
         -k wordsize: the size of a k-mer (<=32; default: 21)
  3. Output
        -o filename:  prefix for the output file name

Test/Sample Data

The test data for Rascaf can be found here: /iplant/home/shared/iplantcollaborative/example_data/rascaf

Test Run

  1. Open Rascaf in DE
  2. Choose the sample input file from /iplant/home/shared/iplantcollaborative/example_data/rascaf/sample.bam
  3. Select the FASTA genome file from /iplant/home/shared/iplantcollaborative/example_data/rascaf/sample.fa
  4. Optionally set the other arguments, such as -ms, -ml, -k
  5. Optionally set the prefix for output file name
  6. Launch the analysis

Test Results: 

There are at least three output files; if multiple BAM input files are provided, the software will generate one file with extension '.out' for each input file. The following is the output when running Rascaf with the data provided and default parameters:

  2. rascaf_scaffold.fa
  3. sample_0.out