Question: PhastCons Scores for Human genes (CDS Regions)
gravatar for jeff.ohana
5.5 years ago by
United States
jeff.ohana10 wrote:

To whom it may concern,

I am looking to calculate the Phastcons scores for the Human coding genes based on FASTA alignments of 27 vertebrate genomes with Human for CDS regions.

Firstly, in my case, how should I type the Phastcons COMMAND ON Unix. I saw the following examples on Phastcons HOWTO:

The program can be run with a command of the form:

    phastCons [OPTIONS] alignment [cons.mod,]noncons.mod > scores.wig

For example,

    phastCons --target-coverage 0.25 --expected-length 12 \
        --msa-format MAF alignment.maf cons.mod,noncons.mod \
        > scores.wig


    phastCons --target-coverage 0.25 --expected-length 12 \
        --rho 0.4 --msa-format MAF alignment.maf noncons.mod \
        > scores.wig

I thought about :

phastCons "my alignment file" "??"  > scores.wig
"??"  = iNSTEAD of "cons.mod,noncons.mod," what should i put? I thought about the ".mod" file that was used to calculate the PhastCons scores based on multiple alignments of 27 vertebrate genomes with Human. But the difference is that it was not only comparing the CDS regions, but the whole genome. Any ideas?

Looking forward to hearing from you soon,

Jeff O.


ADD COMMENTlink modified 5.0 years ago by Biostar ♦♦ 20 • written 5.5 years ago by jeff.ohana10

Did you finally find out how to prepare the input files? If you did, please inform me. Thanks !!

ADD REPLYlink written 4.2 years ago by Bjoux0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1018 users visited in the last hour