ChrM in the Varscan2 copynumber output
0
0
Entering edit mode
5.7 years ago
lhaiyan3 ▴ 60

Hi,

I am using the varscan2 to call the copynumber for the WES. But my output.copynumber file just has the chrM call. I used the following script, can anyone please give me some suggestions? Thanks very much. HY

samtools mpileup -q 1 -f $ref normal.bam tumor.bam | java -jar$VARSCANHOME/varscan.jar copynumber varScan --mpileup 1

java -Xmx64g -jar $VARSCANHOME/varscan.jar copyCaller output.copynumber --output-file output.copynumber.called varscan2 copy number • 1.5k views ADD COMMENT 0 Entering edit mode Do the chromosomes specified in$ref correspond to those in the bam file? Easiest way to check this would be to run samtools idxstats on the bam file and grep '>' on the reference fasta.

0
Entering edit mode

I tried both, I have the chrM to chrY in the bam file and ref file, but the chrM is the first, it looks like the program stopped after running the chrM.

0
Entering edit mode

You mean the chromosomes are not in the same order in both files?

0
Entering edit mode

No, they are in the same order in both files. I mean chrM is the first in both files.