Entering edit mode
2.1 years ago
Chu
•
0
I have a fasta file of 50 sequences that looks like
ID1
> sequence1 : receptor sequence1
ID2
> sequence2 : receptor sequence2
...and so on, 50 entries total.
I want to create a fasta file that plots all 50 sequences to all 50 receptors, ending up with 2500 entries in a fasta file with the appropriate ID names:
ID1_ID1
> sequence1: receptor sequence1
ID1_ID2
> sequence1: receptor sequence2
ID1_ID3
> sequence1: receptor sequence3
ID1_ID4
> sequence1: receptor sequence4
... and so on for 2500 entries total
Thank you!
What do you mean by
you want to paste the sequences together?
I want to make all against all pairs of sequences and receptors.
sequence 1 would pair with all 50 receptors, sequence 2 would do the same, and so on.
Hope that makes sense
Pairwise2to compute alignments between eachThere are almost certainly faster ways to achieve this eg. with
parallel, but this would be quite straightforwardI don't think OP wants to align. This is most definitely an XY problem.
Chu : While trying to reformat your post I see that your examples now looks like
Not sure if they need to look like the following (with two sequences following each other to make a multi-fasta file).
Make changes as necessary (not sure what the word
receptoris doing in that line) and keep using101010button to format ascode.I think "receptor sequence X" is a single entity as is "sequence X". We're supposed to read it as:
OP is likely on an XY quest, it's not worth investing time in IMO.
Chu,
This looks like an XY problem. Please explain to us why you need this solved - in all probability, the task that needs this approach can be tackled in a completely different and probably more efficient manner altogether.