diff --git a/wgs.sh b/wgs.sh
index 177fa76..be9ba9c 100644
--- a/wgs.sh
+++ b/wgs.sh
@@ -24,3 +24,8 @@ grep ">" databases/all_genome.fasta|cut -f 2 -d ">" > results/association.tsv
 
 # Question 3
 
+soft/megahit --k-list 21 --mem-flag 0 -1 EchE_R1.fastq.gz -2 EchE_R2.fastq.gz -o contigs
+
+# Question 4
+
+soft/prodigal -i contigs/final.contigs.fa -f fasta > genes.fna