static int update_bioseq_if_necessary(GtRegionMapping *rm, GtStr *seqid,
GtError *err)
{
int had_err = 0;
gt_error_check(err);
gt_assert(rm && seqid);
if (!rm->sequence_file || gt_str_cmp(rm->sequence_name, seqid)) {
gt_str_delete(rm->sequence_file);
rm->sequence_file = region_mapping_map(rm, gt_str_get(seqid), err);
if (!rm->sequence_file)
had_err = -1;
else {
if (!rm->sequence_name)
rm->sequence_name = gt_str_new();
else
gt_str_reset(rm->sequence_name);
gt_str_append_str(rm->sequence_name, seqid);
gt_bioseq_delete(rm->bioseq);
rm->bioseq = gt_bioseq_new_str(rm->sequence_file, err);
if (!rm->bioseq)
had_err = -1;
}
}
return had_err;
}
static int update_seq_col_if_necessary(GtRegionMapping *rm, GtStr *seqid,
GtError *err)
{
int had_err = 0;
gt_error_check(err);
gt_assert(rm && seqid);
/* for mappings, we need to load the changed sequence, if needed... */
if (rm->mapping) {
if (!rm->sequence_file || (gt_str_cmp(rm->sequence_name, seqid))) {
gt_str_delete(rm->sequence_file);
/* ignore MD5 hashes when using region mappings */
if (gt_md5_seqid_has_prefix(gt_str_get(seqid))) {
rm->sequence_file = region_mapping_map(rm,
gt_str_get(seqid)
+GT_MD5_SEQID_TOTAL_LEN,
err);
} else
rm->sequence_file = region_mapping_map(rm, gt_str_get(seqid), err);
if (!rm->sequence_file)
had_err = -1;
else {
/* load new seqcol */
if (!rm->sequence_filenames)
rm->sequence_filenames = gt_str_array_new();
else
gt_str_array_reset(rm->sequence_filenames);
gt_str_array_add(rm->sequence_filenames, rm->sequence_file);
if (!rm->sequence_name)
rm->sequence_name = gt_str_new();
else
gt_str_reset(rm->sequence_name);
gt_str_append_str(rm->sequence_name, seqid);
gt_seq_col_delete(rm->seq_col);
rm->seq_col = gt_bioseq_col_new(rm->sequence_filenames, err);
if (!rm->seq_col)
had_err = -1;
}
}
} else {
/* ...otherwise, just make sure the seqcol is loaded */
if (!rm->seq_col) {
if (rm->encseq) {
if (!(rm->seq_col = gt_encseq_col_new(rm->encseq, err)))
had_err = -1;
} else {
gt_assert(rm->sequence_filenames);
if (!(rm->seq_col = gt_bioseq_col_new(rm->sequence_filenames, err)))
had_err = -1;
}
}
if (!had_err && rm->usedesc) {
if (rm->seqid2seqnum_mapping)
gt_seqid2seqnum_mapping_delete(rm->seqid2seqnum_mapping);
rm->seqid2seqnum_mapping =
gt_seqid2seqnum_mapping_new_seqcol(rm->seq_col, err);
if (!rm->seqid2seqnum_mapping) {
had_err = -1;
}
}
}
return had_err;
}
