INNER CODE UNIT · C++

seq

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

      seq = std::string(begin(vcf_record.alt) + 1, end(vcf_record.alt));

    std::ostringstream ss;

    if (seqan::length(vcf_record.info) > 0)
      ss << ";";

    ss << "SVTYPE=DEL;SVLEN=" << (-size_diff) << ";SVSIZE=" << (-size_diff)
       << ";END=" << (seq.size() + vcf_record.beginPos + 1 - size_diff);

    if (seq.size() > 0)
      ss << ";SEQ=" << seq;

    vcf_record.ref = vcf_record.ref[0]; // Only first matching base
    vcf_record.alt = "<DEL>";

    seqan::CharString new_info = ss.str();
    seqan::append(vcf_record.info, new_info);

View source record →