INNER CODE UNIT · C++

padding_size

DecodeGenetics/graphtyper · src/graph/constructor.cpp:949

        std::size_t padding_size = EXTRA_SEQUENCE_LENGTH - inv.size();

        /// Read the beginning of the duplicated sequence
        std::copy(inv.begin(), inv.end(), std::back_inserter(inv_begin));
        read_reference_seq(inv_begin, fasta_index, chrom_idx, var.pos + 1, padding_size);
        append_sv_tag_to_node(inv_begin);
        sv.model = "BREAKPOINT1";
        sv.related_sv = static_cast<int>(graph.SVs.size()) + 1; // Next SV
        graph.SVs.push_back(sv);

        /// Read the end of the duplicated sequence
        // Make sure we don't read before the chromosome starts!
        padding_size = std::min(padding_size, static_cast<std::size_t>(var2.pos));
        append_sv_tag_to_node(inv_end);
        read_reference_seq(inv_end, fasta_index, chrom_idx, var2.pos - padding_size + 1, padding_size);

        std::move(inv.begin(), inv.end(), std::back_inserter(inv_end));
        std::copy(sv.ins_seq.begin(), sv.ins_seq.end(), std::back_inserter(inv_end));

View source record →

📰 Research Paper
Loading…
⏳ Fetching content…