Question: varscan copy number error
gravatar for mittu1602
20 months ago by
mittu1602170 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 • 640 views
ADD COMMENTlink written 20 months ago by mittu1602170

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 20 months ago • written 20 months ago by Macspider2.9k
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: 1172 users visited in the last hour