LUYAO REN
5 年之前
共有 1 個檔案被更改，包括 1 行新增 和 1 行删除
+ 1
- 1
tasks/gvcf_homo_ref.wdl
查看文件
| awk '{ if ($5 == "<NON_REF>") { print } }' ${gvcf} | sed s'/:/\t/'g | awk '{ if($16 >= 30) { print }}' | cut -f1,2,8 | sed s'/END=//g' > ${sample}.filtered.bed | awk '{ if ($5 == "<NON_REF>") { print } }' ${gvcf} | sed s'/:/\t/'g | awk '{ if($16 >= 30) { print }}' | cut -f1,2,8 | sed s'/END=//g' > ${sample}.filtered.bed | ||||
| /opt/ccdg/bedtools-2.27.1/bin/bedtools merged -i ${sample}.filtered.bed > ${sample}.homo_ref.bed | |||||
| /opt/ccdg/bedtools-2.27.1/bin/bedtools merge -i ${sample}.filtered.bed > ${sample}.homo_ref.bed | |||||
| >>> | >>> | ||||
Loading…