#Generate 3% divergence
msbar -point 4 -count 16558 toy.fasta > toyheterozygous3percent.fasta
#Cat both files
cat toy.fasta toymutated3percent.fasta > toyheterozygous3percent.fasta
#generated 50X of Illumina paired-end reads
sim_reads --depth 50 toyheterozygous3percent.fasta > toyheterozygous3percentD50.fasta