fn plot<PB: ProtocolBehavior>( input: &str, format: &str, output_prefix: &str, is_multiple: bool, is_tree: bool, ) -> Result<(), Box<dyn Error>>