Hi,
I need to obtain the syntenic net alignments between several mammalian genomes
and I decided to generate them myself.
Given the information provided in the script doBlastzChainNet.pl (in particular
the doSyntenicNet subroutine), I located the
following script:
# filter net for synteny and create syntenic net mafs
netFilter -syn $tDb.$qDb.net.gz \\
| netSplit stdin synNet
chainSplit chain $tDb.$qDb.all.chain.gz
cd ..
mkdir $successDir
foreach f (axtChain/synNet/*.net)
netToAxt \$f axtChain/chain/\$f:t:r.chain \\
$defVars{'SEQ1_DIR'} $defVars{'SEQ2_DIR'} stdout \\
| axtSort stdin stdout \\
| axtToMaf -tPrefix=$tDb. -qPrefix=$qDb. stdin \\
$defVars{SEQ1_LEN} $defVars{SEQ2_LEN} \\
stdout \\
| gzip -c > mafSynNet/\$f:t:r:r:r:r:r.maf.gz
end
rm -fr $runDir/synNet
rm -fr $runDir/chain
Since the $tDb.$qDb.net.gz, $tDb.$qDb.all.chain.gz files can be downloaded from
the goldenPath directory of the download
repository, I figured that I can easily complete the generation procedure
following the script provided. (I do not need to represent
the alignments in MAF form so I stopped at step netToAxt and wrote the output
axt to chr?.$tDb.$qDb.synNet.axt, following UCSC's
naming convention).
In order to ensure that I have successfully generated the synNet.axt files, I
downloaded the following file from UCSC website:
chrY.hg19.rheMac2.synNet.axt.gz and compared it with the one generated by me. A
diff indicated some differences in the two files
for instance, in UCSC's version, such an entry exists which is absent in my
file, though the entries before and after it are identical.
92 chrY 873749 873891 chr3 158314483 158314625 - 9481
It seems that my version has filtered out more entries than UCSC's version. I
am just wondering whether you could kindly offer an
input to this? Thank you very much.
Best,
- Jia
_______________________________________________
Genome maillist - [email protected]
https://lists.soe.ucsc.edu/mailman/listinfo/genome
