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)
}