Provided by: crac_2.5.2+dfsg-6build1_amd64 
NAME
crac v. 2.5.2 - Crac is a tool to analyse RNA-Seq data provided by NGS.
SYNOPSIS
crac [ options ] -i <index_file> -r <reads_file1> [reads_file2] -k <int> -o <output_file>
crac -h|--help
crac -f|--full-help
crac -v|--version
DESCRIPTION
crac CRAC: an integrated approach to the analysis of RNA-seq reads
Whatever the biological questions it addresses, each RNA-seq analysis requires a
computational prediction of either small scale mutations, indels, splice junctions or
fusion RNAs. This prediction is currently performed using complex pipelines involving
multiple tools for mapping, coverage computation, and prediction at distinct steps. We
propose a novel way of analyzing reads that integrates genomic locations and local
coverage, and delivers all above mentioned predictions in a single step. Our program,
CRAC, uses a double k-mer profiling approach to detect candidate mutations, indels, splice
or fusion junctions in each single read. Compared to existing tools, CRAC provides state
of the art sensitivity and improved precision for all types of predictions, yielding high
rates of true positive candidates (99.5% for splice junctions). When applied to four
breast cancer libraries, CRAC recovered 74% of validated fusion RNAs and predicted
reccurrent fusion junctions that were overseen in previous studies. Importantly, CRAC
improves its predictive performance when supplied with e.g. 200 nt reads and should fit
future needs of read analyses.
SPECIAL OPTIONS
crac As a lot of softwares, there are many optional parameters in CRAC but only three are
mandatory. This document is intended to guide users of CRAC to choose the more appropriate
parameters according to their needs:
-h, --help
to print the principal help page of CRAC
-f, --full-help
to print the complete help page of CRAC
-v, --version
to print version of CRAC
USUAL ARGUMENTS
Mandatory options
All these flags must be set.
-i <index_file>
is the name of the index previously built with the crac-index binary file. Note
that crac-index construct the structure <index_file.ssa> with its configuration
<index_file.conf> so only the prefix <index_file> must be specified (without
extension) to consider the structure and the configuration files both in CRAC
-r <reads_file1> [reads_file2]
is the source(s) of the FASTA or FASTQ file(s) containing the reads. Note that the
number of files depends if single or paired-reads. The input file may also be
compressed using gzip
-k <int>
is the length of the k-mer to be used to map the reads on the reference
<index_file>. Note that the condition (k < m) is necessary and reads (or both
paired reads) are ignored if m < k. It must be chosen to ensure (as much as
possible) that a k-mer has a very high probability to occur a single time on the
genome
-o, --sam <output_file>
is the output file in SAM format (see the Documentation of SAM format in CRAC for
more details) or print on STDOUT with "-o -" argument
Optional parameters
--stranded
must be specificied if reads are produced by a stranded protocol of RNA-Seq (not
stranded by default)
--fr/--rf/--ff
set the mates alignement orientation (--rf by default)
-m, --reads-length, -m <int>
must be specified for reads of fixed length. If the read length is fixed, we deeply
recommend you to specify the read length, by using the -m parameter. CRAC will
therefore be much faster. --reads-length <int> is specified for variable or longer
reads, reads shorter are ignored and reads longer are trimmed
--treat-multiple <int>
display alignments with multiple locations (with a fixed limit) rather than a
single alignment per read in the SAM file
--nb-threads <int>
is the number of threads to run crac, computational time is almost divided by the
number of threads (one thread by default)
--max-locs <int>
corresponds to the max number of occurrences retrieved in the index for a given k-
mer: smaller is faster, but with a small value, you may miss some locations that
would help CRAC detecting the right cause
--no-ambiguity <none>
discard biological events (splice, svn, indel, chimera) which have several matches
on the reference index. Indeed, if crac has identified a biological cause in the
read that can match in differents places of the genome we classify this cause as a
biological undetermined event.
Optional output arguments
--gz <none>
all output files specified after this argument are gzipped (included for the sam
file if -o/--sam argument is specified after)
--bam <none>
sam output is encode in binary format(BAM)
--summary <output_file>
save some statistics about mapping and classification
--show-progressbar <none>
show a progress bar for the process times on STDERR
--use-x-in-cigar <none>
use X cigar operator when CRAC identifies a mismatch
Optional output homemade file formats
--all <base_filename>
set output base filename for all causes following. Note that only a base_filename
must be specified. Then, the appropriate file extension is added for each cause
(SNP, chimera, splice, etc) set output base filename for all causes following
--normal <output_file>
save reads that do not contain any break
--almost-normal <output_file>
save reads that do not contain any break but with a variable support
--single <output_file>
save reads which are located in this way: at least --min-percent-single-loc <float>
of k-mers are once located on the reference index
--duplicate <output_file>
save reads which are located in this way: at least --min-percent-duplication-loc
<float> of k-mers are a few times on the reference index (ie. between
--min-duplication <int> and --max-duplication <int> of locations)
--multiple <output_file>
save reads which are located in this way: at least --min-percent-multiple-loc
<float> of k-mers are a many times on the reference index (ie. more than
--max-duplication <int> of locations)
--none <output_file>
save reads which are not located on the reference index
--snv <output_file>
save reads that contain at least a snv
--indel <output_file>
save reads that contain at least a biological indel
--splice <output_file>
save reads that contain at least a splicing junction
--weak-splice <output_file>
save reads that contain at least a low coverage splicing junction
--chimera <output_file>
save reads that contain at least a chimera junction (junction on different
chromosomes, strands or genes)
--paired-end-chimera <output_file>
paired-end-chimera <output_file>= save paired-end reads that contains a chimera in
the non-sequenced part of the original fragment.
--biological <output_file>
save reads that contain a biological cause but for which there is not enough
informations to be more specific. Note that the biological cause is described for
each read
--errors <output_file>
save reads that contain at least a sequence error
--repeat <output_file>
save reads that contain a repeated sequence: at least
--min-percent-repetition-loc <float> percent of k-mers of a given
read are located at least --min-repetition <int> occurrences on the
reference index
--undetermined <output_file>
save reads that contain an undetermined error: some k-mers are not located on the
genome, but the reason for that could not be determined. Note that the error is
described for each read
--nothing <output_file>
save reads that are unclassified
Optional process for specific research
--deep-snv
must be specified to increase sensitivity to find SNVs at the cost of more
computations (only substitution, no indels YET). That process searches for SNV in
border cases reads. Those reads would otherwise be classified in bioundetermined
--stringent-chimera
must be specified to increase accuracy to find chimera junctions in exchange of
sensitivity and computational times
Optional process launcher (once must be selected)
--emt launch an exact matching processing of reads on the index. Either the argument
specified with -k is equal to 0 which means that the entire read is perfectly
mapped on the genome or only a factor of length k per read is mapped (the first one
with a location) and the rest is sofclipped. With this process, reads are not
indexed and it provides a low memory consumption. Note this kind of method is very
useful for DGE reads mapping.
--server
launch a server to query a given read more precisely. That process is useful for
debugging. Note that the output arguments will not be taken into account. Give an
--input-name-server <string> to set the input fifo name (classify.fifo by default)
and give an --output-name-server <string> to set the output fifo name
(classify.out.fifo by default). The server can then be used through a client
crac-client
Additional settings for users
--detailed-sam
more informations are added in SAM output file. See the Documentation of SAM format
in CRAC for more details
--min-percent-single-loc <float>
is, to consider a given read as uniquely mapped, the minimum proportion of k-mers
that are uniquely mapped on the index (0.15 by default)
--min-duplication <int>
is the minimum number of location to consider a duplicated k-mer (2 by default)
--max-duplication <int>
is the maximum number of location to consider a duplicated k-mer (9 by default)
--min-percent-duplication-loc <float>
is, to consider a given read as duplicated, the minimum proportion of k-mers that
are duplicated on the index (0.15 by default)
--min-percent-multiple-loc <float>
is, to consider a given read as “multiple”, the minimum proportion of k-mers that
are multiple mapped on the index (0.50 by default)
--min-repetition <int>
is the minimum number of locations to consider a repeated k-mer (20 by default)
--max-percent-repetition-loc <float>
is, for a given read, the minimum proportion of k-mers that are repeated on the
index to consider a repetition (0.20 by default)
--max-splice-length <int>
is the threshold to consider a splice, ie. a splice is reported if the junction
length is below max-splice-length <int>, a chimera is considered otherwise
(distance by default is 300Kb)
--max-bio-indel <int>
is the threshold to consider a biological indel, ie. an indel is reported if the
gap length is below max-bio-indel, a splice is considered otherwise (distance by
default is 15)
--max-bases-retrieved <int>
is the number of nucleotides to display in outputfile in case of insertion (15 by
default)
--min-support-no-cover <float>
is the minimum coverage to be able to report a biological cause. Note that if a
single read contains a given substitution, it is difficult (if not impossible) to
distinguish a sequence error and a biological cause (1.30 by default)
Additional settings for advanced users
--min-break-length <int>
is the minimal break length (as the percentage of k, the k-mer length) so that a
cause can be reported. Theoretically, for a given cause, the break length is always
>= (kmer_length - 1). Otherwise, the break may be merged with a close enough break,
or the break will be considered as undetermined. (0.5 by default)
--max-bases-randomly-matched <int>
A k-mer overlapping an exon-exon junction, for example, may still match on the
genome if the overlap is at the end of the read (without loss of generality). This
is due to the fact that the nucleotides starting the second exon may be the same as
the nucleotides starting the intron. Theoretically, there is a 0.25 probability
that we have the same nucleotide at the first position of the intron and the exon.
This option specifies how many nucleotides may be matched randomly at most
--max-extension-length <int>
is the maximum number of k-mers extended at each side of a read break. In fact, for
a given break, k-mers with false locations can generate false biological causes, so
the consistency is checked for each side of the break to discard false k-mers and
readjust the good boundaries of the break (10 by default)
--nb-tags-info-stored <int>
is a buffer to store informations for each thread during the computing phase (1000
by default). This value must be increased if threads work below their real
capabilities. With --nb-threads 15, CPU usage must be about 1400%
--reads-index <string>
the reads index data-structure uses by CRAC. Available reads index are: JELLYFISH
and GKARRAYS. (JELLYFISH by default).
--nb-nucleotides-snp-comparison <int>
is the minimum k-mer length tolerated for the deep SNVs search (8 by
default)
--max-number-of-merges <int>
is the maximum number of merges tolerated during the break merge process
for the chimera detection (4 by default)
--min-score-chimera-stringent <float>
is the mimimal score to consider a chimera event
otherwise it is classify as a bioundetermined event (0.6 by default)
SEE ALSO
The full documentation for crac is maintained as a org manual. If the info and crac
programs are properly installed at your site, the command
info crac
should give you access to the complete manual.
AUTHOR
About the crac package.
You can contact Nicolas PHILIPPE, Mikael SALSON, Jerome AUDOUX and Alban MANCHERON by
sending an e-mail to <crac-bugs@lists.gforge.inria.fr>.
Programming:
Nicolas PHILIPPE <nphilippe.research@gmail.com>
Mikaël SALSON <mikael.salson@lifl.fr> Jérome AUDOUX
<jerome.audoux@gmail.com>
with additional contribution for the packaging of:
Alban MANCHERON <alban.mancheron@lirmm.fr>
About the crac publication.
You may cite the following paper if you use our tool:
Gk-arrays: Querying large read collections in main memory: a versatile
data structure
Philippe N., Salson M., Lecroq T., Leonard M., Commes T., Rivals E.
BMC Bioinformatics 2011, 12:242.
Crac: An integrated RNA-Seq read analysis
Philippe N., Salson M., Commes T., Rivals E.
Genome Biology 2013; 14:R30.
2024-04-01 crac(1)