Peter Bazeley

Hi Adam, As the author pointed it was a Perl issue. I installed MOCAT with another account and had no issues. So I would talk to your admin about having the correct version of Perl in your PATH, and use that edit as a last resort. Also the lack of indentation is just my laziness... -Peter

Hi. So in my quest to find abundances for the FMGs in each sample, it appears that there's no good way to associate COGs with the IGC catalog (or is there?). So the alternative appears to be to take the gene.prediction.assembly result sequences, extract the ones that were matched with the FMGs, and align those to the IGC catalog. Does that sound reasonable? Thanks, Peter

Not for taxanomic comparison but as a way to deal with compositional bias. Wouldn't this give a kind of average gene abundane per species, or am I way off?

Abundance scaled by FMG abundances?

Do you have more details of how this is calculated? Is it just a proportion of the gene length that has any read alignment or is it weighted by the number of reads? I'm wondering how to interpret these. Thanks. On Mon, Feb 12, 2018 at 4:21 AM, MOCAT <mocat-pipeline@users.sourceforge.net wrote: No. These are the horizontal gene coverages (what fraction of the gene has coverage), which is not related of horizontal gene transfers. On Mon, Feb 12, 2018 at 05:25 Peter Bazeley theoark@users.sourceforge.net...

Genome size normalization?

horizontal gene coverage

Hi, I ended up just commenting out line 51, and changing line 56 to: " my @line_parts = scalar @line; my @cogs; if (defined($line[$i])) { @cogs = split /[,|]/, $line[$i]; } else { @cogs = '' x $line_parts; } " and this runs fine. HTH, Peter