I have aligned RNA-seq samples in triplicate for two conditions from a colleague (*.BAM files).
I want to examine the log2fc of condition1/condition2 in defined 1kb regions that I have annotated a .BED file for the genome build I am working with (hg19).
Strategies I am thinking about:
bedtools intersect to narrow down the aligned reads (.BAM files) to just the ones aligning to the defined regions of interest (.BED file) --> Deeptools bamCompare to calculate log2fc(condition1/condition2)--> Deeptools computeMatrix and plotheatmap
featureCounts with SAF annotation using the annot.ext argument to define regions of interest then proceeding with countToFPKM to calculate FPKM and then plotting differentially expressed regions in a heatmap
I'm a total novice in RNA-seq analysis, so any feedback/suggestions from experienced bioinformaticians is deeply appreciated! Thank you :)
0 answers
No answers yet.
Log in to answer this question.
-
Deeptools ComputeMatrix returns 0 values for clearly populated regions
written by eglsidlauskaite •
-
peak calling in CUT&RUN(Chipseq) Analysis
written by minoo •
-
Problem with the output of Deeptools PlotProfile: a strange pattern of repetitive summits
written by Alireza •
-
Chip Seq Enrichment Profile
written by Mohammed •
-
ChIP-seq Differential Binding analysis for specified genomic regions of interest? (DiffBind/csaw?)
written by jennxyzhang •
-
deeptools computeMatrix reference-point, want to plot 1kb before TSS and 1kb after TES then use plo…
written by jennxyzhang •
-
deepTools computeMatrix need at least one array to concatenate error
written by leeb4 •
-
Deeptools Normalization Problem
written by silvia.raineri •
-
Extract genomic coverage of specific regions from bam file
written by wsciekly.maciek •
-
Sorting heatmaps by gene expression
written by mmmmcandrew
