diff --git a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/ExtractAlignedFastq.scala b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/ExtractAlignedFastq.scala
index 98f0c704cabe416eb841a38fd7754d5049384b74..d53b1afd0fc6b5ca8e86b4af7e018607afcc1ea6 100644
--- a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/ExtractAlignedFastq.scala
+++ b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/ExtractAlignedFastq.scala
@@ -223,5 +223,17 @@ object ExtractAlignedFastq extends ToolCommand {
     val commandArgs: Args = new OptParser()
       .parse(args, Args())
       .getOrElse(sys.exit(1))
+
+    val memFunc = makeMembershipFunction(
+      iv = makeFeatureFromString(commandArgs.intervals),
+      inAln = commandArgs.inputBam,
+      minMapQ = commandArgs.minMapQ,
+      commonSuffixLength = commandArgs.commonSuffixLength)
+
+    selectFastqReads(memFunc,
+      inputFastq1 = commandArgs.inputFastq1,
+      inputFastq2 = commandArgs.inputFastq2,
+      outputFastq1 = commandArgs.outputFastq1,
+      outputFastq2 = commandArgs.outputFastq2)
   }
 }