Hi, currently, I am trying to build an accurate transcript model of non-model species.
I have performed CAGE-seq, Iso-seq and strand specific RNA-seq on the same sample.
The problem is, the Iso-seq library kit my service provider used did not employ 5' cap selection.Thus, 5' sequences of my Iso-seq data were incomplete.
Now, I am wondering if CAGE-seq and strand-specific RNA-seq data could complement the Iso-seq data.
However, so far I have not found tools which can integrate CAGE-seq and other seq data.
Is there any tool to combine CAGE-seq data, Iso-seq data and RNAseq data?
Any suggestions would be appreciated.
Thanks, freddiejung
•
written
by
freddiejung •
1 answer
I do not have a full solution but in the CAGEr Bioconductor package you can cluster the CAGE signal and extract the coordinates of the dominant peak for each cluster. These dominant peaks might be used to extend or trim the transcript models that have their 5′ end within an arbitrary threshold.
https://bioconductor.org/packages/release/bioc/vignettes/CAGEr/inst/doc/CAGEexp.html#ctss-clustering
Log in to answer this question.
-
Best practices for improving BRAKER annotation with long-read Iso-Seq and short-read RNA-seq transc…
written by Andrei Makhon •
-
GSEA and enrichment analysis of sequence data
written by Bot_COLA •
-
how to choose library type for single-end 454 sequencing?
written by Ruixuan •
-
augustus Bam2hints program creates an empty .gff file
written by freddiejung •
-
Identify transcription start site using RNA-seq data
written by I0110
-
feature count strand explained
written by popayekid55
-
How to analyze CAGE-Seq data?
written by heir_of_isildur88 •
-
Denovo transcript assembly with data from different library types.
written by GouthamAtla
-
How can CAGE and RNA-Seq data complement each other?
written by espop23
-
Ht-Seq Read Count And Strand-Specificity
written by thecuriousbiologist
