busco5.beta not writing fna only faa files
For downstream analysis we were hoping to use the .fna nucleotide versions of the metaeuk output. The 3 directories that hold sequence output contain only the .faa amino acid versions. I noticed that in doing
grep "write_buscos_to_file" GenomeAnalysis.py -n
102: self.hmmer_runner.write_buscos_to_file(self.sequences_aa, self.sequences_nt)
254: self.hmmer_runner.write_buscos_to_file(self.sequences_aa, self.sequences_nt)
358: self.hmmer_runner.write_buscos_to_file(self.sequences_aa)
that the line on 358 does not have the self.sequences_nt
second arg. And this is in class GenomeAnalysisEukaryotesMetaeuk
. I don't suppose that just adding that argument would be all that was necessary since the output is in the codon.fas metaeuk specific file. It would be nice have the fna output for 5beta.
best,
jim henderson