Question: Trying to use Siphy to get pi scores for each position
0
gravatar for Andy Lee
21 months ago by
Andy Lee10
Andy Lee10 wrote:

I am trying to use Siphy (Garber, M. et al, 2009), but am having some difficulties.

I am trying to get the log-odds pi scores for each position in one of the hg18 chromosomes. I used the following command line but the intended output does not seem to be generated:

java -jar bin/scaler_0.5.jar -task 7 -ref hg18 -in chr1.fa -mod primates.mod -minTreeLength 0.5 -out chr1.pi

The command line above results in the following in the terminal:

(((((((hg18:0.006591,)panTro2:0.006639):0.002184,)gorGor1:0.009411):0.009942,)ponAbe2:0.018342):0.014256,)rheMac2:0.036199):0.021496,)calJac1:0.066389):0.056911,)tarSyr1:0.135169):0.011307,)(micMur1:0.091452,)otoGar1:0.128984):0.035463)
base tree total length 0.6507350000000001
Alignment starts at -1 alignment length 247249719

 -3.30858  0.80573  2.63279  0.91909
  0.80573 -4.83099  0.92506  1.90856
  2.63279  0.92506 -6.08850  0.95536
  0.91909  1.90856  0.95536 -2.90928

chr1.pi file is 0 bytes. According to the documentations it's supposed to be a tab limited file with data corresponding to a number of columns.

For your information, I downloaded the hg18 chromosome 1 fasta file from UCSC and the primates.mod from this link (ftp://hgdownload.soe.ucsc.edu/goldenPath/hg18/phastCons44way/primates.mod). I have also tried chromosome 1 MAF file and obtained the same results (i.e. no output).

siphy broad conservation • 549 views
ADD COMMENTlink modified 21 months ago • written 21 months ago by Andy Lee10
0
gravatar for Andy Lee
21 months ago by
Andy Lee10
Andy Lee10 wrote:

Never mind. I figured this out. It turns out the maf file needs to match the mod file.

ADD COMMENTlink written 21 months ago by Andy Lee10
Please log in to add an answer.

Help
Access

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