Skip to content

Commit

Permalink
add warning
Browse files Browse the repository at this point in the history
  • Loading branch information
elileka committed May 1, 2024
1 parent 136c35d commit b2ed745
Showing 1 changed file with 10 additions and 0 deletions.
10 changes: 10 additions & 0 deletions src/exonpredictor/unitesetstofasta.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -385,6 +385,16 @@ int unitesetstofasta(int argn, const char **argv, const Command& command) {

// for the translated result
TranslateNucl translateNucl(static_cast<TranslateNucl::GenCode>(par.translationTable));
if ((par.translationTable != 1) && (par.lenSearchStart > 0)) {
Debug(Debug::WARNING) << "Selected translation table is not canonical and search for start is turned on. Please note that only ATG/atg is considered a start codon for the search!\n";
}
// for now, opting to use only ATG/atg as start codon because:
// In Euks ATG is super dominant, unlike proks
// The other two alternative codons in the standard table: TTG, CTG are NOT translated as M inside the gene
// std::vector<std::string> startCodonsInTable = translateNucl.getStartCodons();
// for (std::vector<std::string>::iterator t=startCodonsInTable.begin(); t!=startCodonsInTable.end(); ++t) {
// std::cout << *t << std::endl;
// }

Debug::Progress progress(predsPerContig.getSize());
#pragma omp parallel
Expand Down

0 comments on commit b2ed745

Please sign in to comment.