func (m MapConfig) MapFromFasta(pdbDb PDBDatabase, seqDb hhsuite.Database, queryFasta string) (*FragmentMap, error) { qseq, err := getOneFastaSequence(queryFasta) if err != nil { return nil, err } queryHHM, err := hhsuite.BuildHHM( hhsuite.HHBlitsDefault, hhsuite.HHMakePseudo, seqDb, queryFasta) if err != nil { return nil, err } return m.computeMap(pdbDb, qseq, queryHHM) }
func main() { inFasta := util.Arg(0) outHHM := util.Arg(1) hhblits := hhsuite.HHBlitsDefault hhmake := hhsuite.HHMakePseudo hhblits.Verbose = !flagQuiet hhmake.Verbose = !flagQuiet HHM, err := hhsuite.BuildHHM( hhblits, hhmake, util.FlagSeqDB, inFasta) util.Assert(err, "Error building HHM") util.Assert(hmm.WriteHHM(util.CreateFile(outHHM), HHM), "Error writing HHM '%s'", outHHM) }