2076 days ago
Samtools commands for bioinformatician !
...ej.fasta sal_sej_sorted.bam.bam > sal_vars.bcf ### calling snp and indels # -c : find snp # -v : outp...e "bcftools" to see current version and commands ### calling snp and indels with no frequency threshold...1603 days ago
1526 days ago
Create random 1000 INDEL in genome !
...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 parame...1211 days ago
Bash script to simulate a genome !
...bgzip -@75 > GCA_003401745.1_ASM340174v1_genomic.fna_upper.fasta.gz # Convert to diploid with approximately 2% heterozygosity rate, max indels=20bp # mutate.sh part of BBT...918 days ago