Question: Help For Identifying Conserved Sequences (Phastcons)
gravatar for vasilislenis
7.1 years ago by
United Kingdom
vasilislenis130 wrote:


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’m starting with the pairwise alignment between human- mouse and human- rat, I’m 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’m 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’m thinking to blast these sequences to the mouse and rat whole genomes. Do you think that this procedure that I’m 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 in advance.

sequence • 2.4k views
ADD COMMENTlink modified 7.1 years ago by Jelena Aleksic910 • written 7.1 years ago by vasilislenis130
gravatar for Jelena Aleksic
7.1 years ago by
Cambridge, UK
Jelena Aleksic910 wrote:


I believe that those particular alignments have already been performed by other people, and should be available through Ensembl Compara ( Unless there's a particular reason you want to do them yourself?

ADD COMMENTlink written 7.1 years ago by Jelena Aleksic910

Thank u for your answer. I know, but I am a 1st year Ph.D and I need to do that because I must find a way to identify the conserved sequences in order to use them as anchors to speed up the alignment process.

ADD REPLYlink written 7.1 years ago by vasilislenis130
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: 919 users visited in the last hour