Gviz Coverage Plots
I have two bam files from single-cell RNA sequencing mapped to the reference genome using CellRanger, I can view them in IGV and I have a particular region where the pattern of reads mapped to the reference genome are different between the two bam files but when I try the following code to plot trackplots with coverage, the resulting plot is not consistent with what I see in IGV. I am not very familiar with this, aren't I supposed to normalize the data as well? please help me.
txdb <-TxDb.Hsapiens.UCSC.hg38.knownGene
txTr <- GeneRegionTrack(txdb, chromosome = "chr1",
start = 153561336, end = 153561564)
dTrack <- DataTrack(range = bam.file, genome = "hg38", type = "l",
name = "Coverage", window = -1,
chromosome = "chr1")
dTrack0 <- DataTrack(range = bam.file.0, genome = "hg38", type = "l",
name = "Coverage", window = -1,
chromosome = "chr1", col="red")
plotTracks(c(dTrack,dTrack0, txTr), from = 153561357, to = 153561585)
• 1,790 views
•
link
updated
by
Ram
•
written
by
Researcher •
•
written
by
Researcher •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
ChIP-seq Error with rtracklayer
written by shamzar •
-
Help manually processing strand in paired-end reads when fishing for lariats in bulk RNA-seq
written by txema.heredia •
-
More efficient way of lifting over and collapsing CNV GenomicRanges using rtracklayer?
written by j.torpy •
-
Tutorial: Methylation Analysis Tutorial in R - Part 1
written by Hamid Ghaedi
-
Can't display protein domain data from UCSC table unipDomain using Gviz
written by paul.jaschke •
-
Detect Transgene Insertion Site using SAM format
written by kspata
-
ChIP-Seq experiment, MACS2 bdgpeakcall 0 fold-change, pvalue and qvalue
written by Bastien Hervé
Do not post on multiple fora just because you want a fast response, it's plain rude. At least point to the other post so we know you're asking on multiple fora.
Post on bioinfo SE: https://bioinformatics.stackexchange.com/questions/21122/gviz-coverage-plots
Sorry, I didn't know that posting it on another forum was rude but I really need help with this.
I understand that but stick to one place at least for a while and always link to other places you've posted (reasonable places, we don't need to see your reddit post if you're not comfortable sharing your reddit username, for example)