attrna / cdna_cupcake

Miscellaneous collection of Python and R scripts for processing Iso-Seq data

cDNA_Cupcake

Last Updated: 06/29/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.4

Python Requirements

• Python >= 2.7
• Biopython

Note: to use scripts in the ToFU suite you will need additional requirements. See wiki for more details.

How to use this repository

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.

Documentation

Please see wiki for the latest maintained list of scripts.

A brief list of currently listed scripts are:

Annotation and Rarefaction

• 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.

Sequence Manipulation

• get_seq_stats.py: Summarize length distribution of a FASTA/FASTQ file.
• rev_comp.py: Reverse complement a sequence from command line.
• fa2fq.py and fq2fa.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

Sequence Simulation

• simulate.py: Simulate error in sequences.

Cupcake ToFU: supporting scripts for Iso Seq after clustering step

• 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.

ToFU2: a beta version of Iso-Seq2

See ToFU2 white paper about improvements and how to install.

Version Changes

2018.06.29 updated to v5.4. collapse,fusion,abundance,demux now works with isoseq3 output.

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.