Question: varscan copy number error
gravatar for mittu1602
15 months ago by
mittu1602150 wrote:

I am running varscan2 and it shows me this error and gives a blank output.copynumber only with the headers and this error on the terminal:

 Min coverage:  10
    Min avg qual:   15
    P-value thresh: 0.01
    Reading input from Normal-SL100.mpileup
    Reading mpileup input...
    Parsing Exception on line:
    chrM    15655   T   0   *   *   0       

I created mpileup file for normal and tumor with the following command:

> samtools mpileup -f ucsc.hg19.fasta -d 8000 -A -E Normal.bam Tumor.bam | awk '{if($4 != 0) print $0}' > Normal-Tumor.mpileup

Thank you

cnv varscan • 508 views
ADD COMMENTlink written 15 months ago by mittu1602150

I suspect that your | awk '{if($4 != 0) print $0}' is removing the header from the raw mpileup file. Try using | awk '{if($4 != 0 || substr($1, 0, 1)=="#") print $0}' instead.

ADD REPLYlink modified 15 months ago • written 15 months ago by Macspider2.8k
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.3.0
Traffic: 1104 users visited in the last hour