Question: Building multiple consensus sequences from multiple fasta sequences
gravatar for miss
6 months ago by
miss0 wrote:

I am new to NGS data analysis and I don't know how to overcome this problem: I have DNA-seq results from 12 different sites stored in one fasta file. Approximately every site has around 1000 reads, so I have 12000 reads in file but I don't know to which site any read belongs. I should build consensus sequences for all 12 sites without having reference sequence or knowing which reads belong to which site. Is it possible to do it and if yes, how?

ADD COMMENTlink modified 6 months ago by Jean-Karim Heriche21k • written 6 months ago by miss0
gravatar for Jean-Karim Heriche
6 months ago by
EMBL Heidelberg, Germany
Jean-Karim Heriche21k wrote:

An approach that doesn't require any prior knowledge is to cluster the sequences based on all pairwise sequence similarities/distances then compute a multiple sequence alignment for each cluster. Refinements can be made based on how much extra information is available.

ADD COMMENTlink written 6 months ago by Jean-Karim Heriche21k

Thank you! Any software tool you can recommend me for this kind of clustering?

ADD REPLYlink written 6 months ago by miss0

You'll need to identify which sequence alignment algorithm is relevant for your sequences, run the similarity computations in parallel then apply a clustering algorithm to the resulting distance or similarity matrix. All this can be done with a bit of scripting in your favourite language.

ADD REPLYlink written 6 months ago by Jean-Karim Heriche21k
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: 2177 users visited in the last hour