Advertisement
Guest User

Untitled

a guest
Jun 27th, 2019
75
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.08 KB | None | 0 0
  1. mouse = "MNGTEGPNFYVPFSNVTGVVRSPFEQPQYYLAEPWQFSMLAAYMFLLIVLGFPINFLTLYV
  2. TVQHKKLRTPLNYILLNLAVADLFMVFGGFTTTLYTSLHGYFVFGPTGCNLEGFFATLGGEIALWSLVVL
  3. AIERYVVVCKPMSNFRFGENHAIMGVVFTWIMALACAAPPLVGWSRYIPEGMQCSCGIDYYTLKPEVNNE
  4. SFVIYMFVVHFTIPMIVIFFCYGQLVFTVKEAAAQQQESATTQKAEKEVTRMVIIMVIFFLICWLPYASV
  5. AFYIFTHQGSNFGPIFMTLPAFFAKSSSIYNPVIYIMLNKQFRNCMLTTLCCGKNPLGDDDASATASKTE
  6. TSQVAPA";
  7. human = "MNGTEGPNFYVPFSNATGVVRSPFEYPQYYLAEPWQFSMLAAYMFLLIVLGFPINFLTLYV
  8. TVQHKKLRTPLNYILLNLAVADLFMVLGGFTSTLYTSLHGYFVFGPTGCNLEGFFATLGGEIALWSLVVL
  9. AIERYVVVCKPMSNFRFGENHAIMGVAFTWVMALACAAPPLAGWSRYIPEGLQCSCGIDYYTLKPEVNNE
  10. SFVIYMFVVHFTIPMIIIFFCYGQLVFTVKEAAAQQQESATTQKAEKEVTRMVIIMVIAFLICWVPYASV
  11. AFYIFTHQGSNFGPIFMTIPAFFAKSAAIYNPVIYIMMNKQFRNCMLTTICCGKNPLGDDEASATVSKTE
  12. TSQVAPA";
  13.  
  14. s = SequenceAlignment[mouse, human]
  15.  
  16. (* {MNGTEGPNFYVPFSN,{V,A},TGVVRSPFE,{Q,Y},
  17. PQYYLAEPWQFSMLAAYMFLLIVLGFPINFLTLYVTVQHKKLRTPLNYILLNLAVADLFMV,
  18. {F,L},GGFT,{T,S},TLYTSLHGYFVFGPTGCNLEGFFATLGGEIALWSLVVLAIERYVVVCKPMSNFRFGENHAIMGV,
  19. {V,A},FTW,{I,V},MALACAAPPL,{V,A},GWSRYIPEG,{M,L},
  20. QCSCGIDYYTLKPEVNNESFVIYMFVVHFTIPMI,{V,I},
  21. IFFCYGQLVFTVKEAAAQQQESATTQKAEKEVTRMVIIMVI,{F,A},
  22. FLICW,{L,V},PYASVAFYIFTHQGSNFGPIFMT,{L,I},PAFFAKS,
  23. {SS,AA},IYNPVIYIM,{L,M},NKQFRNCMLTT,{L,I},CCGKNPLGDD,{D,E},ASAT,{A,V},SKTETSQVAPA} *)
  24.  
  25. Manipulate[
  26. seqLength = Length[s];
  27. shownSize = 5;
  28. Grid[{If[ListQ[#], Column[#, Frame -> All, Spacings -> 1], #] & /@
  29. s[[x ;; x + shownSize]]}, Frame -> All], {x, 1,
  30. seqLength - shownSize, 1}]
  31.  
  32. aln1 = "";
  33. aln2 = "";
  34. ident = "";
  35. Do[{
  36. If[Length[n] == 0,
  37. {
  38. aln1 = aln1 <> n,
  39. aln2 = aln2 <> n,
  40. Do[ident = ident <> "*", {StringLength[n]}]
  41. },
  42. {
  43. While[StringLength[n[[1]]] != StringLength[n[[2]]], {
  44. If[StringLength[n[[1]]] < StringLength[n[[2]]],
  45. n[[1]] = n[[1]] <> "-",
  46. n[[2]] = n[[2]] <> "-"]
  47. }],
  48. Do[ident = ident <> " ", {Max[Map[StringLength, n]]}],
  49. aln1 = aln1 <> n[[1]],
  50. aln2 = aln2 <> n[[2]]
  51. }
  52. ]
  53. }, {n, s}]
  54.  
  55. Manipulate[
  56. seqLength2 = StringLength[ident];
  57. shownSize2 = 40;
  58. Grid[Characters[{ident, aln1, aln2}][[All, n ;; n + shownSize2]]]
  59. , {n, 1, seqLength2 - shownSize2, 1}]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement