Quantification for small RNAseq data
Hi,
I want to quantify smallRNA seq data (single end, mosquito infected with virus). Can you suggest how am I going to do it? I am just new to bioinformatics. Thank you!
quantification
seq
data
smallrna
• 789 views
•
link
updated
by
rfran010 •
•
written
by
ZuelTech •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
RNA seq strandedness (reverse-strand)
written by ZuelTech •Hi. I have an RNA seq data that is reverse-stranded. The library type is ISR. I want to count reads mapping to the genomes. I …
-
Kallisto bustools for scRNA-seq
written by t.ru •Hi, I followed several tutorials about single cell rna-seq analysis for paired end. However, I need to quantify single ends scRNA-seq dataset. I have used …
-
Find HBV integration sites in human genome.
written by Chirag NepalHi there, I want to identify HBV (virus) integration sites in human genome. I have single end CAGE-seq (~30 nucleotide long) data on HBV patients. …
-
True meaning of global transcriptome analysis
written by babasaraki •Dear All, I have two different cancer cell lines (EJ28 and TCCSUP) that are persistently infected with virus and called them EJ28Pi and TCCSUPPi respectively. …
-
annotate and analyse mutation from cosmic
written by sahhomaha •Hi, I want to know the pipline aboout: how can I download and analyse mutations from cosmic cancer cataloge , and how can I get …
-
Any methods available to do QC analysis of Pacbio raw data??
written by karthicHi, I am very new to NGS. Am going to have pacbio raw data in few months. What am looking for is, methods to verify …
-
How to get unquantified reads with featurecounts
written by Chynoki •Hi there, I'm using featurecounts and I'm wondering if their is a way to get all reads which where not quantified by this software (i.e …
-
How to separate htseq-count table to coding RNAs and non-coding RNAs
written by Naresh D J •Hi, I have generated the raw read counts for genes from RNA-seq data using htseq-count. Now I want to separate the this table into coding …
-
Sample-matched RNA-Seq gene expression analysis
written by Nicolas Rosewick<p>HI,</p> <p>I've a bunch of sample-matched (T0 and T1) data of uninfected and infected (virus) RNA-Seq samples. So:</p> <pre> T0 T1 Patient 1 0.4% 40% …
-
Single cell RNAseq - how many cells enough to determine DE?
written by GenosaHello, I am new to this forum. I am interested in SS-RNAseq and would like to know what statistical method should be used to calculate …
You may try looking into our collaborator's pipeline to at least get an idea of the steps, but if you're completely new to bioinformatics, I think there is a lot of ground to cover here.
https://github.com/junchaoshi/sports1.1
PANDORA-seq expands the repertoire of regulatory small RNAs... | Nat. Cell Biol
The accompanying paper focuses on a library construction method, but uses the pipeline, so you can get an idea of what types of questions they ask and answer. I'm not entirely sure how different small RNA analysis would be for mosquito.
I would suggest collecting a few more pipelines/workflows/papers and working through them to understand all the steps more deeply, focusing on the purpose of each step and why that step is important. Of course if you can find a pipeline that matches what you want, or at least a paper, then it's easier to follow that. But if not, then learning about each step can help you properly apply it to your specific question.
In general, you will want to trim reads to remove any readthrough of adapter sequence and then map the small RNA to the reference with minimal mismatches allowed (maybe 1 or 2, if not 0). This may be a problem if your mosquito species is not the same as the reference mosquito species, so I would check on that too.
Trimming reads is important because small RNAs have a higher occurrence of the sequencer reading through the small fragment and into the adapter on the other end (e.g. 3' ends of reads may be enriched with adapter sequence). Then, since reads are so small, multimapping can be more of an issue so it's good to have exact matches only. This includes no "soft clipping", which is a mode where an aligner will potentially ignore the ends of a read if they don't map exactly to a certain loci. Having exact matches and no soft clipping means adapter sequence in your reads can have a larger negative effect on alignments.