For the purpose of displaying the TreeGrafter alignment of the query sequence in the PANTHER tree viewer, we'll need an alignment that fits the same length as the MSA used in the viewer. We can use MAFFT for this and should add the functionality as an option:
- Take the query sequence and use HMMer to find its family hit as usual. This step also produces the alignment that we currently then use in the RAxML step.
- If the new MAFFT option is set, do not use the HMMer alignment from the previous step. Run MAFFT against the family and query sequence.
- Run RAxML using the MAFFT alignment.