https://bioinformaticsonline.com/snippets/view/42592/create-random-1000-indel-in-genome? https://bioinformaticsonline.com/snippets/view/42592/create-random-1000-indel-in-genome Sun, 10 Jan 2021 09:15:26 -0600 https://bioinformaticsonline.com/snippets/view/42592/create-random-1000-indel-in-genome (base) ➜ dupStudy git:(master) ✗ perl ../simuG.pl -refseq simuSNP.simseq.genome.fa -indel_count 1000 -prefix simuINDEL [Sun Jan 10 16:14:00 2021] Starting simuG .. [Sun Jan 10 16:14:00 2021] Check specified options .. Running simuG for SNP/INDEL simulation >> Ignore all options for CNV/inversion/translocation simulation. This simulation use the random seed: 1678440514 The option indel_count has been specified: indel_count = 1000 The option ins_del_ratio has been specified: ins_del_ratio = 1 The option indel_size_powerlaw_alpha has been specified: indel_size_powerlaw_alpha = 2 The option indel_size_powerlaw_constant has been specified: indel_size_powerlaw_constant = 0.5 [Sun Jan 10 16:14:00 2021] Introducing random INDELs based on the following parameters: > indel_count = 1000 > ins_del_ratio = 1 > indel_size_powerlaw_alpha = 2 > indel_size_powerlaw_constant = 0.5 [Sun Jan 10 16:14:00 2021] Simulation completed! :) [Sun Jan 10 16:14:00 2021] Generating output files .. Generating the correspondance map for genomic variants introduced during simulation: simuINDEL.refseq2simseq.map.txt Generating reference-based vcf file for genomic variants introduced during simulation: simuINDEL.refseq2simseq.INDEL.vcf [Sun Jan 10 16:14:00 2021] Done! :)]]> Jit