transcript and fasta sequence

Question

[qudrat]

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

[GenoMax]

Somebody please help me as it is hindering my work.

In what way?

How about this tool that merges assemblies? See more options here: High quality de novo transcriptome assembly rely on merging multiple assembly? Specifically dedupe.sh from BBMap should be very simple to use.

[st.ph.n]

ten thousands of transcripts fasta sequences each with different ids

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

[qudrat]

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

Answer 1

[glihm]

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.

[qudrat]

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.

[glihm]

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