Last Updated: 04/12/2018
cDNA_Cupcake is a miscellaneous collection of Python and R scripts used for analyzing sequencing data. Most of the scripts only require Biopython. For scripts that require additional libraries, it will be specified in documentation.
Current version: 5.3
- Python >= 2.7
- Biopython
Note: to use scripts in the ToFU suite you will need additional requirements. See wiki for more details.
Since most of the scripts are independent (do not depend on each other), you can either clone the whole directory, or, if you are only interested in a specific script, just download that specific script to your local drive.
You can clone the GitHub repository, then add the GitHub repo path to your $PATH
variable. The scripts are organized into different sub-directories (ex: sequence/
, rarefaction/
etc) so you will have to add them individually.
git clone https://github.com/Magdoll/cDNA_Cupcake.git
export PATH=$PATH:<path_to_Cupcake>/sequence/
export PATH=$PATH:<path_to_Cupcake>/rarefaction/
For any issues or bugs, please report to Issues.
Please see wiki for the latest maintained list of scripts.
A brief list of currently listed scripts are:
parse_matchAnnot.py
: Parse matchAnnot results into summary format.make_file_for_subsampling_from_collapsed.py
: Prepare file for running subsampling (rarefaction curve).subsample.py
: Running subsamping. Results can be plotted with Excel graphics and R, etc.
get_seq_stats.py
: Summarize length distribution of a FASTA/FASTQ file.rev_comp.py
: Reverse complement a sequence from command line.fa2fq.py
andfq2fa.py
: Convert between FASTA and FASTQ format.sort_fasta_by_len.py
: sort fasta file by length (increasing or decreasing).get_seqs_from_list.py
: extract list of sequences given a fasta file and a list of IDs.err_correct_w_genome.py
: generate fasta sequences given genom
simulate.py
: Simulate error in sequences.
collapse_isoforms_by_sam.py
: Collapse HQ isoform results to unique isoforms (based on genome alignment).get_abundance_post_collapse.py
: Obtain count information post collapse to unique isoforms.filter_by_count.py
: Filter collapse result by FL count information.filter_away_subset.py
: Filter away 5' degraded isoforms.chain_samples.py
: Chaining together multiple samples.fusion_finder.py
: Finding fusion genes.
See ToFU2 white paper about improvements and how to install.
2018.03.29 updated to v5.3. Update to work with pitchfork SA5.1
2018.03.12 updated to v5.2. Fixed over-collapsing genes in collapse script. Now processing strands separately in correct manner.
2017.11.06 updated to v4.1. pCS merge incorrect in chain_samples.py
. Fixed.
2017.10.31 updated to v4.0. pCS merge incorrect in run_preCluster.py
. Fixed.
2017.10.10 updated to v3.9. Merged pCS branch (--dun_use_partial
) and cdunn's random seed.
2017.09.25 updated to v3.7. Fixed minor printing error in scrubbed.group.txt
for scrub_sample_GFF_junctions.py
.