code: def align(self, out_path=None):
          
          if out_path is None: out_path = self.prefix_path + '.aln'
          sh.muscle38("-in", self.path, "-out", out_path)
          return AlignedFASTA(out_path)
docs: We align the sequences in the fasta file with muscle.