ó M̼]c@s d„ZdS(c Cs¿ddlm}m}m}m}m}m} m} g} g} g} g}xm|D]e}|jƒ}xP|dt |ƒ!D];}d|kryt |jdƒdƒ}| j |ƒqyqyWqSW|j dt t|| ƒƒjƒ| j ƒ| dt |ƒ} xyt|| ƒD]h\}}|jƒ}|djd ƒd}|d }tt |ƒƒ}|| krgPn| j |ƒq Wy=|| ƒ}| |ƒ}dj||ƒƒjd d ƒ}Wnd }d}nX| dt|ƒ}t|ƒt|ƒkrdt|ƒ|}n|dt|ƒkr;dt|ƒ|}ndtt |ƒƒdtt | ƒƒdt| dƒ}dj|d|||ƒjd d ƒdj| ƒ||||gƒdS(Niÿÿÿÿ(t taxidlist2LCAt find_parentstget_ranktfind_rankofparentstname1t find_familyt id_from_namei sNM:t:itkeyit|i t t_tNOMATCH_TAXID_NOT_FOUNDtFAMILY_NOT_FOUNDtNOMATCH_length_below_itNOMATCH_similarity_below_stothits:s_accepted-hits:s_Min-NM:s s (tget_LCA_functionsRRRRRRRtsplittlentinttappendtsorttdicttziptgettstrtjointreplacetfloat(tsamlinestlength_thresholdt id_thresholdtdistanceRRRRRRRtnmst thresholdtidstlseqstlinettexttitnmttaxidtseqtl_seqtlcatfamilytidptstats((sh/mounts/lovelace/aDNA/DNA-Sequence-Analysis/aDNA/software-to-sort/getLCA-master/get_LCA_from_samlines.pytget_LCA_from_samsH4   "      % <N(R0(((sh/mounts/lovelace/aDNA/DNA-Sequence-Analysis/aDNA/software-to-sort/getLCA-master/get_LCA_from_samlines.pyts