Skip to content

Commit

Permalink
update docs for -b and -C options
Browse files Browse the repository at this point in the history
  • Loading branch information
kcotto committed Jan 10, 2023
1 parent 0562e48 commit 3068563
Show file tree
Hide file tree
Showing 3 changed files with 8 additions and 4 deletions.
6 changes: 4 additions & 2 deletions docs/commands/cis-splice-effects-identify.md
Original file line number Diff line number Diff line change
@@ -1,9 +1,7 @@
[csei]: ../images/csei_examples.png


The `cis-splice-effects identify` command is used to identify splicing misregulation events. This command takes in a list of variants in the VCF format and RNAseq alignments produced with a splice-aware aligner in the BAM format. The tool then proceeds to identify non-canonical splicing junctions near the variant sites.


## Usage

`regtools cis-splice-effects identify [options] variants.vcf alignments.bam ref.fa annotations.gtf`
Expand Down Expand Up @@ -33,6 +31,10 @@ The `cis-splice-effects identify` command is used to identify splicing misregula
| -I | Annotate variants in intronic space within a transcript(not to be used with -i). |
| -E | Annotate variants in exonic space within a transcript(not to be used with -e). |
| -S | Don't skip single exon transcripts. |
| -b | The file containing the barcodes of interest for single cell data. |
| -C | Tells cis-splice-effects identify that you want intron-motif method to take priority when assigning strand. i.e. decide strandedness based on the fasta rather than what is encoded in the alignment file. |

**Note** Both junctions extract and cis-splice-effects identify have an intron-motif method that can be used to determine strandedness of junctions extracted from alignment files. Using this method supercedes any strandedness information that might be encoded in the alignment file. To use this method with cis-splice-effects identify, you can the -C option to tell RegTools that you want the intron-motif method to take priority when assigning strandedness.

## Output

Expand Down
4 changes: 3 additions & 1 deletion docs/commands/junctions-extract.md
Original file line number Diff line number Diff line change
Expand Up @@ -2,7 +2,6 @@

The `junctions extract` command can be used to extract exon-exon junctions from an RNAseq BAM file. The output is a BED file in the BED12 format. We have tested this command with alignments from HISAT2, TopHat2, STAR, kallisto, and minimap2 and by comparing the exon-exon junctions with the `junctions.bed` file produced from TopHat.


## Usage

`regtools junctions extract [options] indexed_alignments.bam`
Expand All @@ -24,6 +23,9 @@ The `junctions extract` command can be used to extract exon-exon junctions from
| -r | Region to extract junctions in. This is specified in the format "chr:start-end". If not specified, junctions are extracted from the entire BAM file.|
| -h | Display help message for this command.|
| -s | Strand specificity of RNA library preparation, where the options XS, use XS tags provided by aligner; RF, first-strand; FR, second-strand. This option is required. If your alignments contain XS tags, these will be used in the "unstranded" mode. If you are unsure, we have created this [table](https://rnabio.org/module-09-appendix/0009/12/01/StrandSettings/) to help.
| -b | The file containing the barcodes of interest for single cell data. |

**Note** Both junctions extract and cis-splice-effects identify have an intron-motif method that can be used to determine strandedness of junctions extracted from alignment files. Using this method supercedes any strandedness information that might be encoded in the alignment file. To use this method with junctions extract, you can add a fasta file at the end of your junctions extract command which will tell RegTools that you want the intron-motif method to take priority when assigning strandedness. e.g. `regtools junctions extract [options] indexed_alignments.bam fasta.fa`

## Output

Expand Down
2 changes: 1 addition & 1 deletion docs/index.md
Original file line number Diff line number Diff line change
Expand Up @@ -27,7 +27,7 @@ make
regtools --help
```

If one wishes to test their installation, we include test data under `test_data`.
If one wishes to test their installation, we include test data under `test_data`.

Here's an example command using that data along with the example output.

Expand Down

0 comments on commit 3068563

Please sign in to comment.