diff --git a/nanocaller_src/generate_SNP_pileups.py b/nanocaller_src/generate_SNP_pileups.py index 0df719a..ac09ab3 100644 --- a/nanocaller_src/generate_SNP_pileups.py +++ b/nanocaller_src/generate_SNP_pileups.py @@ -162,6 +162,7 @@ def ex_bed(tree, pos): pos_list=output.keys() current_depth=[] + depth=0 output_pos,output_ref,output_mat,output_dp,output_freq=[],[],[],[],[] if pos_list: @@ -232,5 +233,6 @@ def ex_bed(tree, pos): output_dp=np.array(output_dp) output_freq=np.array(output_freq) + depth=np.mean(current_depth) - return (output_pos,output_ref,output_mat,output_dp,output_freq,np.mean(current_depth)) \ No newline at end of file + return (output_pos,output_ref,output_mat,output_dp,output_freq,depth) \ No newline at end of file