+ ACC=AF086833
+ N=10000
+ mkdir -p refs
+ REF=refs/AF086833.fa
+ BAM=align.bam
+ efetch -db nuccore -format fasta -id AF086833
+ bwa index refs/AF086833.fa
+ samtools faidx refs/AF086833.fa
+ wgsim -N 10000 refs/AF086833.fa r1.fq r2.fq
+ R1=r1.fq
+ R2=r2.fq
+ bwa mem refs/AF086833.fa r1.fq r2.fq
+ samtools sort
+ samtools index align.bam
+ samtools flagstat align.bam