Entering edit mode
12 months 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
Pairwise2
to 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
receptor
is doing in that line) and keep using101010
button 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.