This repository has been archived by the owner on May 3, 2024. It is now read-only.
-
Notifications
You must be signed in to change notification settings - Fork 104
Targeted Iso Seq QC Wiki
Elizabeth Tseng edited this page Oct 2, 2019
·
11 revisions
Last Updated: 10/02/2019
This wiki shows scripts used for targeted Iso-Seq analysis.
No installation is required. You just download the Cupcake repository from GitHub and add the appropriate paths.
$ git clone https://github.com/Magdoll/cDNA_Cupcake.git
$ export PYTHONPATH=$PYTHONPATH:<path_to>/cDNA_Cupcake/sequence
To see that you have successfully gotten BioPython, bx-python, and the sequence/
subdir working, try the following in Python interpreter:
$ python
Python 3.7.4 (default, Aug 13 2019, 20:35:49)
>>> import Bio # testing BioPython
>>> import bx # testing bx-python
>>> import BED # testing Cupcake sequence/
All of the above should work without error messages.
You will need:
- The input fasta file
- The aligned SAM file (see GMAP command below)
- A tab-delimited probe BED file with at least the first three columns and an optional 4th column of
<chrom> <start> <end> <gene_name>
. No headers allowed!
You can generate the aligned SAM file using minimap2:
minimap2 -ax splice -t 30 -uf --secondary=no \
hg38.fa input.fasta > input.fasta.sam
The script usage is:
python calc_probe_hit_from_sam.py [BED] [FASTA] [SAM]
--start_base {0 or 1}
--end_base {0 or 1}
-o [OUTPUT]
Where --start_base, --end_base
indicates whether the start/end index is 0-based or 1-based.
For example:
python <path_to>/cDNA_Cupcake/targeted/calc_probe_hit_from_sam.py \
my-probes.bed input.fasta input.fasta.sam \
--start_base 0 --end_base 0 \
-o input.fasta.sam.probe_hit.txt
The output file will have the following fields:
read_id : the read ID
read_len : length of read
is_fl : by default this is set to 'Y' for now
num_probe : number of distinct probes overlapped
num_base_overlap: total base of the read that overlaps with a probe
loci : mapped start-end of the read on the genome
genes : the 4th column of the probe BED file, if given