This is a test version of Biostars. For the public version, visit https://www.biostars.org.
transcript and fasta sequence

Hello all! I have two file containing ten thousands of transcripts fasta sequences each with different ids and I am interested in finding common sequences between the two files. Somebody please help me as it is hindering my work. Thank you in appreciation

sequence assembly

ten thousands of transcripts fasta sequences each with different ids

Can you post an example? Is this a de novo assembly?

Actually this a de novo assembly produced by using two different softwares to minimizes the false positives

1 answer

Hello qudrat,

  1. If you are interested in IDENTICAL sequences, you can simply write a very short script to extract identical sequences in both files.

  2. You want to apply a "similarity" score, if so I strongly suggest using multi-aligners (BLAST or MUSCLE for instance) and then parse the results to have a global overview of similarity between sequences from your two different files.

  3. EDIT @genomax commentary: Use of assembly merge-tool.

Hi glihm, Actually I was thinking of sort but I do not know script writing. This is a de novo assembly using two different software and I am doing this to minimizes the false positives.

Your request is now clearer. The answer of @genomax is in this case well suited for your issue by using assembly merging.

Log in to answer this question.