ViewBS error in visualizing BSSeq data

Hi Guys,

I did the analysis of Bisulfite sequencing data. Now I want to visualize the results in ViewBS, but I am getting error in this. I am trying with first two samples as I am getting same error if I use all samples. So I am running below commandline for ViewBS

ViewBS MethCoverage --reference /genome --sample sample1_cov.CX_report.txt.gz,sample1 --sample sample2_cov.CX_report.txt.gz,sample2 --outdir coverage

i am getting error

Use of uninitialized value in multiplication (*) at /ViewBS-0.1.6/lib/Meth/Coverage.pm line 103.
Illegal division by zero at /ViewBS-0.1.6/lib/Meth/Coverage.pm line 103.

I have checked the Coverage.pm file but I couldn't understand it clearly why I am getting this error as I am not a perl expert.

Any help would be highly appreciated. Thanks in advance.

bisulfite sequencing analysis software error

• 778 views • link

updated 5.4 years ago by GenoMax

• written 5.4 years ago by abseq •

0 answers

No answers yet.

More posts like this

"Error in h(simpleError(msg, call)) : error in evaluating the argument 'x' in selecting a method fo…
written 4.3 years ago by buffealo •

Hello, When I am trying to use dba.count in DiffBind, I am getting this error; peaks <- dba.count(peaks, bParallel=FALSE) Computing summits... Sample: sample1.bam125 Sample: sample2.bam125 …
ChIP QC
written 5.2 years ago by mropri

I am using ChIPQC to get quality metrics of my reads. I have `.bam` files and `.narrowpeaks` files. I created a sample csv file called …
bcftools isec not producing sites.txt
written 5.3 years ago by steve

I am running a `bcftools` command like this; bcftools isec -p /data --targets-file targets.bed Sample1.vcf.gz Sample2.vcf.gz And I am getting output files like this; 0000.vcf …
snakemake different multiple input paths not in output or rule all
written 6.0 years ago by StartR •

Hi I am new to snakemake - here is my questions: I have different bam files from multiple folders : folder structure is: ``` Folder1/sample1/libx/folder1/alginment.bam …
SPAdes error code -6
written 7.1 years ago by Lina F

Hi all, I recently ran SPAdes on two samples and encountered an error code of -6. I noticed that one of the samples I used …
Question: Error in gmapR: gsnap not functioning
written 7.7 years ago by amer_ghl •

Hey, I'm trying to preform a bisulfite alignment using gmapR package. This is my code: ========================= library("gmapR") ggd<-GmapGenomeDirectory(file.path(getwd(),"indices"), create=T) gg<-GmapGenome(file.path(getwd(),"test.fa"), directory=ggd, name="seq", create=T) cmet_index <- …
error in mapper script in miRDeep2
written 7.8 years ago by nivya.james2016 •

Hi all I'm used the following mapper script for miRNA alignment in miRdeep2. I am getting error and also the reads_vs_genome.arf is empty. I am …
If a string appears twice in column 1, select the higher value in column 2
written 8.1 years ago by drkennetz

Hi all, I subsampled Illumina fastqs twice using 2 different random seeds. I then mapped all of the reads and got a total number of …
Problem in generating GCWinds and MAPPBILITY files for readDepth
written 9.2 years ago by misbahabas •

Hello I used "createCustomAnnotations" scripts from here for non human data: https://xfer.genome.wustl.edu/gxfer1/project/cancer-genomics/readDepth/index.html I have bam files and reference genome, so I convert bam file into …
Perl program: The sequence does not appear to be FASTA format (lacks a descriptor line '>')
written 11.5 years ago by dago

I get the following error when running a perl program: ``` Use of uninitialized value $Bio::DB::NCBIHelper::HOSTBASE in concatenation (.) or string at /usr/share/perl5/Bio/DB/Query/GenBank.pm line 103. …