Fast Reference-Guided Scaffolding of Genome Assembly Contigs. RagTag, the successor to RaGOO, is now available here: https://github.com/malonge/RagTag
RagTag, RaGOO's successor, is now available here! Please transition your work here, if possible, as RaGOO will eventually no longer be supported.
Alonge, Michael, et al. "RaGOO: fast and accurate reference-guided scaffolding of draft genomes." Genome biology 20.1 (2019): 1-17.
RaGOO is a tool for coalescing genome assembly contigs into pseudochromosomes via minimap2 alignments to a closely related reference genome. The focus of this tool is on practicality and therefore has the following features:
RaGOO should install on OSX and most standard flavors of Linux. RaGOO depends on Python3 as well as the following packages:
The first two packages will be installed automatically when installing RaGOO. Minimap2 is straightforward to install following the instructions on its website. Place the minimap2 executable in your path, or specify its location with the -m parameter (see below).
Currently, the only way to install RaGOO is from source. Set up a virtualenv if desired, just be sure to make a python3 environment. Then, enter the following command to install RaGOO:
$ python setup.py install
usage: ragoo.py [-h] [-e ] [-gff ] [-m PATH] [-b] [-R ] [-T sr] [-t 3] [-g 100] [-s] [-i 0.2] [-j ] [-C]
order and orient contigs according to minimap2 alignments to a reference (v1.1)
positional arguments: <contigs.fasta> fasta file with contigs to be ordered and oriented <reference.fasta> reference fasta file
optional arguments: -h, --help show this help message and exit -e <exclude.txt> single column text file of reference headers to ignore -gff <annotations.gff> lift-over gff features to chimera-broken contigs -m PATH path to minimap2 executable -b Break chimeric contigs -R <reads.fasta> Turns on misassembly correction. Align provided reads to the contigs to aid misassembly correction. fastq or fasta allowed. Gzipped files allowed. Turns off '-b'. -T sr Type of reads provided by '-R'. 'sr' and 'corr' accepted for short reads and error corrected long reads respectively. -t 3 Number of threads when running minimap. -g 100 Gap size for padding in pseudomolecules. -s Call structural variants -i 0.2 Minimum grouping confidence score needed to be localized. -j <skip.txt> List of contigs to automatically put in chr0. -C Write unplaced contigs individually instead of making a chr0 </skip.txt></reads.fasta></annotations.gff></exclude.txt></reference.fasta></contigs.fasta></reference.fasta></contigs.fasta></skip.txt></reads.fasta></annotations.gff></exclude.txt>
RaGOO will try to be smart and not redo intermediate analysis already done in previous executions of the pipeline. For example, if the Minimap2 alignment files are already present from previous runs, RaGOO will not recreate them. However, RaGOO is not that smart, so be sure to remove any files that you want to replace. To be safe, one can just remove the entire output directory if a new analysis is desired (see "Output Files" below).
Both the assembly and the reference must be in the current workding directory, so please either copy them or create a symbolic link. For example:
$ cd /path/to/current/working/directory $ ln -s /path/to/contigs.fasta $ ln -s /path/to/reference.fasta $ ragoo.py contigs.fasta reference.fasta
All of the output will be in the "ragoo_output" directory. If breaking chimeric contigs and calling SVs, the contents of this output directory is as follows:
ragoo_output/ ├── ctg_alignments ├── groupings ├── orderings ├── pm_alns └── ragoo.fasta
The final pseudomolecules. Any unlocalized contigs are concatenated and placed in "Chr0_RaGOO".
This directory contains the results from chimeric contig breaking. The most notable file here is the [prefix].intra.chimera.broken.fa, as this is the final corrected assembly used for downstream scaffolding. All of the downstream information, such as confidence scores, refers to this assembly, not the orignal assembly.
There is one file per chromosome listing the contigs assigned to that chromosome and their grouping confidence score. Please note that these contigs are not ordered. Also note that if chimeras were corrected, the headers in these files refer to the broken assembly in "chimera_break", and not the original assembly.
There is one file per chromosome showing the ordering, orientation (second column), location confidence scores (third column), and orientation confidence scores (fourth column).
This directory contains all of the structural variant calling results. The final structural variants can be found in assemblyticsout.Assemblyticsstructural_variants.bed. This bed file can be converted to VCF using SURVIVOR, though the last two columns (overlap with gaps) must be removed first. The alignment used to generate these variant calls are also present in this directory in SAM and delta format (pmcontigsagainst_ref.sam and pmcontigsagainst_ref.sam.delta), and can be used as input for external tools.
Contains the results from misassembly correction. It will contain the corrected contigs in fasta format, as well as an updated gff file if provided.