二代测序数据模拟

dwgsim

source /hwfsxx1/ST_HN/P18Z10200N0423/huty/software/miniconda3/etc/profile.d/conda.sh
conda activate base

conda install dwgsim

模拟数据

 dwgsim -1 100 -2 100 -r 0 -R 0 -X 0 -e 0 -E 0 \
 -N 3000000 -P $base ../00_ref/$i 3m/$base/$base
-1 -2 reads长度
-r 突变率
-R fraction of mutations that are indels
-X probability an indel is extended
-e per base/color/flow error rate of the first read
-E per base/color/flow error rate of the second read
-N reads数
-P 序列前缀
-d outer distance between the two ends for pairs [500]
文件名称前缀

https://davetang.org/wiki/tiki-index.php?page=DWGSIM
https://github.com/nh13/DWGSIM/blob/main/docs/01_Introduction.md

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容