SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
phastCons download from UCSC gene_x Bioinformatics 6 03-08-2013 01:49 PM
PhyloP or PhastCons Avro1986 Bioinformatics 0 12-21-2012 11:30 AM
phyloP or PhastCons Avro1986 Bioinformatics 1 12-18-2012 09:29 AM
phastCons scores behroozt Bioinformatics 1 07-29-2011 05:41 AM
HTSeq output not correlated with Cufflinks output... Help gen2prot Bioinformatics 5 01-31-2011 09:16 AM

Reply
 
Thread Tools
Old 04-16-2013, 09:32 PM   #1
JT Guan
Junior Member
 
Location: China

Join Date: Apr 2013
Posts: 4
Question About phastCons output

Hello Everyone,

I am using phastCons to get conservative elements by inputting human,mouse,rat multiple alignment file, for example chr1.hmr.maf.

And the .bed output is:
chr1 4424 4589 chr1.1-999145.1 43 +
chr1 5434 5511 chr1.1-999145.2 27 +
chr1 6312 6395 chr1.1-999145.3 21 +
chr1 16785 17057 chr1.1-999145.4 104 +
...

The result just specifies every conservative element's location in human chromosome.
But I also want to know every conservative element's location in mouse and rat chromosome.

Does anyone know how to get this kind of result by adding phastCons option or something else?

Thanks in advance for your help

Regards,
Guan
JT Guan is offline   Reply With Quote
Old 08-04-2013, 03:08 AM   #2
Vasilis
Junior Member
 
Location: UK

Join Date: Nov 2012
Posts: 3
Default Help for identifying conserved sequences

Hello,
I am trying to identify conserved regions among Human, mouse and rat.
I am really new to this field and I need your help.
I want to try with the human chromosome 22 as reference.
I知 starting with the pairwise alignment between human- mouse and human- rat, I知 creating the chains and nets, after that I convert the nets to maf files and finally with MULTIZ I am taking the multiple alignment.
I am thinking to use the phastCons with the final maf file but phastCons gives you the probabilities of each nucleotide to be in a conserved sequence.
I知 thinking to find the coordinates of the nucleotides that have probability more than a threshold (e.g 0.97) and extract the sequences from the ref genome (human).
As final step, I知 thinking to blast these sequences to the mouse and rat whole genomes.
Do you think that this procedure that I知 thinking is ok?
Is there a way phastCons to give me the sequences and not only the probabilities?
Also, and more important, which parameters I must use for phastCons?
Thank you very much I advance.
Vasilis is offline   Reply With Quote
Old 08-18-2015, 01:28 PM   #3
milo0615
Member
 
Location: Walnut, California

Join Date: Dec 2012
Posts: 39
Default

Hi JT Guan,

Were you able to use PhastCons to analyze your data? Please let me know. I am currently analyzing a multiple alignment with three plant species.


Thank you,

-Milo

Quote:
Originally Posted by JT Guan View Post
Hello Everyone,

I am using phastCons to get conservative elements by inputting human,mouse,rat multiple alignment file, for example chr1.hmr.maf.

And the .bed output is:
chr1 4424 4589 chr1.1-999145.1 43 +
chr1 5434 5511 chr1.1-999145.2 27 +
chr1 6312 6395 chr1.1-999145.3 21 +
chr1 16785 17057 chr1.1-999145.4 104 +
...

The result just specifies every conservative element's location in human chromosome.
But I also want to know every conservative element's location in mouse and rat chromosome.

Does anyone know how to get this kind of result by adding phastCons option or something else?

Thanks in advance for your help

Regards,
Guan
milo0615 is offline   Reply With Quote
Reply

Tags
conservative element, multiple alignment, phastcons

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 04:00 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO