The RF Map module can process any number of FastQ files, both from single-read or paired-end experiments. Reads are first pre-processed (trimmed and clipped), and mapped to the reference transcriptome.
The resulting SAM/BAM files can be then passed to the RF Count module.
Usage
$ rf-map [options] file1.fastq ... filen.fastq.gz
$ rf-map [options] file1_R1.fastq,file1_R2.fastq ... filen_R1.fastq.gz,filen_R2.fastq.gz
To list the required parameters, simply type:
$ rf-map -h
Parameter | Type | Description |
---|---|---|
-b2 or --bowtie2 | Uses Bowtie v2 for reads mapping (Default: Bowtie v1) | |
-p or --processors | int | Number of processors (threads) to use (Default: 1) |
-wt or --working-threads | int | Number of working threads to use for each instance of SAMTools/Bowtie (Default: 1). Note: RT Counter executes 1 instance of SAMTools/Bowtie for each processor specified by -p . At least -p <processors> * -wt <threads> processors are required. |
-o or --output-dir | string | Output directory for writing mapped reads in SAM/BAM format (Default: rf_map/) |
-ow or --overwrite | Overwrites the output directory if already exists | |
-t or --tmp-dir | string | Path to a directory for temporary files creation (Default: ) Note: If the provided directory does not exist, it will be created |
-nb or --no-bam | Disables conversion of SAM files to BAM format | |
-b or --bowtie | string | Path to bowtie (or bowtie2 ) executable (Default: assumes bowtie /bowtie2 is in PATH) |
-c or --cutadapt | string | Path to cutadapt executable (Default: assumes cutadapt is in PATH) |
-s or --samtools | string | Path to samtools executable (Default: assumes samtools is in PATH) |
-kl or --keep-logs | Disables logs folder deletion (mostly for debugging purposes) | |
Cutadapt options | ||
-ca5 or --cutadapt-5adapter | string | Sequence of 5' adapter to clip (Default: CAAGTCTCAAGATGTCAGGCTGCTAG, Illumina/NEBNext Small RNA 5’ Adapter) Note #1: Sequence of 5' adapter will be automatically reverse-complemented Note #2: Multiple adapter sequences can be provided as a comma-separated list |
-ca3 or --cutadapt-3adapter | string | Sequence of 3' adapter to clip (Default: AGATCGGAAGAGCACACGTCT, NEBNext Small RNA 3’ Adapter) Note: Multiple adapter sequences can be provided as a comma-separated list |
-cq5 or --cutadapt-5quality | int | Quality threshold for trimming bases from read 5'-ends (Phred+33, Default: 0 [no trimming]) Note: 5'-end quality trimming must be avoided when analyzing data from RT-stop-based methods |
-cq3 or --cutadapt-3quality | int | Quality threshold for trimming bases from read 3'-ends (Phred+33, Default: 20) |
-cqo or --cutadapt-quality-only | Disables adapters clipping (only performs quality-based trimming) | |
-cl or --cutadapt-len | int | Minimum length to keep reads after clipping (≥10, Default: 25) |
-cm or --cutadapt-min-align | int | Minimum alignment in nt to adapter’s sequence (>0, Default: 1) |
-ctn or --cutadapt-trim-N | Trims Ns at the end of the read | |
-cmn or --cutadapt-max-N | float | Discards reads with more than this number of Ns (Default: off) Note: if a value between 0 and 1 is provided, this is interpreted as a fraction of read's length |
-cp or --clipped | Assumes that reads have been already clipped | |
Mapping options | ||
-mp or --mapping-params | string | Manually specify additional aligner parameters (e.g. -mp "-n 2 -l 15" )Note: for a complete list of aligner's parameters, please check the aligner's documentation |
-mo or --manual-only | Only uses manually specified aligner's parameters. Any other parameter, except -bi (or --bowtie-index ), will be ignored |
|
-bk or --bowtie-k | int | Reports up to this number of mapping positions for reads (Default: disabled) |
-ba or --bowtie-all | Reports all mapping positions for reads (Default: disabled) | |
-bnr or --bowtie-norc | Maps only to transcript's sense strand (Default: both strands) | |
-b5 or --bowtie-trim5 | int | Number of bases to trim from 5'-end of reads (≥0, Default: 0) |
-b3 or --bowtie-trim3 | int | Number of bases to trim from 3'-end of reads (≥0, Default: 0) |
-bi or --bowtie-index | string | Path to transcriptome reference index (see rf-index ) |
Bowtie v1 options | ||
-bl or --bowtie-seedlen | int | Seed length (≥5, Default: 28) |
-bn or --bowtie-n | int | Use Bowtie mapper in -n mode (0-3, Default: 2) Note: in -n mode, Bowtie admits no more than -bn mismatches in the seed |
-bv or --bowtie-v | int | Use Bowtie mapper in -v mode (0-3, Default: disabled) Note: in -v mode, Bowtie admits no more than -bv mismatches in the entire read (Phred quality ignored) |
-bm or --bowtie-max | int | Discard read if more than this number of alignments exist (Default: 1) |
-bc or --bowtie-chunkmbs | int | Maximum MB of RAM for best-first search frames (Default: 128) |
Bowtie v2 options | ||
-bl or --bowtie-seedlen | int | Seed length (3 ≤ l ≤ 32, Default: 22) |
-bN or --bowtie-N | int | Bowtie seed mismatches (0-1, Default: 0) |
-bD or --bowtie-D | int | Maximum number of seed extension attempts (≥0, Default: 15) |
-bR or --bowtie-R | int | Maximum number of re-seeding attempts for reads with repetitive seeds (≥0, Default: 2) |
-bmp or --bowtie-mp | int[,int] | Maximum and minimum mismatch penalities (≥0, Default: 6,2) |
-bdp or --bowtie-dpad | int | Number of extra reference bases included on sides of the DP table (≥0, Default: 15) |
-bdg or --bowtie-rdg | int[,int] | Read's gap open and extend penalities (≥0, Default: 5,3) |
-bfg or --bowtie-rfg | int[,int] | Reference's gap open and extend penalities (≥0, Default: 5,3) |
-bs or --bowtie-softclip | Enables local alignment mode (Default: entire read must align) | |
-bma or --bowtie-ma | int | Match bonus in local alignment mode (Default: 2) |
-bd or --bowtie-dovetail | Paired-end dovetailed reads are allowed (mates extend past each other) |
Important
When using Bowtie v1, Bowtie's --best
and --strata
parameters are automatically added. Please check Bowtie v1 documentation for additional information.
Important
When using Bowtie v2 with paired-end reads, Bowtie's --no-mixed
parameter is automatically added to discard those reads for which only one of the two mates can be mapped. Please check Bowtie v2 documentation for additional information.