-
Notifications
You must be signed in to change notification settings - Fork 25
Description
Dear authors,
Thank you very much for creating this great tool.
I have encountered the following problem.
gene_pairs = gpf.find_all(n='', align_thr=0.2, n_genes=100)
I used gpffind_all to search for gene pairs enriched in pairs of cell type. In my understanding, this step returns homologous gene pairs that are specifically expressed in cell types similar to the one I am interested in. However, when I used the following code to test, I found that all the gene pairs output had a blast value of 0, indicating that they were not homologous genes. I have no idea how to solve this result. Can you provide me with some suggestions.
blast_score = sm.query_gene_pair('', '') print(blast_score)
part of the results:
-
gf_gfas1.m1.2879;sp_Spis_XP_022788977_1
BLAST: 0.0
Correlation: 0.0017043123726123018 -
gf_gfas1.m1.18621;sp_Spis_XP_022788977_1
BLAST: 0.0
Correlation: 0.04213641462550106 -
gf_gfas1.m1.10632;sp_Spis_XP_022786627_1
BLAST: 0.0
Correlation: 0.00658905127332004