# --missing-to-ref表示缺失的GT表示为0/0,-merge操作多等位基因位点,这里表示不产生多等位基因位点。 # 当然对于header使用--use-header指定,info列的合并使用--info-rules指定规则。 bcftools merge--missing-to-ref--merge none--output'$workdir'/test/combined.vcf.gz--output-type z--threads30 5. 上...
这里做下修改ls *.raw.snv_indel.vcf > all.vcf.lst while read vcf ;do bcftools view -i "INFO/DP>10" -o ${vcf}.gz -Oz ${vcf} bcftools index ${vcf}.gz done < all.vcf.lst bcftools merge -Ov -o all_sample.merge.vcf *.gz 参考: http://www.htslib.org/doc/1.0/bcftools.html#...
#define MERGE_MAX 1 #define MERGE_AVG 2 static double get_min_score(args_t *args, int iskip) { int i; double score, min_score = HUGE_VAL; for (i=0; i<args->nfold; i++) { if ( i==iskip ) continue; score = som_get_score(args->som[i], args->vals, args->bmu_th);...
to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the following conditions: The above copyright notice and this permission notice shall be included in all copies or...
[ci skip] Merge PR 21878 … Verified 6d5b92d Contributor Author dariober commented May 2, 2020 Thanks for the pointer John! Setting bcftools >= 1.9 seems to have fixed it. I wonder why the solver chose version 1.8 rather than the latest or any other... But nevermind. dariober close...
(the "Software"), to deal in the Software without restriction, including without limitation the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the ...