How to find genes that do not have any homologous in the OrthoMCL output?
1
0
Entering edit mode
5.6 years ago
2351688559 ▴ 20

Hi!

Is there anyone using the OrthoMCL programs?

Recently, I met a problem, that is, I could not find the genome-specific genes that do not have any homologous(even paralogous) in the final output file of OrthoMCL.

I think those genes are of great importance and should not be ignored in the comparative genomic analysis.

I really want to find them out, but I don't have any idea how.

Many thanks!

genome alignment • 1.3k views
ADD COMMENT
0
Entering edit mode

What databases have you used?

ADD REPLY
0
Entering edit mode

MySQL, as proposed by the user guide.

ADD REPLY
0
Entering edit mode
5.6 years ago
5heikki 9.8k

For example with pyfasta:

pyfasta extract -help
Usage: extract some sequences from a fasta file. e.g.:
               pyfasta extract --fasta some.fasta --header at2g26540 at3g45640

Options:
  -h, --help     show this help message and exit
  --fasta=FASTA  path to the fasta file
  --header       include headers
  --exclude      extract all sequences EXCEPT those listed
  --file         if this flag is used, the sequences to extract are read from
                 the file specified in args
  --space        use the fasta identifier only up to the space as the key
ADD COMMENT

Login before adding your answer.

Traffic: 860 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6