I have a lot of unassembled sequencing data and I want to build phylogeny for some genes from these data. I can retrieved the sequene by mapping reads to a reference sequence but if the two species are too divergent, reads will not match and I won't get the full CDS. Is there a better way to deal with that? Maybe I can do that to retrieve the well conserved part and then extend this region with reads that overlap (as if it was gap-filling)? What do you think of this approach? Do you know a software for such extension process? Thanks
Add a comment
|
1 Answer
$\begingroup$
$\endgroup$
2
That's a good question. I remember seeing a poster about making trees directly from reads, but unfortunately, I can not find any other record than this tweet. Perhaps you can try to email David, if he would be keen to share the implementation.
While searching for any information about David's work, I run into a GitHub repository called reads2trees. It's quite old and it's not super popular, so hard to say if it can be useful. But the name matches the problem description quite well.
I also remember there were some attempts (An example paper) to get phylogeny directly from kmers, but I have never tried that on my own. I also think these approached were focused on microbial phylogienies. So again, not sure if it's going to be useful for you.
That's all I got.
answered Apr 24, 2020 at 13:53
-
$\begingroup$ Thank you very much for your suggestions! I'm going to check that, however the ideal for me would be to retrieve the sequences not just the phylogenetic tree, for example to do selection analysis. $\endgroup$– LauraRApr 25, 2020 at 8:36
-
$\begingroup$ The "reads2trees" techniques should be able to give you sequences too (it's just a fancy name for a special case of reference-based assembly anyway). The kmer-based methods certainly won't. Good luck $\endgroup$ Apr 27, 2020 at 9:22