Entering edit mode
8.8 years ago
tejaswikoganti ▴ 70
I am trying to run Varscan somatic command for the first time. I have single tumor-normal mpileup file with reads only for certain regions on chr 17 and with a tumor purity of 50%. My command looks like this -
java -jar Varscan.jar somatic file.mpileup OUTPUT --mpileup 1 --tumor-purity 0.5
I encountered an error saying
Error: Invalid format or not enough samples in mpileup: ï¿ï¿½BCï¿½
I am not sure if this is because my file is really small. Has anybody else faced this problem?
Can varscan somatic call be stored in vcf format. Because i could get output only in .snp and .indel.
Since some other downstream tools use vcf format. What tools did you use after u called somatic mutations.
I just manually convert the file format for annotation in ANNOVAR.
I am only aware of the mpileup commands (mpileup2snp, mpileup2indel, mpileup2cns) having a parameter to produce .vcf files. However, I didn't develop the program, so it is possible that I missed something.
Thanks Charles !
I wanted to extract the read count that support the reference/variant call for each Normal/Tumor pair, to see if the variant is occurring in all or few libraries. There is readCount script from Varscan
But the problem is that I have assembled all Normal/Tumor for 5 patients in the same file. Do you know some easy way on how to get read count supporting reference/variant frequency for each patient. If not, I might just have to redo mpileup for each normal/tumor pair.
I would recommend asking your question on a VarScan discussion group.
I think the separate
.pileupfiles will work, but I'm guessing that this function won't really save you much time (since the pileup files already contain a coverage column). However, I honestly haven't tried this particular function.