(bitscore,length,orfQid,orfSid) = uqkey dpcpacbpcol.add_edge(nodeQ,nodeS,wt=bitscore) ############################################################################ if verbose: print pacbpcol, "bitscores:", print [ pacbporf.bitscore for pacbporf in dpcpacbpcol.pacbps.values() ] ############################################################################ # do some transformations on the pacbpcol pacbpcol.remove_low_connectivity_nodes(min_connectivity=cbg.node_count()-1) splittedCBGs = pacbpcol.find_fully_connected_subgraphs( edges=cbg.node_count()-1 , max_missing_edges=0 ) # convert to list of CBGs and do some transformations cbgList = ListOfCodingBlockGraphs(splittedCBGs,input={},crossdata={}) cbgList.remove_all_but_cbgs() cbgList.remove_cbgs_with_lt_nodes(cbg.node_count()) cbgList.harvest_pacbps_from_pacbpcollection(dpcpacbpcol) cbgList.remove_cbgs_without_omsr() cbgList.update_edge_weights_by_minimal_spanning_range() cbgList.order_graphlist_by_total_weight_and_identity() ############################################################################ if verbose: print "FScbgs (%s)" % len(cbgList) for fscbg in cbgList: print fscbg ############################################################################ if not cbgList: # no (better) frameshifted CBG return None