@@ -845,51 +845,73 @@ filter_read_support <- function(vcf_df, config = NULL) {
845845 if (" DP_TUMOR" %in% colnames(vcf_df )){
846846 if (! any(is.na(vcf_df $ DP_TUMOR )) &
847847 ! is.null(config $ somatic_snv $ allelic_support $ tumor_dp_min )) {
848+ # cat("DP_TUMOR - Number of rows before filtering - ", nrow(vcf_df), "\n")
849+ # cat("DP_TUMOR - Filtering criteria used: ", config$somatic_snv$allelic_support$tumor_dp_min, "\n")
848850 vcf_df <- dplyr :: filter(
849851 vcf_df ,
850852 .data $ DP_TUMOR > = config $ somatic_snv $ allelic_support $ tumor_dp_min )
853+ # cat("DP_TUMOR - Number of rows after filtering: ", nrow(vcf_df), "\n")
851854 }
852855 }
853856 if (" VAF_TUMOR" %in% colnames(vcf_df )){
854857 if (! any(is.na(vcf_df $ VAF_TUMOR )) &
855858 ! is.null(config $ somatic_snv $ allelic_support $ tumor_af_min )) {
859+ # cat("VAF_TUMOR - Number of rows before filtering - ", nrow(vcf_df), "\n")
860+ # cat("VAF_TUMOR - Filtering criteria used: ", config$somatic_snv$allelic_support$tumor_af_min, "\n")
856861 vcf_df <- dplyr :: filter(
857862 vcf_df ,
858863 .data $ VAF_TUMOR > = config $ somatic_snv $ allelic_support $ tumor_af_min )
864+ # cat("VAF_TUMOR - Number of rows after filtering: ", nrow(vcf_df), "\n")
859865 }
860866 }
861867 if (" AD_TUMOR" %in% colnames(vcf_df )){
868+
862869 if (! any(is.na(vcf_df $ AD_TUMOR )) &
863870 ! is.null(config $ somatic_snv $ allelic_support $ tumor_ad_min )) {
871+ # cat("AD_TUMOR - Number of rows before filtering - ", nrow(vcf_df), "\n")
872+ # cat("AD_TUMOR - Filtering criteria used: ", config$somatic_snv$allelic_support$tumor_ad_min, "\n")
864873 vcf_df <- dplyr :: filter(
865874 vcf_df ,
866875 .data $ AD_TUMOR > = config $ somatic_snv $ allelic_support $ tumor_ad_min )
876+ # cat("AD_TUMOR - Number of rows after filtering: ", nrow(vcf_df), "\n")
867877 }
878+
868879 }
869880
870881 if (" VAF_CONTROL" %in% colnames(vcf_df )){
882+
871883 if (! any(is.na(vcf_df $ VAF_CONTROL )) &
872884 ! is.null(config $ somatic_snv $ allelic_support $ control_af_max )) {
885+ # cat("VAF_CONTROL - Number of rows before filtering - ", nrow(vcf_df), "\n")
886+ # cat("VAF_CONTROL - Filtering criteria used: ", config$somatic_snv$allelic_support$control_af_max, "\n")
873887 vcf_df <- dplyr :: filter(
874888 vcf_df ,
875889 .data $ VAF_CONTROL < = config $ somatic_snv $ allelic_support $ control_af_max )
890+ # cat("VAF_CONTROL - Number of rows after filtering: ", nrow(vcf_df), "\n")
876891 }
892+
877893 }
878894 if (" DP_CONTROL" %in% colnames(vcf_df )){
879895 if (! any(is.na(vcf_df $ DP_CONTROL )) &
880896 ! is.null(config $ somatic_snv $ allelic_support $ control_dp_min )) {
897+ # cat("Number of rows before filtering - ", nrow(vcf_df), "\n")
898+ # cat("DP_CONTROL - Filtering criteria used: ", config$somatic_snv$allelic_support$control_dp_min, "\n")
881899 vcf_df <- dplyr :: filter(
882900 vcf_df ,
883901 .data $ DP_CONTROL > = config $ somatic_snv $ allelic_support $ control_dp_min )
902+ # cat("DP_CONTROL - Number of rows after filtering: ", nrow(vcf_df), "\n")
884903 }
885904 }
886905
887906 if (" AD_CONTROL" %in% colnames(vcf_df )){
888907 if (! any(is.na(vcf_df $ AD_CONTROL )) &
889908 ! is.null(config $ somatic_snv $ allelic_support $ control_ad_max )) {
909+ # cat("AD_CONTROL - Number of rows before filtering - ", nrow(vcf_df), "\n")
910+ # cat("AD_CONTROL - Filtering criteria used: ", config$somatic_snv$allelic_support$control_ad_max, "\n")
890911 vcf_df <- dplyr :: filter(
891912 vcf_df ,
892913 .data $ AD_CONTROL < = config $ somatic_snv $ allelic_support $ control_ad_max )
914+ # cat("AD_CONTROL - Number of rows after filtering: ", nrow(vcf_df), "\n")
893915 }
894916 }
895917
0 commit comments