This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Match trimmed sequences with full length reference set

I have a set of sequences that have been trimmed and that only consist of the ITS2 subregion. Is there a way I could match this set against the untrimmed dataset to retrieve a set of full length reads for every sequence in the trimmed dataset?

sequence dna

Are these fasta sequences? Are the headers in the trimmed dataset same as in the untrimmed dataset? If yes, you can use these headers to get the original sequences. Otherwise you can "grep" the trimmed sequences to identify the original sequences in the untrimmed file. But this will also fetch you extra sequences in case your trimmed sequences has more than one matches in the untrimmed file.

Yes they are in fasta format and the headers are the same. How would I go about to do that?

0 answers

No answers yet.

Log in to answer this question.