Index of /public/dataset/redwood2020/isoseq
Name Last modified Size Description
Parent Directory -
Final-MappedTranscripts/ 2021-01-14 07:52 -
Final-UnmappedTranscripts/ 2021-01-08 14:28 -
Intermediate-FullLengthReads/ 2021-01-08 14:27 -
README.txt 2021-01-13 21:46 7.4K
README (Last Updated 01/13/2020)
Edited by: Elizabeth Tseng (email@example.com)
IMPORTANT: Please note that this release of Iso-Seq data maps to the
Hifiasm v12 version of the genome at
This README file describes the contents in this directory.
This dataset contains intermediate and processed files for a
Redwood Iso-Seq® dataset. The library was sequenced on the Sequel® II System
and processed using SMRTLink v10.1 followed by community tool analysis.
For moreinformation on Iso-Seq® methods, bioinformatics analysis,
see the PacBio® Iso-Seq GitHub and additional references below.
Needles were collected from the same redwood tree as used for the genome sequencing,
flash frozen and stored at -80C.
Iso-Seq® Express Template Preparation for Sequel® and Sequel® II Systems
Sequel II System with Sequel II Binding Kit 2.1 and Sequel II Sequencing Kit 2.0
SMRT Link v10.1 "IsoSeq" protocol, followed by mapping to PacBio's redwood
reference genome of the same species
and collapsed into non-redundant transcript set using Cupcake.
The mapping was done using gmapl which can handle large complex genomes:
gmapl -n 0 -t 60 --max-intronlength-ends 200000 --max-intronlength-middle 200000 --cross-species -z sense_force
After mapping, the Cupcake collapse script was run using a cutoff of
95% alignment coverage and 90% alignment identity. The finalfiltered Iso-Seq
HQ transcripts above this threshold are in the "Final-MappedTranscripts"
subdirectory in the data release.
Iso-Seq transcripts that were unmapped or poorly mapped are additionally run
through Cogent analysis and put in the "Final-UnMmappedTranscripts" subdirectory.
Additionally, BLASTN is run for both mapped and unmapped transcripts
against NT database using E-value cutoff of 0.1.
Mapped transcripts were further phased using IsoPhase as described in  and .
The `run.sh` script in the subdirectories contain detailed run parameters.
Final Mapped Transcripts: Mapped, Collapsed, and Phased
Users wishing to immediately make use the processed, mapped,
results should use the following files:
│ ├── redwood_isoseq.IsoPhase.by_loci_info.txt
│ ├── redwood_isoseq.IsoPhase.by_loci_results.tar.gz
│ ├── redwood_isoseq.IsoPhase.summarized_loci_results.txt
│ ├── redwood_isoseq.IsoPhase.use_partial.cleaned.vcf
│ ├── redwood_isoseq.IsoPhase.use_partial.pre_cleaned.vcf
│ ├── run_phasing_in_dir.sh
│ └── run.sh
Where the companion reference genome is:
(Internal version number: redwood_v12)
Phasing is run using IsoPhase as described in  and , with a minimum required
coverage per locus of 40 FLNC reads. Only SNPs are called (no indel calls).
The "cleaning" (or error correction) step IsoPhase can potentially reduce the number of true
alleles called given the high ploidy (n=6) and uneven allelic isoform expressions of this species.
As such, we have provided both the pre-cleaned and cleaned VCF file. A lack of SNPs called
for a locus does not necessarily indicate lack of variants, rather could also be due to lack of
Final Unmapped Transcripts: Additional Cogent Analysis
For potentially missing or poorly assembled genes
Users wishing to look at unmapped (or poorly mapped) transcripts
that were excluded from the "Final-MappedTranscripts" bin, can use
the files in the directory below.
The file "redwood_isoseq.unmapped_or_badmapped.fasta" are the
individual Iso-Seq HQ transcripts.
The file "redwood_isoseq.unmapped_or_badmapped.cogent_reconstructed_contigs.fasta"
are the result of further running Cogent analysis  where gene family information
between the HQ transcripts were inferred based on k-mer similarities and transcripts
thought to be from the same gene family and were reconstructed using Cogent. Not all
unmapped HQ transcripts have a Cogent output. For those that do, the reconstructed
contigs can be used to visualize alternative splicing of the isoforms. Further, these
reconstructed contigs can be used to supplement the reference genome for missing genes
or mis-assembled genes.
For an example of Cogent applied to ref genome QC, see Warr et al. 
Intermediate: Full-Length, Non-Concatemer (FLNC) Reads
We DO NOT recommend most users re-analyzing from intermediate (FLNC) data.
 PacBio Iso-Seq Landing Page: https://www.pacb.com/applications/rna-sequencing/
 PacBio Iso-Seq GitHub Wiki: https://github.com/PacificBiosciences/IsoSeq_SA3nUP
 Community Tool Cupcake: https://github.com/Magdoll/cDNA_Cupcake
 Community Tool Cogent: https://github.com/Magdoll/Cogent
 Community Tool IsoPhase: https://github.com/Magdoll/cDNA_Cupcake/wiki/IsoPhase:-Haplotyping-using-Iso-Seq-data
 Wang et al. "Variant phasing and haplotypic expression from long-read sequencing in maize", Communications Biology (2020)
 Warr et al. "An improved pig reference genome sequence to enable pig genetics and genomics research ", GigaScience (2020)
For Research Use Only. Not for use in diagnostic procedures. Copyright 2021,
Pacific Biosciences of California, Inc. All rights reserved. The data provided in
these files is subject to change without notice and Pacific Biosciences assumes no
responsibility for any errors or omissions. Certain notices, terms, conditions and/or
use restrictions may pertain to your use of Pacific Biosciences data, products and/or
third party products. Please refer to the applicable Pacific Biosciences Terms and
Conditions of Sale and to the applicable license terms at