Question: Plotting DataTrack for non-model organism.
gravatar for anupriyaverma1408
11 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 11 months ago • written 11 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: 172 users visited in the last hour