Open IAlibay opened 7 months ago
This is a reasonably large change - pinging @MDAnalysis/coredevs
fasta2select()
https://github.com/MDAnalysis/mdanalysis/blob/b9f321c1d27efe5562421edf33d78ea4fdf82734/package/MDAnalysis/analysis/align.py#L1102 is a reasonably interesting function. We could probably just replace the clustalw call with Bio.Align.PairwiseAligner
.
I would say remove the clustalW wrapper, unsure how that intersects with @orbeckst's suggestions
I think we can remove clustalw functionality and just use Bio.Align.PairwiseAligner
internally. If we don't get regressions then I'd be ok just noting that we internally changed the alignment algorithm from clustalw to Biopython with no changes in our tests.
The alignment commandline wrappers via Biopython are deprecated and will be removed in a future version.
This means that the use of the clustalw wrapper in
align.py
will no longer work once this has been acheived.Options are: