WES

run gistic2 with sequenza segmentation output from whole exome sequencing

Gistic was designed for SNP6 array data. I saw many papers use it for whole exome sequencing data as well.\ I have the segment files from sequenza and want…

bwa-mem multi-mapping reads

I was using TEQC to do quality control of my WES bam files aligned by bwa-mem. My data are paired end, so a function reads2pairs is called to make the…