fix: tweak mpileup params

Danderson123 · Jan 23, 2025 · f23bfa2 · f23bfa2
1 parent 440096e
commit f23bfa2
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/amira/result_operations.py b/amira/result_operations.py
@@ -767,7 +767,7 @@ def genotype_promoters(
 
 def get_mean_read_depth_per_contig(bam_file, core_genes=None):
     # Run samtools depth command and capture output
-    command = ["samtools", "mpileup", "-aa", "-Q", "0", "--ff", "UNMAP", bam_file]  # -aa ensures all positions are included
+    command = ["samtools", "mpileup", "-aa", "-Q", "0", "--ff", "UNMAP,QCFAIL,DUP", bam_file]
     # Run the command and capture the output
     result = subprocess.run(command, capture_output=True, text=True, check=True)
     # Parse the mpileup output