Add a feature where if user specifies somatic_variation_id, readcounts will be done from the normal and tumor BAMs associated with the som_var_model