diff --git a/processbam.c b/processbam.c
index f251eaa..dd04544 100644
--- a/processbam.c
+++ b/processbam.c
@@ -75,6 +75,7 @@ void load_bam( parameters *params, configuration *cfg, bam_info* in_bam, char* p
 		( in_bam->libraries)[i]->fastq1 = NULL;
 		( in_bam->libraries)[i]->fastq2 = NULL;
 		( in_bam->libraries)[i]->divet = NULL;
+		( in_bam->libraries)[i]->is_valid = 0;
 		( in_bam->libraries)[i]->mappings_discordant = ( discordantMapping **) getMem( (NHASH + 1) * sizeof( discordantMapping *));
 		for( j = 0; j < NHASH; j++)
 			( in_bam->libraries)[i]->mappings_discordant[j] = NULL;
@@ -245,6 +246,7 @@ void load_bam( parameters *params, configuration *cfg, bam_info* in_bam, char* p
 		fprintf( stderr, "\nLibrary %s\n\tMean: %f\n\tStdev: %f\n\n", ( in_bam->libraries)[i]->libname, ( in_bam->libraries)[i]->frag_avg, ( in_bam->libraries)[i]->frag_std);
 		fprintf( logFile, "\nLibrary %s\n\tMean: %f\n\tStdev: %f\n\tRead Length: %d\n\n", ( in_bam->libraries)[i]->libname, ( in_bam->libraries)[i]->frag_avg, ( in_bam->libraries)[i]->frag_std, ( in_bam->libraries)[i]->read_length);
 		set_library_min_max( ( in_bam->libraries)[i]);
+		is_library_valid(( in_bam->libraries)[i]);
 		plot_histogram( params, cfg, in_bam->sample_name, in_bam->libraries[i]->libname, second_test_pass[i], second_pass_fragments[i], in_bam->libraries[i]->frag_avg, in_bam->libraries[i]->frag_std);
 	}