2079 days ago
Corona variant calling steps !
...bam bcftools mpileup -O b -o results/bcf/$basename.bcf -f data/ref_genome/sequences.fasta results/bam/$basename.sorted.bam bcftools call --ploidy 1 -m -v -o results/bcf/$basen...1049 days ago
Bash script to simulate a genome !
...CA_003401745.1_ASM340174v1_genomic.fna_upper.diploid.vcf.gz \ out=GCA_003401745.1_ASM340174v1_genomic.fna_upper.diploid.fasta.gz \ ziplevel=6 \ ploidy=2 \ subrate=0.0192 \ indelr...922 days ago