Question: Aligning genome wide coding sequences of two species for Ka/Ks calculation
gravatar for shreyasibiswas88
4.6 years ago by
United States
shreyasibiswas8830 wrote:


I have two fasta files each containing ~300,000 sequences of coding exons for one species. Both species are mapped to the same reference and hence, have same co-ordinates (synteny cannot be detected) example: File 1: >Rnor_chr1:268298874-268299113 TTACGCACACGGGGCACAGCCGCACTTGGTGGGCTTCT.....


I want to align each of the entries and then calculate Ka/Ks using PAML. Basically, pairwise alignment of 1000s of pairs. Looking at Muscle or Prank, I get the idea that they require an input file which contains both sequences to be aligned from the two species. That would mean, I need 300,000 input fasta files? Can anybody give me an idea how this can be solved?

Thank you.

next-gen alignment • 1.2k views
ADD COMMENTlink modified 4.5 years ago by Biostar ♦♦ 20 • written 4.6 years ago by shreyasibiswas8830
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: 1006 users visited in the last hour