This is a test version of Biostars. For the public version, visit https://www.biostars.org.
CNVRs visualization with karyoploteR

Hi,

I am using karyoploteR to visualize CNV. I am getting following error. But I do have the genomic region in a bed file. Could anyone provide a suggestion what might have gone wrong?

Error in plotKaryotype(zoom = "Chr01:3647-12170735", main = "BedGraphs", : You are trying to set the zoom to a region not part of the current genome. In addition: Warning message: In .Seqinfo.mergexy(x, y) : The 2 combined objects have no sequence levels in common. (Use suppressWarnings() to suppress this warning.)

cnv
• 949 views • link
written by rthapa • rthapa

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • summarizeOverlaps warning A. thaliana differential expression
    written by eva.lucarelli • eva.lucarelli
    Hello, I am currently working on a differential expression tutorial regarding gene expression in Arabidopsis thaliana. At some point in the exercise, we use the …
  • KaryoploteR for copy number
    written by rthapa • rthapa
    Hi, I am trying to plot a genome-wide Copy number variation using KaryoploteR. I am getting an error saying duplicate chromosome names in the genome. …
  • Error in Gviz BiomartGeneRegionTrack
    written by codezy • codezy
    Dear friends, I would like to run the "BiomartGeneRegionTrack" in package Gviz to get the annotation file. But it seems I can not use the …
  • Error while annotating with Chipseeker
    written by hamaor • hamaor
    I followd the Diffbind pipeline for a ChipSeq data: peaksets2 = dba(sampleSheet="SampleSheet.csv",config = data.frame(RunParallel=TRUE, AnalysisMethod=DBA_DESEQ2, bCorPlot=FALSE, bUsePval=TRUE, minQCth=10),minOverlap = 2) peaksets2 <- dba.count(peaksets2) peaksets2 <- …
  • Error using Tabix
    written by AS • AS
    Hi, I am trying to use vcf files downloaded from public database (NCBI SRA Bioproject:PRJNA486011). I am using only R to do that because I …
  • CHIPseeker error now
    written by lzhizhao9179 • lzhizhao9179
    tagMatrix_1 <- getTagMatrix(peak_1, weightCol =NULL, windows = promoter_1) > Error in split.default(1:length(windows), as.factor(seqnames(windows))) : > group length is 0 but data length > 0 > …
  • RNA Seq Workflow
    written by jadepinket • jadepinket
    I am using the RNA-seq workflow: gene-level exploratory analysis and differential expression. I have summerizeoverlap output 0 so <- summarizeOverlaps(ebg, "Epigenomics02-96.alignments.bam", mode="Union", + singleEnd=FALSE, + …
  • Annotate DMR information
    written by yassimaryam • yassimaryam
    We have used “GRCh37/hg19” for reference genome in our project. We convert output file to GenomicRange object that has been 3 columns included that: (Seqnames, …
  • intansv - problem at the annotation step
    written by rajini28m • rajini28m
    Hi Everyone, I am doing structural variant detection using different programs such as DELLY, Softsearch, Svseq. After detection I am integrating the SV calls from …
  • problem with methylKit of the annotation
    written by hua.peng1314 hua.peng1314
    I am dealing with RRBS data follow [the manual][1]. I met some problem when doing annotation.The problem showed below. ``` > read.transcript.features("mm9.refseq.genes.bed",remove.unsual=TRUE,up.flank=1500,down.flank=500,unique.prom=TRUE) > diffAnn <- …