lizhan96

Dose the batch file is required? I don't have such file.

So could you give me reference code for using MoChA with array data? I can't find that in the link https://github.com/freeseek/mocha/tree/master/wdl#illumina-example you gave me

`bcftools +mocha \ --genome $assembly \ --input-stats $tsv \ --no-version \ --output - \ --output-type b \ --variants ^$dir/$pfx.xcl.bcf \ --calls $dir/$pfx.calls.tsv \ --stats $dir/$pfx.stats.tsv \ --ucsc-bed $dir/$pfx.ucsc.bed \ --cnp...

got it. I have solved the problem. Thanks for your help!

If I want to get the tumor somatic CNV, could you give me some suggestions to process the result of the separate sample?

``` > kegg.dn = enrichKEGG(gene.dn2$ENTREZID, organism = 'mmu', keyType = 'kegg', pvalueCutoff = 0.05, qvalueCutoff = 0.2) --> No gene can be mapped.... --> Expected input gene ID: --> return...