find . -maxdepth 1 -type f -name "*.g.vcf.gz" ! -name "*interval*" -exec...
1、群體區(qū)分 將每個(gè)群體中所包含的個(gè)體一行一個(gè),生成popA,popB文件 2、vcftools進(jìn)行對(duì)snp文件進(jìn)行劃窗計(jì)算fst和pi 窗口大...
gstacks -I aligned/ -M ../../00.extracted/pop_newall1.txt -t 32 -O gstac...
vcf2treemix.sh $FILE.vcf.gz $FILE.clust .clust 格式是三列,第一列和第二列相同均為樣品名稱,第三列...
使用user_modify.py 建立新用戶之后, 使用sacctmgr add user User=xxx Account=normal 命令...
ls *.fq.gz | while read f; do b="${f%.fq.gz}.bam"; if [ ! -f "$b" ]; the...
1. bam2fasta -o xx.ccs.bam 2.hifiasm -o xjg04 -t 50 --hg-size 1.5g xjg.h...
宏基因組二代測(cè)序結(jié)果binning之后的序列可能存在重復(fù),首先使用prokka將binning后的序列轉(zhuǎn)化為蛋白序列, 再使用CD-hit 命令...
for file in */umap.svg; do mv "$file" "$(dirname "$file")/$(basename $(d...