rasusa --input control_lambda_7_2_2021.fastq.gz --coverage 30 --genomes-size 50000 --output lambda_subsample.fastq.gz