Question: Plotting DataTrack for non-model organism.
gravatar for anupriyaverma1408
14 months ago by
anupriyaverma14080 wrote:

Hi all!

Today I was trying the package called 'Gviz' for plotting chip-seq peaks as histogram across the whole genome of Mycobacterium tuberculosis. But I was unsuccessful as it keep giving this error: 

Error in .normargGenome(value, seqnames(x)) : 
  supplied 'genome' must be a vector or factor

 These were the steps I was running:


txdb <- makeTxDbFromGFF("/home/m.tub1.gff3",circ_seqs=character(), organism="Mycobacterium tuberculosis")
annoData <- genes(txdb)
annoData1 <-
a <- annoData1[which(annoData1$end == max(annoData1$end)),]

a$label= "+"
a$label1= "-"

bed_file1 = read.csv("/home/peak_annotations_forward_strand.txt", header = T, sep = "\t")
bed_file2 = read.csv("/home/peak_annotations_reverse_strand.txt", header = T, sep = "\t") 

b <- bed_file1$pileup
c <- -(bed_file2$pileup)

gtrack <- GenomeAxisTrack(range = makeGRangesFromDataFrame(data.frame(seqnames=as.character(a$seqname), start=1, end=a$end)))

lim <- c(0, a$end)
dtrack <- DataTrack(data = c, start = bed_file1$start, end = bed_file1$end, name = "Peaks (+)")
dtrack1 <- DataTrack(data = d, start = bed_file2$start, end = bed_file2$end, name = "Peak (-)")

plotTracks(list(dtrack, dtrack1), from = lim[1], to = lim[2], type = "histogram")

Can Gviz handle any non-model organism by taking FASTA or any other method ? Please help me out .

Thank you in advance.

ADD COMMENTlink modified 14 months ago • written 14 months ago by anupriyaverma14080
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 368 users visited in the last hour