Bootstrapping of unaligned sequeces
1
1
Entering edit mode
6.6 years ago
arit.mhp ▴ 10

Hi everyone, I have developed a new alignment free method to generate the phylogenetic tree from the Mitochondrial Genome sequences. I need to check the bootstrap supports of my newly generated trees. But I have doubt how can I bootstrap unaligned sequences and how to compute the final bootstrap scores for each clade of the tree. So in this regard, if you would suggest what should my approach to compute the confidence scores of our proposed method. Thanking you,

bootstrapping genome sequence • 2.0k views
ADD COMMENT
0
Entering edit mode

If your alignment-free approach is based on k-mers you can resample the original k-mers with replacement.

ADD REPLY
0
Entering edit mode

no its based on graphical method, representing a genome sequence graphically and computes the distance matrix.

ADD REPLY
0
Entering edit mode
6.6 years ago
kloetzl ★ 1.1k

Bootstrapping requires an alignments. However, there are other approaches to compute support values. See the following paper, for example: Support Values for Genome Phylogenies.

ADD COMMENT
0
Entering edit mode

Hi kloetzl, thank u for sharing this very interesting paper. I got a different view on bootstrapping. But I have a confusion, as I have developed a method which computes the distance matrix by taking the unaligned sequences as input, how this method can be applicable for my case. Waiting for your suggestions.

ADD REPLY
0
Entering edit mode

If your method produces a distance matrix in the PHYLIP file format, you are half way there. Download afra (mentioned in the above paper) and it will give you a tree with "bootstrap" values.

ADD REPLY
0
Entering edit mode

Basically, we have developed a alignment-free method using mitochondrial genome sequences which computes three distance matrices using three different methods. From each distance matrix, we get a tree by applying UPGMA. Final tree is formed applying the consensus based method. Is afra is applicable here?

ADD REPLY
0
Entering edit mode

afra should be close to what you want, except it does neighbor-joining, not UPGMA. So you might have to fix that in afra or write your own tool.

ADD REPLY
0
Entering edit mode

Ok... Thanks a lot.....

ADD REPLY

Login before adding your answer.

Traffic: 1807 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