![](https://img.haomeiwen.com/i25274977/dec9d7c81ace475a.png)
环境
##Assemblytics
git clone https://github.com/MariaNattestad/Assemblytics.git
###需要Assemblytics_uniq_anchor.py
ln -s /public/home/fengting/work/6.25111/Assemblytics/scripts/* ~/miniconda3/bin
chmod 754
测试
1构建泛基因组
python /public/home/fengting/demo/psvcp/test/psvcp_v1.01-main/1Genome_construct_Pangenome.py -t 12 /public/home/fengting/demo/psvcp/test/psvcp_v1.01-main/example7/genome_gff_dir_example /public/home/fengting/demo/psvcp/test/psvcp_v1.01-main/example7/genome_list
###把主程序(.sh)numcer 中的-t 32 参数去掉
![](https://img.haomeiwen.com/i25274977/9a36796d761299e1.png)
2 reads比对到泛基因组
python ~/demo/psvcp/test/psvcp_v1.01-main/Construct_pan_and_Call_sv.py example7/genome_gff_dir_example/ example7/genome_list -fqd example7/fq_dir_example/ -o pop
## -o 具体参数暂时不知
![](https://img.haomeiwen.com/i25274977/ff5735198062c33e.png)
![](https://img.haomeiwen.com/i25274977/9da2f71c2c09de92.png)
![](https://img.haomeiwen.com/i25274977/6cd0388656e871e8.png)
![](https://img.haomeiwen.com/i25274977/fd07008fd88fd2c0.png)
![](https://img.haomeiwen.com/i25274977/78b4fa64bc3f7737.png)
![](https://img.haomeiwen.com/i25274977/7e97be42017ef6d1.png)
网友评论