二代测序数据模拟

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