Add Letters in a sequence

1 vue (au cours des 30 derniers jours)
A C
A C le 15 Nov 2022
Commenté : David Hill le 15 Nov 2022
I want to add Letters in a sequence
a b c d e f g h i j k l m n o p q r s t u v w x y z
a b c - - - ghi - - - - n o p q - - - - - - x y -
I got it to work with one gap but the sequence can have more than one gap.
I know the location n = [4 5 6 10 11 12 13 18 19 20 21 22 23 26] and the letters gap = [d e f j k l m r s t u v w z]
I do not know how to make it run for more in gap.
clc;clear;
% sequence we want to add because the sequence have gaps
% We have the amino acids however the position is unknow
% So the Letters are not include
% This a problem when there is over 300 amino acids
sequence = 'ABCGHIJ';
n = [4,5,6];
gap = 'DEF';
% loop
store = [];
% We want the gaps to be filled
for K = 1
% removing the letters where the gap needs
% to go
store = [sequence(n(K)):sequence(end)]
num = [n(end)+1:n(end)+1+numel(gap)]
for L = 1:numel(n)
% filling the gap
sequence(n(L)) = gap(L)
end
for M = 1:numel(num)
% Replacing the reoved letters
sequence(num(M)) = store(M)
end
end
sequence

Réponse acceptée

David Hill
David Hill le 15 Nov 2022
Modifié(e) : David Hill le 15 Nov 2022
l='abcdefghijklmnopqrstuvwxyz';
gap='defjklmrstuvwz'-96;
l(gap)='-'
l = 'abc---ghi----nopq------xy-'
  2 commentaires
A C
A C le 15 Nov 2022
That is neat. I did not explain it right.
We start with
l = 'abcghinopqxy'
We want to filled it with
gap = 'defjklmrstuvwz'
We know the locations
n = [4 5 6 10 11 12 13 18 19 20 21 22 23 26]
Results is
l = 'abcdefghijklmnopqrstuvwxyz'
David Hill
David Hill le 15 Nov 2022
l = 'abcghinopqxy';
gap = 'defjklmrstuvwz';
n = [4 5 6 10 11 12 13 18 19 20 21 22 23 26];
L=length(l)+length(gap);
r=zeros(1,L);
r(n)=gap;
r(~ismember(1:L,n))=l;
r=char(r)
r = 'abcdefghijklmnopqrstuvwxyz'

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

En savoir plus sur Logical dans Help Center et File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by