08简化测序流程2-无参ddrad

主流程

没有参考基因组,直接根据序列聚类来做。

配置文件:   
perl /share/nas6/xul/program/ddrad/create_profile_for_ddrad.pl -n -i 邮件里路径/2.clean_data/
注意老师给的数据命名,自动生成的名字是文件名字前缀,老师可能另外要求了编号,-n 表示无参,-i 输入所有测序数据所在的所有目录。
如果老师重新命名,可以自己编写了对照表,第一列是老的名字,第二列是新的名字,一一对应,然后 -l 参数修改生成配置文件中的名字。
  
主流程: 
$ nohup perl /share/nas6/zhouxy/pipline/stacks-seq/current/stacks.pip.v2.pl -i ddrad_config.yaml -o analysis &
$ nohup perl /share/nas1/yuj/pipline/stacks-seq/v1.0/stacks.pip.v2.pl -i ddrad_config.yaml -o analysis & # 修改了resub循环数,暂时还没验证好不好用
   
整理结果:
$ perl /share/nas6/xul/program/ddrad/get_result_noref.pl  analysis complete_dir
无分组.不用管treemix结果

$ cp /share/nas6/xul/program/ddrad/html_for_ddrad_noref/ddrad_noref_report.cfg . && realpath ddrad_noref_report.cfg

生成报告:
$ perl /share/nas1/yuj/program/ddrad/html_for_ddrad_noref/ddrad_noref_report.pl -id  complete_dir -cfg ddrad_noref_report.cfg # 测序平台改成novaseq 修改电话号码
# perl /share/nas6/xul/program/ddrad/html_for_ddrad_noref/ddrad_noref_report.pl -id  complete_dir -cfg ddrad_noref_report.cfg

子代与亲本相似度计算(无参简化)

cd xxxx
 python3 /share/nas1/yuj/script/Denovo/calculate_similarity.py -h

  -i1 samples.pop.snp.recode.vcf
  -i2 tags.consensus.fa
  -g relationship.txt    三列,f1 parent1 parent2
  -o similarity_result.txt