diff --git a/src/SNP_calling.config b/src/SNP_calling.config
index 1f8c6670d329058afcf64922302e346061f531ac..34660e06afb9bcaff9d1855155369ee51a5ea4cf 100644
--- a/src/SNP_calling.config
+++ b/src/SNP_calling.config
@@ -34,10 +34,10 @@ profiles {
         container = "samtools:1.7"
       }
       withName: samtools_SNP_tumor {
-        container = "samtools:1.7"
+        container = "bcftools:1.7"
       }
       withName: samtools_SNP_norm {
-        container = "samtools:1.7"
+        container = "bcftools:1.7"
       }
       withName: vcf_to_csv_tumor {
         container = "gatk:4.0.8.1"
diff --git a/src/SNP_calling.nf b/src/SNP_calling.nf
index 0de5eb46d463ee825ed7fd80674ab4df1948ddaa..2c9cdfcfb977dae64ff17a842f05c0ad223e60c1 100644
--- a/src/SNP_calling.nf
+++ b/src/SNP_calling.nf
@@ -355,7 +355,8 @@ process samtools_SNP_tumor {
 
   script:
 """
-samtools mpileup -AE -uf ${fasta} ${bam_tumor} | \
+bcftools mpileup -AE -f ${fasta} ${bam_tumor} --output-type v \
+-a FORMAT/AD,FORMAT/ADF,FORMAT/ADR,FORMAT/DP,FORMAT/SP,INFO/AD,INFO/ADF,INFO/ADR | \
 bcftools call -mv --output-type v > ${file_id_tumor}_raw.vcf
 bcftools filter -s LowQual -e '%QUAL<20 || DP>100' ${file_id_tumor}_raw.vcf \
 > ${file_id_tumor}_filtered.vcf
@@ -379,7 +380,8 @@ process samtools_SNP_norm {
 
   script:
 """
-samtools mpileup -AE -uf ${fasta} ${bam_norm} | \
+bcftools mpileup -AE -f ${fasta} ${bam_norm} --output-type v \
+-a FORMAT/AD,FORMAT/ADF,FORMAT/ADR,FORMAT/DP,FORMAT/SP,INFO/AD,INFO/ADF,INFO/ADR | \
 bcftools call -mv --output-type v  > ${file_id_norm}_raw.vcf
 bcftools filter -s LowQual -e '%QUAL<20 || DP>100' ${file_id_norm}_raw.vcf \
 > ${file_id_norm}_filtered.vcf
@@ -400,10 +402,10 @@ process vcf_to_csv_tumor {
   script:
 """
 gatk VariantsToTable -V ${file_id_tumor}_raw.vcf \
--F CHROM -F POS -F TYPE -GF GT -GF AD -GF AF \
+-F CHROM -F POS -F TYPE -GF GT -GF AD -F AD -F DP \
 -O ${file_id_tumor}_raw.csv
 gatk VariantsToTable -V ${file_id_tumor}_filtered.vcf \
--F CHROM -F POS -F TYPE -GF GT -GF AD -GF AF \
+-F CHROM -F POS -F TYPE -GF GT -GF AD -F AD -F DP \
 -O ${file_id_tumor}_filtered.csv
 """
 }
@@ -422,10 +424,10 @@ process vcf_to_csv_norm {
   script:
 """
 gatk VariantsToTable -V ${file_id_norm}_raw.vcf \
--F CHROM -F POS -F TYPE -GF GT -GF AD -GF AF \
+-F CHROM -F POS -F TYPE -GF GT -GF AD -F AD -F DP \
 -O ${file_id_norm}_raw.csv
 gatk VariantsToTable -V ${file_id_norm}_filtered.vcf \
--F CHROM -F POS -F TYPE -GF GT -GF AD -GF AF \
+-F CHROM -F POS -F TYPE -GF GT -GF AD -F AD -F DP \
 -O ${file_id_norm}_filtered.csv
 """
 }