Available software for checking all species have a 4d site for a particular codon in alignment?
1
1
Entering edit mode
8.6 years ago
RoryC ▴ 30

Hi guys,

I have a whole genome alignment for many species and a file containing coding sequence coordinates for one focal species, which I will use to extract 4d sites from the alignment and produce a phylogeny. What I would like to check is if what is a 4d site in the focal species actually corresponds in the alignment to 4d sites for each of the other species i.e. for each codon in the focal species which contains a 4d site, check each aligned codon for each species individually to ensure that they also have a 4d site at that position. Does anybody know of any existing software which can check this?

Cheers, Rory

genome gene alignment • 2.3k views
ADD COMMENT
0
Entering edit mode

Can't you just extract all the 4d sites into a new alignment and then check it column by column with something like sort?

ADD REPLY
0
Entering edit mode

Hi, thanks for the suggestion! I could see a way of writing a code for this, but I imagine it would be more complicated than that as you'd have to compare to a codon table to make sure the codon for each specific species has a 4d site in the 3rd position.

ADD REPLY
0
Entering edit mode

Hi RoryC, I have the same question as you. Have you solved this problem? Could you share the script or software to me? Many thanks in advance.

ADD REPLY
0
Entering edit mode
7.4 years ago
mtollis ▴ 30

If your focal species is a reference species in your alignment as in one constructed with MULTIZ, and it has available CDS annotations, you should try msa_view from the PHAST package.

http://compgen.cshl.edu/phast/

ADD COMMENT

Login before adding your answer.

Traffic: 985 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6