Loss of orthogroups after query_pangenes function?

jtlovell / GENESPACE

Other

180 stars 24 forks source link

Loss of orthogroups after query_pangenes function? #164

Closed FSU-pjtb closed 1 month ago

FSU-pjtb commented 1 month ago

In the REFERENCE_pangenes.txt output file, we have 27000 unique orthogroups ('og' column, syntenic and non-syntenic orthologs).

However after running query_pangenes(gsParam, refGenome = "REFERENCE"), we get an output with only 21000 unique orthogroups.

The query_pangenes table otherwise seems to be exactly what we need, so my question is what might be happening to the other 6000 orthogroups, and whether their loss is important.

Thank you!

FSU-pjtb commented 1 month ago

Nevermind - apparently the lost OGs are resulting from if (!showUnPlacedPgs) {out <- subset(out, !is.na(interpOrd))}, where showUnPlacedPgs defaults to FALSE.