python extract_gene_sequence.py -i fasta -l list -o extract.fa

sh run_blast.sh -d <database_file> -q <query_file> [-o <output_prefix>]

Rscript EDEGR_RNA_seq.R -i sample.list -n 3 -p 0.01 -f 0.05 -l 1 -c 10 -g comparisons.txt

python rename_gff.py -i input.gff -o output.gff

python convert_blast_results.py QTL_A02.fasta.ECD10.7.out QTL_A02.fasta.ECD10.csv -chr A02 -e 0.0001 -b 10000