Commit e38417b7 authored by Carsten Kemena's avatar Carsten Kemena

correct longest orf setting

parent bff9b2f6
Pipeline #562 passed with stage
in 1 minute and 4 seconds
......@@ -290,7 +290,7 @@ main(int argc, char const *argv[])
BSDL::SequenceSet<BSDL::Sequence<> > newSeqSet;
for (auto &seq : seqSet)
{
auto pair = BSDL::longestOrf(seq, {"ATG"}, {"TAA", "TAG", "TGA"}, 0);
auto pair = BSDL::longestOrf(seq, {"ATG"}, {"TAA", "TAG", "TGA"}, 1);
if (pair.first != pair.second)
{
if (pair.second < 0)
......@@ -349,8 +349,8 @@ main(int argc, char const *argv[])
{
auto domInfo = infoSet.find(domain.accession());
out << std::left << setw(nameLength) << seqSet[i].name() << " " << std::right << setw(6) << domain.start() << " " << setw(6) << domain.end() \
<< " " << std::right << setw(6) << domain.start() << " " << setw(6) << domain.end() << " " \
<< domain.accession() << std::left << " " << setw(15) << domInfo->second.name << " " << setw(11) << domInfo->second.type << " 0 0 0 0 0 0 " << domInfo->second.clan << "\n";
<< " " << std::right << setw(6) << domain.start() << " " << setw(6) << domain.end() << " " \
<< domain.accession() << std::left << " " << setw(15) << domInfo->second.name << " " << setw(11) << domInfo->second.type << " 0 0 0 0 0 0 " << domInfo->second.clan << "\n";
}
}
else
......@@ -358,8 +358,8 @@ main(int argc, char const *argv[])
for (auto domain : da)
{
auto domInfo = infoSet.find(domain.accession());
out << std::left << setw(nameLength) << seqSet[i].name() << " " << std::right << setw(6) << domain.start() << " " << setw(6) << domain.end() << " " <<
domain.accession() << std::left << " " << setw(15) << domInfo->second.name << " " << setw(11) << domInfo->second.type << " " << domInfo->second.clan << "\n";
out << std::left << setw(nameLength) << seqSet[i].name() << " " << std::right << setw(6) << domain.start() << " " << setw(6) << domain.end() << " " <<
domain.accession() << std::left << " " << setw(15) << domInfo->second.name << " " << setw(11) << domInfo->second.type << " " << domInfo->second.clan << "\n";
}
}
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment