More posts like this
-
quantifying splicing transcripts for nanopore RNAseq
written by jshen 0Hi all, I'm looking for some help in quantifying the number of transcripts in a virus genome from nanopore RNAseq data. I've currently been mapping …
-
Is there a way to distinguish doublets from phagocytes performing efferocytosis in scRNA-seq data?
written by m.zarodniuk 0Hello everyone! Here is an interesting conundrum. It has [been reported][1] that phagocytes may contain so-called “passenger” transcripts that originate from engulfed apoptotic cells, which …
-
Finding assembled transcript sequences from StringTie output
written by Jen 10Hi All, I am currently using StringTie to identify novel transcripts (MSTRG ID transcripts) from RNAseq data in mouse. I have the gene.abund.tab output files …
-
Denovo RNA seq
written by kanwarjag 120I am trying to find which tool or work flow is best/ commonly used for de novo assembly of RNAseq PE data? I want to …
-
Sequential contig mapping
written by rh1234 0Hi, I am using short read data to map against assemblies as part of a bulk segregant analysis. The result varies significantly for certain contigs …
-
What transcript quantification tool would you recommend for 3' prime RNA-seq?
written by roy.granit 89I've searched around and could not find (possibly missed) advice on which transcript quantification tool is well suited for 3' prime RNA-seq lib? I've considered …
-
RNAseq transcript quantification: how to deal with different number of genes in host organism?
written by Lina F 20Hi all, I have RNAseq datasets from two organisms: one is the basic E. coli organism and the other one is the same E. coli, …
-
Aligning reads to genome/transcriptome
written by EVR 61Hi all, I have a RNAseq data and I am interested in only those reads which maps to the exactly between two scaffolds or transcripts. …
-
Transcriptome reconstruction from both short reads and long sequences
written by Prakki Rama 271Hi all, Could I please know if there is any tool to reconstruct transcriptome from reference genome using both short read and long sequences from …
-
merged transcripts from RNA de novo assembly to create a reference transcriptome
written by teixe005 3Hello, I have created de novo assemblies from RNAseq reads using velvet/oases for different subjects at several time points. For every subject I have a …
Can you elaborate on what you mean by the 'reference transcript HMM profiles'?
The idea is that I would not map to the reference genome but tot a set of HMM profiles (describing the nucleotide sequences).
Why would you want to do that, as opposed to mapping it to the actual transcriptome?
Yes, it is a good question and I cannot really go into the details but in general, I would not know the exact transcriptome but I would only have the HMM profiles based on the homologous sequences that were used to create the HMM profile.
Thanks
K
Do you have access to the homologous sequences themselves? It'd be simple enough to just align to them. With an HMM I suspect you'll have to invent a scoring method for assessing the quality of an "alignment" (really, the probability of a given HMM producing an observed sequence compared to all other HMMs producing said sequence).
No I do not have access to that.
How about searching for the HMM profiles in the reads directly?
These are 10s of millions short reads RNA-seq data and HMM profiles are several 100s bp long
Would that work?
Thanks all.
RNAseq data would probably everything from 35bp to 250bp reads and the HMM profiles are predominantly more than 1kb, thus direct mapping does not work.
I was thinking about doing something like that; first using the RNAseq data do de novo transcripts assembly and then using HMMER to identify the transcripts which contain particular HMM profiles.
K
Why would direct alignment not work?
As I said, I would not have a reference sequence to map to, only a HMM profile of the sequences.