malacology/RGBEPP
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

fix: skip genes

Browse source
This commit is contained in:
kuoi 2024-09-19 13:48:13 +10:00
parent 56c0293057
commit 1227280755
1 changed files with 6 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

9
RGBEPP.d
View file

@ -473,9 +473,12 @@ void processTrimming(string[] ARG_G, string DirAlign, string DirTrim, string Pat
string inputFastaAA = DirAA_out ~ "/" ~ gene ~ ".fasta";
string inputBackTransNT = DirNT_out ~ "/" ~ gene ~ ".fasta";
string outputFastaNT = DirTrimNT ~ "/" ~ gene ~ ".fasta";
string[] cmdTrim = [PathTrimal, "-in", inputFastaAA, "-backtrans", inputBackTransNT, "-out", outputFastaNT, "-automated1"];
executeCommand(cmdTrim);
if (exists(inputFastaAA) && exists(inputBackTransNT)) {
string[] cmdTrim = [PathTrimal, "-in", inputFastaAA, "-backtrans", inputBackTransNT, "-out", outputFastaNT, "-automated1"];
executeCommand(cmdTrim);
} else {
writeln("Skipping gene: ", gene, " as files are missing.");
}
}
writeln("Trimming::End");