Use when simulating RNA-seq or DNA-seq reads from a reference genome and GTF annotation file, optionally incorporating SNP variants and controlling expression profiles.
hisat2_simulate_reads.py [genome_file] [gtf_file] [snp_file] [base_fname]/home/vimalinx/miniforge3/envs/bio/bin/hisat2_simulate_reads.pyreferences/help.md for complete options and usage detailshisat2_simulate_reads.py when you need synthetic HISAT2-style benchmark reads from a genome, annotation, and SNP set.-d.# Default paired-end RNA-seq simulation
hisat2_simulate_reads.py genome.fa annotation.gtf variants.snp sim
# Single-end RNA-seq with shorter reads
hisat2_simulate_reads.py genome.fa annotation.gtf variants.snp sim --single-end -r 75 -n 100000
# DNA-seq simulation with explicit random seed
hisat2_simulate_reads.py genome.fa annotation.gtf variants.snp sim -d --random-seed 42
# Run internal sanity checks and print extra statistics
hisat2_simulate_reads.py genome.fa annotation.gtf variants.snp sim --sanity-check -v
hisat2_simulate_reads.py with required positional arguments and desired options (e.g., -n for fragment count, -e for expression profile)--sanity-check and review statistics with -v if needed--single-end and/or -d for alternative modes--random-seed for reproducible simulationsbase_fname, including base_fname.sam and at least base_fname_1.fa (plus base_fname_2.fa for paired-end mode)-h prints a non-fatal Python SyntaxWarning before the usage text; treat that as a script warning, not a failed invocation