Package: reallyrun
ClaretBio's REALLY Library Processing SoftwareThis software is for the basic informatic processing of sequencing data generated using ClaretBio'sREALLYlibrary prep kit with or without using unique molecular identifiers (UMIs).InstallationThis software can be installed as a python package using the commandpip install reallyrunUsageThe basic analysis can be run withreally runsampleswhen running on standard libraries orreally runsamples --umiwhen running on libraries with UMIs. The software takes in raw fastqs and trims adapters, aligns to a user-specified reference transcriptome, and marks duplicates. For UMI aware demltiplexing of REALLY libraries please use our SRSLYumi python package (more info athttps://github.com/claretbio/SRSLYumi)In order to run, this software requires an installation of conda. For speed, we recommendmambawhich is best installed frommambaforge. If you prefer to use standard conda, installation instructions can be foundhere.Required Arguments--starindex : a path to the STAR index if one exists

OR

--reference: a path to the reference genome you wish to have converted to a STAR transcriptome for alignment

--gtf: a path to the GTF of genes that corresponds to the reference you are aligning to

--refflat: a path to the refFlat of genes that correspond to the reference you are aligning to

--ribosomal: a path to the rRNA interval list that corresponds to the reference you are aligning to

--libraries or --libfile: the library IDs you would like analzed in comma separated format or the path to a file with one ID per line, repsectivelyOptional Arguments--fastqdir : a path to the directory containing the raw fastqs you wish to process (if not specified, defaults to current working directory)

--resultsdir : a path to the directory you would like the output to be in (if not specified, defaults to current working directory)

--indexdir: a path to the directory where you would like the STAR index to be created (if not specified, defalts to current working directory). Use one of --starindex or --indexdir, not bothHelpful information about argumentsThe library IDs provided should match the beginning of the fastq files. For example, the library ID for the fastq files namedlib1_R1.fastq.gzandlib1_R2.fastq.gzwould belib1. This can be provided directly on the command line with a comma separated list:--libraries lib1,lib2or as a file that lists one library ID per line:--libfile libfile.txt.Details on GTF format can be foundherePicard'sBedToIntervalList toolcan be used to generate the ribosomal interval list. An example interval list can be foundhere.For an example refFlat file for GRCh38, seerefFlat.txt.gzathttps://hgdownload.cse.ucsc.edu/goldenPath/hg38/database/.Example Commandreally runsamples --fastqdir /home/user/fastqfiles \
--resultsdir /home/user/amazing-results \
--reference /home/user/data/hg38.fa \
--gtf /home/user/data/hg38.gtf \
--refflat /home/user/data/hg38_refflat.txt \
--ribosomal /home/user/data/hg38_rrna.interval_list \
--libraries lib1,lib2,lib3For reproducibility's sake and to ensure appropriate versions we use snakemake wrappers for many of the tools in this pipeline, which are often slow to create the first time they are used. As a result, your first time running the software may take a long time - don't worry, this is totally normal!