bioroebe 0.10.80 → 0.12.12

Sign up to get free protection for your applications and to get access to all the features.

Potentially problematic release.


This version of bioroebe might be problematic. Click here for more details.

Files changed (242) hide show
  1. checksums.yaml +4 -4
  2. data/README.md +3612 -2781
  3. data/bin/bioroebe +7 -1
  4. data/bin/bioroebe_hash +7 -0
  5. data/bin/codon_to_aminoacid +1 -0
  6. data/bioroebe.gemspec +3 -3
  7. data/doc/README.gen +3612 -2742
  8. data/doc/quality_control/commandline_applications.md +3 -3
  9. data/doc/todo/bioroebe_java_todo.md +22 -0
  10. data/doc/todo/bioroebe_todo.md +2059 -2615
  11. data/lib/bioroebe/aminoacids/aminoacid_substitution.rb +1 -9
  12. data/lib/bioroebe/aminoacids/codon_percentage.rb +1 -9
  13. data/lib/bioroebe/aminoacids/deduce_aminoacid_sequence.rb +1 -9
  14. data/lib/bioroebe/aminoacids/display_aminoacid_table.rb +1 -0
  15. data/lib/bioroebe/aminoacids/show_hydrophobicity.rb +1 -6
  16. data/lib/bioroebe/base/colours_for_base/colours_for_base.rb +18 -8
  17. data/lib/bioroebe/base/commandline_application/commandline_arguments.rb +15 -11
  18. data/lib/bioroebe/base/commandline_application/misc.rb +66 -49
  19. data/lib/bioroebe/base/commandline_application/opn.rb +8 -8
  20. data/lib/bioroebe/base/commandline_application/reset.rb +3 -2
  21. data/lib/bioroebe/base/misc.rb +35 -0
  22. data/lib/bioroebe/base/prototype/misc.rb +11 -1
  23. data/lib/bioroebe/codon_tables/frequencies/parse_frequency_table.rb +2 -10
  24. data/lib/bioroebe/codons/codons.rb +1 -1
  25. data/lib/bioroebe/codons/convert_this_codon_to_that_aminoacid.rb +192 -58
  26. data/lib/bioroebe/codons/possible_codons_for_this_aminoacid.rb +1 -9
  27. data/lib/bioroebe/codons/show_codon_tables.rb +6 -2
  28. data/lib/bioroebe/codons/show_codon_usage.rb +15 -4
  29. data/lib/bioroebe/colours/rev.rb +4 -1
  30. data/lib/bioroebe/constants/aminoacids_and_proteins.rb +1 -0
  31. data/lib/bioroebe/constants/database_constants.rb +1 -1
  32. data/lib/bioroebe/constants/files_and_directories.rb +31 -4
  33. data/lib/bioroebe/constants/misc.rb +20 -0
  34. data/lib/bioroebe/conversions/dna_to_aminoacid_sequence.rb +58 -24
  35. data/lib/bioroebe/count/count_amount_of_aminoacids.rb +3 -2
  36. data/lib/bioroebe/count/count_amount_of_nucleotides.rb +3 -0
  37. data/lib/bioroebe/crystal/README.md +2 -0
  38. data/lib/bioroebe/crystal/to_rna.cr +19 -0
  39. data/lib/bioroebe/data/README.md +11 -8
  40. data/lib/bioroebe/data/electron_microscopy/pos_example.pos +396 -0
  41. data/lib/bioroebe/data/electron_microscopy/test_particles.star +36 -0
  42. data/lib/bioroebe/data/fasta/human/Homo_sapiens_hemoglobin_subunit_alpha_HBB_mRNA.fasta +9 -0
  43. data/lib/bioroebe/data/fasta/human/Homo_sapiens_hemoglobin_subunit_beta_HBB_mRNA.fasta +8 -0
  44. data/lib/bioroebe/data/fasta/human/README.md +2 -0
  45. data/lib/bioroebe/electron_microscopy/coordinate_analyzer.rb +15 -18
  46. data/lib/bioroebe/{fasta_and_fastq/parse_fasta/run.rb → electron_microscopy/electron_microscopy_module.rb} +16 -8
  47. data/lib/bioroebe/electron_microscopy/fix_pos_file.rb +1 -9
  48. data/lib/bioroebe/electron_microscopy/flipy.rb +83 -0
  49. data/lib/bioroebe/electron_microscopy/parse_coordinates.rb +2 -10
  50. data/lib/bioroebe/electron_microscopy/read_file_xmd.rb +1 -9
  51. data/lib/bioroebe/electron_microscopy/simple_star_file_generator.rb +4 -9
  52. data/lib/bioroebe/enzymes/has_this_restriction_enzyme.rb +10 -3
  53. data/lib/bioroebe/enzymes/restriction_enzyme.rb +23 -1
  54. data/lib/bioroebe/enzymes/restriction_enzymes/statistics.rb +65 -0
  55. data/lib/bioroebe/fasta_and_fastq/autocorrect_the_name_of_this_fasta_file.rb +1 -9
  56. data/lib/bioroebe/fasta_and_fastq/compact_fasta_file/compact_fasta_file.rb +7 -9
  57. data/lib/bioroebe/fasta_and_fastq/fasta_defline/fasta_defline.rb +1 -5
  58. data/lib/bioroebe/fasta_and_fastq/fasta_to_yaml/fasta_to_yaml.rb +81 -0
  59. data/lib/bioroebe/fasta_and_fastq/parse_fasta/parse_fasta.rb +1465 -7
  60. data/lib/bioroebe/fasta_and_fastq/return_fasta_subsection_of_this_file.rb +11 -2
  61. data/lib/bioroebe/fasta_and_fastq/show_fasta_headers.rb +27 -12
  62. data/lib/bioroebe/fasta_and_fastq/simplify_fasta_header/simplify_fasta_header.rb +1 -5
  63. data/lib/bioroebe/fasta_and_fastq/split_this_fasta_file_into_chromosomes/constants.rb +0 -5
  64. data/lib/bioroebe/genome/README.md +4 -0
  65. data/lib/bioroebe/genome/genome.rb +67 -0
  66. data/lib/bioroebe/genomes/genome_pattern.rb +3 -9
  67. data/lib/bioroebe/gui/gtk +1 -0
  68. data/lib/bioroebe/gui/gtk3/alignment/alignment.rb +73 -128
  69. data/lib/bioroebe/gui/gtk3/controller/controller.rb +45 -27
  70. data/lib/bioroebe/gui/gtk3/dna_to_aminoacid_widget/dna_to_aminoacid_widget.rb +76 -50
  71. data/lib/bioroebe/gui/gtk3/fasta_table_widget/fasta_table_widget.rb +99 -21
  72. data/lib/bioroebe/gui/gtk3/hamming_distance/hamming_distance.rb +42 -28
  73. data/lib/bioroebe/gui/gtk3/nucleotide_analyser/nucleotide_analyser.rb +119 -71
  74. data/lib/bioroebe/gui/gtk3/protein_to_DNA/protein_to_DNA.rb +18 -18
  75. data/lib/bioroebe/gui/gtk3/random_sequence/random_sequence.rb +19 -11
  76. data/lib/bioroebe/gui/jruby/alignment/alignment.rb +165 -0
  77. data/lib/bioroebe/gui/libui/alignment/alignment.rb +3 -1
  78. data/lib/bioroebe/gui/libui/controller/controller.rb +116 -0
  79. data/lib/bioroebe/gui/libui/random_sequence/random_sequence.rb +18 -2
  80. data/lib/bioroebe/gui/libui/show_codon_table/show_codon_table.rb +2 -0
  81. data/lib/bioroebe/gui/libui/three_to_one/three_to_one.rb +8 -6
  82. data/lib/bioroebe/gui/shared_code/alignment/alignment_module.rb +102 -0
  83. data/lib/bioroebe/gui/shared_code/levensthein_distance/levensthein_distance_module.rb +18 -16
  84. data/lib/bioroebe/gui/shared_code/protein_to_DNA/protein_to_DNA_module.rb +14 -14
  85. data/lib/bioroebe/gui/swing/three_to_one/ThreeToOne$1.class +0 -0
  86. data/lib/bioroebe/gui/swing/three_to_one/ThreeToOne$CloseListener.class +0 -0
  87. data/lib/bioroebe/gui/swing/three_to_one/ThreeToOne.class +0 -0
  88. data/lib/bioroebe/gui/swing/three_to_one/ThreeToOne.java +141 -0
  89. data/lib/bioroebe/images/FORWARD_PRIMER.png +0 -0
  90. data/lib/bioroebe/images/REVERSE_PRIMER.png +0 -0
  91. data/lib/bioroebe/java/README.md +4 -0
  92. data/lib/bioroebe/java/bioroebe/Sequence.java +25 -1
  93. data/lib/bioroebe/java/bioroebe/src/main/java/bioroebe/Base.class +0 -0
  94. data/lib/bioroebe/java/bioroebe/{Base.java → src/main/java/bioroebe/Base.java} +15 -2
  95. data/lib/bioroebe/java/bioroebe/{BisulfiteTreatment.class → src/main/java/bioroebe/BisulfiteTreatment.class} +0 -0
  96. data/lib/bioroebe/java/bioroebe/{Codons.class → src/main/java/bioroebe/Codons.class} +0 -0
  97. data/lib/bioroebe/java/bioroebe/src/main/java/bioroebe/Codons.java +34 -0
  98. data/lib/bioroebe/java/bioroebe/src/main/java/bioroebe/Commandline.class +0 -0
  99. data/lib/bioroebe/java/bioroebe/src/main/java/bioroebe/Commandline.java +101 -0
  100. data/lib/bioroebe/java/bioroebe/{GenerateRandomDnaSequence.class → src/main/java/bioroebe/GenerateRandomDnaSequence.class} +0 -0
  101. data/lib/bioroebe/java/bioroebe/{GenerateRandomDnaSequence.java → src/main/java/bioroebe/GenerateRandomDnaSequence.java} +8 -2
  102. data/lib/bioroebe/java/bioroebe/{IsPalindrome.class → src/main/java/bioroebe/IsPalindrome.class} +0 -0
  103. data/lib/bioroebe/java/bioroebe/{IsPalindrome.java → src/main/java/bioroebe/IsPalindrome.java} +5 -1
  104. data/lib/bioroebe/java/bioroebe/src/main/java/bioroebe/PartnerNucleotide.class +0 -0
  105. data/lib/bioroebe/java/bioroebe/src/main/java/bioroebe/PartnerNucleotide.java +56 -0
  106. data/lib/bioroebe/java/bioroebe/{RemoveFile.java → src/main/java/bioroebe/RemoveFile.java} +10 -4
  107. data/lib/bioroebe/java/bioroebe/{RemoveNumbers.class → src/main/java/bioroebe/RemoveNumbers.class} +0 -0
  108. data/lib/bioroebe/java/bioroebe/{RemoveNumbers.java → src/main/java/bioroebe/RemoveNumbers.java} +1 -0
  109. data/lib/bioroebe/java/bioroebe/{ToCamelcase.class → src/main/java/bioroebe/ToCamelcase.class} +0 -0
  110. data/lib/bioroebe/java/bioroebe/{ToCamelcase.java → src/main/java/bioroebe/ToCamelcase.java} +3 -3
  111. data/lib/bioroebe/java/bioroebe/src/main/java/bioroebe/ToRNA.class +0 -0
  112. data/lib/bioroebe/java/bioroebe/src/main/java/bioroebe/ToRNA.java +42 -0
  113. data/lib/bioroebe/java/bioroebe/src/main/java/bioroebe/toplevel_methods/BaseComposition.class +0 -0
  114. data/lib/bioroebe/java/bioroebe/src/main/java/bioroebe/toplevel_methods/BaseComposition.java +75 -0
  115. data/lib/bioroebe/misc/ruler.rb +11 -2
  116. data/lib/bioroebe/nucleotides/most_likely_nucleotide_sequence_for_this_aminoacid_sequence.rb +1 -9
  117. data/lib/bioroebe/nucleotides/show_nucleotide_sequence.rb +7 -7
  118. data/lib/bioroebe/parsers/genbank_parser.rb +347 -26
  119. data/lib/bioroebe/parsers/gff.rb +1 -9
  120. data/lib/bioroebe/patterns/scan_for_repeat.rb +1 -5
  121. data/lib/bioroebe/pdb/fetch_fasta_sequence_from_pdb.rb +1 -9
  122. data/lib/bioroebe/pdb/parse_mmCIF_file.rb +1 -9
  123. data/lib/bioroebe/pdb/parse_pdb_file.rb +4 -10
  124. data/lib/bioroebe/project/project.rb +1 -1
  125. data/lib/bioroebe/python/README.md +1 -0
  126. data/lib/bioroebe/python/__pycache__/mymodule.cpython-39.pyc +0 -0
  127. data/lib/bioroebe/python/gui/gtk3/all_in_one.css +4 -0
  128. data/lib/bioroebe/python/gui/gtk3/all_in_one.py +59 -0
  129. data/lib/bioroebe/python/gui/gtk3/widget1.py +20 -0
  130. data/lib/bioroebe/python/gui/tkinter/all_in_one.py +91 -0
  131. data/lib/bioroebe/python/mymodule.py +8 -0
  132. data/lib/bioroebe/python/protein_to_dna.py +33 -0
  133. data/lib/bioroebe/python/shell/shell.py +19 -0
  134. data/lib/bioroebe/python/to_rna.py +14 -0
  135. data/lib/bioroebe/python/toplevel_methods/esystem.py +12 -0
  136. data/lib/bioroebe/python/toplevel_methods/open_in_browser.py +20 -0
  137. data/lib/bioroebe/python/toplevel_methods/palindromes.py +42 -0
  138. data/lib/bioroebe/python/toplevel_methods/rds.py +13 -0
  139. data/lib/bioroebe/python/toplevel_methods/shuffleseq.py +23 -0
  140. data/lib/bioroebe/python/toplevel_methods/three_delimiter.py +37 -0
  141. data/lib/bioroebe/python/toplevel_methods/time_and_date.py +43 -0
  142. data/lib/bioroebe/python/toplevel_methods/to_camelcase.py +21 -0
  143. data/lib/bioroebe/requires/require_the_bioroebe_project.rb +3 -1
  144. data/lib/bioroebe/sequence/alignment.rb +14 -4
  145. data/lib/bioroebe/sequence/dna.rb +1 -0
  146. data/lib/bioroebe/sequence/nucleotide_module/nucleotide_module.rb +28 -25
  147. data/lib/bioroebe/sequence/protein.rb +105 -3
  148. data/lib/bioroebe/sequence/sequence.rb +87 -21
  149. data/lib/bioroebe/shell/menu.rb +3829 -3714
  150. data/lib/bioroebe/shell/misc.rb +59 -4307
  151. data/lib/bioroebe/shell/readline/readline.rb +1 -1
  152. data/lib/bioroebe/shell/shell.rb +11255 -28
  153. data/lib/bioroebe/siRNA/siRNA.rb +81 -1
  154. data/lib/bioroebe/string_matching/find_longest_substring.rb +3 -2
  155. data/lib/bioroebe/string_matching/hamming_distance.rb +1 -9
  156. data/lib/bioroebe/taxonomy/class_methods.rb +3 -8
  157. data/lib/bioroebe/taxonomy/constants.rb +4 -3
  158. data/lib/bioroebe/taxonomy/edit.rb +2 -1
  159. data/lib/bioroebe/taxonomy/help/help.rb +10 -10
  160. data/lib/bioroebe/taxonomy/help/helpline.rb +2 -2
  161. data/lib/bioroebe/taxonomy/info/check_available.rb +15 -9
  162. data/lib/bioroebe/taxonomy/info/info.rb +18 -11
  163. data/lib/bioroebe/taxonomy/info/is_dna.rb +46 -36
  164. data/lib/bioroebe/taxonomy/interactive.rb +140 -104
  165. data/lib/bioroebe/taxonomy/menu.rb +27 -18
  166. data/lib/bioroebe/taxonomy/parse_fasta.rb +3 -1
  167. data/lib/bioroebe/taxonomy/shared.rb +1 -0
  168. data/lib/bioroebe/taxonomy/taxonomy.rb +1 -0
  169. data/lib/bioroebe/toplevel_methods/aminoacids_and_proteins.rb +31 -24
  170. data/lib/bioroebe/toplevel_methods/colourize_related_methods.rb +164 -0
  171. data/lib/bioroebe/toplevel_methods/databases.rb +1 -1
  172. data/lib/bioroebe/toplevel_methods/digest.rb +18 -8
  173. data/lib/bioroebe/toplevel_methods/fasta_and_fastq.rb +107 -63
  174. data/lib/bioroebe/toplevel_methods/file_and_directory_related_actions.rb +14 -2
  175. data/lib/bioroebe/toplevel_methods/frequencies.rb +8 -1
  176. data/lib/bioroebe/toplevel_methods/misc.rb +142 -12
  177. data/lib/bioroebe/toplevel_methods/nucleotides.rb +118 -46
  178. data/lib/bioroebe/toplevel_methods/open_in_browser.rb +2 -0
  179. data/lib/bioroebe/toplevel_methods/palindromes.rb +1 -2
  180. data/lib/bioroebe/toplevel_methods/taxonomy.rb +2 -2
  181. data/lib/bioroebe/toplevel_methods/to_camelcase.rb +5 -0
  182. data/lib/bioroebe/utility_scripts/align_open_reading_frames.rb +1 -9
  183. data/lib/bioroebe/utility_scripts/check_for_mismatches/check_for_mismatches.rb +1 -9
  184. data/lib/bioroebe/utility_scripts/compacter.rb +1 -9
  185. data/lib/bioroebe/utility_scripts/compseq/compseq.rb +1 -9
  186. data/lib/bioroebe/utility_scripts/consensus_sequence.rb +6 -6
  187. data/lib/bioroebe/utility_scripts/create_batch_entrez_file.rb +1 -9
  188. data/lib/bioroebe/utility_scripts/dot_alignment.rb +1 -9
  189. data/lib/bioroebe/utility_scripts/move_file_to_its_correct_location.rb +1 -4
  190. data/lib/bioroebe/utility_scripts/parse_taxonomy.rb +2 -2
  191. data/lib/bioroebe/utility_scripts/showorf/constants.rb +0 -5
  192. data/lib/bioroebe/utility_scripts/showorf/reset.rb +1 -4
  193. data/lib/bioroebe/version/version.rb +2 -2
  194. data/lib/bioroebe/www/embeddable_interface.rb +103 -54
  195. data/lib/bioroebe/www/sinatra/sinatra.rb +186 -70
  196. data/lib/bioroebe/yaml/aminoacids/amino_acids_long_name_to_one_letter.yml +2 -2
  197. data/lib/bioroebe/yaml/configuration/browser.yml +1 -1
  198. data/lib/bioroebe/yaml/configuration/temp_dir.yml +1 -1
  199. data/lib/bioroebe/yaml/genomes/README.md +3 -4
  200. data/lib/bioroebe/yaml/restriction_enzymes/restriction_enzymes.yml +27 -27
  201. metadata +81 -64
  202. data/doc/setup.rb +0 -1655
  203. data/lib/bioroebe/fasta_and_fastq/parse_fasta/constants.rb +0 -50
  204. data/lib/bioroebe/fasta_and_fastq/parse_fasta/initialize.rb +0 -86
  205. data/lib/bioroebe/fasta_and_fastq/parse_fasta/menu.rb +0 -117
  206. data/lib/bioroebe/fasta_and_fastq/parse_fasta/misc.rb +0 -981
  207. data/lib/bioroebe/fasta_and_fastq/parse_fasta/report.rb +0 -156
  208. data/lib/bioroebe/fasta_and_fastq/parse_fasta/reset.rb +0 -128
  209. data/lib/bioroebe/genbank/genbank_parser.rb +0 -291
  210. data/lib/bioroebe/java/bioroebe/Base.class +0 -0
  211. data/lib/bioroebe/java/bioroebe/Codons.java +0 -22
  212. data/lib/bioroebe/java/bioroebe/PartnerNucleotide.class +0 -0
  213. data/lib/bioroebe/java/bioroebe/PartnerNucleotide.java +0 -19
  214. data/lib/bioroebe/java/bioroebe.jar +0 -0
  215. data/lib/bioroebe/shell/add.rb +0 -108
  216. data/lib/bioroebe/shell/assign.rb +0 -360
  217. data/lib/bioroebe/shell/chop_and_cut.rb +0 -281
  218. data/lib/bioroebe/shell/constants.rb +0 -166
  219. data/lib/bioroebe/shell/download.rb +0 -335
  220. data/lib/bioroebe/shell/enable_and_disable.rb +0 -158
  221. data/lib/bioroebe/shell/enzymes.rb +0 -310
  222. data/lib/bioroebe/shell/fasta.rb +0 -345
  223. data/lib/bioroebe/shell/gtk.rb +0 -76
  224. data/lib/bioroebe/shell/history.rb +0 -132
  225. data/lib/bioroebe/shell/initialize.rb +0 -217
  226. data/lib/bioroebe/shell/loop.rb +0 -74
  227. data/lib/bioroebe/shell/prompt.rb +0 -107
  228. data/lib/bioroebe/shell/random.rb +0 -289
  229. data/lib/bioroebe/shell/reset.rb +0 -335
  230. data/lib/bioroebe/shell/scan_and_parse.rb +0 -135
  231. data/lib/bioroebe/shell/search.rb +0 -337
  232. data/lib/bioroebe/shell/sequences.rb +0 -200
  233. data/lib/bioroebe/shell/show_report_and_display.rb +0 -2901
  234. data/lib/bioroebe/shell/startup.rb +0 -127
  235. data/lib/bioroebe/shell/taxonomy.rb +0 -14
  236. data/lib/bioroebe/shell/tk.rb +0 -23
  237. data/lib/bioroebe/shell/user_input.rb +0 -88
  238. data/lib/bioroebe/shell/xorg.rb +0 -45
  239. /data/lib/bioroebe/java/bioroebe/{BisulfiteTreatment.java → src/main/java/bioroebe/BisulfiteTreatment.java} +0 -0
  240. /data/lib/bioroebe/java/bioroebe/{Esystem.class → src/main/java/bioroebe/Esystem.class} +0 -0
  241. /data/lib/bioroebe/java/bioroebe/{Esystem.java → src/main/java/bioroebe/Esystem.java} +0 -0
  242. /data/lib/bioroebe/java/bioroebe/{RemoveFile.class → src/main/java/bioroebe/RemoveFile.class} +0 -0
@@ -1,2823 +1,2267 @@
1
- -------------------------------------------------------------------------------
2
- (1) → https://biopython.org/DIST/docs/tutorial/Tutorial.html#sec15
3
-
4
- ^^^ work through the above, also integrate it + write docs
5
-
6
- https://raw.githubusercontent.com/biopython/biopython/master/Doc/examples/ls_orchid.fasta
7
-
8
- -------------------------------------------------------------------------------
1
+ --------------------------------------------------------------------------------
2
+ Bioroebe::RestrictionEnzymes::Statistics.show
3
+ ^^^ improve these
4
+ and then add it to the documentation.
5
+ --------------------------------------------------------------------------------
6
+ use glimmer + nebula for widgets
7
+ ^^^
8
+ improve the nucleotide sequence analyser
9
+
10
+ --------------------------------------------------------------------------------
11
+ add the possibility to show what the effect of enzymes
12
+ are
13
+ AND inhibitors of enzymes. perhaps bioroebe can be
14
+ used in system biology one day
15
+ --------------------------------------------------------------------------------
16
+ (1) Bioroebe::Sequence.new('AGCTTAGCGTACAGCTACGACGTAGTCTGACGA').cut_with? :AluI
17
+ ^^^ support this API and document it too
18
+ --------------------------------------------------------------------------------
9
19
  (2) → integrate electrno microscopy slowly and also add documentation
10
20
  about this AS YOU GO!!!!!
11
21
  ^^^ yup add more of it
12
- -------------------------------------------------------------------------------
22
+ --------------------------------------------------------------------------------
13
23
  (3) → Add save session support
14
24
  to reload our last activity completely ...
15
- hmmm..
16
- This has to be well designed...
17
- Perhaps before we do so, we will add some
18
- class that anylizes what we have
19
- call it:
20
-
21
- class AnalyseLocalDataset
22
-
23
- And it is called when the bioshell is
24
- called. Can be enabled and disabled.
25
- AND document it then.
26
- The idea is to provide additional information
27
- upon startup of the bioroebe shell.
28
- This is in preparation for save-session support.
29
-
30
-
31
- -------------------------------------------------------------------------------
32
- (5) → Lys-Asp-Glu-Leu
33
-
34
- if i.include?('-') and Bioroebe.is_in_the_three_letter_code?(i)
35
- end
36
-
37
- - Lys-Asp-Glu-Leu-COO-
38
-
39
- Lys-Asp-Glu-Leu
40
-
41
- ^^^ this is "retention in lumen of ER"
42
- find this too
43
-
44
- BUT!
45
-
46
- we must verify it
47
-
48
- ^^ yep this is also called KDEL
49
- https://en.wikipedia.org/wiki/KDEL_(amino_acid_sequence)
50
- -------------------------------------------------------------------------------
51
- (6) → Add "orthologs". this shall show us the top 25 orthologs or
25
+ hmmm..
26
+ This has to be well designed...
27
+ Perhaps before we do so, we will add some
28
+ class that anylizes what we have
29
+ call it:
30
+ class AnalyseLocalDataset
31
+ And it is called when the bioshell is
32
+ called. Can be enabled and disabled.
33
+ AND document it then.
34
+ The idea is to provide additional information
35
+ upon startup of the bioroebe shell.
36
+ This is in preparation for save-session support.
37
+ --------------------------------------------------------------------------------
38
+ (4) Lys-Asp-Glu-Leu
39
+ if i.include?('-') and Bioroebe.is_in_the_three_letter_code?(i)
40
+ end
41
+ - Lys-Asp-Glu-Leu-COO-
42
+ Lys-Asp-Glu-Leu
43
+ ^^^ this is "retention in lumen of ER"
44
+ find this too
45
+ BUT!
46
+ we must verify it
47
+ ^^ yep this is also called KDEL
48
+ https://en.wikipedia.org/wiki/KDEL_(amino_acid_sequence)
49
+ --------------------------------------------------------------------------------
50
+ (5) → Add "orthologs". this shall show us the top 25 orthologs or
52
51
  something. In the bioshell? Hmm. Not sure yet.
53
- -------------------------------------------------------------------------------
54
- (7) → clone the functionality of this:
55
-
56
- http://www.kazusa.or.jp/codon/cgi-bin/countcodon.cgi
57
- http://www.kazusa.or.jp/codon/countcodon.html
58
-
59
- In other words, create a class that can generate such an output.
60
- ^^^ This is now done.
61
- Then add this to a GUI as well as the www output.
62
- ^^^ This still has to be done, though. We will use a ruby-gtk3
63
- widget first. And sinatra output too.
64
- AND document it as well
65
-
66
- -------------------------------------------------------------------------------
67
- (8) → SARS genom analyisere in bioroebe
52
+ --------------------------------------------------------------------------------
53
+ (6) → clone the functionality of this:
54
+ http://www.kazusa.or.jp/codon/cgi-bin/countcodon.cgi
55
+ http://www.kazusa.or.jp/codon/countcodon.html
56
+ In other words, create a class that can generate such an output.
57
+ ^^^ This is now done.
58
+ Then add this to a GUI as well as the www output.
59
+ ^^^ This still has to be done, though. We will use a ruby-gtk3
60
+ widget first. And sinatra output too.
61
+ AND document it as well
62
+ --------------------------------------------------------------------------------
63
+ (7) SARS genom analyisere in bioroebe
68
64
  eventuell auch graphisch
69
-
70
- Gibt es neue GUIs die wir kombinieren könnten? Hmmm.
71
- -------------------------------------------------------------------------------
72
- (9) → In bioroebe, generate that .ps thingy graphical thing from the
65
+ Gibt es neue GUIs die wir kombinieren könnten? Hmmm.
66
+ --------------------------------------------------------------------------------
67
+ (8) → In bioroebe, generate that .ps thingy graphical thing from the
73
68
  vienna RNA tutorial. Hmmm.
74
-
75
- https://www.tbi.univie.ac.at/RNA/tutorial/
76
- -------------------------------------------------------------------------------
77
- (1) → get insulin squence frmo NCBI
78
- human
69
+ https://www.tbi.univie.ac.at/RNA/tutorial/
70
+ --------------------------------------------------------------------------------
71
+ (9) → get insulin squence frmo NCBI
72
+ human
79
73
  then apply trypsin onto it
80
74
  and try it like this:
81
-
82
- trypsin --insulin
83
- ^^^
84
- also document it then. well .....
85
-
86
- Also add:
87
- insulin?
88
- ^^^ to show it
89
- Hmm. Perhaps also auto-download or something.
90
-
91
- -------------------------------------------------------------------------------
92
- (1) → in bioroebe: UAG?
75
+ trypsin --insulin
76
+ ^^^
77
+ also document it then. well .....
78
+ Also add:
79
+ insulin?
80
+ ^^^ to show it
81
+ Hmm. Perhaps also auto-download or something.
82
+ --------------------------------------------------------------------------------
83
+ (10) in bioroebe: UAG?
93
84
  ^^^ show all stop codons with that in the bioshell
94
85
  all UAG sequences... hmm. and TAG?
95
86
  Finish that.
96
- ..........................................................................
97
- (1) → The position of a symbol in a string is the total number of
87
+ --------------------------------------------------------------------------------
88
+ (11) → The position of a symbol in a string is the total number of
98
89
  symbols found to its left, including itself (e.g., the positions
99
90
  of all occurrences of 'U' in "AUGCUUCAGAAAGGUCUUACG" are 2, 5,
100
91
  6, 15, 17, and 18). The symbol at position i
101
92
  of s is denoted by s[i].
102
-
103
- ^^^ add a solution there, a toplevel API
104
- !!!!!
105
- -------------------------------------------------------------------------------
106
- (1) → http://bioruby.org/rdoc/Bio/Blast.html
93
+ ^^^ add a solution there, a toplevel API
94
+ !!!!!
95
+ --------------------------------------------------------------------------------
96
+ (12) → http://bioruby.org/rdoc/Bio/Blast.html
107
97
  ^^^ add support for BLAST
108
- ..........................................................................
109
- (1) → add: parse_pdb()
98
+ --------------------------------------------------------------------------------
99
+ (13) → add: parse_pdb()
110
100
  With this we shall just show some info, about a given
111
101
  .pdb file at hand.
112
102
  Also make it commandline based too + bioshell variant
113
103
  here, and a sinatra interface once this all works.
114
104
  Don't forget to document it!!!!!
115
105
  ^^^ and google a bit how others do that
116
- ..........................................................................
117
- (2) → pdb 1a6m
106
+ --------------------------------------------------------------------------------
107
+ (14) → pdb 1a6m
118
108
  ^^^ download this when that is used in the bioshell; we also have
119
- to use the download directory for this, so make sure that
120
- we do.
109
+ to use the download directory for this, so make sure that
110
+ we do.
121
111
  ^^^ And then, also document this clearly.
122
- -------------------------------------------------------------------------------
123
- (3) show_string
124
- ^^^ slowly port this ... find out differences
125
- then unify into one method. right now we used
126
- two or something.
127
- -------------------------------------------------------------------------------
128
- (4) → Try to see if we can integrate this into our GUI:
129
-
130
- https://cdn.snapgene.com/assets/7.6.11/assets/images/snapgene/homepage/homepage-hero.png
131
- -------------------------------------------------------------------------------
132
- (5) → Scan for leucine zipper!
133
-
112
+ --------------------------------------------------------------------------------
113
+ (15) show_string
114
+ ^^^ slowly port this ... find out differences
115
+ then unify into one method. right now we used
116
+ two or something.
117
+ --------------------------------------------------------------------------------
118
+ (16) → Try to see if we can integrate this into our GUI:
119
+ https://cdn.snapgene.com/assets/7.6.11/assets/images/snapgene/homepage/homepage-hero.png
120
+ --------------------------------------------------------------------------------
121
+ (17) → Scan for leucine zipper!
134
122
  This is ~25% implemented. We need to double-check what
135
123
  exactly is a leucine zipper.
136
- ..........................................................................
137
- (6) → Extend the sinatra-interface for the Rosalind task,
124
+ --------------------------------------------------------------------------------
125
+ (18) → Extend the sinatra-interface for the Rosalind task,
138
126
  perhaps add a sub-link to show which parts are solved
139
127
  as-is. Hmm. I am not continuing on this though.
140
- ^^^^
141
- well - make rosalind anew again or something.
142
-
143
- ...........................................................................
144
- (7) - Add a blast interface; both via the web-interface, GUI,
128
+ ^^^^
129
+ well - make rosalind anew again or something.
130
+ --------------------------------------------------------------------------------
131
+ (19) → Add a blast interface; both via the web-interface, GUI,
145
132
  and also from the commandline.
146
- -------------------------------------------------------------------------------
147
- (8) - Write a tutorial about primer design.
133
+ --------------------------------------------------------------------------------
134
+ (20) Write a tutorial about primer design.
148
135
  also make sure that the GUI has support for this.
149
- ..........................................................................
150
- (9) - In the documentation examples, show some exampls for how to work
136
+ --------------------------------------------------------------------------------
137
+ (21) In the documentation examples, show some exampls for how to work
151
138
  with different organisms.
152
- ..........................................................................
153
- (10) - In the bioshell, if "stop?" is issued, then the colouring isn't
154
- correct. It currently does not show any result. This has to
155
- be fixed.
156
- ..........................................................................
157
- (11) https://www.rubydoc.info/gems/biomart
158
- ^^^ integrate biomart
159
-
160
- p biomart.list_datasets
161
- p biomart.datasets?
162
- -------------------------------------------------------------------------------
163
- (12) Add Trypsin und Trypsinogen sequences, both as FASTA
164
- but also as shortcut via the commandline such as:
139
+ --------------------------------------------------------------------------------
140
+ (22) In the bioshell, if "stop?" is issued, then the colouring isn't
141
+ correct. It currently does not show any result. This has to
142
+ be fixed.
143
+ --------------------------------------------------------------------------------
144
+ (23) https://www.rubydoc.info/gems/biomart
145
+ ^^^ integrate biomart
146
+ p biomart.list_datasets
147
+ p biomart.datasets?
148
+ --------------------------------------------------------------------------------
149
+ (24) → Add Trypsin und Trypsinogen sequences, both as FASTA
150
+ but also as shortcut via the commandline such as:
165
151
  show_orf :trypsine
166
152
  show_orf :trypsin
167
- Or something like this; and document it as well.
168
- -------------------------------------------------------------------------------
169
- (13) → 1..60
170
-
171
- setdna 57
172
- append stop
173
- 1..60
174
-
175
- Next showing the nucleotides 1370 to 1462 (including 1370 and 1462).
176
- The length of the fragment will be 93 nucleotides.
177
- 5' - ATGTGCAGTCAGGTGAATTTATTGAAAAATTTGAGGCTCCTGGTGGTGCAAATCAAAGAACTGCTCCTCAGTGGATGTTGCCTTTACTTCTAG - 3'
178
- ^^^ hier beim colourize, wenn das letzte codon ein STOP codon ist
179
- dann colourizen wir das auch.
180
- -------------------------------------------------------------------------------
181
- (14) → MG1655
153
+ Or something like this; and document it as well.
154
+ --------------------------------------------------------------------------------
155
+ (25) → 1..60
156
+ setdna 57
157
+ append stop
158
+ 1..60
159
+ Next showing the nucleotides 1370 to 1462 (including 1370 and 1462).
160
+ The length of the fragment will be 93 nucleotides.
161
+ 5' - ATGTGCAGTCAGGTGAATTTATTGAAAAATTTGAGGCTCCTGGTGGTGCAAATCAAAGAACTGCTCCTCAGTGGATGTTGCCTTTACTTCTAG - 3'
162
+ ^^^ hier beim colourize, wenn das letzte codon ein STOP codon ist
163
+ dann colourizen wir das auch.
164
+ --------------------------------------------------------------------------------
165
+ (26) MG1655
182
166
  ^^^ input this to download the sequence. Also show it to the user.
183
- -------------------------------------------------------------------------------
184
- (15) → extend virus-information into the bioroebe project.
185
- -------------------------------------------------------------------------------
186
- (16) → Add a way to analyse the chemical structure of all
167
+ --------------------------------------------------------------------------------
168
+ (27) → extend virus-information into the bioroebe project.
169
+ --------------------------------------------------------------------------------
170
+ (28) → Add a way to analyse the chemical structure of all
187
171
  aminoacids. We wish to show the chemical formula.
188
-
189
172
  E. g. if we input:
190
-
191
173
  "phenylalanin"
192
-
193
- Then the C9N should be shown, of its -R part.
194
- ^^^ H wird aber rausgelöscht und O ebenso.
195
- wtf?
196
- I don't understand why it removes H and 0 so perhaps
197
- dont remove that part. But still show the -R.
198
-
199
- -------------------------------------------------------------------------------
200
- (17) FIX THE COLOURIZATION BUG; THIS ONE TRIGGERED THE WHOLE
201
- REWRITE AFTER ALL!
202
- -------------------------------------------------------------------------------
203
- (18) FIX TAXONOMY related-problems AS WELL
204
- ^^^^^^ AND DOCUMENT THIS related-problems.
205
- -------------------------------------------------------------------------------
206
- (19) Do note that z will then be a String, not a sequence object anymore.
207
- (This may be subject to change in the future, but for now, aka
208
- **February 2020**, it is that way.)
209
- ^^^^
210
- -------------------------------------------------------------------------------
211
- (20) ^^^ colours are appended. That should not be the case!
212
- ADD SOMETHING NEW ... some todo entries
213
- and some python tool
214
- -------------------------------------------------------------------------------
215
- (21) → rewrite the whole project anew
174
+ Then the C9N should be shown, of its -R part.
175
+ ^^^ H wird aber rausgelöscht und O ebenso.
176
+ wtf?
177
+ I don't understand why it removes H and 0 so perhaps
178
+ dont remove that part. But still show the -R.
179
+ --------------------------------------------------------------------------------
180
+ (29) → FIX THE COLOURIZATION BUG; THIS ONE TRIGGERED THE WHOLE
181
+ REWRITE AFTER ALL!
182
+ --------------------------------------------------------------------------------
183
+ (30) FIX TAXONOMY related-problems AS WELL
184
+ ^^^^^^ AND DOCUMENT THIS related-problems.
185
+ --------------------------------------------------------------------------------
186
+ (31) Do note that z will then be a String, not a sequence object anymore.
187
+ (This may be subject to change in the future, but for now, aka
188
+ **February 2020**, it is that way.)
189
+ ^^^^
190
+ --------------------------------------------------------------------------------
191
+ (32) → ^^^ colours are appended. That should not be the case!
192
+ ADD SOMETHING NEW ... some todo entries
193
+ and some python tool
194
+ --------------------------------------------------------------------------------
195
+ (33) rewrite the whole project anew
216
196
  - improve the documentation
217
- - focus on class Protein first and add
218
- all_dna_combinations or somethingl ike
219
- that, as well as:
220
- .backtrans
221
- .reverse_translate
222
- -------------------------------------------------------------------------------
223
- (22) → AND THEN test on windows as well.
224
- ^^^^^^^^^^^^^^
225
- -------------------------------------------------------------------------------
226
- (23)
227
- Reduced alphabets for proteins | [not implemented yet]
228
- ^^^ check this as well
229
-
230
- require 'bioroebe/base/commandline_application/aminoacids.rb'
231
- ^^^ verify whether we need this really in the commandline
232
- application part. Or perhaps it should be part of
233
- Base.
234
-
235
- - in bioroebe, test sequel for taxonomy...
236
-
237
-
238
- - LEARN FUCKING JAVA; combine it with bioroebe though. work through
239
- bioroebe and as you go, also write related-problems into java.
240
- Add restriction thingy complete in java AND bioroebe
241
- show how many will be cut.
242
- Improve this also in bioroebe at the same time
243
- add a GUI in swing too, also in ruby-gtk.
244
- document publish IMPROVE
245
- for now just show how many segments we will
246
- generate.
247
-
248
- First focus on bioroebe.
249
- ^^^^^^^^^^^^^^^
250
-
251
- -
252
- efetch "https://www.ncbi.nlm.nih.gov/gene/744779"
253
- ^^^ test this. again
254
-
255
- -------------------------------------------------------------------------------
256
- (25) fix tk-levensthein
257
- -------------------------------------------------------------------------------
258
- (26) → rewrite the whole project anew
197
+ - focus on class Protein first and add
198
+ all_dna_combinations or somethingl ike
199
+ that, as well as:
200
+ .backtrans
201
+ .reverse_translate
202
+ --------------------------------------------------------------------------------
203
+ (34) →
204
+ Reduced alphabets for proteins | [not implemented yet]
205
+ ^^^ check this as well
206
+ require 'bioroebe/base/commandline_application/aminoacids.rb'
207
+ ^^^ verify whether we need this really in the commandline
208
+ application part. Or perhaps it should be part of
209
+ Base.
210
+ - in bioroebe, test sequel for taxonomy...
211
+ - LEARN FUCKING JAVA; combine it with bioroebe though. work through
212
+ bioroebe and as you go, also write related-problems into java.
213
+ Add restriction thingy complete in java AND bioroebe
214
+ show how many will be cut.
215
+ Improve this also in bioroebe at the same time
216
+ add a GUI in swing too, also in ruby-gtk.
217
+ document publish IMPROVE
218
+ for now just show how many segments we will
219
+ generate.
220
+ First focus on bioroebe.
221
+ ^^^^^^^^^^^^^^^
222
+ -
223
+ efetch "https://www.ncbi.nlm.nih.gov/gene/744779"
224
+ ^^^ test this. again
225
+ --------------------------------------------------------------------------------
226
+ (35) → fix tk-levensthein
227
+ --------------------------------------------------------------------------------
228
+ (36) rewrite the whole project anew
259
229
  - improve the documentation
260
- - rework the WHOLE tutorial as well
261
- - focus on class Protein first and add
262
- all_dna_combinations or somethingl ike
263
- that
264
- .backtrans
265
- .reverse_translate
266
- -------------------------------------------------------------------------------
267
- (27) → analyze /Depot/Temp/Bioroebe/1CEZ.pdb
268
-
269
- ^^^
270
- support this. Already works half-way, we started writing a pdb parser.
271
- this should work in general, for .fasta files as well.
272
- ..........................................................................
273
- (28) → SINATRA STUFF:
230
+ - rework the WHOLE tutorial as well
231
+ - focus on class Protein first and add
232
+ all_dna_combinations or somethingl ike
233
+ that
234
+ .backtrans
235
+ .reverse_translate
236
+ --------------------------------------------------------------------------------
237
+ (37) → analyze /Depot/Temp/Bioroebe/1CEZ.pdb
238
+ ^^^
239
+ support this. Already works half-way, we started writing a pdb parser.
240
+ this should work in general, for .fasta files as well.
241
+ --------------------------------------------------------------------------------
242
+ (38) → SINATRA STUFF:
274
243
  FIX AND EXTEND SINATRA IN BIOROEBE.
275
244
  extend it too.
276
-
277
245
  http://localhost:4567/random_aminoacids
278
246
  ^^^ add a form there
279
-
280
247
  add emboss show orf or so
281
248
  and special-dispaly on sinatra kaa
282
249
  where the nucleotide sequence has numbers
283
250
  ^^^
284
- -------------------------------------------------------------------------------
285
- (29) pick any virus and begin to amass tons of data; and then when done
286
- also connect this into a GUI for use therein.
287
-
288
- See:
251
+ --------------------------------------------------------------------------------
252
+ (39) pick any virus and begin to amass tons of data; and then when done
253
+ also connect this into a GUI for use therein.
254
+ See:
289
255
  https://raw.githubusercontent.com/labsquare/fastQt/master/screenshot.gif
290
- ^^^^^^^
291
-
292
- begin with a circovirus
293
-
294
- ^^^^^^^
295
-
296
- DOCUMENT THIS AS YOU GO
297
- research about circovirus too
298
-
299
- https://www.ncbi.nlm.nih.gov/nuccore/NC_038391.1
300
-
301
-
302
-
303
-
304
-
305
- -------------------------------------------------------------------------------
306
- (1) → Fix:
307
-
308
- require 'bioroebe/toplevel_methods/open_reading_frames.rb'
309
-
310
- Something is wrong; it returns regions that contain
311
- a stop codon, which can not be true.
312
-
313
- -------------------------------------------------------------------------------
314
- (3) → Fix: extend glycovirology parts
315
- seek stuff in viral genomes
316
- -------------------------------------------------------------------------------
317
- (4) →
318
-
319
- seq = Bio::Sequence::NA.new("atgcatgcaaaaaaa")
320
- puts seq
321
- puts seq.complement
322
- puts seq.subseq(3,8)
323
- puts seq.subseq(3,8).complement #wont work
324
- p seq.gc_percent
325
- p (100 - seq.gc_percent) # at_percent
326
- p seq.composition
327
- puts seq.translate
328
- puts seq.translate(2)
329
- puts seq.translate(1,11)
330
- puts seq.translate.codes
331
- puts seq.translate.names
332
- puts seq.translate.composition
333
- puts seq.translate.molecular_weight
334
- puts seq.complement.translate
335
- ^^^ make sure this works
336
- seq = Bioroebe::Sequence.new("atgcatgcaaaaaaa")
337
- puts seq
338
- puts seq.complement
339
- -------------------------------------------------------------------------------
340
- (5) →
341
- make sure we have a good fasta-showing widget
342
- show how many nucleotides are
343
- AND add support to modify this as-is
344
- ^^^^
345
- -------------------------------------------------------------------------------
346
- (6) → In BioRoebe:
347
-
348
- Add a table showing how compatible bioroebe is compared to the other
349
- bio-projects, staring with biophp.
350
-
351
- Also show the status how much is complete in each,
352
- including Bio (ruby-bio) the main ruby project here.
353
- And add a table which functionality is implemented
354
- in Java already.
355
- -------------------------------------------------------------------------------
356
- (7) →
357
- ********************************************************************************
358
- Was passiert wenn wir das Lambda-Genom mit EcoRI behandeln?
359
- ********************************************************************************
360
-
361
- Es entstehen 3 chromosomale Fragmente.
362
- ^^^ dies testen
363
- also zerst lambda genom herunterladen.
364
- download lambda
365
- download lambda_genome
366
- ^^^
367
- dann:
368
-
369
- Bioroebe.digest_this_dna(:lambda_genome, with: :EcoRI)
370
- Bioroebe.digest_this_dna("/root/Bioroebe/fasta/NC_001416.1_Enterobacteria_phage_lambda_complete_genome.fasta", with: :EcoRI)
371
-
372
-
373
- ^^^ test this API and document it as well.
374
- ^^^ and say how many fragments will be created in this CIRCULAR
375
- DNA.
376
- ^^^ this now works kind of ... but it must be better
377
- documented and we must test this with more data.
378
- -------------------------------------------------------------------------------
379
- (8) → add the bioroebe logo to sinatra, but as appropriate size,
256
+ ^^^^^^^
257
+ begin with a circovirus
258
+ ^^^^^^^
259
+ DOCUMENT THIS AS YOU GO
260
+ research about circovirus too
261
+ https://www.ncbi.nlm.nih.gov/nuccore/NC_038391.1
262
+ --------------------------------------------------------------------------------
263
+ (40) → Fix:
264
+ require 'bioroebe/toplevel_methods/open_reading_frames.rb'
265
+ Something is wrong; it returns regions that contain
266
+ a stop codon, which can not be true.
267
+ --------------------------------------------------------------------------------
268
+ (41) → Fix: extend glycovirology parts
269
+ seek stuff in viral genomes
270
+ --------------------------------------------------------------------------------
271
+ (42) →
272
+ seq = Bio::Sequence::NA.new("atgcatgcaaaaaaa")
273
+ puts seq
274
+ puts seq.complement
275
+ puts seq.subseq(3,8)
276
+ puts seq.subseq(3,8).complement #wont work
277
+ p seq.gc_percent
278
+ p (100 - seq.gc_percent) # at_percent
279
+ p seq.composition
280
+ puts seq.translate
281
+ puts seq.translate(2)
282
+ puts seq.translate(1,11)
283
+ puts seq.translate.codes
284
+ puts seq.translate.names
285
+ puts seq.translate.composition
286
+ puts seq.translate.molecular_weight
287
+ puts seq.complement.translate
288
+ ^^^ make sure this works
289
+ seq = Bioroebe::Sequence.new("atgcatgcaaaaaaa")
290
+ puts seq
291
+ puts seq.complement
292
+ --------------------------------------------------------------------------------
293
+ (43) → In BioRoebe:
294
+ Add a table showing how compatible bioroebe is compared to the other
295
+ bio-projects, staring with biophp.
296
+ Also show the status how much is complete in each,
297
+ including Bio (ruby-bio) the main ruby project here.
298
+ And add a table which functionality is implemented
299
+ in Java already.
300
+ --------------------------------------------------------------------------------
301
+ (44)
302
+ ********************************************************************************
303
+ Was passiert wenn wir das Lambda-Genom mit EcoRI behandeln?
304
+ ********************************************************************************
305
+ Es entstehen 3 chromosomale Fragmente.
306
+ ^^^ dies testen
307
+ also zerst lambda genom herunterladen.
308
+ download lambda
309
+ download lambda_genome
310
+ ^^^
311
+ dann:
312
+ Bioroebe.digest_this_dna(:lambda_genome, with: :EcoRI)
313
+ Bioroebe.digest_this_dna("/root/Bioroebe/fasta/NC_001416.1_Enterobacteria_phage_lambda_complete_genome.fasta", with: :EcoRI)
314
+ ^^^ test this API and document it as well.
315
+ ^^^ and say how many fragments will be created in this CIRCULAR
316
+ DNA.
317
+ ^^^ this now works kind of ... but it must be better
318
+ documented and we must test this with more data.
319
+ --------------------------------------------------------------------------------
320
+ (45) add the bioroebe logo to sinatra, but as appropriate size,
380
321
  via base64. perhaps width 50 or so. need to determine
381
322
  which size fits here.
382
- -------------------------------------------------------------------------------
383
- (9) → Integrate http://nc2.neb.com/NEBcutter2/cutshow.php?name=ffe1d68e-
384
-
323
+ --------------------------------------------------------------------------------
324
+ (46) → Integrate http://nc2.neb.com/NEBcutter2/cutshow.php?name=ffe1d68e-
385
325
  in particular the visual part.
386
- -------------------------------------------------------------------------------
387
- (10) → https://international.neb.com/products/r0196-ncii#Product%20Information
388
- ^^^ autogenerate such an image, aka restriction cutting enzyme
389
- to indicate the target sequence.
390
- -------------------------------------------------------------------------------
391
- (6) → how to do codon optimiation in e.coli? bioroebe must support this!
392
-
393
- we must first get a display which codon is very commonly used in
394
- E. coli, from some remote site ... japanese site I think.
395
-
396
- then, we analyse all possibilities.
397
-
398
- and then we look which codons may be improvable - display
399
- them on the commandline
400
-
401
- class: OptimizeCodons.new(of_this_sequence)
402
- -------------------------------------------------------------------------------
403
- (7) → Molekulare Grösse von "Ubiquitin"? "8.5 kd".
326
+ --------------------------------------------------------------------------------
327
+ (47) → https://international.neb.com/products/r0196-ncii#Product%20Information
328
+ ^^^ autogenerate such an image, aka restriction cutting enzyme
329
+ to indicate the target sequence.
330
+ --------------------------------------------------------------------------------
331
+ (48) → how to do codon optimiation in e.coli? bioroebe must support this!
332
+ we must first get a display which codon is very commonly used in
333
+ E. coli, from some remote site ... japanese site I think.
334
+ then, we analyse all possibilities.
335
+ and then we look which codons may be improvable - display
336
+ them on the commandline
337
+ class: OptimizeCodons.new(of_this_sequence)
338
+ --------------------------------------------------------------------------------
339
+ (49) Molekulare Grösse von "Ubiquitin"? "8.5 kd".
404
340
  ^^^ das sollte automatisch ausgerechnet werden
405
- -------------------------------------------------------------------------------
406
- (8) → taxonomy !!!!!!!!!!!!!!!!!!
407
- -------------------------------------------------------------------------------
408
- (9) → Given a list of gene names that I would like to get chromosome/position
341
+ --------------------------------------------------------------------------------
342
+ (50) → taxonomy !!!!!!!!!!!!!!!!!!
343
+ --------------------------------------------------------------------------------
344
+ (51) → Given a list of gene names that I would like to get chromosome/position
409
345
  information for (in mm10). Is there some service online where I can
410
346
  paste this list? ^^^ enable this
411
- -------------------------------------------------------------------------------
412
- (10) → Show the frequency of codons in different tables
413
-
414
- This works quite ok, but right now the approach is to store
415
- this in a .yml file which is not ideal.
416
-
417
- Thus, we have to add two things:
418
- - The ability to store this into a SQL database
419
- - The ability to batch-download all of these codons,
420
- which first requires that we have a way to obtain all
421
- taxonomic ids.
422
- -------------------------------------------------------------------------------
423
- (11) → Add a way in bioroebe to store a gene into a yaml file
424
- or so, and to also load it up again. Perhaps simplify
425
- this automatically. Need some ways to describe that.
426
- -------------------------------------------------------------------------------
427
- (12) → Make bioroebe very useful from the www, no matter if via sinatra
347
+ --------------------------------------------------------------------------------
348
+ (52) → Make bioroebe very useful from the www, no matter if via sinatra
428
349
  or rails. It should be a tool-set project on the www as well.
429
- -------------------------------------------------------------------------------
430
- (13) → Suppose you have a GenBank file which you want to turn into a
431
- Fasta file. For example, lets consider the file cor6_6.gb
432
- which is included in the Biopython unit tests under the
350
+ --------------------------------------------------------------------------------
351
+ (53) → Suppose you have a GenBank file which you want to turn into a
352
+ Fasta file. For example, lets consider the file cor6_6.gb
353
+ which is included in the Biopython unit tests under the
433
354
  GenBank directory.
434
-
435
-
436
- need to check that this is equivalent, think about the API
437
- document it and then remove this entry.
438
-
439
- ^^^ also build a GUI for this.
440
- call it format-converter or so
441
- the GUI works somewhat but needs to be polished up.
442
- THEN THIS CAN BE REMOVED!!!!!!!
443
-
444
- -------------------------------------------------------------------------------
445
- (14) → Wir brauchen eine table wo wir die starken promotoren verschiedener
355
+ need to check that this is equivalent, think about the API
356
+ document it and then remove this entry.
357
+ ^^^ also build a GUI for this.
358
+ call it format-converter or so
359
+ the GUI works somewhat but needs to be polished up.
360
+ THEN THIS CAN BE REMOVED!!!!!!!
361
+ --------------------------------------------------------------------------------
362
+ (54) Wir brauchen eine table wo wir die starken promotoren verschiedener
446
363
  Organismen zusammenstellen und vergleichen können.
447
-
448
- strong_promoters.yml
449
- -------------------------------------------------------------------------------
450
- (15) add:
451
- start position of exons
452
- and show the sequence based on that file
453
-
454
- Normally there's a "gene" entry for each gene, so:
455
- awk 'BEGIN{FS="\t"; OFS="\t"}{if($3 == "gene") print $1, $4, $5}' foo.gtf
456
-
457
- -------------------------------------------------------------------------------
458
- (16) → also add 30-33 to aminoacids hmmm difficult.
459
- -------------------------------------------------------------------------------
460
- (17) → http://bioinformatics.oxfordjournals.org/content/18/8/1135
364
+ strong_promoters.yml
365
+ --------------------------------------------------------------------------------
366
+ (55) → add:
367
+ start position of exons
368
+ and show the sequence based on that file
369
+ Normally there's a "gene" entry for each gene, so:
370
+ awk 'BEGIN{FS="\t"; OFS="\t"}{if($3 == "gene") print $1, $4, $5}' foo.gtf
371
+ --------------------------------------------------------------------------------
372
+ (56) also add 30-33 to aminoacids hmmm difficult.
373
+ --------------------------------------------------------------------------------
374
+ (57) → http://bioinformatics.oxfordjournals.org/content/18/8/1135
461
375
  "TFBS: Computational framework for transcription factor
462
- binding site analysis"
463
- study the above and see if it can be included
464
- into bioroebe
465
-
466
- http://tfbs.genereg.net/
467
- -------------------------------------------------------------------------------
468
- (18) → They include trypsin, chymotrypsin, thrombin, plasmin, papain and factor Xa.
376
+ binding site analysis"
377
+ study the above and see if it can be included
378
+ into bioroebe
379
+ http://tfbs.genereg.net/
380
+ --------------------------------------------------------------------------------
381
+ (58) → They include trypsin, chymotrypsin, thrombin, plasmin, papain and factor Xa.
469
382
  ^^^ provide means to identify where they cut,
470
- and show this then by simualting a digest.
471
- return an array with the starting aminoacids.
472
- also document this on bioroebe todo
383
+ and show this then by simualting a digest.
384
+ return an array with the starting aminoacids.
385
+ also document this on bioroebe todo
473
386
  this is done via digestion/digestions
474
387
  but it's not quite perfect yet.
475
- -------------------------------------------------------------------------------
476
- (19) → a) add a commandline way to generate a random protein
477
- with a specified length and then display it on the
388
+ --------------------------------------------------------------------------------
389
+ (59) → a) add a commandline way to generate a random protein
390
+ with a specified length and then display it on the
478
391
  commandline [DONE] !!!
479
-
480
- bioroebe --random-aminoacids=33
481
- bioroebe --n-aminoacids=33
482
- sinatra:
483
-
484
- random_aminoacids/33 [Implemented! 23.09.2019]
485
-
486
- also added the gtk-GUI code here; needs to be
487
- documented briefly, then this part is completelty
488
- done. contiu on random_aminoacids: in particular
489
- add a gtk_entry that specifie, no, a spin button
490
- that states how many no... an entry or so
491
- to state how many aminoacids to generate
492
- randomly
493
-
494
- b) add a way to generate a cDNA sequence from such a
495
- protein and view all possible sequences from that
496
- sequence. ^^^
497
-
498
- Enable this BOTH from the commandline AND from the
499
- interactive variant and from sinatra! Hmmmm.
500
-
501
- -------------------------------------------------------------------------------
502
- (1) → add an option to design a
503
-
504
- degenerate primer
505
- -------------------------------------------------------------------------------
506
- (2) Add upcase to sequences and ensure that it works; also document it
507
- internally and in the .pdf tutorial
508
- what does that mean? upcase as method? hmmm.
509
-
510
- ..........................................................................
511
- (1) → http://www.biomart.org/other/user-docs.pdf
512
- ^^^ work through this
513
- ^^^ integrate the old .cgi part and improve as you go
514
- ..........................................................................
515
- (1) → Access geninfo numbers easily.
392
+ bioroebe --random-aminoacids=33
393
+ bioroebe --n-aminoacids=33
394
+ sinatra:
395
+ random_aminoacids/33 [Implemented! 23.09.2019]
396
+ also added the gtk-GUI code here; needs to be
397
+ documented briefly, then this part is completelty
398
+ done. contiu on random_aminoacids: in particular
399
+ add a gtk_entry that specifie, no, a spin button
400
+ that states how many no... an entry or so
401
+ to state how many aminoacids to generate
402
+ randomly
403
+ b) add a way to generate a cDNA sequence from such a
404
+ protein and view all possible sequences from that
405
+ sequence. ^^^
406
+ Enable this BOTH from the commandline AND from the
407
+ interactive variant and from sinatra! Hmmmm.
408
+ --------------------------------------------------------------------------------
409
+ (60) → add an option to design a
410
+ degenerate primer
411
+ --------------------------------------------------------------------------------
412
+ (61) Add upcase to sequences and ensure that it works; also document it
413
+ internally and in the .pdf tutorial
414
+ what does that mean? upcase as method? hmmm.
415
+ --------------------------------------------------------------------------------
416
+ (62) → http://www.biomart.org/other/user-docs.pdf
417
+ ^^^ work through this
418
+ ^^^ integrate the old .cgi part and improve as you go
419
+ --------------------------------------------------------------------------------
420
+ (63) Access geninfo numbers easily.
516
421
  Die suchen und runterladen.
517
- -------------------------------------------------------------------------------
518
- - Add all of bioruby into bioroebe:
519
-
520
- continous project
521
- https://github.com/biopython/biopython
522
- https://github.com/bioruby/bioruby/tree/master/lib/bio
523
- -------------------------------------------------------------------------------
524
- (3) https://github.com/bioruby/bioruby/issues/134
525
- ^^^ check this, for restriction enzymes
526
- http://rebase.neb.com/rebase/enz/MboII.html
527
-
528
- Bio::RestrictionEnzyme.cut(seq, 'MboII').primary rescue [seq]
529
- => ["agaagattaggatt", "gatgat"]
530
- > seq = seq.reverse_complement
531
- > Bio::RestrictionEnzyme.cut(seq, 'MboII').primary rescue [seq]
532
- => ["atcatcaatcctaatcttct"]
533
- -------------------------------------------------------------------------------
534
- (4) → Document how an ORF is defined for the bioroebe project.
535
- ..........................................................................
536
- (5) Continue with biojava in bioroebe.
537
-
538
- → We need to make some table that tells us what is implemented
422
+ --------------------------------------------------------------------------------
423
+ (64) Add all of bioruby into bioroebe:
424
+ continous project
425
+ https://github.com/biopython/biopython
426
+ https://github.com/bioruby/bioruby/tree/master/lib/bio
427
+ --------------------------------------------------------------------------------
428
+ (65) → https://github.com/bioruby/bioruby/issues/134
429
+ ^^^ check this, for restriction enzymes
430
+ http://rebase.neb.com/rebase/enz/MboII.html
431
+ Bio::RestrictionEnzyme.cut(seq, 'MboII').primary rescue [seq]
432
+ => ["agaagattaggatt", "gatgat"]
433
+ > seq = seq.reverse_complement
434
+ > Bio::RestrictionEnzyme.cut(seq, 'MboII').primary rescue [seq]
435
+ => ["atcatcaatcctaatcttct"]
436
+ --------------------------------------------------------------------------------
437
+ (66) → Document how an ORF is defined for the bioroebe project.
438
+ --------------------------------------------------------------------------------
439
+ (67) → Continue with biojava in bioroebe.
440
+ → We need to make some table that tells us what is implemented
539
441
  in java.
540
- → Make it possible to randomly generate aminoacids, and then,
442
+ → Make it possible to randomly generate aminoacids, and then,
541
443
  based on that, design degenarate DNA that matches to it) ←
542
444
  this also must work standalone, and be documented.
543
-
544
- document on bioroebe.cgi as well
545
-
546
- We can generate degenerate primers now:
547
-
548
- dprimer M-T-T-Y-Y-T-A-A-A-STOP
549
-
550
- ..........................................................................
551
- (1) → The codon tables:
552
- → In January we added a codon-table GUI to ruby-gtk3.
553
-
445
+ document on bioroebe.cgi as well
446
+ We can generate degenerate primers now:
447
+ dprimer M-T-T-Y-Y-T-A-A-A-STOP
448
+ --------------------------------------------------------------------------------
449
+ (68) → The codon tables:
450
+ → In January we added a codon-table GUI to ruby-gtk3.
554
451
  also enable an inverse table.
555
- Ala/A GCT, GCC, GCA, GCG GCN Leu/L TTA, TTG, CTT, CTC, CTA, CTG YTR, CTN
556
- Arg/R CGT, CGC, CGA, CGG, AGA, AGG CGN, MGR Lys/K AAA, AAG AAR
557
- Asn/N AAT, AAC AAY Met/M ATG
558
- Asp/D GAT, GAC GAY Phe/F TTT, TTC TTY
559
- Cys/C TGT, TGC TGY Pro/P CCT, CCC, CCA, CCG CCN
560
- Gln/Q CAA, CAG CAR Ser/S TCT, TCC, TCA, TCG, AGT, AGC TCN, AGY
561
- Glu/E GAA, GAG GAR Thr/T ACT, ACC, ACA, ACG ACN
562
- Gly/G GGT, GGC, GGA, GGG GGN Trp/W TGG
563
- His/H CAT, CAC CAY Tyr/Y TAT, TAC TAY
564
- Ile/I ATT, ATC, ATA ATH Val/V GTT, GTC, GTA, GTG GTN
565
- START ATG STOP TAA, TGA, TAG TAR, TRA
566
-
567
- I think this is already stored in:
568
- inverse_rna_codon_table.yml
569
-
570
- table = Bio::CodonTable[1]
571
- ^^^^ this is quite a useful feature of bioruby. We need to
572
- add this as well. Then document it.
573
-
574
- ^^^ document this better too
452
+ Ala/A GCT, GCC, GCA, GCG GCN Leu/L TTA, TTG, CTT, CTC, CTA, CTG YTR, CTN
453
+ Arg/R CGT, CGC, CGA, CGG, AGA, AGG CGN, MGR Lys/K AAA, AAG AAR
454
+ Asn/N AAT, AAC AAY Met/M ATG
455
+ Asp/D GAT, GAC GAY Phe/F TTT, TTC TTY
456
+ Cys/C TGT, TGC TGY Pro/P CCT, CCC, CCA, CCG CCN
457
+ Gln/Q CAA, CAG CAR Ser/S TCT, TCC, TCA, TCG, AGT, AGC TCN, AGY
458
+ Glu/E GAA, GAG GAR Thr/T ACT, ACC, ACA, ACG ACN
459
+ Gly/G GGT, GGC, GGA, GGG GGN Trp/W TGG
460
+ His/H CAT, CAC CAY Tyr/Y TAT, TAC TAY
461
+ Ile/I ATT, ATC, ATA ATH Val/V GTT, GTC, GTA, GTG GTN
462
+ START ATG STOP TAA, TGA, TAG TAR, TRA
463
+ I think this is already stored in:
464
+ inverse_rna_codon_table.yml
465
+ table = Bio::CodonTable[1]
466
+ ^^^^ this is quite a useful feature of bioruby. We need to
467
+ add this as well. Then document it.
468
+ ^^^ document this better too
575
469
  that we can now display all the different codon tables.
576
-
577
- This now sorta works semi-ok.
578
-
579
- -------------------------------------------------------------------------------
580
- (1) In the bioroebe-shell, enable input such as:
581
-
582
- NC_000011.10
583
-
584
- This shall quickly download this sequence into the
585
- local file, and also rename it properly.
586
- -------------------------------------------------------------------------------
587
- clone all of bioruby
588
- -------------------------------------------------------------------------------
589
- (1) → bioinf bücher udrhclesen und zeug inkludiere !!!
590
- ^^^^^ mehr bilderchen hinzufügen ... auchv on den GUIs eventuell.
591
- Und auch biopython durcharbeiten und alles wichtige nach
592
- bioroebe übertragen.
593
- -------------------------------------------------------------------------------
594
- - Add: DetectMotif
595
-
596
- This class shall be used for detecting subsequences.
597
- -------------------------------------------------------------------------------
598
- - Neue funktionälit rein
599
- -------------------------------------------------------------------------------
600
- - mehr doku!
601
- -------------------------------------------------------------------------------
602
- - continue on bioroebe, and when it is done, write to the guy.
603
- -------------------------------------------------------------------------------
604
- - Rewrite bioroebe completely - add some tests, too or so, to
605
- test this. ^^^
606
- That way we learn how to write tests.
607
- AND ... we will actually start with the taxonomy project
608
- so that it finally works again.
609
- continue work on bioroebe
610
- MAKE BIOROEBE EPIC because this is what I will make money with.
611
- CONTINUE THE BIOROEBE PORT !
612
- ^^^^
613
- require 'bioroebe/constants/remote_urls.rb
614
- ^^^
615
- ncbi taxonomy databse: move this into this file.
616
- # require 'bioroebe/constants/aminoacid_families.rb'
617
- ^^^ also ... fix this here.
618
- also continue bioroebe port...
619
- hmm. and perhaps add something else, like the option to have
620
- multiple genes and multiple proteis
621
- and define workspaces.
622
- but start with taxonomy first.
623
- ^^^
624
- Also during that rewrite, make sure that the quality of the
625
- documentation improves. That way the whole project serves
626
- as advertisement too.
627
- wait with rewrite though......
628
-
629
- clone bioroebe in java, STEP BY STEP, as means of learning java too,
630
- for preparation at the TU Wien.
631
- ^^^ but first read up some java tutorial because I dont know related-problems about
632
- java
633
-
634
- CPK: international colour scheme
635
- add document in bioroebe + registration
636
- https://proteopedia.org/wiki/index.php/CPK
637
- ^^^
638
- bioroebe erweitern... auch rosalind
639
- ^^^
640
- improve bioroebe so that it is supper
641
- and add C++
642
- extend bioroebe sinatra interface
643
- also add a footer to show which entries are available or so
644
- in bioroebe, mach das die postgresql datenbank wieder funktioniert ...
645
-
646
-
647
-
648
- ..........................................................................
649
-
650
- ^^^ improve this whole project a lot
651
-
652
- before uploading then send email
653
-
654
-
655
- - 1fat.pdb
656
-
657
- ^^^ download this, also via bioshell
658
- download 1fat
659
- ^^^ notify the user about this
660
- but put it into the dir of bioshell
661
-
662
- → add:
663
-
664
- set_dna :insulin
665
- set_dna insulin
666
-
667
- This shall allow us to use the sequence of human insulin
668
- here. Also document this. Shall just make things more
669
- convenient for us.
670
-
671
- http://www.ncbi.nlm.nih.gov/gene/3630
672
-
673
- insulin = 'ncbi_gene: 3630'
674
- → becomes: http://www.ncbi.nlm.nih.gov/gene/3630
675
-
676
- wtf ... better to learn how NCBI uworks
677
- -------------------------------------------------------------------------------
678
- - Add a seuqence table int obioroebe for GFP, YFP etc
679
- and mae this show in both the interactio bioshell but
680
- also the main README.md
681
- -------------------------------------------------------------------------------
682
- - stop_frame1?
683
- ^^^ add support for this
684
- and stop_frame2?
685
- etcc
686
- to show stop-codons in this colour
687
- THEN UPLOAD!
688
- ^^^ this works now but is not documented
689
-
690
-
691
- -------------------------------------------------------------------------------
692
-
693
- - chop to first ATG
694
-
695
- chop :ATG
696
-
697
- ^^^^ enable this, to chop towards the first ATG
698
- sequence in the string
699
-
700
- -------------------------------------------------------------------------------
701
- → http://www.biophp.org/stats/describe_data/demo.php?show=formula
702
-
703
- ^^^ should also add documentation like this, also via www interface
704
- -------------------------------------------------------------------------------
705
- → add mouse chromsoome URL, also in the bioshell
706
- and the main README, to be of help for the
707
- user. add a mouse subsection.
708
- ..........................................................................
709
- → fix the taxonomy stuff...
710
- ..........................................................................
711
- (1) → add 2nd_orf
712
- → this shall scan for the 2nd orf
713
- → and third ORF as well, then, and document it.
714
- ..........................................................................
715
- (2) → Add a "cutter-range example" in restriction enzymes +
470
+ This now sorta works semi-ok.
471
+ --------------------------------------------------------------------------------
472
+ (69) → In the bioroebe-shell, enable input such as:
473
+ NC_000011.10
474
+ This shall quickly download this sequence into the
475
+ local file, and also rename it properly.
476
+ --------------------------------------------------------------------------------
477
+ (70) → clone all of bioruby
478
+ --------------------------------------------------------------------------------
479
+ (71) bioinf bücher udrhclesen und zeug inkludiere !!!
480
+ ^^^^^ mehr bilderchen hinzufügen ... auchv on den GUIs eventuell.
481
+ Und auch biopython durcharbeiten und alles wichtige nach
482
+ bioroebe übertragen.
483
+ --------------------------------------------------------------------------------
484
+ (72) Add: DetectMotif
485
+ This class shall be used for detecting subsequences.
486
+ --------------------------------------------------------------------------------
487
+ (73) → Neue funktionälit rein
488
+ --------------------------------------------------------------------------------
489
+ (74) → mehr doku!!!
490
+ --------------------------------------------------------------------------------
491
+ (75) → Rewrite bioroebe completely - add some tests, too or so, to
492
+ test this. ^^^
493
+ That way we learn how to write tests.
494
+ AND ... we will actually start with the taxonomy project
495
+ so that it finally works again.
496
+ continue work on bioroebe
497
+ MAKE BIOROEBE EPIC because this is what I will make money with.
498
+ CONTINUE THE BIOROEBE PORT !
499
+ ^^^^
500
+ require 'bioroebe/constants/remote_urls.rb
501
+ ^^^
502
+ ncbi taxonomy databse: move this into this file.
503
+ # require 'bioroebe/constants/aminoacid_families.rb'
504
+ ^^^ also ... fix this here.
505
+ also continue bioroebe port...
506
+ hmm. and perhaps add something else, like the option to have
507
+ multiple genes and multiple proteis
508
+ and define workspaces.
509
+ but start with taxonomy first.
510
+ ^^^
511
+ Also during that rewrite, make sure that the quality of the
512
+ documentation improves. That way the whole project serves
513
+ as advertisement too.
514
+ wait with rewrite though......
515
+ clone bioroebe in java, STEP BY STEP, as means of learning java too,
516
+ for preparation at the TU Wien.
517
+ ^^^ but first read up some java tutorial because I dont know related-problems about
518
+ java
519
+ CPK: international colour scheme
520
+ add document in bioroebe + registration
521
+ https://proteopedia.org/wiki/index.php/CPK
522
+ ^^^
523
+ bioroebe erweitern... auch rosalind
524
+ ^^^
525
+ improve bioroebe so that it is supper
526
+ and add C++
527
+ extend bioroebe sinatra interface
528
+ also add a footer to show which entries are available or so
529
+ in bioroebe, mach das die postgresql datenbank wieder funktioniert ...
530
+ --------------------------------------------------------------------------------
531
+ (76) → ^^^ improve this whole project a lot
532
+ before uploading then send email
533
+ → add:
534
+ set_dna :insulin
535
+ set_dna insulin
536
+ This shall allow us to use the sequence of human insulin
537
+ here. Also document this. Shall just make things more
538
+ convenient for us.
539
+ http://www.ncbi.nlm.nih.gov/gene/3630
540
+ insulin = 'ncbi_gene: 3630'
541
+ → becomes: http://www.ncbi.nlm.nih.gov/gene/3630
542
+ wtf ... better to learn how NCBI uworks
543
+ --------------------------------------------------------------------------------
544
+ (77) Add a seuqence table into bioroebe for GFP, YFP etc
545
+ and mae this show in both the interactio bioshell but
546
+ also the main README.md
547
+ --------------------------------------------------------------------------------
548
+ (78) → http://www.biophp.org/stats/describe_data/demo.php?show=formula
549
+ ^^^ should also add documentation like this, also via www interface
550
+ --------------------------------------------------------------------------------
551
+ (79) Add a "cutter-range example" in restriction enzymes +
716
552
  table + examples + tutorial
717
-
718
- one example each in this overview.
719
-
553
+ one example each in this overview.
720
554
  Also, add in the documentation where this
721
555
  can be found.
722
- ..........................................................................
723
- (3) → Add aaruler, similar to "ruler"; in the bioshell.
724
- But we want to do this on the dna-sequence rather
725
- than the aminoacid sequence.
726
- This works but the display is not ideal.
727
- ..........................................................................
728
- (4) → Add some codon-usage analyzer. What shall it show? It
556
+ --------------------------------------------------------------------------------
557
+ (80) → Add some codon-usage analyzer. What shall it show? It
729
558
  should show how many codons are used, frequencies etc...
730
559
  by an organism, and compare that to other data.
731
- ..........................................................................
732
- (5) → Implement a GPCR interface.
733
-
560
+ --------------------------------------------------------------------------------
561
+ (81) → Implement a GPCR interface.
734
562
  This is for "G-protein coupled receptors."
735
563
  Denote which variants exist and so forth. Document it as well.
736
- ..........................................................................
737
- (6) → alu?
738
-
564
+ --------------------------------------------------------------------------------
565
+ (82) → alu?
739
566
  Will read from the file `/Programs/Ruby/2.3.0/lib/ruby/site_ruby/2.3.0/bioroebe/yaml/alu_elements.yml`.
740
567
  Bioroebe::ParseFasta: This sequence is assumed to be a protein.
741
568
  This sequence has 1317 aminoacids.
742
-
743
569
  We have identified a total of 1 entries in this fasta dataset.
744
570
  The ALU sequence in humans may be:
745
-
746
- GC ...
747
-
748
- ^^^ das stimmt aber net ... hmmm.
749
-
750
- (3) → The .pdb file that used to be distributed via bioroebe was way
751
- too large. Perhaps add a way to download it instead if needed
752
- e. g.:
753
-
754
- common_downloads:
755
-
756
- ^^^ add this and document it or something like that
757
- And perhaps add a small protein as an example how to
758
- work with .pdb files instead.
759
- -------------------------------------------------------------------------------
760
- (4) → Extend bioroebe to allow download
761
-
762
- PDB files
763
-
764
- id 3030
765
-
766
- and then display nice thingies to the user.
767
-
768
- http://www.pdb.org/pdb/download/downloadFile.do?fileFormat=pdb&compression=NO&structureId=3030
769
- http://www.pdb.org/pdb/download/downloadFile.do?fileFormat=pdb&compression=NO&structureId=2VEZ
770
-
771
- in 3EML 2VTP 2VEZ
772
- do
773
- ..........................................................................
774
- (1) → Fully integrate electron microscopy then remove the old entry.
571
+ GC ...
572
+ ^^^ das stimmt aber net ... hmmm.
573
+ (3) → The .pdb file that used to be distributed via bioroebe was way
574
+ too large. Perhaps add a way to download it instead if needed
575
+ e. g.:
576
+ common_downloads:
577
+ ^^^ add this and document it or something like that
578
+ And perhaps add a small protein as an example how to
579
+ work with .pdb files instead.
580
+ --------------------------------------------------------------------------------
581
+ (83) → Extend bioroebe to allow download
582
+ PDB files
583
+ id 3030
584
+ and then display nice thingies to the user.
585
+ http://www.pdb.org/pdb/download/downloadFile.do?fileFormat=pdb&compression=NO&structureId=3030
586
+ http://www.pdb.org/pdb/download/downloadFile.do?fileFormat=pdb&compression=NO&structureId=2VEZ
587
+ in 3EML 2VTP 2VEZ
588
+ do
589
+ --------------------------------------------------------------------------------
590
+ (84) → Fully integrate electron microscopy then remove the old entry.
775
591
  Test it though.
776
592
  Hmm... but ... we will first polish the main bioroebe
777
593
  gem AND the taxonomy gem and THEN AFTERWARDS
778
594
  integate elctron microsopcy.
779
- ..........................................................................
780
- (1) → ORF Finder:
781
-
595
+ --------------------------------------------------------------------------------
596
+ (85) → ORF Finder:
782
597
  We must add an ORF finder for the bioroebe project,
783
598
  similar to the NCBI ORF Finder.
784
-
785
599
  This works partially... start_stop works but we do not
786
600
  yet find all subsequences.
787
-
788
- ..........................................................................
789
- (1) → must change determine whether we have protein or nucleotide or
601
+ --------------------------------------------------------------------------------
602
+ (86) → must change determine whether we have protein or nucleotide or
790
603
  so via a topelvel method!
791
- ..........................................................................
792
- (1) → there is a talens module.
793
- we have to improve on it for a while
794
- better docu
795
- more testing
796
- then we can get rid of this entry here
797
- ..........................................................................
798
- (1) → 33.44
604
+ --------------------------------------------------------------------------------
605
+ (87) → there is a talens module.
606
+ we have to improve on it for a while
607
+ better docu
608
+ more testing
609
+ then we can get rid of this entry here
610
+ --------------------------------------------------------------------------------
611
+ (88) → 33.44
799
612
  Next showing the nucleotides 33 to 44 (including 33 and 44).
800
- The length of the fragment will be 12 nucleotides.
801
- 5' - 2;70;130;180 - 3'
802
- ^^^ there is some problem; we somehow embed the colour codes,
803
- which should not happen.
804
- ..........................................................................
805
- (1) → set_aa DTLCIGYHAN NSTDTVDTVL EKNVTVTHSV NLLEDKHNGK LCKLRGVAPL HLGKCNIAGW ILGNPECESL STASSWSYIV ETSNSDNGTC YPGDFINYEE LREQLSSVSS FERFEIFPKT SSWPNHDNKG VTAACPHAGA KSFYKNLIWL VKKGNSYPKL NQSYINDKGK EVLVLWGIHH PSTTADQQSL YQNADAYVFV GTSRYSKKFK PEIATRPKVR DQEGRMNYYW TLVEPGDKIT FEATGNLVVP RYAFMERNAG SGIIISDTPV HDCNTTCQTP EGAINTSLPF QNIHPITIGK CPKYVKSTKL RLATGLRNVP SIQSRGLFGA IAGFIEGGWT GMVDGWYGYH HQNEQGSGYA ADLKSTQNAI DKITNKVNSV IKMNTQFTAV GKEFNHLEKR IENLNKKVDD GFLDIWTYNA ELLVLLENER TLDYHDSNVK NLYEKVRNQL KNNAKEIGNG CFEFYHKCDN TCMESVKNGT YDYPKYSEEA KLNREKIDGV KLESTRIYHH HHHH
806
-
807
- ^^^ enable copy/pasting,
808
- then reverse_sequence
809
- dna_sequence?
810
-
811
- @type=:dna>
812
-
813
- BIO SHELL> aaseq?
613
+ The length of the fragment will be 12 nucleotides.
614
+ 5' - 2;70;130;180 - 3'
615
+ ^^^ there is some problem; we somehow embed the colour codes,
616
+ which should not happen.
617
+ --------------------------------------------------------------------------------
618
+ (89) → set_aa DTLCIGYHAN NSTDTVDTVL EKNVTVTHSV NLLEDKHNGK LCKLRGVAPL HLGKCNIAGW ILGNPECESL STASSWSYIV ETSNSDNGTC YPGDFINYEE LREQLSSVSS FERFEIFPKT SSWPNHDNKG VTAACPHAGA KSFYKNLIWL VKKGNSYPKL NQSYINDKGK EVLVLWGIHH PSTTADQQSL YQNADAYVFV GTSRYSKKFK PEIATRPKVR DQEGRMNYYW TLVEPGDKIT FEATGNLVVP RYAFMERNAG SGIIISDTPV HDCNTTCQTP EGAINTSLPF QNIHPITIGK CPKYVKSTKL RLATGLRNVP SIQSRGLFGA IAGFIEGGWT GMVDGWYGYH HQNEQGSGYA ADLKSTQNAI DKITNKVNSV IKMNTQFTAV GKEFNHLEKR IENLNKKVDD GFLDIWTYNA ELLVLLENER TLDYHDSNVK NLYEKVRNQL KNNAKEIGNG CFEFYHKCDN TCMESVKNGT YDYPKYSEEA KLNREKIDGV KLESTRIYHH HHHH
619
+ ^^^ enable copy/pasting,
620
+ then reverse_sequence
621
+ dna_sequence?
622
+ @type=:dna>
623
+ BIO SHELL> aaseq?
814
624
  DTLCIGYHANNSTDTVDTVLEKNVTVTHSVNLLEDKHNGKLCKLRGVAPLHLGKCNIAGWILGNPECESLSTASSWSYIVETSNSDNGTCYPGDFINYEELREQLSSVSSFERFEIFPKTSSWPNHDNKGVTAACPHAGAKSFYKNLIWLVKKGNSYPKLNQSYINDKGKEVLVLWGIHHPSTTADQQSLYQNADAYVFVGTSRYSKKFKPEIATRPKVRDQEGRMNYYWTLVEPGDKITFEATGNLVVPRYAFMERNAGSGIIISDTPVHDCNTTCQTPEGAINTSLPFQNIHPITIGKCPKYVKSTKLRLATGLRNVPSIQSRGLFGAIAGFIEGGWTGMVDGWYGYHHQNEQGSGYAADLKSTQNAIDKITNKVNSVIKMNTQFTAVGKEFNHLEKRIENLNKKVDDGFLDIWTYNAELLVLLENERTLDYHDSNVKNLYEKVRNQLKNNAKEIGNGCFEFYHKCDNTCMESVKNGTYDYPKYSEEAKLNREKIDGVKLESTRIYHHHHHH
815
- BIO SHELL> aasize?
816
- This sequence has 50 aminoacids.
817
- ^^^ das stimmt net.
818
-
819
- ..........................................................................
820
- (1) → add this functionality:
821
-
625
+ BIO SHELL> aasize?
626
+ This sequence has 50 aminoacids.
627
+ ^^^ das stimmt net.
628
+ --------------------------------------------------------------------------------
629
+ (90) → add this functionality:
822
630
  meting temper
823
631
  melting temper
824
632
  melting_temperature?
825
-
826
- ^^^ for short primers
827
-
828
- 4°C for each G/C base pair
829
- → 2°C for each A/T base pair
830
-
831
- Also add an explanation as can be seen here:
832
-
833
- http://comments.gmane.org/gmane.comp.lang.ruby.bio/1182
834
-
835
- "I discovered the above discussion by accident, about two
836
- years lateron. :)
837
-
838
- I am not using email-discussions or usegroups/newsgroups
839
- in general, largely because I have never been able to
840
- keep up to date with them and deal with emails properly;
841
- I am more a casual emails user myself, growing up in a
842
- www-world where phpBB really made it convenient to
843
- communicate with other people. So probably a bit after the
844
- emails-people use emails.
845
-
846
- At any rate, when I noticed it, I decided on my todo-list
847
- that I will improve the melting-temperature calculation
848
- of BioRoebe.
849
-
850
- Note that NCBI Blast and several other sites also already
851
- have very good algorithms in this regards, so the prime
852
- use case for BioRoebe is to explain a bit the algorithms
853
- and also provide a commandline-way to calculate them,
854
- using ruby. The latter may be useful and rather easy for
855
- scripted use.
856
- ..........................................................................
857
- (1) → show insulin
633
+ ^^^ for short primers
634
+ 4°C for each G/C base pair
635
+ → 2°C for each A/T base pair
636
+ Also add an explanation as can be seen here:
637
+ http://comments.gmane.org/gmane.comp.lang.ruby.bio/1182
638
+ "I discovered the above discussion by accident, about two
639
+ years lateron. :)
640
+ I am not using email-discussions or usegroups/newsgroups
641
+ in general, largely because I have never been able to
642
+ keep up to date with them and deal with emails properly;
643
+ I am more a casual emails user myself, growing up in a
644
+ www-world where phpBB really made it convenient to
645
+ communicate with other people. So probably a bit after the
646
+ emails-people use emails.
647
+ At any rate, when I noticed it, I decided on my todo-list
648
+ that I will improve the melting-temperature calculation
649
+ of BioRoebe.
650
+ Note that NCBI Blast and several other sites also already
651
+ have very good algorithms in this regards, so the prime
652
+ use case for BioRoebe is to explain a bit the algorithms
653
+ and also provide a commandline-way to calculate them,
654
+ using ruby. The latter may be useful and rather easy for
655
+ scripted use.
656
+ --------------------------------------------------------------------------------
657
+ (91) → show insulin
858
658
  ^^^ to show the insulin structure
859
- how to find it? no idea...
860
- but we should have these structures already made available somewhere.
861
- ..........................................................................
862
- (1) → Todo: find family of enzymes, based on sequence structure
659
+ how to find it? no idea...
660
+ but we should have these structures already made available somewhere.
661
+ --------------------------------------------------------------------------------
662
+ (92) → Todo: find family of enzymes, based on sequence structure
863
663
  alone.
864
- ..........................................................................
865
- (1) → https://pubchem.ncbi.nlm.nih.gov/compound/16131099#section=Top
866
-
867
- ^^^ this website is quite interesting; try to use components
868
- from it.
869
- -------------------------------------------------------------------------------
870
- (1) → Add some option to show the aminoacid sequence, at the least
871
- store it; and optionally show it.
872
-
873
- possibly always report how many aminoacids are
874
- part of that file; and optionally also show
875
- the whole sequence.
876
- -------------------------------------------------------------------------------
877
- (1) → WORK THROUGH the PROTOCOL AT BOKU. THEN WORK THROUGH THE VARIOUST
664
+ --------------------------------------------------------------------------------
665
+ (93) → WORK THROUGH the PROTOCOL AT BOKU. THEN WORK THROUGH THE VARIOUST
878
666
  TIDBIDS AT UNI WIEN STARTING WITH HEIKO.
879
- ^^^ da sind wir nun.
667
+ ^^^ da sind wir nun.
880
668
  wir sind an beginn von 1b ... hmmmm, also zerst mal das an der
881
669
  BOKU durchgehen. Dann das löschen.
882
- -------------------------------------------------------------------------------
883
- (1) → Begin tk-bindings for bioroebe, following the gtk stuff.
884
- -------------------------------------------------------------------------------
885
- (2) → frame_value = position_of_the_stop_codon - position_of_the_start_codon
670
+ --------------------------------------------------------------------------------
671
+ (94) → Begin tk-bindings for bioroebe, following the gtk stuff.
672
+ --------------------------------------------------------------------------------
673
+ (95) → frame_value = position_of_the_stop_codon - position_of_the_start_codon
886
674
  ^^^ continue on this ...
887
- -------------------------------------------------------------------------------
888
- (1) → improve both the gtk-apps parts, and the sinatra web-interface,
675
+ --------------------------------------------------------------------------------
676
+ (96) → improve both the gtk-apps parts, and the sinatra web-interface,
889
677
  and other GUI-like elements. The idea is to make this software
890
678
  more useful for people around the world, which should help
891
679
  increase its adoption rate.
892
- -------------------------------------------------------------------------------
893
- (2) → Look to integrate this:
894
-
895
- http://www.ncbi.nlm.nih.gov/nuccore/NM_007315.3?report=fasta&log$=seqview&format=text
680
+ --------------------------------------------------------------------------------
681
+ (97) → Look to integrate this:
682
+ http://www.ncbi.nlm.nih.gov/nuccore/NM_007315.3?report=fasta&log$=seqview&format=text
896
683
  ^^^
897
- -------------------------------------------------------------------------------
898
- (1) → Clone and document the getorf functionality properly.
899
-
900
- See: http://emboss.sourceforge.net/apps/cvs/emboss/apps/getorf.html
901
- -------------------------------------------------------------------------------
902
- (2) → set_dna_sequence alu
903
-
904
- ^^^ fetch random alu
905
-
906
- ^^^ alu sequence
907
- Ok we started this now adding more details, but we
908
- need to become better at searching for this
909
- sequence.
910
- -------------------------------------------------------------------------------
911
- (3) → We need to make available the ... thingy magick
684
+ --------------------------------------------------------------------------------
685
+ (98) → We need to make available the ... thingy magick
912
686
  emboss functionality. that may seem useful
913
687
  but also feel free to extend these parts for
914
688
  bioroebe as necessary.
915
- -------------------------------------------------------------------------------
916
- (4) → integrate electron_microscopy fully
917
- This will take more time, so first we finish with the
689
+ --------------------------------------------------------------------------------
690
+ (99) → integrate electron_microscopy fully
691
+ This will take more time, so first we finish with the
918
692
  taxonomy module instead.
919
- -------------------------------------------------------------------------------
920
- (5) → Improve support for BLAST up until
921
-
693
+ --------------------------------------------------------------------------------
694
+ (100) → Improve support for BLAST up until
922
695
  middle of 2015 so that I am better prepared
923
696
  for work-related stuff. In order for this
924
- to succed, we first have to understand
697
+ to succed, we first have to understand
925
698
  BLAST very well.
926
-
927
699
  So, work on BLAST tutorial at bioinf page:
928
-
929
- bl bioinf; rf bioinf
930
- -------------------------------------------------------------------------------
931
- (3) → integrate a "codon usage database", whatever this means.
700
+ bl bioinf; rf bioinf
701
+ --------------------------------------------------------------------------------
702
+ (101) → integrate a "codon usage database", whatever this means.
932
703
  It is a cool database anyway. Then document this.
933
704
  First, create a codon-usage analyze on a per-FASTA
934
705
  site basis. Meaning we download a fasta sequence
935
706
  and calculate the codon usage from there.
936
-
937
- ^^^ and add some GUI to this. hmmm
938
- ..........................................................................
939
- (4) → Input sequence:
940
-
707
+ ^^^ and add some GUI to this. hmmm
708
+ --------------------------------------------------------------------------------
709
+ (102) → Input sequence:
941
710
  MFLMVSPTAYHQNKDECFLP
942
711
  TAYHQNKDECMVSPTAYHQN
943
712
  KDECFLPTAYHQMVSPTAYH
944
713
  QNKDECFLPTAYHQ
945
714
  Reverse Translated sequence
946
-
947
715
  ATG TTY YTNATGGTNWSNCCNACNGCNTAYCAYCARAAYAARGAYGARTGYTTYYTNCCN
948
716
  ACNGCNTAYCAYCARAAYAARGAYGARTGYATGGTNWSNCCNACNGCNTAYCAYCARAAY
949
717
  AARGAYGARTGYTTYYTNCCNACNGCNTAYCAYCARATGGTNWSNCCNACNGCNTAYCAY
950
718
  CARAAYAARGAYGARTGYTTYYTNCCNACNGCNTAYCAYCAR
951
-
952
719
  ^^^ we should also show this on the commandline AND the
953
- www ... hmmm.
954
- ..........................................................................
955
- (5) → enable a graphical layer so that we can find out which
720
+ www ... hmmm.
721
+ --------------------------------------------------------------------------------
722
+ (103) → enable a graphical layer so that we can find out which
956
723
  transcription factor activates which gene(s). This
957
724
  should show e. g. a transcription factor highlighting
958
725
  a target genetic area.
959
- ..........................................................................
960
- (2) → We should add more screenshots, make them available on imgur
726
+ --------------------------------------------------------------------------------
727
+ (104) → We should add more screenshots, make them available on imgur
961
728
  as well, after storing them locally. Start with the more
962
729
  important functionality.
963
-
964
- ..........................................................................
965
- (2) → clone serial cloner or whatever the name was, that GUI,
966
- so that we can offer the same functionality.
967
- ..........................................................................
968
- (1) →
969
-
970
- # * searching for PubMed IDs given a query string:
971
- # * Bio::PubMed#esearch (recommended)
972
- # * Bio::PubMed#search (only retrieves top 20 hits; will be deprecated)
973
- ^^^ implement this
974
-
975
-
976
- ..........................................................................
977
- (3) → Aufgabe 16 in bioroebe lösen könnnen
978
- ..........................................................................
979
- (4) The taxonomy part should be fully integrated, without it
980
- being a standalone part anymore.
981
- continue on the taxonomy stuff.
982
- ne day this will work again *shake fist*
983
- -------------------------------------------------------------------------------
984
- (5) → re1 = Bio::RestrictionEnzyme::DoubleStranded.new(enzyme1)
985
-
986
- ^^^ add this? hmmmm
987
- ^^^ from here.
988
- -------------------------------------------------------------------------------
989
- (1) → Colourize exon/intron boundaries.
990
- -------------------------------------------------------------------------------
991
- (2) → In bioroebe: enhance phylogeny stuff and perhaps automatically
730
+ --------------------------------------------------------------------------------
731
+ (105) → clone serial cloner or whatever the name was, that GUI,
732
+ so that we can offer the same functionality.
733
+ --------------------------------------------------------------------------------
734
+ (106) →
735
+ # * searching for PubMed IDs given a query string:
736
+ # * Bio::PubMed#esearch (recommended)
737
+ # * Bio::PubMed#search (only retrieves top 20 hits; will be deprecated)
738
+ ^^^ implement this
739
+ --------------------------------------------------------------------------------
740
+ (107) → Aufgabe 16 in bioroebe lösen könnnen
741
+ --------------------------------------------------------------------------------
742
+ (108) → re1 = Bio::RestrictionEnzyme::DoubleStranded.new(enzyme1)
743
+ ^^^ add this? hmmmm
744
+ ^^^ from here.
745
+ --------------------------------------------------------------------------------
746
+ (109) Colourize exon/intron boundaries.
747
+ --------------------------------------------------------------------------------
748
+ (110) In bioroebe: enhance phylogeny stuff and perhaps automatically
992
749
  generate pictures here.
993
- -------------------------------------------------------------------------------
994
- (1) → In sinatra: add a backtranseq entry point, perhaps
750
+ --------------------------------------------------------------------------------
751
+ (111) → In sinatra: add a backtranseq entry point, perhaps
995
752
  alias it as well.
996
-
997
- ^^ sync this to ruby-gtk3? hmm
998
-
999
- bioroebe --protein-to-dna
1000
-
1001
- ^^^ this shall start the GTK3 variant
1002
-
1003
- -------------------------------------------------------------------------------
1004
- (1) require 'rubygems/text'
1005
- include Gem::Text
1006
- levenshtein_distance 'shevy', 'chevy' # => 1
1007
-
1008
- ^^^ add some class that outpus, on the commandline
1009
- the levensthein distance ont he commandline
1010
- and also in the interactive bioshell
1011
-
1012
- https://github.com/rubygems/rubygems/blob/master/lib/rubygems/text.rb
1013
- ^^^ actually move that part into bioroebe itself...
1014
-
1015
- -------------------------------------------------------------------------------
1016
- (1) → add _source to all APIs in sinatra there. Ensure that this works
753
+ ^^ sync this to ruby-gtk3? hmm
754
+ bioroebe --protein-to-dna
755
+ ^^^ this shall start the GTK3 variant
756
+ --------------------------------------------------------------------------------
757
+ (112) → require 'rubygems/text'
758
+ include Gem::Text
759
+ levenshtein_distance 'shevy', 'chevy' # => 1
760
+ ^^^ add some class that outpus, on the commandline
761
+ the levensthein distance ont he commandline
762
+ and also in the interactive bioshell
763
+ https://github.com/rubygems/rubygems/blob/master/lib/rubygems/text.rb
764
+ ^^^ actually move that part into bioroebe itself...
765
+ --------------------------------------------------------------------------------
766
+ (113) add _source to all APIs in sinatra there. Ensure that this works
1017
767
  too. The user should be able to view the source code.
1018
768
  ^^^ it has been added for 2 methods so far in sinatra; we need
1019
- to add it for the remaining ones too. Then we can remove
1020
- this entry point.
1021
- -------------------------------------------------------------------------------
1022
- (2) → Check out expasy
1023
- peptidcutter
1024
- also offer this functionality, through commandline, GUI
769
+ to add it for the remaining ones too. Then we can remove
770
+ this entry point.
771
+ --------------------------------------------------------------------------------
772
+ (114) → Check out expasy
773
+ peptidcutter
774
+ also offer this functionality, through commandline, GUI
1025
775
  and sinatra.
1026
- https://web.expasy.org/peptide_cutter/
1027
- We now have added trypsin but we should add more here; and
776
+ https://web.expasy.org/peptide_cutter/
777
+ We now have added trypsin but we should add more here; and
1028
778
  still have to add support for sinatra here.
1029
- -------------------------------------------------------------------------------
1030
- (3) → melting temperature subsection
1031
-
779
+ --------------------------------------------------------------------------------
780
+ (115) → melting temperature subsection
1032
781
  hmmm .... molecular weight calculation works now ... but
1033
782
  ... is it correct for a ssDNA string? hmm...
1034
- -------------------------------------------------------------------------------
1035
- (3) → Add useful formulas for bioshell.
1036
-
1037
-
1038
- ...........................................................................
1039
- (1) → Degenerate Primers
1040
-
1041
- You can try to determine the degenerate primers via the Shell
1042
- component. Issue the following instructions:
1043
-
1044
- degenerate_primer
1045
-
1046
- ^^^ epxnad that subsection
1047
- more explanations and examples
1048
-
1049
- -------------------------------------------------------------------------------
1050
- (1) → Copy the functionality of plotorf:
1051
-
783
+ --------------------------------------------------------------------------------
784
+ (116) → Degenerate Primers
785
+ You can try to determine the degenerate primers via the Shell
786
+ component. Issue the following instructions:
787
+ degenerate_primer
788
+ ^^^ epxnad that subsection
789
+ more explanations and examples
790
+ --------------------------------------------------------------------------------
791
+ (117) Copy the functionality of plotorf:
1052
792
  See:
1053
-
1054
- http://www.bioinformatics.nl/cgi-bin/emboss/plotorf
1055
-
793
+ http://www.bioinformatics.nl/cgi-bin/emboss/plotorf
1056
794
  Also extend emboss info on the main homepage.
1057
795
  For plotorf we also need to be able to generate images.
1058
796
  We also need a simpler toplevel API here, something like
1059
797
  Bioroebe.return_all_open_reading_frames(of_this_sequence, use_these_as_start_codons = start_codons?, use_these_as_stop_codons = stop_codons?)
1060
798
  ^^^
1061
799
  Bioroebe.return_all_ORFS
1062
-
1063
-
1064
-
1065
- -------------------------------------------------------------------------------
1066
- (2) Start nucleotide position is at: 142
1067
-
1068
- See the following example:
1069
-
1070
- BIO SHELL> highlight AAA
1071
- 5' - GTAACTGTTAAACTGTCAGGCAGGCGCTCAGGTGTACGTTTGATGCTCAGTAGTATTCCATTCTCGCGAGGGTCACGATACCCAAGATCTCCATGGCTTTCTGTTAGACGCAGCCGTGGACGACTAGAGCGTTTTTTTTTGGAAAGTATATGACCAGCACTCTACATCCTAACTAGAAGGTCTTCTAGGCGTACCAATATTAACGAATAGTGAGTGGTTACCCGTACCCGTCATGACGTCTATCATTAATT - 3'
1072
- BIO SHELL>
800
+ --------------------------------------------------------------------------------
801
+ (118) → Start nucleotide position is at: 142
802
+ See the following example:
803
+ BIO SHELL> highlight AAA
804
+ 5' - GTAACTGTTAAACTGTCAGGCAGGCGCTCAGGTGTACGTTTGATGCTCAGTAGTATTCCATTCTCGCGAGGGTCACGATACCCAAGATCTCCATGGCTTTCTGTTAGACGCAGCCGTGGACGACTAGAGCGTTTTTTTTTGGAAAGTATATGACCAGCACTCTACATCCTAACTAGAAGGTCTTCTAGGCGTACCAATATTAACGAATAGTGAGTGGTTACCCGTACCCGTCATGACGTCTATCATTAATT - 3'
805
+ BIO SHELL>
1073
806
  ^^^ this does not work; nothing is highlighted.
1074
-
1075
- -------------------------------------------------------------------------------
1076
- (2) → Add a myristoylierung-signal
1077
-
1078
- Met-Gly-Xaa-Xaa-YXaa-Ser/Thr-Lys-Lys
1079
-
1080
- 1^^ but check first.
1081
-
1082
- -------------------------------------------------------------------------------
1083
- (3) → integrate the bioroebe_tutorial.cgi into the .md file completely.
1084
-
1085
- -------------------------------------------------------------------------------
1086
- (4) → Integrate everything from the biopython tutorial, if it makes
807
+ --------------------------------------------------------------------------------
808
+ (119) → Add a myristoylierung-signal
809
+ Met-Gly-Xaa-Xaa-YXaa-Ser/Thr-Lys-Lys
810
+ 1^^ but check first.
811
+ --------------------------------------------------------------------------------
812
+ (120) → integrate the bioroebe_tutorial.cgi into the .md file completely.
813
+ --------------------------------------------------------------------------------
814
+ (121) → Integrate everything from the biopython tutorial, if it makes
1087
815
  sense.
1088
-
1089
- -------------------------------------------------------------------------------
1090
- (5) → Improve the codon-optimizer in Bioroebe, including the
816
+ --------------------------------------------------------------------------------
817
+ (122) → Improve the codon-optimizer in Bioroebe, including the
1091
818
  documentation. We need to make this really useful.
1092
- -------------------------------------------------------------------------------
1093
- (6) →
1094
- 5'- TACACGGCACAT -3'
1095
- 3'- ATGTGCCGTGTA -5'
1096
-
1097
- Imperfect DNA mirror repeats (IMRs) are less than 100% symmetrical.
1098
-
1099
- ^^^ integrate mirror repeats creation
1100
- and searching for them. Hmmm.
1101
- -------------------------------------------------------------------------------
1102
- (7) continue porting bioroebe/taxonomy
1103
-
1104
- ^^^^^^^^^^
1105
- It has been 5 years ...
1106
-
1107
- ^^^ taxonomy/colours/colours wird integriert
819
+ --------------------------------------------------------------------------------
820
+ (123) →
821
+ 5'- TACACGGCACAT -3'
822
+ 3'- ATGTGCCGTGTA -5'
823
+ Imperfect DNA mirror repeats (IMRs) are less than 100% symmetrical.
824
+ ^^^ integrate mirror repeats creation
825
+ and searching for them. Hmmm.
826
+ --------------------------------------------------------------------------------
827
+ (124) continue porting bioroebe/taxonomy
828
+ ^^^^^^^^^^
829
+ It has been 5 years ...
830
+ ^^^ taxonomy/colours/colours wird integriert
1108
831
  ^^^ das ist der nächste schritt, so das
1109
- wir das nit mehr benötigen.
1110
-
1111
- -------------------------------------------------------------------------------
1112
- (8) → find out which bacteria all contain the needle complex; find out
832
+ wir das nit mehr benötigen.
833
+ --------------------------------------------------------------------------------
834
+ (125) → find out which bacteria all contain the needle complex; find out
1113
835
  the sequence for the needle complex as well and study it;
1114
836
  find the positions of the genes responsible.
1115
-
1116
- -------------------------------------------------------------------------------
1117
- (9) → Add trypsin_digest, also in the shell, but possibly
837
+ --------------------------------------------------------------------------------
838
+ (126) → Add trypsin_digest, also in the shell, but possibly
1118
839
  on toplevel as well (if the input is a protein sequence.
1119
-
1120
840
  Also, more generally in the shell, add this:
1121
-
1122
- digest trypsin
1123
-
841
+ digest trypsin
1124
842
  ^^^ onto the aminoacid
1125
-
1126
- "This is the result when digsting with trypsine."
1127
- And document it; but do not digest if a prolin
1128
- follows !!!
1129
- ^^^ document this too into .md
1130
-
1131
- -------------------------------------------------------------------------------
1132
- (10) → in bioroebe, add a commassie check... do we include
1133
- arginine or not.
1134
-
1135
- ..........................................................................
1136
- (11) add codon usage in bioroebe
1137
- -------------------------------------------------------------------------------
1138
- (12) → Clone the following functionality.
1139
-
1140
- http://www.bioinformatics.nl/cgi-bin/emboss/help/sirna
1141
- -------------------------------------------------------------------------------
1142
- (13) Improve the "find and scan" subsection. We must be able to find
1143
- subsequences; check for "matches" as well, including the bioshell.
1144
- -------------------------------------------------------------------------------
1145
- (14) → Clone the CLUSTAL format aligment.
1146
- -------------------------------------------------------------------------------
1147
- (15) We need to be able to load up a whole geneome into bioroebe,
1148
- and then be able to manipulate it.
1149
-
1150
- ^^^ perhaps test this with some example
1151
- data or so...
1152
- -------------------------------------------------------------------------------
1153
- (16) → Restriction enzymes:
1154
-
1155
- Add a subsection about restritction enzymes including
1156
- examples, and also explain how to use this in bioroebe.
1157
- Minute by minute...
1158
-
1159
- AND add resources to useful sites.
1160
- ^^^ start it... have to expand on it
1161
-
1162
- Also, improve the part about restriction enzymes in
1163
- general, so that we can reproduce and verify the
1164
- information there.
1165
-
1166
- -------------------------------------------------------------------------------
1167
- (18) → clone pepinfo
1168
-
1169
- The program "pepinfo" plots various amino acid properties in
1170
- parallel for an input protein sequence.
1171
-
1172
- The types of plot available are:
1173
-
1174
- i. Hydrophobicity plots using the method of Kyte & Doolittle, the
1175
- optimal matching hydrophobicity scale (OHM) of Sweet & Eisenberg,
1176
- or consensus parameters (Eisenberg et al).
1177
-
1178
- ii. Histogram of the presence of residues with the physico-chemical
1179
- properties: Tiny, Small, Aliphatic, Aromatic, Non-polar, Polar,
1180
- Charged, Positive, Negative.
1181
-
1182
- The data are also written out to an output file.
1183
-
1184
- -------------------------------------------------------------------------------
1185
- (19) gff?
1186
-
1187
- There are 6 .gff3 files in the current directory.
1188
- We will simply pass the first entry there into class Bioroebe::Parser::GFF.
1189
- The accession id is `NZ_CP011602.1`.
1190
- Bioroebe::Parser::GFF: We are instructed to split into standalone files, but we
1191
- Bioroebe::Parser::GFF: can not do so, as there is not more than one accession id
1192
- Bioroebe::Parser::GFF: in this file.
1193
-
1194
- ^^^ we need an analyze-mode as well.
1195
-
1196
- ..........................................................................
1197
- (20) → ^^^^ add the ability to
1198
- show a ruler AND highlighting as well
1199
- ^^^ then document it.
1200
- ..........................................................................
1201
- (21) → https://github.com/bioperl/bioperl-live
1202
- Look what we can take from ^^^.
1203
-
1204
- https://github.com/bioperl/bioperl-live/tree/master/examples
1205
-
1206
- ..........................................................................
1207
- (23) → continue biojava, and bioroebe a bit
1208
-
1209
- Ideally we should have biojava o a working point.
1210
- ..........................................................................
1211
- (24) → Clone all of Emboss. :)
1212
- ..........................................................................
1213
- (25) → clone the functionality found at https://web.expasy.org/protparam/
1214
-
1215
- https://web.expasy.org/cgi-bin/protparam/protparam
1216
- ^^^ this is halfway done...
1217
-
1218
- ^^^^ also add pI calculation:
1219
-
1220
- Theoretical pI: 5.78
1221
-
1222
- -------------------------------------------------------------------------------
1223
- (27) → NP_417539.1
1224
-
843
+ "This is the result when digsting with trypsine."
844
+ And document it; but do not digest if a prolin
845
+ follows !!!
846
+ ^^^ document this too into .md
847
+ --------------------------------------------------------------------------------
848
+ (127) → add codon usage in bioroebe
849
+ --------------------------------------------------------------------------------
850
+ (128) → Clone the following functionality.
851
+ http://www.bioinformatics.nl/cgi-bin/emboss/help/sirna
852
+ --------------------------------------------------------------------------------
853
+ (129) → Improve the "find and scan" subsection. We must be able to find
854
+ subsequences; check for "matches" as well, including the bioshell.
855
+ --------------------------------------------------------------------------------
856
+ (130) → Clone the CLUSTAL format aligment.
857
+ --------------------------------------------------------------------------------
858
+ (131) → We need to be able to load up a whole geneome into bioroebe,
859
+ and then be able to manipulate it.
860
+ ^^^ perhaps test this with some example
861
+ data or so...
862
+ --------------------------------------------------------------------------------
863
+ (132) → Restriction enzymes:
864
+ Add a subsection about restritction enzymes including
865
+ examples, and also explain how to use this in bioroebe.
866
+ Minute by minute...
867
+ AND add resources to useful sites.
868
+ ^^^ start it... have to expand on it
869
+ Also, improve the part about restriction enzymes in
870
+ general, so that we can reproduce and verify the
871
+ information there.
872
+ --------------------------------------------------------------------------------
873
+ (133) clone pepinfo
874
+ The program "pepinfo" plots various amino acid properties in
875
+ parallel for an input protein sequence.
876
+ The types of plot available are:
877
+ i. Hydrophobicity plots using the method of Kyte & Doolittle, the
878
+ optimal matching hydrophobicity scale (OHM) of Sweet & Eisenberg,
879
+ or consensus parameters (Eisenberg et al).
880
+ ii. Histogram of the presence of residues with the physico-chemical
881
+ properties: Tiny, Small, Aliphatic, Aromatic, Non-polar, Polar,
882
+ Charged, Positive, Negative.
883
+ The data are also written out to an output file.
884
+ --------------------------------------------------------------------------------
885
+ (134) → gff?
886
+ There are 6 .gff3 files in the current directory.
887
+ We will simply pass the first entry there into class Bioroebe::Parser::GFF.
888
+ The accession id is `NZ_CP011602.1`.
889
+ Bioroebe::Parser::GFF: We are instructed to split into standalone files, but we
890
+ Bioroebe::Parser::GFF: can not do so, as there is not more than one accession id
891
+ Bioroebe::Parser::GFF: in this file.
892
+ ^^^ we need an analyze-mode as well.
893
+ --------------------------------------------------------------------------------
894
+ (135) → ^^^^ add the ability to
895
+ show a ruler AND highlighting as well
896
+ ^^^ then document it.
897
+ --------------------------------------------------------------------------------
898
+ (136) → https://github.com/bioperl/bioperl-live
899
+ Look what we can take from ^^^.
900
+ https://github.com/bioperl/bioperl-live/tree/master/examples
901
+ --------------------------------------------------------------------------------
902
+ (137) → continue biojava, and bioroebe a bit
903
+ Ideally we should have biojava o a working point.
904
+ --------------------------------------------------------------------------------
905
+ (138) → clone the functionality found at https://web.expasy.org/protparam/
906
+ https://web.expasy.org/cgi-bin/protparam/protparam
907
+ ^^^ this is halfway done...
908
+ ^^^^ also add pI calculation:
909
+ Theoretical pI: 5.78
910
+ --------------------------------------------------------------------------------
911
+ (139) → NP_417539.1
1225
912
  https://www.ncbi.nlm.nih.gov/protein/NP_417539.1
1226
913
  https://www.ncbi.nlm.nih.gov/protein/NP_417539.1?report=fasta
1227
-
1228
- ^^^ if the input is exactly like the above, on the first line,
914
+ ^^^ if the input is exactly like the above, on the first line,
1229
915
  download the sequence.
1230
- -------------------------------------------------------------------------------
1231
- (28) → Integrate these nice GUI parts parts:
1232
-
1233
- https://dev.to/kojix2/introduction-to-gr-rb-data-visualization-with-ruby-2c39
1234
-
1235
- -------------------------------------------------------------------------------
1236
- (29) → http://insilico.ehu.es/
1237
-
1238
- ^^^ check if we have all of this incorporated
1239
- -------------------------------------------------------------------------------
1240
- (30) → http://www.biostars.org/
1241
-
1242
- ^^^ regularly work through this
1243
- and try to help
1244
- and extend bioruby at the same time.
1245
- -------------------------------------------------------------------------------
1246
- (31) → The taxonomy-submodule should work one day, and be properly
1247
- documented as well. Perhaps integrate the parts of Taxonomy
1248
- that can be included into the toplevel domain.
1249
- -------------------------------------------------------------------------------
1250
- (32) → Enable:
1251
-
1252
- Bioroebe.set_genetic_code()
1253
- Bioroebe.set_genetic_code(to: 'Vertebrate Mitochondrial')
1254
-
916
+ --------------------------------------------------------------------------------
917
+ (140) → http://www.biostars.org/
918
+ ^^^ regularly work through this
919
+ and try to help
920
+ and extend bioruby at the same time.
921
+ --------------------------------------------------------------------------------
922
+ (141) → The taxonomy-submodule should work one day, and be properly
923
+ documented as well. Perhaps integrate the parts of Taxonomy
924
+ that can be included into the toplevel domain.
925
+ --------------------------------------------------------------------------------
926
+ (142) → Enable:
927
+ Bioroebe.set_genetic_code()
928
+ Bioroebe.set_genetic_code(to: 'Vertebrate Mitochondrial')
1255
929
  ^^^ enable this
1256
-
1257
- Seq("ATGGCCATTGTAATGGGCCGCTGAAAGGGTGCCCGATAG")
1258
- coding_dna.translate(table=2)
1259
- coding_dna.translate(table="Vertebrate Mitochondrial")
1260
- Seq('MAIVMGRWKGAR*', HasStopCodon(IUPACProtein(), '*'))
1261
- Seq('MAIVMGRWKGAR*')
1262
-
930
+ Seq("ATGGCCATTGTAATGGGCCGCTGAAAGGGTGCCCGATAG")
931
+ coding_dna.translate(table=2)
932
+ coding_dna.translate(table="Vertebrate Mitochondrial")
933
+ Seq('MAIVMGRWKGAR*', HasStopCodon(IUPACProtein(), '*'))
934
+ Seq('MAIVMGRWKGAR*')
1263
935
  ^^^ enable this as well; extent documentation too.
1264
-
1265
- -------------------------------------------------------------------------------
1266
- (34) We have found a restriction enzyme called NheI.
1267
-
1268
- The sequence this 6-cutter relates to is: `5' - GCTAGC - 3'`
1269
-
1270
- This restriction enzyme will produce a blunt overhang.
1271
-
1272
- ^^^ nope das ist falsch
1273
- -------------------------------------------------------------------------------
1274
- (35) → Sau3A?
1275
- ^^^ enable this restriction site
1276
-
1277
- -------------------------------------------------------------------------------
1278
- (37) → Add matplotlib support.
1279
-
1280
- try_to_use_matplotlib
1281
- -------------------------------------------------------------------------------
1282
- (38) → https://www.ncbi.nlm.nih.gov/CBBresearch/Lu/Demo/tmTools/RESTfulAPIs.html
1283
- -------------------------------------------------------------------------------
1284
- (39) → The following input:
1285
-
1286
- downcase; orf?; seq?
1287
-
1288
- leads to strange display. Something is wrong here, must be checked.
1289
- -------------------------------------------------------------------------------
1290
- (40) → Continue with rosalind problems.
1291
-
1292
- These challenges can be found here:
1293
-
1294
- http://rosalind.info/problems/sign/
1295
-
1296
- Also integrate these rosalind-quizzes into bioroebe
1297
- when possible.
1298
- -------------------------------------------------------------------------------
1299
- (41) → https://web.expasy.org/cgi-bin/peptide_mass/peptide-mass.pl
1300
-
1301
- ^^^ make the above usable in sinaitra as well
1302
- -------------------------------------------------------------------------------
1303
- (42) → Integrate a way to search for commonly known promoters:
1304
-
1305
- promoters?
1306
- ^^^ this functionality
1307
- ^^^ this has to be expanded
1308
- and ...
1309
- -------------------------------------------------------------------------------
1310
- (43) → Integrate:
1311
-
1312
- http://biotools.nubic.northwestern.edu/OligoCalc.html
1313
- -------------------------------------------------------------------------------
1314
- (44) Extend the Java part of BioRoebe systematically..
1315
-
1316
- What should come next? Let's make a list.
1317
-
1318
- remove_numbers [DONE]
1319
- -------------------------------------------------------------------------------
1320
- (46) → Study gnuplot; one day we have to draw graphs.
1321
-
1322
- -------------------------------------------------------------------------------
1323
- (47) Add a genome browser, both ascii without GUI and also
1324
- with. In ruby-gtk.
1325
- -------------------------------------------------------------------------------
1326
- (48) Clone the functionality of:
1327
-
1328
- http://www.biophp.org/minitools/restriction_digest/demo.php
1329
- -------------------------------------------------------------------------------
1330
- (50) Add the loxP sequence to readme [DONE] and explain this
1331
- better on the main readme; and perhaps also assign
1332
- the sequence via the bioshell.
1333
- -------------------------------------------------------------------------------
1334
- (51) 33. Cephalodiscidae Mitochondrial UAA-Tyr Code (transl_table=33)
1335
-
1336
- AAs = FFLLSSSSYYY*CCWWLLLLPPPPHHQQRRRRIIIMTTTTNNKKSSSKVVVVAAAADDEEGGGG
1337
- Starts = ---M-------*-------M---------------M---------------M------------
1338
- Base1 = TTTTTTTTTTTTTTTTCCCCCCCCCCCCCCCCAAAAAAAAAAAAAAAAGGGGGGGGGGGGGGGG
1339
- Base2 = TTTTCCCCAAAAGGGGTTTTCCCCAAAAGGGGTTTTCCCCAAAAGGGGTTTTCCCCAAAAGGGG
1340
- Base3 = TCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAG
1341
-
1342
- ^^^ add a parser, and document it, that can take this input
1343
- and output the corresponding code, in a valid .yml file.
1344
- -------------------------------------------------------------------------------
1345
- (52) → Add to bioroebe the ability to add cloning vectors
1346
- and molecular_weight calcuation
1347
- for this
1348
-
1349
- and also to show the sequence of a vector
1350
-
1351
- and how to add them to bioroebe
1352
-
1353
- also download sequence data for a vector - this
1354
- should probably be some interactive table or
1355
- something of the sort here.
1356
-
1357
- ^^^ this works a tiny bit now, must be documented still
1358
- via seq2? and so forth
1359
-
1360
- "pBR322 contains the genes for resistance to ampicillin and
1361
- tetracycline, and can be amplified with chloramphenicol.
1362
- The molecular weight is 2.83 x 106 daltons."
1363
-
1364
- ^^^ we also need a way to find out what resistance genes
1365
- are carried there.
1366
- -------------------------------------------------------------------------------
1367
- (53) → In the lambda genome sequence there are 10 EcoB and
1368
- 5 EcoK sites.
1369
- ^^^ verify this too, as an example as well
1370
- -------------------------------------------------------------------------------
1371
- (54) show restriction sites, composable and compatible with
1372
- serial clone ... hmm
1373
- -------------------------------------------------------------------------------
1374
- (55) → enable:
1375
- BIOROEBE_USE_COLOURS:
1376
- can be 0 or 1
1377
- what is this?
1378
- -------------------------------------------------------------------------------
1379
- (56) → Burrows-Wheeler-Transform (BWT)
1380
-
1381
- ^^^ add some method here
1382
- Bioroebe.burrows_wheeler_transform
1383
- ^^^ if no '$' char is in the input, then append it
1384
-
1385
- then, output the array of a SORTED BWT transform,
1386
- frmo the given input string.
1387
- document this properly as well
1388
-
1389
- also test this against my paper-result
1390
- with input being: "GATAG$".
1391
- -------------------------------------------------------------------------------
1392
- (56) Enable working with several genes... hmm and store that somewhere.
1393
- Something like a per-project workspace thingy.
1394
- -------------------------------------------------------------------------------
1395
- (57) → Add:
1396
-
1397
- http://nar.oxfordjournals.org/content/35/suppl_2/W71.long
1398
-
1399
- -------------------------------------------------------------------------------
1400
- (58) → Now, you may want to translate the nucleotides up to
1401
- the first in frame stop codon, and then stop (as
1402
- happens in nature):
1403
-
1404
- coding_dna.translate()
1405
- Seq('MAIVMGR*KGAR*', HasStopCodon(IUPACProtein(), '*'))
1406
- >>> coding_dna.translate(to_stop=True)
1407
- Seq('MAIVMGR', IUPACProtein())
1408
- ^^^^ support this hmmm.
1409
-
1410
- Then continue from here:
1411
-
1412
- https://people.duke.edu/~ccc14/pcfb/biopython/BiopythonSequences.html
1413
- -------------------------------------------------------------------------------
1414
- (59) → Add:
1415
-
1416
- set_dna :Ubiquitin
1417
- set_dna :ubiquitin
1418
-
1419
- ^^^ we want to obtain the ubuiqitin sequence
1420
- -------------------------------------------------------------------------------
1421
- (59) → Telomers
1422
-
1423
- Telomeres are listed from 5' to 3'.
1424
-
1425
- Example for the human telomeres would be:
1426
- 5'-TTAGGG-3
1427
-
1428
- ^^^ stimmt das?
1429
-
1430
- add:
1431
- doc_telomeres
1432
-
1433
- ^^^ add this to say the human telomere sequence
1434
- -------------------------------------------------------------------------------
1435
- (60) → ORF_positions?
936
+ --------------------------------------------------------------------------------
937
+ (143) → We have found a restriction enzyme called NheI.
938
+ The sequence this 6-cutter relates to is: `5' - GCTAGC - 3'`
939
+ This restriction enzyme will produce a blunt overhang.
940
+ ^^^ nope das ist falsch
941
+ --------------------------------------------------------------------------------
942
+ (144) Sau3A?
943
+ ^^^ enable this restriction site
944
+ --------------------------------------------------------------------------------
945
+ (145) → Add matplotlib support.
946
+ try_to_use_matplotlib
947
+ --------------------------------------------------------------------------------
948
+ (146) → https://www.ncbi.nlm.nih.gov/CBBresearch/Lu/Demo/tmTools/RESTfulAPIs.html
949
+ --------------------------------------------------------------------------------
950
+ (147) → The following input:
951
+ downcase; orf?; seq?
952
+ leads to strange display. Something is wrong here, must be checked.
953
+ --------------------------------------------------------------------------------
954
+ (148) → Continue with rosalind problems.
955
+ These challenges can be found here:
956
+ http://rosalind.info/problems/sign/
957
+ Also integrate these rosalind-quizzes into bioroebe
958
+ when possible.
959
+ --------------------------------------------------------------------------------
960
+ (149) https://web.expasy.org/cgi-bin/peptide_mass/peptide-mass.pl
961
+ ^^^ make the above usable in sinaitra as well
962
+ --------------------------------------------------------------------------------
963
+ (150) → Integrate a way to search for commonly known promoters:
964
+ promoters?
965
+ ^^^ this functionality
966
+ ^^^ this has to be expanded
967
+ and ...
968
+ --------------------------------------------------------------------------------
969
+ (151) → Integrate:
970
+ http://biotools.nubic.northwestern.edu/OligoCalc.html
971
+ --------------------------------------------------------------------------------
972
+ (152) → Extend the Java part of BioRoebe systematically..
973
+ What should come next? Let's make a list.
974
+ → remove_numbers [DONE]
975
+ --------------------------------------------------------------------------------
976
+ (153) → Study gnuplot; one day we have to draw graphs.
977
+ --------------------------------------------------------------------------------
978
+ (154) Add a genome browser, both ascii without GUI and also
979
+ with. In ruby-gtk.
980
+ --------------------------------------------------------------------------------
981
+ (155) → Clone the functionality of:
982
+ http://www.biophp.org/minitools/restriction_digest/demo.php
983
+ --------------------------------------------------------------------------------
984
+ (156) → Add the loxP sequence to readme [DONE] and explain this
985
+ better on the main readme; and perhaps also assign
986
+ the sequence via the bioshell.
987
+ --------------------------------------------------------------------------------
988
+ (157) 33. Cephalodiscidae Mitochondrial UAA-Tyr Code (transl_table=33)
989
+ AAs = FFLLSSSSYYY*CCWWLLLLPPPPHHQQRRRRIIIMTTTTNNKKSSSKVVVVAAAADDEEGGGG
990
+ Starts = ---M-------*-------M---------------M---------------M------------
991
+ Base1 = TTTTTTTTTTTTTTTTCCCCCCCCCCCCCCCCAAAAAAAAAAAAAAAAGGGGGGGGGGGGGGGG
992
+ Base2 = TTTTCCCCAAAAGGGGTTTTCCCCAAAAGGGGTTTTCCCCAAAAGGGGTTTTCCCCAAAAGGGG
993
+ Base3 = TCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAG
994
+ ^^^ add a parser, and document it, that can take this input
995
+ and output the corresponding code, in a valid .yml file.
996
+ --------------------------------------------------------------------------------
997
+ (158) → Add to bioroebe the ability to add cloning vectors
998
+ and molecular_weight calcuation
999
+ for this
1000
+ and also to show the sequence of a vector
1001
+ and how to add them to bioroebe
1002
+ also download sequence data for a vector - this
1003
+ should probably be some interactive table or
1004
+ something of the sort here.
1005
+ ^^^ this works a tiny bit now, must be documented still
1006
+ via seq2? and so forth
1007
+ "pBR322 contains the genes for resistance to ampicillin and
1008
+ tetracycline, and can be amplified with chloramphenicol.
1009
+ The molecular weight is 2.83 x 106 daltons."
1010
+ ^^^ we also need a way to find out what resistance genes
1011
+ are carried there.
1012
+ --------------------------------------------------------------------------------
1013
+ (159) → In the lambda genome sequence there are 10 EcoB and
1014
+ 5 EcoK sites.
1015
+ ^^^ verify this too, as an example as well
1016
+ --------------------------------------------------------------------------------
1017
+ (160) → show restriction sites, composable and compatible with
1018
+ serial clone ... hmm
1019
+ --------------------------------------------------------------------------------
1020
+ (161) → enable:
1021
+ BIOROEBE_USE_COLOURS:
1022
+ can be 0 or 1
1023
+ what is this?
1024
+ --------------------------------------------------------------------------------
1025
+ (162) Burrows-Wheeler-Transform (BWT)
1026
+ ^^^ add some method here
1027
+ Bioroebe.burrows_wheeler_transform
1028
+ ^^^ if no '$' char is in the input, then append it
1029
+ then, output the array of a SORTED BWT transform,
1030
+ frmo the given input string.
1031
+ document this properly as well
1032
+ also test this against my paper-result
1033
+ with input being: "GATAG$".
1034
+ --------------------------------------------------------------------------------
1035
+ (163) → Enable working with several genes... hmm and store that somewhere.
1036
+ Something like a per-project workspace thingy.
1037
+ --------------------------------------------------------------------------------
1038
+ (164) → Add:
1039
+ http://nar.oxfordjournals.org/content/35/suppl_2/W71.long
1040
+ --------------------------------------------------------------------------------
1041
+ (165) Now, you may want to translate the nucleotides up to
1042
+ the first in frame stop codon, and then stop (as
1043
+ happens in nature):
1044
+ coding_dna.translate()
1045
+ Seq('MAIVMGR*KGAR*', HasStopCodon(IUPACProtein(), '*'))
1046
+ >>> coding_dna.translate(to_stop=True)
1047
+ Seq('MAIVMGR', IUPACProtein())
1048
+ ^^^^ support this hmmm.
1049
+ Then continue from here:
1050
+ https://people.duke.edu/~ccc14/pcfb/biopython/BiopythonSequences.html
1051
+ --------------------------------------------------------------------------------
1052
+ (166) → Add:
1053
+ set_dna :Ubiquitin
1054
+ set_dna :ubiquitin
1055
+ ^^^ we want to obtain the ubuiqitin sequence
1056
+ --------------------------------------------------------------------------------
1057
+ (167) Telomers
1058
+ Telomeres are listed from 5' to 3'.
1059
+ Example for the human telomeres would be:
1060
+ 5'-TTAGGG-3
1061
+ ^^^ stimmt das?
1062
+ add:
1063
+ doc_telomeres
1064
+ ^^^ add this to say the human telomere sequence
1065
+ --------------------------------------------------------------------------------
1066
+ (168) → ORF_positions?
1436
1067
  ^^^ change this a bit, to actually show the positions
1437
- of the various ORFs with the start-position.
1438
- -------------------------------------------------------------------------------
1439
- (62) → add:
1440
-
1441
- setgene2
1442
- add_dna2
1443
- dna2
1444
- dna? <--- this one is not a setter but a query.
1445
- -------------------------------------------------------------------------------
1446
- (63) improve the TM calculation. must be better, must have more
1447
- documentation, and a small tutorial.
1448
- -------------------------------------------------------------------------------
1449
- (64) → Compare bioroebe to:
1450
-
1451
- https://www.ncbi.nlm.nih.gov/orffinder
1452
-
1453
- whether both return the same
1454
- also possibly add a web-gui
1455
- -------------------------------------------------------------------------------
1456
- (65) Find out ratios from:
1457
-
1458
- Doolittle RF. 1989. Redundancies in protein sequences. I
1459
-
1460
- http://onlinelibrary.wiley.com/doi/10.1110/ps.9.6.1203/pdf
1461
- ^^^ that table perhaps
1462
-
1463
- (1) require 'bioroebe'
1464
-
1465
- NameError (uninitialized constant Bioroebe::Blosum)
1466
- irb(main):003:0> module LibSSW
1467
- irb(main):004:1> BLOSUM50 = [
1468
- irb(main):005:2*
1469
-
1470
-
1471
- ^^^ also enable Bioroebe::Blosum.matrix?,
1472
- Bioroebe::Blosum.matrix?
1473
- #^^^ add this
1068
+ of the various ORFs with the start-position.
1069
+ --------------------------------------------------------------------------------
1070
+ (169) → add:
1071
+ setgene2
1072
+ add_dna2
1073
+ dna2
1074
+ dna? <--- this one is not a setter but a query.
1075
+ --------------------------------------------------------------------------------
1076
+ (170) → improve the TM calculation. must be better, must have more
1077
+ documentation, and a small tutorial.
1078
+ --------------------------------------------------------------------------------
1079
+ (171) → Compare bioroebe to:
1080
+ https://www.ncbi.nlm.nih.gov/orffinder
1081
+ whether both return the same
1082
+ also possibly add a web-gui
1083
+ --------------------------------------------------------------------------------
1084
+ (172) Find out ratios from:
1085
+ Doolittle RF. 1989. Redundancies in protein sequences. I
1086
+ http://onlinelibrary.wiley.com/doi/10.1110/ps.9.6.1203/pdf
1087
+ ^^^ that table perhaps
1088
+ (1) require 'bioroebe'
1089
+ NameError (uninitialized constant Bioroebe::Blosum)
1090
+ irb(main):003:0> module LibSSW
1091
+ irb(main):004:1> BLOSUM50 = [
1092
+ irb(main):005:2*
1093
+ ^^^ also enable Bioroebe::Blosum.matrix?,
1094
+ Bioroebe::Blosum.matrix?
1095
+ #^^^ add this
1474
1096
  Bioroebe::Blosum[50]
1475
1097
  ^^^ add this, and show an error if the file does not exist.
1476
- .show_matrix
1477
- ^^^ and so forth, also:
1478
- Bioroebe::Blosum[50] as an API.
1479
- and document it in general.
1480
-
1481
- ..........................................................................
1482
- (65) http://www.biomart.org/other/user-docs.pdf
1483
- ^^^ work through this
1484
- -------------------------------------------------------------------------------
1485
- (66) → add:
1486
-
1487
- class Cell
1488
- ^^^ simulate a cell
1098
+ .show_matrix
1099
+ ^^^ and so forth, also:
1100
+ Bioroebe::Blosum[50] as an API.
1101
+ and document it in general.
1102
+ --------------------------------------------------------------------------------
1103
+ (173) → http://www.biomart.org/other/user-docs.pdf
1104
+ ^^^ work through this
1105
+ --------------------------------------------------------------------------------
1106
+ (174) → add:
1107
+ class Cell
1108
+ ^^^ simulate a cell
1489
1109
  Hmmm. Needs specific components ... and needs a better plan.
1490
- -------------------------------------------------------------------------------
1491
- (68) → class Protein:
1492
-
1493
- add glycosyslation patteren
1494
- .glycosylated? yes no
1495
- + glycoslated?
1496
- need to somehow add the modiication type
1497
-
1498
- https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5358406/
1499
- -------------------------------------------------------------------------------
1500
- (69) → In the BioShell we must be able to do probes - completementary
1501
- to amino acids.
1502
- -------------------------------------------------------------------------------
1503
- (70) Add www-related functionality to bioroebe eventually make use
1504
- of rails, but start with sinatra possibly. In the long run,
1505
- make it flexible to work with as many different frameworks
1506
- as possible, though.
1507
- -------------------------------------------------------------------------------
1508
- (71) → Spaltstellen anzeigen zum beispiel lambda-DNA verdau
1509
- BgI II.
1510
- -------------------------------------------------------------------------------
1511
- (72) → dnaanalyze
1512
-
1513
- In the DNA string `TCCGTCGCAACACATCGCCTCAACAAACCGACCGGGATATGCAATACCGGAATCCGATCCTTTAGAAGCTGCATTCCAAACGCTTGCAATAACACCCACTCGACTATTCAGCATTGGCAAAGGGTACGAATTCGACGAAGGGAGGGTGCTATATTTTCCAAGTTGCTCGCCGATTGATACGGAGCCTGTGGAAAGATTTCGCGGCTCTAGTCTTTAGCTTTGATGTCACCCCTGAGTAGTAACCCGGCGTGGTAGCTTTCATTAGACTTCTCGGAGAGAGTATTAAGCAAAGGTGGAGGTCCCAGGGGTCCAGTGAGCTGTATCGCACTAAAAGCATGCCTACGGGCAATGCTATTTTGCTCACAGGAACTTTGGGGGAGCCACAAACTCTCGAAGCCGGATTGTTGTGGCGGCTAACTTTCCAAAGGCGACCATTCATGGTCTGAATGGGCCCTCACCAGAAGAACGTTTTCGACGGGCATTCTTCCCCGGGGTTTCGAAGGCAAGGGTCAGCACGGCGCGGAAAAGTACGCGACGCATACCGGACTAGTCATGCAACTCCCTCGGAACTGGCGATTCCCACCCAAGAGACGCACGCTGATCATTGCCCATGCCGACTGGAGATGCTGAATTTGGTATGCGGGTCTGTTGCCAGCGCTGACATTATCGGACATTGTGGGGAGAACCGTGTGATTGATTGAGCTGGCGCATTTGTCCGCATGCTCTCCTCATGTGGACACCTTCGCAGGTTCTTTCCGCGGCCACAGTGTCGGGATCTACCCCTGGTGCGTCGCCGCGAGTACAGGTGGGGTTTCGCGCATGAGAACCAATGTTGCACGCCTCAAAACATGGCTGTAACATATTAGCGCCAATAAAAATTTTTGGCAACAAAGAAACAAGGCCAACCGAAGTGCTAAGCCGCGATCATGAAGGGGCGATGCCAGAATGGGAGTCTGCCTTTCCTGTGTGGACGTGAGATTGTACCTAGACAGAGAACGCC` we found these Nucleotides:
1514
- ================================================================================
1515
- Adenines: 244 | 24.40 %
1516
- Guanines: 273 | 27.30 %
1517
- Cytosines: 255 | 25.50 %
1518
- Thymine: 228 | 22.80 %
1519
-
1520
-
1521
- ^^^ created balanced composition
1522
-
1110
+ --------------------------------------------------------------------------------
1111
+ (175) → class Protein:
1112
+ add glycosyslation patteren
1113
+ .glycosylated? yes no
1114
+ + glycoslated?
1115
+ need to somehow add the modiication type
1116
+ https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5358406/
1117
+ --------------------------------------------------------------------------------
1118
+ (176) → In the BioShell we must be able to do probes - completementary
1119
+ to amino acids.
1120
+ --------------------------------------------------------------------------------
1121
+ (177) → Add www-related functionality to bioroebe eventually make use
1122
+ of rails, but start with sinatra possibly. In the long run,
1123
+ make it flexible to work with as many different frameworks
1124
+ as possible, though.
1125
+ --------------------------------------------------------------------------------
1126
+ (178) Spaltstellen anzeigen zum beispiel lambda-DNA verdau
1127
+ BgI II.
1128
+ --------------------------------------------------------------------------------
1129
+ (179) → dnaanalyze
1130
+ In the DNA string `TCCGTCGCAACACATCGCCTCAACAAACCGACCGGGATATGCAATACCGGAATCCGATCCTTTAGAAGCTGCATTCCAAACGCTTGCAATAACACCCACTCGACTATTCAGCATTGGCAAAGGGTACGAATTCGACGAAGGGAGGGTGCTATATTTTCCAAGTTGCTCGCCGATTGATACGGAGCCTGTGGAAAGATTTCGCGGCTCTAGTCTTTAGCTTTGATGTCACCCCTGAGTAGTAACCCGGCGTGGTAGCTTTCATTAGACTTCTCGGAGAGAGTATTAAGCAAAGGTGGAGGTCCCAGGGGTCCAGTGAGCTGTATCGCACTAAAAGCATGCCTACGGGCAATGCTATTTTGCTCACAGGAACTTTGGGGGAGCCACAAACTCTCGAAGCCGGATTGTTGTGGCGGCTAACTTTCCAAAGGCGACCATTCATGGTCTGAATGGGCCCTCACCAGAAGAACGTTTTCGACGGGCATTCTTCCCCGGGGTTTCGAAGGCAAGGGTCAGCACGGCGCGGAAAAGTACGCGACGCATACCGGACTAGTCATGCAACTCCCTCGGAACTGGCGATTCCCACCCAAGAGACGCACGCTGATCATTGCCCATGCCGACTGGAGATGCTGAATTTGGTATGCGGGTCTGTTGCCAGCGCTGACATTATCGGACATTGTGGGGAGAACCGTGTGATTGATTGAGCTGGCGCATTTGTCCGCATGCTCTCCTCATGTGGACACCTTCGCAGGTTCTTTCCGCGGCCACAGTGTCGGGATCTACCCCTGGTGCGTCGCCGCGAGTACAGGTGGGGTTTCGCGCATGAGAACCAATGTTGCACGCCTCAAAACATGGCTGTAACATATTAGCGCCAATAAAAATTTTTGGCAACAAAGAAACAAGGCCAACCGAAGTGCTAAGCCGCGATCATGAAGGGGCGATGCCAGAATGGGAGTCTGCCTTTCCTGTGTGGACGTGAGATTGTACCTAGACAGAGAACGCC` we found these Nucleotides:
1131
+ ================================================================================
1132
+ Adenines: 244 | 24.40 %
1133
+ Guanines: 273 | 27.30 %
1134
+ Cytosines: 255 | 25.50 %
1135
+ Thymine: 228 | 22.80 %
1136
+ ^^^ created balanced composition
1523
1137
  "Enter the percentage"
1524
-
1525
1138
  "interactive_string"
1526
1139
  ^^^ here we ask user input
1527
-
1528
1140
  otherwise, we assume it to be 25%
1529
1141
  Ok 25% works...
1530
- The other part does not yet work. I am so lazy...
1531
- ^^^ add a GUI part too. The GUI has been added;
1532
- we need to make it so that an input sequence
1533
- can be assigned, and dnaanalyse --GUI should
1534
- start it too. ALSO document it once this works.
1535
- -------------------------------------------------------------------------------
1536
- (73) → go through the individual components slowly and improve them,
1537
- step by step, including the documentation. Then eventually
1538
- remove this todo-entry here.
1539
- -------------------------------------------------------------------------------
1540
- (74) → Add a consensus sequence for:
1541
-
1542
- Asn-X-Ser/Thr-Conesnsus
1543
-
1544
- first in a yaml file; also documented, then also add
1142
+ The other part does not yet work. I am so lazy...
1143
+ ^^^ add a GUI part too. The GUI has been added;
1144
+ we need to make it so that an input sequence
1145
+ can be assigned, and dnaanalyse --GUI should
1146
+ start it too. ALSO document it once this works.
1147
+ --------------------------------------------------------------------------------
1148
+ (180) → go through the individual components slowly and improve them,
1149
+ step by step, including the documentation. Then eventually
1150
+ remove this todo-entry here.
1151
+ --------------------------------------------------------------------------------
1152
+ (181) → Add a consensus sequence for:
1153
+ Asn-X-Ser/Thr-Conesnsus
1154
+ first in a yaml file; also documented, then also add
1545
1155
  a way to scan for these something like:
1546
- N-Glycosylation?
1547
- NGlycosylation?
1548
- NGlyc
1549
- /N-?Glyc/i
1550
- ^^^ use that regex
1551
- -------------------------------------------------------------------------------
1552
- (74) → make sure that newly generated files respect the
1553
- default chmod value on the system. from bioroebe.
1554
- right now we default to 755 which I assume is
1555
- hardcoded but perhaps this is wrong.
1556
-
1557
- -------------------------------------------------------------------------------
1558
- (75) require 'bio'
1559
-
1560
- # creating a Bio::Sequence::NA object containing ambiguous alphabets
1561
- ambiguous_seq = Bio::Sequence::NA.new("atgcyrwskmbdhvn")
1562
-
1563
- # show the contents and class of the DNA sequence object
1564
- p ambiguous_seq # => "atgcyrwskmbdhvn"
1565
- p ambiguous_seq.class # => Bio::Sequence::NA
1566
-
1567
- # convert the sequence to a Regexp object
1568
- p ambiguous_seq.to_re # => /atgc[tc][ag][at][gc][tg][ac][tgc][atg][atc][agc][atgc]/
1569
- p ambiguous_seq.to_re.class # => Regexp
1570
-
1571
- ^^^ add .to_re to this. it must generate a regexp object.
1572
-
1573
-
1574
-
1575
- - Also add some restirctione nzymes example
1576
- on the bioroebe readme... and bio.cgi
1577
-
1578
- perhaps add the restriction enyme
1579
- part nto a standalone file
1580
- so taht it can be used by both the .cgi and
1581
- well rdoc...
1582
- -------------------------------------------------------------------------------
1583
- - Add more protein-specific thingies to bioroebe.
1584
- -------------------------------------------------------------------------------
1585
- - Die bioshell vorantreiben und durch std_biology.rb abarbeiten.
1586
- Vielleicht können wir ja etwas davon auslagern in eine Klasse
1587
- oder so.
1588
-
1589
- Das ganze sollte auch mit Webmin (biomin) verknüpft werden, so das
1590
- wir die Bioshell auch elegant über das www verwenden können!
1591
- -------------------------------------------------------------------------------
1592
- - ^^^ when we find restriction enzyme sites in a DNA
1156
+ N-Glycosylation?
1157
+ NGlycosylation?
1158
+ NGlyc
1159
+ /N-?Glyc/i
1160
+ ^^^ use that regex
1161
+ --------------------------------------------------------------------------------
1162
+ (182) → require 'bio'
1163
+ # creating a Bio::Sequence::NA object containing ambiguous alphabets
1164
+ ambiguous_seq = Bio::Sequence::NA.new("atgcyrwskmbdhvn")
1165
+ # show the contents and class of the DNA sequence object
1166
+ p ambiguous_seq # => "atgcyrwskmbdhvn"
1167
+ p ambiguous_seq.class # => Bio::Sequence::NA
1168
+ # convert the sequence to a Regexp object
1169
+ p ambiguous_seq.to_re # => /atgc[tc][ag][at][gc][tg][ac][tgc][atg][atc][agc][atgc]/
1170
+ p ambiguous_seq.to_re.class # => Regexp
1171
+ ^^^ add .to_re to this. it must generate a regexp object.
1172
+ - Also add some restirctione nzymes example
1173
+ on the bioroebe readme... and bio.cgi
1174
+ perhaps add the restriction enyme
1175
+ part nto a standalone file
1176
+ so taht it can be used by both the .cgi and
1177
+ well rdoc...
1178
+ --------------------------------------------------------------------------------
1179
+ (183) Add more protein-specific thingies to bioroebe.
1180
+ --------------------------------------------------------------------------------
1181
+ (184) Die bioshell vorantreiben und durch std_biology.rb abarbeiten.
1182
+ Vielleicht können wir ja etwas davon auslagern in eine Klasse
1183
+ oder so.
1184
+ Das ganze sollte auch mit Webmin (biomin) verknüpft werden, so das
1185
+ wir die Bioshell auch elegant über das www verwenden können!
1186
+ --------------------------------------------------------------------------------
1187
+ (185) → ^^^ when we find restriction enzyme sites in a DNA
1593
1188
  string, colourize them RED.
1594
-
1595
- also set it to
1596
- set_restriction_size()
1597
- -------------------------------------------------------------------------------
1598
- - also ... while learning C++ we extend the project here...
1599
- Useful C++ things will be combined.
1600
- -------------------------------------------------------------------------------
1601
- - As of April 2003, there were 176,890 total taxa represented.
1602
-
1603
- ^^^ we need a way to also output how many entries we
1604
- have there.
1605
- -------------------------------------------------------------------------------
1606
- - Replace bioruby with bioroebe completely!
1607
- In order for this to work, we first need to find out
1608
- what bioruby is able to do. :P
1609
- -------------------------------------------------------------------------------
1610
- - append 33
1611
- # ^^^ in the bioshell
1612
- Only numbers were given: Adding 33 random nucleotides to the main string next.
1613
- Traceback (most recent call last):
1614
- 10: from /usr/bin/bioshell:26:in `<main>'
1615
- 9: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/misc.rb:4121:in `shell'
1616
- 8: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/misc.rb:4121:in `new'
1617
- 7: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/initialize.rb:168:in `initialize'
1618
- 6: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/loop.rb:18:in `enter_main_loop'
1619
- 5: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/loop.rb:18:in `loop'
1620
- 4: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/loop.rb:30:in `block in enter_main_loop'
1621
- 3: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/loop.rb:30:in `each'
1622
- 2: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/loop.rb:31:in `block (2 levels) in enter_main_loop'
1623
- 1: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/menu.rb:3565:in `menu'
1624
- /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/misc.rb:3979:in `append': undefined method `return_sequence_that_is_cut_via_restriction_enzyme' for Bioroebe:Module (NoMethodError)
1625
- Did you mean? return_random_codon_sequence_for_this_aminoacid_sequence
1626
-
1627
-
1628
- ^^^^^ BUG!
1629
- -------------------------------------------------------------------------------
1630
- > rest?
1631
-
1632
- We found these restriction sites within the sequence `TTCAGAACTCAACGCCTGGTTGGCCGTCCAGTAAGCTGACTAAGTAAGTCTATGCCCGCGATAACCAGGATACAGATATCGTGAAACCTGGTTTATCTCCTTCTATAAGAGTCTGCACATCTAGC`:
1633
-
1634
- AccII → CGCG ( 1 times found)
1635
- AluIAGCT ( 1 times found)
1636
- BfaI CTAG ( 1 times found)
1637
- BshI GGCC ( 1 times found)
1638
- Bsh1236I CGCG ( 1 times found)
1639
- BshFI GGCC ( 1 times found)
1640
- BstFNI CGCG ( 1 times found)
1641
- BstUI CGCG ( 1 times found)
1642
- BsuRI -> GGCC ( 1 times found)
1643
- CviRI -> TGCA ( 1 times found)
1644
- Eco32I -> GATATC ( 1 times found)
1645
- EcoRV -> GATATC ( 1 times found)
1646
- FnuDII -> CGCG ( 1 times found)
1647
- HaeIII -> GGCC ( 1 times found)
1648
- HpyCH4V -> TGCA ( 1 times found)
1649
- MaeI -> CTAG ( 1 times found)
1650
- MvnI -> CGCG ( 1 times found)
1651
- PalI -> GGCC ( 1 times found)
1652
- SelI -> CGCG ( 1 times found)
1653
- ThaI -> CGCG ( 1 times found)
1654
- XspI CTAG ( 1 times found)
1655
-
1656
- ^^^^ also show the position
1657
-
1658
-
1659
- -------------------------------------------------------------------------------
1660
-
1661
- PMID entries are:
1662
-
1663
- x = 'Goldman, JM, Melo JV 2003 NEJM 349:1451 14534339
1664
- Lewis GD 1993 Cancer Immunol Immun other 37: 255 8102322
1665
- McShane LM 2009 Clin Canc Res 15: 1898 19276274
1666
- Fox JL 2007 Nature Biotech 25: 489 17483821
1667
- Bodin L 2005 Blood 106: 135 15790782'
1668
-
1669
- require 'rubygems'
1670
- require 'bio'
1671
-
1672
- my_file = File.new(ARGV[0])
1673
- refs = my_file.readlines
1674
- ids = []
1675
-
1676
- refs.each do |line|
1189
+ also set it to
1190
+ set_restriction_size()
1191
+ --------------------------------------------------------------------------------
1192
+ (186) → also ... while learning C++ we extend the project here...
1193
+ Useful C++ things will be combined.
1194
+ --------------------------------------------------------------------------------
1195
+ (187) → As of April 2003, there were 176,890 total taxa represented.
1196
+ ^^^ we need a way to also output how many entries we
1197
+ have there.
1198
+ --------------------------------------------------------------------------------
1199
+ (188) → Replace bioruby with bioroebe completely!
1200
+ In order for this to work, we first need to find out
1201
+ what bioruby is able to do. :P
1202
+ --------------------------------------------------------------------------------
1203
+ (189) append 33
1204
+ # ^^^ in the bioshell
1205
+ Only numbers were given: Adding 33 random nucleotides to the main string next.
1206
+ Traceback (most recent call last):
1207
+ 10: from /usr/bin/bioshell:26:in `<main>'
1208
+ 9: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/misc.rb:4121:in `shell'
1209
+ 8: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/misc.rb:4121:in `new'
1210
+ 7: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/initialize.rb:168:in `initialize'
1211
+ 6: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/loop.rb:18:in `enter_main_loop'
1212
+ 5: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/loop.rb:18:in `loop'
1213
+ 4: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/loop.rb:30:in `block in enter_main_loop'
1214
+ 3: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/loop.rb:30:in `each'
1215
+ 2: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/loop.rb:31:in `block (2 levels) in enter_main_loop'
1216
+ 1: from /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/menu.rb:3565:in `menu'
1217
+ /home/Programs/Ruby/2.7.1/lib/ruby/site_ruby/2.7.0/bioroebe/shell/misc.rb:3979:in `append': undefined method `return_sequence_that_is_cut_via_restriction_enzyme' for Bioroebe:Module (NoMethodError)
1218
+ Did you mean? return_random_codon_sequence_for_this_aminoacid_sequence
1219
+ ^^^^^ BUG!
1220
+ --------------------------------------------------------------------------------
1221
+ (190) → > rest?
1222
+ We found these restriction sites within the sequence `TTCAGAACTCAACGCCTGGTTGGCCGTCCAGTAAGCTGACTAAGTAAGTCTATGCCCGCGATAACCAGGATACAGATATCGTGAAACCTGGTTTATCTCCTTCTATAAGAGTCTGCACATCTAGC`:
1223
+ AccII CGCG ( 1 times found)
1224
+ AluI → AGCT ( 1 times found)
1225
+ BfaI → CTAG ( 1 times found)
1226
+ BshI → GGCC ( 1 times found)
1227
+ Bsh1236I CGCG ( 1 times found)
1228
+ BshFI → GGCC ( 1 times found)
1229
+ BstFNI → CGCG ( 1 times found)
1230
+ BstUICGCG ( 1 times found)
1231
+ BsuRI -> GGCC ( 1 times found)
1232
+ CviRI -> TGCA ( 1 times found)
1233
+ Eco32I -> GATATC ( 1 times found)
1234
+ EcoRV -> GATATC ( 1 times found)
1235
+ FnuDII -> CGCG ( 1 times found)
1236
+ HaeIII -> GGCC ( 1 times found)
1237
+ HpyCH4V -> TGCA ( 1 times found)
1238
+ MaeI -> CTAG ( 1 times found)
1239
+ MvnI -> CGCG ( 1 times found)
1240
+ PalI -> GGCC ( 1 times found)
1241
+ SelI -> CGCG ( 1 times found)
1242
+ ThaI -> CGCG ( 1 times found)
1243
+ XspI CTAG ( 1 times found)
1244
+ ^^^^ also show the position
1245
+ --------------------------------------------------------------------------------
1246
+ (191) PMID entries are:
1247
+ x = 'Goldman, JM, Melo JV 2003 NEJM 349:1451 14534339
1248
+ Lewis GD 1993 Cancer Immunol Immun other 37: 255 8102322
1249
+ McShane LM 2009 Clin Canc Res 15: 1898 19276274
1250
+ Fox JL 2007 Nature Biotech 25: 489 17483821
1251
+ Bodin L 2005 Blood 106: 135 15790782'
1252
+ require 'rubygems'
1253
+ require 'bio'
1254
+ my_file = File.new(ARGV[0])
1255
+ refs = my_file.readlines
1256
+ ids = []
1257
+ refs.each do |line|
1677
1258
  pmid = line.strip().split("\t")
1678
1259
  ids.push(pmid[2])
1679
- end
1680
-
1681
- ids.each { |id|
1260
+ end
1261
+ ids.each { |id|
1682
1262
  entry = Bio::PubMed.query(id)
1683
1263
  medline = Bio::MEDLINE.new(entry)
1684
1264
  reference = medline.reference
1685
1265
  puts reference.endnote
1686
- }
1687
- - Clone, in bioroebe, the get_ORF functionality.
1688
-
1689
- - in bioroebe: Find out how to count where a restriction
1690
- enzyme was found + add it into a table, also WHEN it
1691
- was found and by WHO.
1692
-
1693
- We should make a good reference table there, so that
1694
- people can reproduce where the information is kept
1695
- or obtained from.
1696
-
1697
- - Enable:
1698
- download ecoli
1699
- first fix a bug
1700
- and then make it so that we can download the ecoli
1701
- sequence from that file.... yaml file
1702
- - Rewrite Bioroebe (Pending since as of September 2019 ...)
1703
- Start with the functionality that shall decode something,
1704
- and put this back ... as a first class citizen with
1705
- descriptions how to work with this.
1706
- - improve quality
1707
- - also fix the taxonomy project as you go...
1708
- - Add C++ wrapper stuff, starting with qt-widgets specifically
1709
- but also extend the base part of bioroebe as far as C++ is
1710
- concerned. Perhaps at a later time also add/embed mruby
1711
- into the project.
1712
-
1713
- - Need to overlay 2 exceptionally large DNA data sets and
1714
- analyze the overlap.
1715
-
1716
- Want to determine the frequency of one event within the DNA
1717
- reads of the other data set. Prefer to discuss the details
1718
- of the data sets one-on-one.
1719
- ^^^ we need to compare two DNA data sets and analyze the overlap.
1720
-
1721
- - Make sure that the sinatra interface has the parts of
1722
- emboss available, too. List the progress here.
1723
- - interactivefasta interactive_fasta
1724
- ^^^ enable this... but I am not sure what is meant with
1725
- that. :\
1726
-
1727
-
1728
-
1729
-
1730
-
1731
- ..........................................................................
1732
- Bei der Datenbanksuche werden die gemessenen Massen mit den Peptidmassen
1733
- aller Proteine bzw. Gene in einer Datenbank (NCBI, Uniprot) verglichen. DNA-
1734
- Sequenzen werden dazu in Proteinsequenzen übersetzt und in silico mit der beim
1735
- Verdau benutzten Protease geschnitten.
1736
- ^^^ enable digestions
1737
-
1738
-
1739
-
1740
-
1741
- ..........................................................................
1742
- Complexity of libraries:
1743
- How many independent clones are necessary to represent a genome (plant,
1744
- animal/fungus) or how many such clones have to be screened to have realistic
1745
- chance of finding the gene of interest?
1746
- This can be calculated by the formula:
1747
- P = 1 - (1 - F/G) N
1748
- N = ln(1 - P) / ln(1 - F/G)
1749
- P.... Propability that a certain insert is present in library consisting of N clones
1750
- F.... average insert size (kb)
1751
- G... Genome size (kb) of the organism from which a library should be constructed
1752
- N.... Number of clones in the library
1753
- ^^^^ add this formula + documentation
1754
-
1755
- Example: 16 kb average insert size in a replacement vector
1756
- Genome sizes:
1757
- Yeast:
1758
- 16 Mb = 16 000 kb (1000 clones with 16 kb = 1 genome equivalent)
1759
- F/G=0.001.
1760
- ln(1-F/G)= - 0.0010005
1761
- 95%: ln0.05= - 2.99573
1762
- 2,99/0.001 3000
1763
- Wheat: 16 000 Mb
1764
- 3 .10 6 clones
1765
- How many plaques can be screended on one 9 cm petri dish or filter?
1766
- Plaques close to confluency: about 10 4 pfu per plate.
1767
- While in case of yeast the whole library can be screened on 1 plate,
1768
- 300 plates would be needed for wheat - which is impracticable!
1769
- (BAC "bacterial artificial chromosome" libraries with 150 -300 kb inserts are used!)
1770
- Most plants have reasonable genome size (e.g. tomato about 800 Mb) - 15 filters
1771
- have to be hybridized.
1772
-
1773
-
1774
-
1775
-
1776
- ..........................................................................
1777
-
1778
- BIO SHELL> BglI?
1779
-
1780
- We have found a restriction enzyme called BglI.
1781
-
1782
- The sequence this 11-cutter relates to is: `5' - GCCNNNNNGGC - 3'`
1783
-
1784
- It will specifically cut between: 5' - GCCNNNN|NGGC - 3'
1785
- 3' - CGG(A/T/G/C)(A/T/G/C)(A/T/G/C)(A/T/G/C)(A/T/|G/C)CCG - 5'
1786
-
1787
- ^
1788
- also add a line to say
1789
- "This is a blunt-end cutter."
1790
- or
1791
- "This enzyme cre
1792
- ates a staggered cut."
1793
-
1794
-
1795
- http://biopython.org/DIST/docs/api/Bio.Restriction.Restriction.Blunt-class.html
1796
- ^^^ also add something like this... as a query thingy.
1797
-
1798
- catalyse(cls, dna, linear=True)
1799
- List the sequence fragments after cutting dna with enzyme. source code
1800
-
1801
- catalyze(cls, dna, linear=True)
1802
- List the sequence fragments after cutting dna with enzyme. source code
1803
-
1804
- is_blunt(cls)
1805
- Return if the enzyme produces blunt ends. source code
1806
-
1807
- is_5overhang(cls)
1808
- Return if the enzymes produces 5' overhanging ends. source code
1809
-
1810
- is_3overhang(cls)
1811
- Return if the enzyme produces 3' overhanging ends. source code
1812
-
1813
- overhang(cls)
1814
- Return the type of the enzyme's overhang as string. source code
1815
-
1816
- compatible_end(cls, batch=None)
1817
- List all enzymes that produce compatible ends for the enzyme.
1818
- http://biopython.org/DIST/docs/api/Bio.Restriction.Restriction.Blunt-class.html
1819
-
1820
-
1821
- ..........................................................................
1822
- https://www.reddit.com/r/bioinformatics/comments/5o3kn8/bioinformatics_contest_2017_jan_23rd29th_solve_as/
1823
- ..........................................................................
1824
- (1) Finish all of biophp integration into bioroebe.
1825
- http://www.biophp.org/
1826
- -------------------------------------------------------------------------------
1827
-
1828
- locate oriC here:
1829
-
1830
- ttcgttaagtaacttcactgcccgtagtgtaccggcattcgctagcaagagtctttctg
1831
- ggcaagcttcacttgtgatcgcggcctgtgcccccggaatgaaacaaccacgtccctgct
1832
- aacaacgacgggaaaagggaagtgatccgtcggcagacccagactagtgcccttctccgg
1833
- cttccaacaccaacgagtcggaccgaattgagcactcgaatgcacggcgctttttgccgg
1834
- ccgaaacggcgcctccgcattgatcgacgcacggcctcttttggctacagcgcatggctt
1835
- tacactcggcatgcatttccagtgctaatcaaacagaattccttgtaaagtccttcaacc
1836
- gtgacagactatcgctaaggagcctttccagtcgtgcctgcaatcactcgcgaaatcaac
1837
- aaatctacatctaagcacgctcgtggttcggagtcccgccctcatgtggaccatagccgg
1838
- ttcgcccgagtcctaggcacgatcagaggacctatctttcgccactcaactcttctgagt
1839
- gaaacaatatcgaccgaaaccttgctcggttttgtccacaacaacgtcaggcccataagc
1840
- agacgacattagtccgctgttgtcgcgggcgtcccatagccgtacgatgtcccgtcgga
1841
-
1842
- ori?
1843
-
1844
- ^^^ this shall give us all ORI in a sequence.
1845
- DnaA protein binds to DNA-Box in an ori.
1846
-
1847
- '9cut'
1848
- '8cut'
1849
- '7cut'
1850
- ^^^ these give us slices
1851
-
1852
- But I do not know how to locate ORIs.
1853
-
1854
-
1855
-
1856
-
1857
-
1858
-
1859
-
1860
-
1861
- -------------------------------------------------------------------------------
1862
- ^^^ also integrate git into bioroebe.
1863
- -------------------------------------------------------------------------------
1864
- WIR MÜSSEN DAS HIER EXTREM VERBESSERN.
1865
-
1866
- DANN UPLOADEN UND ALS BASIS FÜR APPLICATIONS NUTZEN.
1867
- -------------------------------------------------------------------------------
1868
-
1869
- Study MetaCyc
1870
- ^^^ study metabolic pathways.
1871
-
1872
- http://metacyc.org/
1873
-
1874
- → Create KuroMetaCyc, in Analogy towards Metabolic Cycle.
1875
-
1876
- -------------------------------------------------------------------------------
1877
-
1878
- Welcome to BioShell May 2012. Type "help" to get some help.
1879
-
1880
- Hello and welcome to the Bio Shell Version, last updated: May 2012
1881
-
1882
- BIO SHELL> IPEYVDWRQKGAVTPVKNQGSCGSCWAFSAVVTIEGIIKIRTGNLNQYSEQELLDCDRRSYGCNGGYPWSALQLVAQYGI
1883
- BIO SHELL> HYRNTYPYEGVQRYCRSREKGPYAAKTDGVRQVQPYNQGALLYSIANQPVSVVLQAAGKDFQLYRGGIFVGPCGNKVDHA
1884
- BIO SHELL> VAAVGYGPNYILIKNSWGTGWGENGYIRIKRGTGNSYGVCGLYTSSFYPVKN
1885
- BIO SHELL>
1886
- BIO SHELL> input?
1887
- BIO SHELL> pdb
1888
- BIO SHELL>
1889
-
1890
- ^^^ (1)
1891
- Add a pdb submodule
1892
- When we type this, we then ask:
1893
- "Please input your FASTA format now:"
1894
-
1895
-
1896
-
1897
-
1898
- -------------------------------------------------------------------------------
1899
-
1900
- http://biopython.org/DIST/docs/cookbook/Restriction.html#mozTocId101269
1901
-
1902
- ^^^ support this also:
1903
-
1904
- >>> from Bio import Restriction
1905
- >>> dir()
1906
- ['Restriction', '__builtins__', '__doc__', '__name__']
1907
- >>> Restriction.EcoRI
1908
- EcoRI
1909
- >>> Restriction.EcoRI.site
1910
- 'GAATTC'
1911
- >>>
1912
-
1913
- and document it somewhere; perhaps in a new .cgi page.
1914
-
1915
- The above will return the exact site, without verbosity.
1916
-
1917
-
1918
- Restriction.EcoRI
1919
- Restriction.EcoRI.site
1920
-
1921
- result = Bioroebe.restriction_enzyme 'EcoRI.site'
1266
+ }
1267
+ - Clone, in bioroebe, the get_ORF functionality.
1268
+ - in bioroebe: Find out how to count where a restriction
1269
+ enzyme was found + add it into a table, also WHEN it
1270
+ was found and by WHO.
1271
+ We should make a good reference table there, so that
1272
+ people can reproduce where the information is kept
1273
+ or obtained from.
1274
+ - Enable:
1275
+ download ecoli
1276
+ first fix a bug
1277
+ and then make it so that we can download the ecoli
1278
+ sequence from that file.... yaml file
1279
+ - Rewrite Bioroebe (Pending since as of September 2019 ...)
1280
+ Start with the functionality that shall decode something,
1281
+ and put this back ... as a first class citizen with
1282
+ descriptions how to work with this.
1283
+ - improve quality
1284
+ - also fix the taxonomy project as you go...
1285
+ - Add C++ wrapper stuff, starting with qt-widgets specifically
1286
+ but also extend the base part of bioroebe as far as C++ is
1287
+ concerned. Perhaps at a later time also add/embed mruby
1288
+ into the project.
1289
+ - Need to overlay 2 exceptionally large DNA data sets and
1290
+ analyze the overlap.
1291
+ Want to determine the frequency of one event within the DNA
1292
+ reads of the other data set. Prefer to discuss the details
1293
+ of the data sets one-on-one.
1294
+ ^^^ we need to compare two DNA data sets and analyze the overlap.
1295
+ - Make sure that the sinatra interface has the parts of
1296
+ emboss available, too. List the progress here.
1297
+ - interactivefasta interactive_fasta
1298
+ ^^^ enable this... but I am not sure what is meant with
1299
+ that. :\
1300
+ --------------------------------------------------------------------------------
1301
+ (192) Bei der Datenbanksuche werden die gemessenen Massen mit den Peptidmassen
1302
+ aller Proteine bzw. Gene in einer Datenbank (NCBI, Uniprot) verglichen. DNA-
1303
+ Sequenzen werden dazu in Proteinsequenzen übersetzt und in silico mit der beim
1304
+ Verdau benutzten Protease geschnitten.
1305
+ ^^^ enable digestions
1306
+ --------------------------------------------------------------------------------
1307
+ (193) → Complexity of libraries:
1308
+ How many independent clones are necessary to represent a genome (plant,
1309
+ animal/fungus) or how many such clones have to be screened to have realistic
1310
+ chance of finding the gene of interest?
1311
+ This can be calculated by the formula:
1312
+ P = 1 - (1 - F/G) N
1313
+ N = ln(1 - P) / ln(1 - F/G)
1314
+ P.... Propability that a certain insert is present in library consisting of N clones
1315
+ F.... average insert size (kb)
1316
+ G... Genome size (kb) of the organism from which a library should be constructed
1317
+ N.... Number of clones in the library
1318
+ ^^^^ add this formula + documentation
1319
+ Example: 16 kb average insert size in a replacement vector
1320
+ Genome sizes:
1321
+ Yeast:
1322
+ 16 Mb = 16 000 kb (1000 clones with 16 kb = 1 genome equivalent)
1323
+ F/G=0.001.
1324
+ ln(1-F/G)= - 0.0010005
1325
+ 95%: ln0.05= - 2.99573
1326
+ 2,99/0.001 3000
1327
+ Wheat: 16 000 Mb
1328
+ 3 .10 6 clones
1329
+ How many plaques can be screended on one 9 cm petri dish or filter?
1330
+ Plaques close to confluency: about 10 4 pfu per plate.
1331
+ While in case of yeast the whole library can be screened on 1 plate,
1332
+ 300 plates would be needed for wheat - which is impracticable!
1333
+ (BAC "bacterial artificial chromosome" libraries with 150 -300 kb inserts are used!)
1334
+ Most plants have reasonable genome size (e.g. tomato about 800 Mb) - 15 filters
1335
+ have to be hybridized.
1336
+ --------------------------------------------------------------------------------
1337
+ (194) → BIO SHELL> BglI?
1338
+ We have found a restriction enzyme called BglI.
1339
+ The sequence this 11-cutter relates to is: `5' - GCCNNNNNGGC - 3'`
1340
+ It will specifically cut between: 5' - GCCNNNN|NGGC - 3'
1341
+ 3' - CGG(A/T/G/C)(A/T/G/C)(A/T/G/C)(A/T/G/C)(A/T/|G/C)CCG - 5'
1342
+ ^
1343
+ also add a line to say
1344
+ "This is a blunt-end cutter."
1345
+ or
1346
+ "This enzyme cre
1347
+ ates a staggered cut."
1348
+ http://biopython.org/DIST/docs/api/Bio.Restriction.Restriction.Blunt-class.html
1349
+ ^^^ also add something like this... as a query thingy.
1350
+ catalyse(cls, dna, linear=True)
1351
+ List the sequence fragments after cutting dna with enzyme. source code
1352
+ catalyze(cls, dna, linear=True)
1353
+ List the sequence fragments after cutting dna with enzyme. source code
1354
+ is_blunt(cls)
1355
+ Return if the enzyme produces blunt ends. source code
1356
+ is_5overhang(cls)
1357
+ Return if the enzymes produces 5' overhanging ends. source code
1358
+ is_3overhang(cls)
1359
+ Return if the enzyme produces 3' overhanging ends. source code
1360
+ overhang(cls)
1361
+ Return the type of the enzyme's overhang as string. source code
1362
+ compatible_end(cls, batch=None)
1363
+ List all enzymes that produce compatible ends for the enzyme.
1364
+ http://biopython.org/DIST/docs/api/Bio.Restriction.Restriction.Blunt-class.html
1365
+ --------------------------------------------------------------------------------
1366
+ (195) → https://www.reddit.com/r/bioinformatics/comments/5o3kn8/bioinformatics_contest_2017_jan_23rd29th_solve_as/
1367
+ --------------------------------------------------------------------------------
1368
+ (196) Finish all of biophp integration into bioroebe.
1369
+ http://www.biophp.org/
1370
+ --------------------------------------------------------------------------------
1371
+ (197) locate oriC here:
1372
+ ttcgttaagtaacttcactgcccgtagtgtaccggcattcgctagcaagagtctttctg
1373
+ ggcaagcttcacttgtgatcgcggcctgtgcccccggaatgaaacaaccacgtccctgct
1374
+ aacaacgacgggaaaagggaagtgatccgtcggcagacccagactagtgcccttctccgg
1375
+ cttccaacaccaacgagtcggaccgaattgagcactcgaatgcacggcgctttttgccgg
1376
+ ccgaaacggcgcctccgcattgatcgacgcacggcctcttttggctacagcgcatggctt
1377
+ tacactcggcatgcatttccagtgctaatcaaacagaattccttgtaaagtccttcaacc
1378
+ gtgacagactatcgctaaggagcctttccagtcgtgcctgcaatcactcgcgaaatcaac
1379
+ aaatctacatctaagcacgctcgtggttcggagtcccgccctcatgtggaccatagccgg
1380
+ ttcgcccgagtcctaggcacgatcagaggacctatctttcgccactcaactcttctgagt
1381
+ gaaacaatatcgaccgaaaccttgctcggttttgtccacaacaacgtcaggcccataagc
1382
+ agacgacattagtccgctgttgtcgcgggcgtcccatagccgtacgatgtcccgtcgga
1383
+ ori?
1384
+ ^^^ this shall give us all ORI in a sequence.
1385
+ DnaA protein binds to DNA-Box in an ori.
1386
+ '9cut'
1387
+ '8cut'
1388
+ '7cut'
1389
+ ^^^ these give us slices
1390
+ But I do not know how to locate ORIs.
1391
+ --------------------------------------------------------------------------------
1392
+ (198) → ^^^ also integrate git into bioroebe.
1393
+ --------------------------------------------------------------------------------
1394
+ (199) WIR MÜSSEN DAS HIER EXTREM VERBESSERN.
1395
+ DANN UPLOADEN UND ALS BASIS FÜR APPLICATIONS NUTZEN.
1396
+ --------------------------------------------------------------------------------
1397
+ (200) Study MetaCyc
1398
+ ^^^ study metabolic pathways.
1399
+ http://metacyc.org/
1400
+ → Create KuroMetaCyc, in Analogy towards Metabolic Cycle.
1401
+ --------------------------------------------------------------------------------
1402
+ (201) → Welcome to BioShell May 2012. Type "help" to get some help.
1403
+ Hello and welcome to the Bio Shell Version, last updated: May 2012
1404
+ BIO SHELL> IPEYVDWRQKGAVTPVKNQGSCGSCWAFSAVVTIEGIIKIRTGNLNQYSEQELLDCDRRSYGCNGGYPWSALQLVAQYGI
1405
+ BIO SHELL> HYRNTYPYEGVQRYCRSREKGPYAAKTDGVRQVQPYNQGALLYSIANQPVSVVLQAAGKDFQLYRGGIFVGPCGNKVDHA
1406
+ BIO SHELL> VAAVGYGPNYILIKNSWGTGWGENGYIRIKRGTGNSYGVCGLYTSSFYPVKN
1407
+ BIO SHELL>
1408
+ BIO SHELL> input?
1409
+ BIO SHELL> pdb
1410
+ BIO SHELL>
1411
+ ^^^ (1)
1412
+ Add a pdb submodule
1413
+ When we type this, we then ask:
1414
+ "Please input your FASTA format now:"
1415
+ --------------------------------------------------------------------------------
1416
+ (202) → http://biopython.org/DIST/docs/cookbook/Restriction.html#mozTocId101269
1417
+ ^^^ support this also:
1418
+ >>> from Bio import Restriction
1419
+ >>> dir()
1420
+ ['Restriction', '__builtins__', '__doc__', '__name__']
1421
+ >>> Restriction.EcoRI
1422
+ EcoRI
1423
+ >>> Restriction.EcoRI.site
1424
+ 'GAATTC'
1425
+ >>>
1426
+ and document it somewhere; perhaps in a new .cgi page.
1427
+ The above will return the exact site, without verbosity.
1428
+ Restriction.EcoRI
1429
+ Restriction.EcoRI.site
1430
+ result = Bioroebe.restriction_enzyme 'EcoRI.site'
1922
1431
  # => "GAATTC"
1923
-
1924
- ^^^ funktioniert bereits Teilweise, aber noch nit
1925
- ausreichend.
1926
-
1927
-
1928
-
1929
- - Clone the following functionality from Bio:
1930
-
1931
- require 'bio'
1932
- quality_threshold = 60
1933
- Bio::FlatFile.open('sample.fastq').each {|entry|
1432
+ ^^^ funktioniert bereits Teilweise, aber noch nit
1433
+ ausreichend.
1434
+ - Clone the following functionality from Bio:
1435
+ require 'bio'
1436
+ quality_threshold = 60
1437
+ Bio::FlatFile.open('sample.fastq').each {|entry|
1934
1438
  hq_seq = entry.mask(quality_threshold)
1935
1439
  puts hq_seq.output_fasta(entry.entry_id)
1936
- }
1937
-
1938
-
1939
- - Document the workflow of all scripts in a reproducible
1940
- manner, e. g. so that others can use it...
1941
-
1942
-
1943
- it must also allow for different tables to be
1944
- used! check this... so that we can search in
1945
- standard ORF but also in different ORFs
1946
-
1947
- und die länge angeben, zumindest vom längsten ORF
1948
- start + stop... also so das das ergebnis auch
1949
- passt
1950
-
1951
-
1952
-
1953
- - Before we start to use rails for bioroebe, let's polish the GUI
1954
- components more.
1955
- But we really should use rails for this project too, at the
1956
- least optionally.
1957
- Ideally every small class has a tiny widget that can be
1958
- interconnected. Perhaps do this via sinatra first and think
1959
- of ways how to generalize on it.
1960
-
1961
- We should probably be systematic about this and go through
1962
- each class, then write a GUI for it:
1963
- ^^^ first GUI
1964
- ^^^ then rails
1965
-
1966
- Ok, first task - write a lot of GUIs.
1967
- gui/hamming_distance.rb - works ok-ish but it needs an in-widget
1968
- notification, that is, we need to somehow colourize this thing.
1969
-
1970
-
1971
-
1972
- - fix batch generation of .star files
1973
-
1974
- - add script that will generate those weird
1975
- files ... from the tch shell script
1976
-
1977
-
1978
- - in mohitstar --help auch examples hinzufügen
1979
- --examples
1980
- we also need the Z coordinates into the .star file
1981
- test.xmd
1982
-
1983
- /home/kumar/Desktop/test_exit/InputMicrographs/
1984
-
1985
-
1986
-
1987
-
1988
- ..........................................................................
1989
- BioTodo - GENESIS, science fiction.
1990
-
1991
- - create virus(:which_one, :amount) # Note the difference to the below
1992
- - create hydra(:amount)
1993
- - create bread
1994
- ..........................................................................
1995
- → both
1996
- ^ should work, does not work right now.
1997
- ..........................................................................
1998
- → Taxonomy is now integrated into bioroebe. This is good but we need more
1999
- documentation, some more tests, a rethinking of the layout and the
2000
- structures, and a fixing of the query-part of the database.
2001
-
2002
- Also, make sure that it does the main functions.
2003
-
2004
- rewrite for taxonomy in bioroebe
2005
- and while doing this, also continue with the
2006
- protokoll in bioinformatik
2007
- so that we can finish both related-problemsters
2008
- at about the same time \o/
2009
- AND document this related-problems too
2010
- Integrate this some other day...
2011
- ..........................................................................
2012
- - http://www.restrictionmapper.org/cgi-bin/sitefind3.pl
2013
-
2014
- ^^^ Das sollte man integrieren, die Funktionalität, so das
1440
+ }
1441
+ - Document the workflow of all scripts in a reproducible
1442
+ manner, e. g. so that others can use it...
1443
+ it must also allow for different tables to be
1444
+ used! check this... so that we can search in
1445
+ standard ORF but also in different ORFs
1446
+ und die länge angeben, zumindest vom längsten ORF
1447
+ start + stop... also so das das ergebnis auch
1448
+ passt
1449
+ - Before we start to use rails for bioroebe, let's polish the GUI
1450
+ components more.
1451
+ But we really should use rails for this project too, at the
1452
+ least optionally.
1453
+ Ideally every small class has a tiny widget that can be
1454
+ interconnected. Perhaps do this via sinatra first and think
1455
+ of ways how to generalize on it.
1456
+ We should probably be systematic about this and go through
1457
+ each class, then write a GUI for it:
1458
+ ^^^ first GUI
1459
+ ^^^ then rails
1460
+ Ok, first task - write a lot of GUIs.
1461
+ gui/hamming_distance.rb - works ok-ish but it needs an in-widget
1462
+ notification, that is, we need to somehow colourize this thing.
1463
+ - fix batch generation of .star files
1464
+ - add script that will generate those weird
1465
+ files ... from the tch shell script
1466
+ - in mohitstar --help auch examples hinzufügen
1467
+ --examples
1468
+ we also need the Z coordinates into the .star file
1469
+ test.xmd
1470
+ /home/kumar/Desktop/test_exit/InputMicrographs/
1471
+ --------------------------------------------------------------------------------
1472
+ (203) BioTodo - GENESIS, science fiction.
1473
+ - create virus(:which_one, :amount) # Note the difference to the below
1474
+ - create hydra(:amount)
1475
+ - create bread
1476
+ --------------------------------------------------------------------------------
1477
+ (204) → both
1478
+ ^ should work, does not work right now.
1479
+ --------------------------------------------------------------------------------
1480
+ (205) → Taxonomy is now integrated into bioroebe. This is good but we need more
1481
+ documentation, some more tests, a rethinking of the layout and the
1482
+ structures, and a fixing of the query-part of the database.
1483
+ Also, make sure that it does the main functions.
1484
+ rewrite for taxonomy in bioroebe
1485
+ and while doing this, also continue with the
1486
+ protokoll in bioinformatik
1487
+ so that we can finish both related-problemsters
1488
+ at about the same time \o/
1489
+ AND document this related-problems too
1490
+ Integrate this some other day...
1491
+ --------------------------------------------------------------------------------
1492
+ (206) → http://www.restrictionmapper.org/cgi-bin/sitefind3.pl
1493
+ ^^^ Das sollte man integrieren, die Funktionalität, so das
2015
1494
  man ALLE Restriktion-Enzymes ausprobiert ausgehend von
2016
1495
  einer bestimmten Sequenz.
2017
- ..........................................................................
2018
- A search is essentially substring search across a database of strings
2019
- (albeit with a smaller alphabet). Some common use cases: one,
2020
- scientists will search for certain genes that they've used in engineered
2021
- plasmids. Two, since multiple codons can translate to the same amino
2022
- acid, a process called "codon optimization" might replace codons with
2023
- equivalent ones that work better given a certain kind of organism -
2024
- your CAA (glutamine) might become CAG (still glutamine) instead. You
2025
- might want to search by amino acid ("glutamine") so that you find
2026
- it regardless of whether CAA or CAG was used.
2027
-
2028
- ^^^^ yeah enable this too.
2029
- Perhaps start with a "codon optimizer". This shall, for a given
2030
- organism, replace the given codons with the "optimal" ones.
2031
- Then also add this to sinatra interafce.
2032
-
2033
- Bioroebe::DetermineOptimalCodons
2034
- ^^^ this is currently incomplete.
2035
-
2036
- ..........................................................................
2037
- Redo restrictions enzymes completely.
2038
- And polish this a LOT.
2039
- This may take some days. But we want this to be REALLY good and
2040
- lasting for a long time.
2041
- Need to keep on working at that!
2042
- ..........................................................................
2043
- Add: average_aminoacid_weight?
2044
-
2045
-
2046
- === LV-Nummer 300214 UE Übung III B Sequenzanalysen in der Molekularbiologie
2047
- Pubmed
2048
- Finding sequences
2049
- Sequence homology search (Blast, FASTA)
2050
- Pairwise sequence alignment
2051
- DNA analysis, translation
2052
- Gene finding in genomes
2053
- ^^^ find genes
2054
- Plasmid Cloning
2055
- Primer for PCR
2056
- Protein analysis
2057
- PROSITE
2058
- EMBOSS
2059
- ^^^ das nochmals in Ruhe durchgehen.
2060
-
2061
- in bioroebe viel C++ schreiben, möglichst optimieren,
2062
- und dann irgendwann später ruby bindings dazu liefern
2063
- add the ability to compare two FASTA files
2064
- can probably do so via the two sequences.
2065
- integrate "grundlagen der bioinfo" slides
2066
- gleichzeitig während wir dafür lernen.
2067
- ntSeq.each_entry do |e|
2068
- # pep = Bio::Sequence.new(e.naseq.reverse_complement!.translate)
2069
- pep = Bio::Sequence.new(e.naseq.translate)
2070
- if strand == 0
2071
- pep = Bio::Sequence.new(e.naseq.reverse_complement!.translate)
2072
- end
2073
- puts pep.output_fasta(e.definition,20)
2074
- end
2075
- ^^^ we must show this in n characters per line:
2076
- see https://raw.githubusercontent.com/zorino/bioruby-scripts/master/transeq.rb
2077
- We must be able to align not only nucleotides but also aminoacids.
2078
- But where is the alignment comparer? perhaps hamming distance?
2079
- hmm we have to see.
2080
- ..........................................................................
2081
- /Programs/Ruby/2.3.1/lib/ruby/site_ruby/2.3.0/bioroebe/bioshell/menu.rb:311:in `menu': undefined method `upcase' for ["EcoRI"]:Array (NoMethodError)
2082
- from /Programs/Ruby/2.3.1/lib/ruby/site_ruby/2.3.0/bioroebe/bioshell/user_input.rb:31:in `block in enter_main_loop'
2083
- from /Programs/Ruby/2.3.1/lib/ruby/site_ruby/2.3.0/bioroebe/bioshell/user_input.rb:12:in `loop'
2084
- from /Programs/Ruby/2.3.1/lib/ruby/site_ruby/2.3.0/bioroebe/bioshell/user_input.rb:12:in `enter_main_loop'
2085
- from /Programs/Ruby/2.3.1/lib/ruby/site_ruby/2.3.0/bioroebe/bioshell/initialize.rb:42:in `initialize'
2086
- from /Programs/Ruby/2.3.1/lib/ruby/site_ruby/2.3.0/bioroebe/bioshell/bioshell.rb:52:in `new'
2087
- from /Programs/Ruby/2.3.1/lib/ruby/site_ruby/2.3.0/bioroebe/bioshell/bioshell.rb:52:in `shell'
2088
- from /System/Executables/bioshell:6:in `<main>'
2089
-
2090
- ^^^ also fix this stupid bug.
2091
- perhaps redo the whole restriction enzyme stuff.
2092
-
2093
- → Taxonomy components:
2094
-
2095
- (1) Wir wollen nur einen query machen statt vieler kleiner
1496
+ --------------------------------------------------------------------------------
1497
+ (207) A search is essentially substring search across a database of strings
1498
+ (albeit with a smaller alphabet). Some common use cases: one,
1499
+ scientists will search for certain genes that they've used in engineered
1500
+ plasmids. Two, since multiple codons can translate to the same amino
1501
+ acid, a process called "codon optimization" might replace codons with
1502
+ equivalent ones that work better given a certain kind of organism -
1503
+ your CAA (glutamine) might become CAG (still glutamine) instead. You
1504
+ might want to search by amino acid ("glutamine") so that you find
1505
+ it regardless of whether CAA or CAG was used.
1506
+ ^^^^ yeah enable this too.
1507
+ Perhaps start with a "codon optimizer". This shall, for a given
1508
+ organism, replace the given codons with the "optimal" ones.
1509
+ Then also add this to sinatra interafce.
1510
+ Bioroebe::DetermineOptimalCodons
1511
+ ^^^ this is currently incomplete.
1512
+ --------------------------------------------------------------------------------
1513
+ (208) → Redo restrictions enzymes completely.
1514
+ And polish this a LOT.
1515
+ This may take some days. But we want this to be REALLY good and
1516
+ lasting for a long time.
1517
+ Need to keep on working at that!
1518
+ --------------------------------------------------------------------------------
1519
+ (209) → Add: average_aminoacid_weight?
1520
+ === LV-Nummer 300214 UE Übung III B Sequenzanalysen in der Molekularbiologie
1521
+ → Pubmed
1522
+ Finding sequences
1523
+ → Sequence homology search (Blast, FASTA)
1524
+ → Pairwise sequence alignment
1525
+ DNA analysis, translation
1526
+ Gene finding in genomes
1527
+ ^^^ find genes
1528
+ Plasmid Cloning
1529
+ Primer for PCR
1530
+ Protein analysis
1531
+ PROSITE
1532
+ EMBOSS
1533
+ ^^^ das nochmals in Ruhe durchgehen.
1534
+ in bioroebe viel C++ schreiben, möglichst optimieren,
1535
+ und dann irgendwann später ruby bindings dazu liefern
1536
+ add the ability to compare two FASTA files
1537
+ can probably do so via the two sequences.
1538
+ integrate "grundlagen der bioinfo" slides
1539
+ gleichzeitig während wir dafür lernen.
1540
+ ntSeq.each_entry do |e|
1541
+ # pep = Bio::Sequence.new(e.naseq.reverse_complement!.translate)
1542
+ pep = Bio::Sequence.new(e.naseq.translate)
1543
+ if strand == 0
1544
+ pep = Bio::Sequence.new(e.naseq.reverse_complement!.translate)
1545
+ end
1546
+ puts pep.output_fasta(e.definition,20)
1547
+ end
1548
+ ^^^ we must show this in n characters per line:
1549
+ see https://raw.githubusercontent.com/zorino/bioruby-scripts/master/transeq.rb
1550
+ We must be able to align not only nucleotides but also aminoacids.
1551
+ But where is the alignment comparer? perhaps hamming distance?
1552
+ hmm we have to see.
1553
+ --------------------------------------------------------------------------------
1554
+ (210) → /Programs/Ruby/2.3.1/lib/ruby/site_ruby/2.3.0/bioroebe/bioshell/menu.rb:311:in `menu': undefined method `upcase' for ["EcoRI"]:Array (NoMethodError)
1555
+ from /Programs/Ruby/2.3.1/lib/ruby/site_ruby/2.3.0/bioroebe/bioshell/user_input.rb:31:in `block in enter_main_loop'
1556
+ from /Programs/Ruby/2.3.1/lib/ruby/site_ruby/2.3.0/bioroebe/bioshell/user_input.rb:12:in `loop'
1557
+ from /Programs/Ruby/2.3.1/lib/ruby/site_ruby/2.3.0/bioroebe/bioshell/user_input.rb:12:in `enter_main_loop'
1558
+ from /Programs/Ruby/2.3.1/lib/ruby/site_ruby/2.3.0/bioroebe/bioshell/initialize.rb:42:in `initialize'
1559
+ from /Programs/Ruby/2.3.1/lib/ruby/site_ruby/2.3.0/bioroebe/bioshell/bioshell.rb:52:in `new'
1560
+ from /Programs/Ruby/2.3.1/lib/ruby/site_ruby/2.3.0/bioroebe/bioshell/bioshell.rb:52:in `shell'
1561
+ from /System/Executables/bioshell:6:in `<main>'
1562
+ ^^^ also fix this stupid bug.
1563
+ perhaps redo the whole restriction enzyme stuff.
1564
+ Taxonomy components:
1565
+ Wir wollen nur einen query machen statt vieler kleiner
2096
1566
  queries.
2097
-
2098
- Taxonomy component
2099
-
2100
- (2) Send email, also enable disable_email notification
1567
+ Taxonomy component
1568
+ (2) Send email, also enable disable_email notification
2101
1569
  this should be simple if we use the email part
2102
1570
  of cyberweb-project.
2103
1571
  We send an email when everything has finished.
2104
-
2105
- data = 'Last update... job is now finished,
2106
- at this date.'
2107
- SendEmail.new to: Roebe.email?, data
2108
-
2109
- ..........................................................................
2110
-
2111
-
2112
- Document which parts of emboss have already been copied.
2113
- EMBOSS.md
2114
- ..........................................................................
2115
-
2116
-
2117
-
2118
-
2119
-
2120
- - Trametes_versicolor_FP-101664_SS1_pyranose_2-oxidase_partial_mRNA_XM_008046051.1.fasta
2121
-
2122
-
2123
-
2124
-
2125
- Bioroebe::Shell: Now loading from `Trametes_versicolor_FP-101664_SS1_pyranose_2-oxidase_partial_mRNA_XM_008046051.1.fasta`.
2126
- Bioroebe::ParseFasta: Will read from the file `Trametes_versicolor_FP-101664_SS1_pyranose_2-oxidase_partial_mRNA_XM_008046051.1.fasta`.
2127
- This sequence is assumed to be DNA or RNA.
2128
- The GC content of "XM_008046051.1 Trametes versicolor FP-101664 SS1 pyranose 2-oxidase partial mRNA" is:
2129
- 60.41667 %
2130
- This sequence has 1872 nucleotides.
2131
- We have identified a total of 1 entry in this fasta dataset.
2132
- Setting DNA sequence to (1872 nucleotides):
2133
- 5'- ATGTCTACCAGCTCGAGCGACCCGTTCTTCAACTTCACGAAGTCGAGCTTCAGGAGCGCGGCGGCGCAGAAGGCCTCGGCGACTTCTCTGCCGCCGCTGCCTGGTCCCGACAAGAAAGTCCCTGGAATGGACATCAAGTACGACGTTGTCATAGTAGGCTCCGGACCGATTGGATGCACCTATGCCCGTGAGCTCGTCGAAGCCGGTTACAAGGTCGCGATGTTCGACATCGGAGAGATCGACTCCGGCCTGAAGATCGGTGCCCACAAGAAGAACACTGTCGAATACCAGAAGAACATTGACAAGTTTGTGAACGTCATTCAGGGTCAATTGATGTCTGTTTCCGTTCCCGTCAATACCCTCGTGATCGATACGCTCAGCCCGACGTCTTGGCAAGCTTCATCGTTCTTCGTCCGTAACGGCTCGAACCCAGAGCAGGACCCGCTTCGTAACCTCAGTGGTCAGGCGGTCACACGCGTCGTCGGGGGAATGTCCACGCATTGGACTTGCGCGACACCCCGCTTTGACCGCGAGCAGCGCCCACTGCTCGTGAAGGATGACACGGACGCCGACGACGCCGAGTGGGACCGGCTGTACACCAAGGCCGAGTCGTACTTCAAGACCGGGACGGACCAGTTCAAGGAGTCGATCCGCCACAACCTCGTGCTCAACAAGCTCGCGGAGGAATACAAAGGTCAGCGCGACTTCCAGCAGATCCCGCTGGCGGCAACGCGCCGCAGTCCGACCTTCGTCGAGTGGAGCTCGGCGAACACCGTGTTCGACCTCCAGAACAGGCCGAACACGGACGCGCCGAATGAGCGCTTCAACCTCTTCCCCGCGGTCGCATGTGAGCGCGTCGTGCGCAACACGTCGAACTCCGAGATCGAGAGTCTGCACATCCACGACCTCATCTCAGGCGACCGCTTCGAAATCAAGGCAGACGTGTTTGTTCTCACAGCCGGGGCGGTCCACAACGCGCAGCTTCTCGTGAACTCTGGCTTTGGACAGCTGGGCCGGCCGGACCCCGCGAACCCGCCGCGGTTGCTGCCTTCCCTGGGAAGCTACATCACCGAGCAGTCGCTCGTCTTCTGCCAGACCGTGATGAGCACCGAGCTCATCGACAGCGTCAAGTCCGACATGATCATCAGGGGCAACCCTGGTGATCCGGGGTATAGCGTCACGTACACGCCGGGCGCGTCGACCAACAAGCACCCGGACTGGTGGAACGAGAAGGTGAAGAACCACATGATGCAGCACCAGGAGGACCCGCTCCCGATCCCGTTCGAGGACCCCGAGCCCCAGGTCACCACCCTGTTCCAGCCGTCGCACCCGTGGCACACCCAGATCCACCGCGACGCTTTCAGTTACGGTGCGGTGCAGCAAAGCATCGACTCGCGTCTCATCGTCGACTGGCGCTTTTTCGGAAGGACGGAGCCCAAGGAGGAAAACAAGCTCTGGTTCTCGGACAAGATCACCGACACGTACAACATGCCGCAGCCGACGTTCGACTTCCGCTTCCCGGCAGGCCGCACGAGCAAGGAGGCGGAGGACATGATGACCGACATGTGCGTCATGTCGGCGAAGATTGGTGGCTTCCTGCCCGGCTCTCTCCCGCAATTCATGGAGCCCGGTCTTGTCCTTCACCTCGGTGGTACGCACCGCATGGGCTTCGATGAGCAGGAGGACAAGTGCTGCGTCAACACGGACTCCCGCGTGTTCGGCTTTAAGAACCTTTTCCTCGGCGGCTGCGGCAACATTCCCACCGCGTACGGCGCGAACCCGACGCTCACCGCAATGTCGCTCGCGATCAAGAGTTGCGAGTACATCAAGAACAACTTCACACCGAGCCCTTTCACAGATCAGGCTCAGTGA - 3'
2134
- BIO SHELL> gc?
2135
- Traceback (most recent call last):
2136
- 12: from /System/Index/bin/bioshell:27:in `<main>'
2137
- 11: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/shell.rb:109:in `shell'
2138
- 10: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/shell.rb:109:in `new'
2139
- 9: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/initialize.rb:152:in `initialize'
2140
- 8: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/user_input.rb:18:in `enter_main_loop'
2141
- 7: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/user_input.rb:18:in `loop'
2142
- 6: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/user_input.rb:41:in `block in enter_main_loop'
2143
- 5: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/menu.rb:997:in `menu'
2144
- 4: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/shell.rb:2605:in `calculcate_gc_content'
2145
- 3: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/shell.rb:2605:in `new'
2146
- 2: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/calculate/calculate_gc_content.rb:41:in `initialize'
2147
- 1: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/calculate/calculate_gc_content.rb:71:in `set_data'
2148
- /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/calculate/calculate_gc_content.rb:71:in `exist?': can't convert Bioroebe::Sequence to String (Bioroebe::Sequence#to_str gives Bioroebe::Sequence) (TypeError)
2149
-
2150
-
2151
-
2152
- - require 'bio'
2153
-
2154
- keywords = ARGV.join(' ')
2155
- options = {
2156
- 'retmax' => 1
2157
- }
2158
- entries = Bio::PubMed.esearch(keywords, options)
2159
- Bio::PubMed.efetch(entries).each do |entry|
2160
- medline = Bio::MEDLINE.new(entry)
2161
- reference = medline.reference
2162
- puts reference.bibtex
2163
- end
2164
- ^^^ enable BioPubMed access, similar to bioruby,
2165
- then docment it as well.
2166
-
2167
- - Learn from:
2168
-
2169
- http://www.snapgene.com/products/snapgene_viewer/
2170
-
2171
- -------------------------------------------------------------------------------
2172
- (1) → Wir sollten GFP tagging unterstützen, also wie das
1572
+ data = 'Last update... job is now finished,
1573
+ at this date.'
1574
+ SendEmail.new to: Roebe.email?, data
1575
+ --------------------------------------------------------------------------------
1576
+ (211) → Document which parts of emboss have already been copied.
1577
+ → EMBOSS.md
1578
+ --------------------------------------------------------------------------------
1579
+ (212) → Trametes_versicolor_FP-101664_SS1_pyranose_2-oxidase_partial_mRNA_XM_008046051.1.fasta
1580
+ Bioroebe::Shell: Now loading from `Trametes_versicolor_FP-101664_SS1_pyranose_2-oxidase_partial_mRNA_XM_008046051.1.fasta`.
1581
+ Bioroebe::ParseFasta: Will read from the file `Trametes_versicolor_FP-101664_SS1_pyranose_2-oxidase_partial_mRNA_XM_008046051.1.fasta`.
1582
+ This sequence is assumed to be DNA or RNA.
1583
+ The GC content of "XM_008046051.1 Trametes versicolor FP-101664 SS1 pyranose 2-oxidase partial mRNA" is:
1584
+ 60.41667 %
1585
+ This sequence has 1872 nucleotides.
1586
+ We have identified a total of 1 entry in this fasta dataset.
1587
+ Setting DNA sequence to (1872 nucleotides):
1588
+ 5'- ATGTCTACCAGCTCGAGCGACCCGTTCTTCAACTTCACGAAGTCGAGCTTCAGGAGCGCGGCGGCGCAGAAGGCCTCGGCGACTTCTCTGCCGCCGCTGCCTGGTCCCGACAAGAAAGTCCCTGGAATGGACATCAAGTACGACGTTGTCATAGTAGGCTCCGGACCGATTGGATGCACCTATGCCCGTGAGCTCGTCGAAGCCGGTTACAAGGTCGCGATGTTCGACATCGGAGAGATCGACTCCGGCCTGAAGATCGGTGCCCACAAGAAGAACACTGTCGAATACCAGAAGAACATTGACAAGTTTGTGAACGTCATTCAGGGTCAATTGATGTCTGTTTCCGTTCCCGTCAATACCCTCGTGATCGATACGCTCAGCCCGACGTCTTGGCAAGCTTCATCGTTCTTCGTCCGTAACGGCTCGAACCCAGAGCAGGACCCGCTTCGTAACCTCAGTGGTCAGGCGGTCACACGCGTCGTCGGGGGAATGTCCACGCATTGGACTTGCGCGACACCCCGCTTTGACCGCGAGCAGCGCCCACTGCTCGTGAAGGATGACACGGACGCCGACGACGCCGAGTGGGACCGGCTGTACACCAAGGCCGAGTCGTACTTCAAGACCGGGACGGACCAGTTCAAGGAGTCGATCCGCCACAACCTCGTGCTCAACAAGCTCGCGGAGGAATACAAAGGTCAGCGCGACTTCCAGCAGATCCCGCTGGCGGCAACGCGCCGCAGTCCGACCTTCGTCGAGTGGAGCTCGGCGAACACCGTGTTCGACCTCCAGAACAGGCCGAACACGGACGCGCCGAATGAGCGCTTCAACCTCTTCCCCGCGGTCGCATGTGAGCGCGTCGTGCGCAACACGTCGAACTCCGAGATCGAGAGTCTGCACATCCACGACCTCATCTCAGGCGACCGCTTCGAAATCAAGGCAGACGTGTTTGTTCTCACAGCCGGGGCGGTCCACAACGCGCAGCTTCTCGTGAACTCTGGCTTTGGACAGCTGGGCCGGCCGGACCCCGCGAACCCGCCGCGGTTGCTGCCTTCCCTGGGAAGCTACATCACCGAGCAGTCGCTCGTCTTCTGCCAGACCGTGATGAGCACCGAGCTCATCGACAGCGTCAAGTCCGACATGATCATCAGGGGCAACCCTGGTGATCCGGGGTATAGCGTCACGTACACGCCGGGCGCGTCGACCAACAAGCACCCGGACTGGTGGAACGAGAAGGTGAAGAACCACATGATGCAGCACCAGGAGGACCCGCTCCCGATCCCGTTCGAGGACCCCGAGCCCCAGGTCACCACCCTGTTCCAGCCGTCGCACCCGTGGCACACCCAGATCCACCGCGACGCTTTCAGTTACGGTGCGGTGCAGCAAAGCATCGACTCGCGTCTCATCGTCGACTGGCGCTTTTTCGGAAGGACGGAGCCCAAGGAGGAAAACAAGCTCTGGTTCTCGGACAAGATCACCGACACGTACAACATGCCGCAGCCGACGTTCGACTTCCGCTTCCCGGCAGGCCGCACGAGCAAGGAGGCGGAGGACATGATGACCGACATGTGCGTCATGTCGGCGAAGATTGGTGGCTTCCTGCCCGGCTCTCTCCCGCAATTCATGGAGCCCGGTCTTGTCCTTCACCTCGGTGGTACGCACCGCATGGGCTTCGATGAGCAGGAGGACAAGTGCTGCGTCAACACGGACTCCCGCGTGTTCGGCTTTAAGAACCTTTTCCTCGGCGGCTGCGGCAACATTCCCACCGCGTACGGCGCGAACCCGACGCTCACCGCAATGTCGCTCGCGATCAAGAGTTGCGAGTACATCAAGAACAACTTCACACCGAGCCCTTTCACAGATCAGGCTCAGTGA - 3'
1589
+ BIO SHELL> gc?
1590
+ Traceback (most recent call last):
1591
+ 12: from /System/Index/bin/bioshell:27:in `<main>'
1592
+ 11: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/shell.rb:109:in `shell'
1593
+ 10: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/shell.rb:109:in `new'
1594
+ 9: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/initialize.rb:152:in `initialize'
1595
+ 8: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/user_input.rb:18:in `enter_main_loop'
1596
+ 7: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/user_input.rb:18:in `loop'
1597
+ 6: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/user_input.rb:41:in `block in enter_main_loop'
1598
+ 5: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/menu.rb:997:in `menu'
1599
+ 4: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/shell.rb:2605:in `calculcate_gc_content'
1600
+ 3: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/shell/shell.rb:2605:in `new'
1601
+ 2: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/calculate/calculate_gc_content.rb:41:in `initialize'
1602
+ 1: from /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/calculate/calculate_gc_content.rb:71:in `set_data'
1603
+ /Programs/Ruby/2.6.4/lib/ruby/site_ruby/2.6.0/bioroebe/calculate/calculate_gc_content.rb:71:in `exist?': can't convert Bioroebe::Sequence to String (Bioroebe::Sequence#to_str gives Bioroebe::Sequence) (TypeError)
1604
+ - require 'bio'
1605
+ keywords = ARGV.join(' ')
1606
+ options = {
1607
+ 'retmax' => 1
1608
+ }
1609
+ entries = Bio::PubMed.esearch(keywords, options)
1610
+ Bio::PubMed.efetch(entries).each do |entry|
1611
+ medline = Bio::MEDLINE.new(entry)
1612
+ reference = medline.reference
1613
+ puts reference.bibtex
1614
+ end
1615
+ ^^^ enable BioPubMed access, similar to bioruby,
1616
+ then docment it as well.
1617
+ - Learn from:
1618
+ http://www.snapgene.com/products/snapgene_viewer/
1619
+ --------------------------------------------------------------------------------
1620
+ (213) Wir sollten GFP tagging unterstützen, also wie das
2173
1621
  Protein-Konstrukt aussehen soll und so weiter.
2174
1622
  Das geht teilweise...
2175
1623
  GFP? zeigt die Sequenz an.
2176
- assign :GFP
2177
- fügt die sequence asl main dna sequenz ein.
2178
- Was fehlt? Hmmmm... eventuell noch mehr an
2179
- dokumentation.
2180
- -------------------------------------------------------------------------------
2181
-
2182
- - in bioroebe, create subsequences for siRNA, then scan for
2183
- submatcher + report where these are. Should be fast too.
2184
- -------------------------------------------------------------------------------
2185
- - Reverse complement now works quite well, also via the sinatra
2186
- interface. We still should have a way to show 5' and
2187
- 3', both on the commandline, and via sinatra.
2188
- Perhaps via --fancy commandline flag or so.
2189
- -------------------------------------------------------------------------------
2190
- - Cn3D files?
2191
- ^^^ add support for these; research what they are, too.
2192
- -------------------------------------------------------------------------------
2193
- - Consider adding graphviz, perhaps to the taxonomy project
2194
- where we make graphs towards different nodes or so...
2195
- -------------------------------------------------------------------------------
2196
- - in parse fasta
2197
- @colourize_sequence = false
2198
- ^^^ change this lateron...
2199
- perhaps create a toplevel method
2200
- this method now exists, but we still have to make
2201
- the check better whether it is a protein or a DNA/RNA
2202
- add a toplevel method for this.
2203
- -------------------------------------------------------------------------------
2204
- - clone the BLast ident matcher functionality for aminacids into
2205
- Bioroebe.
2206
-
2207
- - fasta_download AAC76198.2
2208
-
2209
- ^^^ enable the above in the bioshell, and perhaps also outside
1624
+ assign :GFP
1625
+ fügt die sequence asl main dna sequenz ein.
1626
+ Was fehlt? Hmmmm... eventuell noch mehr an
1627
+ dokumentation.
1628
+ --------------------------------------------------------------------------------
1629
+ (214) → in bioroebe, create subsequences for siRNA, then scan for
1630
+ submatcher + report where these are. Should be fast too.
1631
+ --------------------------------------------------------------------------------
1632
+ (215) → Reverse complement now works quite well, also via the sinatra
1633
+ interface. We still should have a way to show 5' and
1634
+ 3', both on the commandline, and via sinatra.
1635
+ Perhaps via --fancy commandline flag or so.
1636
+ --------------------------------------------------------------------------------
1637
+ (216) → Cn3D files?
1638
+ ^^^ add support for these; research what they are, too.
1639
+ --------------------------------------------------------------------------------
1640
+ (217) → Consider adding graphviz, perhaps to the taxonomy project
1641
+ where we make graphs towards different nodes or so...
1642
+ --------------------------------------------------------------------------------
1643
+ (218) → in parse fasta
1644
+ @colourize_sequence = false
1645
+ ^^^ change this lateron...
1646
+ perhaps create a toplevel method
1647
+ this method now exists, but we still have to make
1648
+ the check better whether it is a protein or a DNA/RNA
1649
+ add a toplevel method for this.
1650
+ --------------------------------------------------------------------------------
1651
+ (219) → clone the BLast ident matcher functionality for aminacids into
1652
+ Bioroebe.
1653
+ - fasta_download AAC76198.2
1654
+ ^^^ enable the above in the bioshell, and perhaps also outside
2210
1655
  of the bioshell.
2211
-
2212
- http://www.ncbi.nlm.nih.gov/protein/145693187?report=fasta
2213
-
2214
- ^^^ shall use something such as the above
2215
-
2216
-
2217
-
2218
- -------------------------------------------------------------------------------
2219
- - Be able to mark exon/intron boundaries.
2220
-
2221
- - Add "taxid?" to tell us the name of the organism. This works now.
2222
-
2223
- ^^^ should also work with a local database. ← we integrate this
1656
+ http://www.ncbi.nlm.nih.gov/protein/145693187?report=fasta
1657
+ ^^^ shall use something such as the above
1658
+ --------------------------------------------------------------------------------
1659
+ (220) Be able to mark exon/intron boundaries.
1660
+ - Add "taxid?" to tell us the name of the organism. This works now.
1661
+ ^^^ should also work with a local database. ← we integrate this
2224
1662
  at a later point.
2225
-
2226
- - We have identified a total of 24 entries in this fasta dataset.
2227
-
2228
- There is a total of 4038 letters/nucleotides stored in total.
2229
- Setting DNA sequence to (1582 nucleotides):
2230
- ^^^^ hmm should enable:
2231
-
2232
- @seq1
2233
- @seq2
2234
-
2235
- and so forth
2236
-
2237
- - SUMOylation
2238
- "small ubiquitin modifier"
2239
- chemistry
2240
- SUMO proteins are small
2241
- 100 aa / 12 kD.
2242
-
2243
- Similar structural fold as ubiqitin
2244
-
2245
- Most SUMO-modified proteins contain the tetrapeptide consensus motif.
2246
- phi is a hydrophobic residue
2247
- kappa is the lysine conjugated to SUMO
2248
- x is any aa
2249
- D or E is an acidic residue
2250
- SomethingHydrophobic-K-x-D/E
2251
-
2252
- Prediction programmes, e.g. SUMOplot http://www.abgent.com/sumoplot
2253
-
2254
- >MYB44 LENGTH=305
2255
- MADRIKGPWSPEEDEQLRRLVVKYGPRNWTVISKSIPGRSGKSCRLRWCNQLSPQVEHRPFSAEEDETIARAHAQFGNKWATI
2256
- ARLLNGRTDNAVKNHWNSTLKRKCGGYDHRGYDGSEDHRPVKRSVSAGSPPVVTGLYMSPGSPTGSDVSDSSTIPILPSVELF
2257
- KPVPRPGAVVLPLPIETSSSSDDPPTSLSLSLPGADVSEESNRSHESTNINNTTSSRHNHNNTVSFMPFSGGFRGAIEEMGKS
2258
- FPGNGGEFMAVVQEMIKAEVRSYMTEMQRNNGGGFVGGFIDNGMIPMSQIGVGRIE
2259
-
2260
- ^^^
2261
- study sumoplot ...
2262
- -------------------------------------------------------------------------------
2263
- - http://a-little-book-of-r-for-bioinformatics.readthedocs.io/en/latest/src/chapter7.html
2264
- -------------------------------------------------------------------------------
2265
- - http://biopython.org/DIST/docs/tutorial/Tutorial.html#htoc22
2266
- ^^^ continue here; "You can also specify the table using the
2267
- NCBI table number which is shorter, and often included in
2268
- the feature annotation of GenBank files:"
2269
-
2270
- ^^^ work through this and see if it is good.
2271
-
2272
- -------------------------------------------------------------------------------
2273
-
2274
- - Clone ALL of biophp, if it us useful.
2275
-
2276
- Then state so too, then get rid of this entry here.
2277
-
2278
- But remember, we must also be able to do so via a webinterface!
2279
-
2280
- Oligos now work. Hmm.
2281
-
2282
- http://www.biophp.org/
2283
-
2284
- Let's use a table for now to show which variants we have enabled:
2285
-
2286
- DNA to protein | Bioroebe.to_aa
2287
- Protein to DNA | [IMPLEMENTED FULLY] → http://www.biophp.org/minitools/protein_to_dna/
2288
- | GTK GUI bindings now exist in a simple manner.
2289
- Restriction digest of DNA | http://www.biophp.org/minitools/restriction_digest/demo.php
2290
- Find Palindromic Sequences | [IMPLEMENTED FULLY]
2291
- Sequence manipulation and data |
2292
- Melting Temperature (Tm) Calculator |
2293
- PCR Amplification |
2294
- Microsatellite Repeats Finder |
2295
- Alignment of DNA/Protein sequences | http://www.biophp.org/minitools/seq_alignment/demo.php
2296
-
2297
- Microarray analysis: adaptive quantification |
2298
- Protein sequence information | http://www.biophp.org/minitools/protein_properties/demo.php
2299
- Reduced alphabets for proteins | http://www.biophp.org/minitools/reduce_protein_alphabet/ (started; implemented the first one for now)
2300
- Chaos Game Representation |
2301
- GC-, AT-, KETO- and oligo-skews generator |
2302
- Oligonucleotide Frequency | [IMPLEMENTED PARTIALLY]
2303
- ^^^ we need a way to obtain a hash with the frequencies;
1663
+ - We have identified a total of 24 entries in this fasta dataset.
1664
+ There is a total of 4038 letters/nucleotides stored in total.
1665
+ Setting DNA sequence to (1582 nucleotides):
1666
+ ^^^^ hmm should enable:
1667
+ @seq1
1668
+ @seq2
1669
+ and so forth
1670
+ - SUMOylation
1671
+ "small ubiquitin modifier"
1672
+ chemistry
1673
+ SUMO proteins are small
1674
+ 100 aa / 12 kD.
1675
+ Similar structural fold as ubiqitin
1676
+ Most SUMO-modified proteins contain the tetrapeptide consensus motif.
1677
+ phi is a hydrophobic residue
1678
+ kappa is the lysine conjugated to SUMO
1679
+ x is any aa
1680
+ D or E is an acidic residue
1681
+ SomethingHydrophobic-K-x-D/E
1682
+ Prediction programmes, e.g. SUMOplot http://www.abgent.com/sumoplot
1683
+ >MYB44 LENGTH=305
1684
+ MADRIKGPWSPEEDEQLRRLVVKYGPRNWTVISKSIPGRSGKSCRLRWCNQLSPQVEHRPFSAEEDETIARAHAQFGNKWATI
1685
+ ARLLNGRTDNAVKNHWNSTLKRKCGGYDHRGYDGSEDHRPVKRSVSAGSPPVVTGLYMSPGSPTGSDVSDSSTIPILPSVELF
1686
+ KPVPRPGAVVLPLPIETSSSSDDPPTSLSLSLPGADVSEESNRSHESTNINNTTSSRHNHNNTVSFMPFSGGFRGAIEEMGKS
1687
+ FPGNGGEFMAVVQEMIKAEVRSYMTEMQRNNGGGFVGGFIDNGMIPMSQIGVGRIE
1688
+ ^^^
1689
+ study sumoplot ...
1690
+ --------------------------------------------------------------------------------
1691
+ (221) → http://a-little-book-of-r-for-bioinformatics.readthedocs.io/en/latest/src/chapter7.html
1692
+ --------------------------------------------------------------------------------
1693
+ (222) → http://biopython.org/DIST/docs/tutorial/Tutorial.html#htoc22
1694
+ ^^^ continue here; "You can also specify the table using the
1695
+ NCBI table number which is shorter, and often included in
1696
+ the feature annotation of GenBank files:"
1697
+ ^^^ work through this and see if it is good.
1698
+ --------------------------------------------------------------------------------
1699
+ (223) Clone ALL of biophp, if it us useful.
1700
+ Then state so too, then get rid of this entry here.
1701
+ But remember, we must also be able to do so via a webinterface!
1702
+ Oligos now work. Hmm.
1703
+ http://www.biophp.org/
1704
+ Let's use a table for now to show which variants we have enabled:
1705
+ DNA to protein | Bioroebe.to_aa
1706
+ Protein to DNA | [IMPLEMENTED FULLY] → http://www.biophp.org/minitools/protein_to_dna/
1707
+ | GTK GUI bindings now exist in a simple manner.
1708
+ Restriction digest of DNA | http://www.biophp.org/minitools/restriction_digest/demo.php
1709
+ Find Palindromic Sequences | [IMPLEMENTED FULLY]
1710
+ Sequence manipulation and data |
1711
+ Melting Temperature (Tm) Calculator |
1712
+ PCR Amplification |
1713
+ Microsatellite Repeats Finder |
1714
+ Alignment of DNA/Protein sequences | http://www.biophp.org/minitools/seq_alignment/demo.php
1715
+ Microarray analysis: adaptive quantification |
1716
+ Protein sequence information | http://www.biophp.org/minitools/protein_properties/demo.php
1717
+ Reduced alphabets for proteins | http://www.biophp.org/minitools/reduce_protein_alphabet/ (started; implemented the first one for now)
1718
+ Chaos Game Representation |
1719
+ GC-, AT-, KETO- and oligo-skews generator |
1720
+ Oligonucleotide Frequency | [IMPLEMENTED PARTIALLY]
1721
+ ^^^ we need a way to obtain a hash with the frequencies;
2304
1722
  we also need to extend this to 3 or 4 etc... oligos
2305
- Oligonucleotides for distance among sequences |
2306
- Random sequences |
2307
- Useful formulas
2308
-
2309
- rf biophp
2310
-
2311
-
2312
- Palindromic sequences finder
2313
-
2314
- ^^^ enable this next.
2315
-
2316
- We should also put this poart into doc/ subsection
2317
- to keep track of what is missing and what is not.
2318
-
2319
- -------------------------------------------------------------------------------
2320
- (1) → sizeseq
2321
-
2322
- ^^^ clone this functionality and describe it in detail.
1723
+ Oligonucleotides for distance among sequences |
1724
+ Random sequences |
1725
+ Useful formulas
1726
+ rf biophp
1727
+ Palindromic sequences finder
1728
+ ^^^ enable this next.
1729
+ We should also put this poart into doc/ subsection
1730
+ to keep track of what is missing and what is not.
1731
+ --------------------------------------------------------------------------------
1732
+ (224) sizeseq
1733
+ ^^^ clone this functionality and describe it in detail.
2323
1734
  also for the www. Hmmm. Need to add this for the
2324
1735
  www.
2325
-
2326
1736
  http://www.bioinformatics.nl/cgi-bin/emboss/sizeseq
2327
- ^^^ hmm implemented that I think.
2328
- GTTGTTGCAAGATACAATCTGGTGTGTACTAGA
2329
- AGCTAACTCCAGACCGATACAT
2330
- CGGACTCGGCC
2331
- AATACCAGCGTAGGCTGTGAGCTCGCGGCTGACAAAC
2332
- GGAAACGTTTCCTATGTCGGGATTC
2333
-
2334
- Output file outseq
2335
-
2336
- >four
2337
- ATGC
2338
- >two
2339
- GTTGTTGCAAGATACAATCTGGTGTGTACTAGACCGATACATCGGACTCGGCCAATACCA
2340
- GCGTAGGCTGTGAGCTCGCGGCTGACAAACGGAAACGTTTCCTATGTCGGGAT
2341
- >one
2342
- GTTGTTGCAAGATACAATCTGGTGTGTACTAGAAGCTAACTCCAGACCGATACATCGGAC
2343
- TCGGCCAATACCAGCGTAGGCTGTGAGCTCGCGGCTGACAAACGGAAACGTTTCCTATGT
2344
- CGGGATTC
2345
- >three
2346
- GTTGTTGCAAGATACAATCTGGTGTGTACTAGAAGCTAACTCCAGACGTTGTTGCAAGAT
2347
- ACAATCTGGTGTGTACTAGACGATACATCGGGTTGTTGCAAGATACAATCTGGTGTGTAC
2348
- TAGAACTCGGCCAATACCAGCGTAGGCTGTGAGCTCGCGGCTGACAAACGGAAACGTTTC
2349
- CTATGTCGGGATTC
2350
-
2351
- foobar.fasta
2352
- ^^^ demonstrate via foobar.fasta
2353
-
2354
- ALSO ADD A GUI; sizeseq.rb was added in February 2021.
2355
-
2356
- -------------------------------------------------------------------------------
2357
- - In the sinatra-web-interface for Bioroebe:
2358
- continue quiz in rosalind !!!
2359
- also, at to_dna: default to RNA
2360
- And improve the general quality.
2361
- also add the ability to tchange the codon table
2362
- via URL through the sinatra interface
2363
- Example:
2364
- codon_table/1,2,3
2365
- view it, change it, document it
2366
- also add:
2367
- view_codon_table
2368
- ^^^ shall display it on-line
2369
- and give a formatted-view
2370
- output-view numbering
2371
- Something like:
2372
- → formatted_view
2373
- 111^^^^ in ncbi format
2374
- and document all of this.
2375
- ..........................................................................
2376
- -------------------------------------------------------------------------------
2377
- - Add a ruby-GUI stuff, probably the old biology/ subsection
2378
- will be moved into the project.
2379
-
2380
- Also tell how to start or get this GUI stuff to run, then add
2381
- components that can be a part of bioroebe into it.
2382
-
2383
- ^^^ we should push this before asking for a job in the summer
1737
+ ^^^ hmm implemented that I think.
1738
+ GTTGTTGCAAGATACAATCTGGTGTGTACTAGA
1739
+ AGCTAACTCCAGACCGATACAT
1740
+ CGGACTCGGCC
1741
+ AATACCAGCGTAGGCTGTGAGCTCGCGGCTGACAAAC
1742
+ GGAAACGTTTCCTATGTCGGGATTC
1743
+ Output file outseq
1744
+ >four
1745
+ ATGC
1746
+ >two
1747
+ GTTGTTGCAAGATACAATCTGGTGTGTACTAGACCGATACATCGGACTCGGCCAATACCA
1748
+ GCGTAGGCTGTGAGCTCGCGGCTGACAAACGGAAACGTTTCCTATGTCGGGAT
1749
+ >one
1750
+ GTTGTTGCAAGATACAATCTGGTGTGTACTAGAAGCTAACTCCAGACCGATACATCGGAC
1751
+ TCGGCCAATACCAGCGTAGGCTGTGAGCTCGCGGCTGACAAACGGAAACGTTTCCTATGT
1752
+ CGGGATTC
1753
+ >three
1754
+ GTTGTTGCAAGATACAATCTGGTGTGTACTAGAAGCTAACTCCAGACGTTGTTGCAAGAT
1755
+ ACAATCTGGTGTGTACTAGACGATACATCGGGTTGTTGCAAGATACAATCTGGTGTGTAC
1756
+ TAGAACTCGGCCAATACCAGCGTAGGCTGTGAGCTCGCGGCTGACAAACGGAAACGTTTC
1757
+ CTATGTCGGGATTC
1758
+ foobar.fasta
1759
+ ^^^ demonstrate via foobar.fasta
1760
+ ALSO ADD A GUI; sizeseq.rb was added in February 2021.
1761
+ --------------------------------------------------------------------------------
1762
+ (225) In the sinatra-web-interface for Bioroebe:
1763
+ continue quiz in rosalind !!!
1764
+ also, at to_dna: default to RNA
1765
+ And improve the general quality.
1766
+ → also add the ability to tchange the codon table
1767
+ via URL through the sinatra interface
1768
+ Example:
1769
+ codon_table/1,2,3
1770
+ view it, change it, document it
1771
+ also add:
1772
+ view_codon_table
1773
+ ^^^ shall display it on-line
1774
+ and give a formatted-view
1775
+ output-view numbering
1776
+ Something like:
1777
+ → formatted_view
1778
+ 111^^^^ in ncbi format
1779
+ and document all of this.
1780
+ --------------------------------------------------------------------------------
1781
+ (226)
1782
+ --------------------------------------------------------------------------------
1783
+ (227) Add a ruby-GUI stuff, probably the old biology/ subsection
1784
+ will be moved into the project.
1785
+ Also tell how to start or get this GUI stuff to run, then add
1786
+ components that can be a part of bioroebe into it.
1787
+ ^^^ we should push this before asking for a job in the summer
2384
1788
  months.
2385
-
2386
-
2387
- - http://www.biophp.org/minitools/seq_alignment/demo.php
2388
-
2389
- ^^^ implement smith waterman alignment
2390
-
2391
- swalign AAGGGGAGGACGATGCGGATGTTC AGGGAGGACGATGCGG
2392
- ---
2393
- Query: cmdline (16 nt)
2394
- Ref : cmdline (24 nt)
2395
-
2396
- Query: 1 A-GGGAGGACGATGCGG 16
2397
- | |||||||||||||||
2398
- Ref : 2 AGGGGAGGACGATGCGG 18
2399
-
2400
- Score: 31
2401
- Matches: 16 (94.1%)
2402
- Mismatches: 1
2403
- CIGAR: 1M1D15M
2404
-
2405
- ---
2406
-
2407
- ^^^^ also add this commandline tool
2408
-
2409
-
2410
- bin/swalin
2411
-
2412
- should have same output
2413
- - Gene finding in genomes
2414
-
2415
- ^^^ find genes;
1789
+ - http://www.biophp.org/minitools/seq_alignment/demo.php
1790
+ ^^^ implement smith waterman alignment
1791
+ swalign AAGGGGAGGACGATGCGGATGTTC AGGGAGGACGATGCGG
1792
+ --------------------------------------------------------------------------------
1793
+ (228) Query: cmdline (16 nt)
1794
+ Ref : cmdline (24 nt)
1795
+ Query: 1 A-GGGAGGACGATGCGG 16
1796
+ | |||||||||||||||
1797
+ Ref : 2 AGGGGAGGACGATGCGG 18
1798
+ Score: 31
1799
+ Matches: 16 (94.1%)
1800
+ Mismatches: 1
1801
+ CIGAR: 1M1D15M
1802
+ --------------------------------------------------------------------------------
1803
+ (229) → ^^^^ also add this commandline tool
1804
+ bin/swalin
1805
+ should have same output
1806
+ - Gene finding in genomes
1807
+ ^^^ find genes;
2416
1808
  and add: return all genes from the ORFs
2417
-
2418
-
2419
-
2420
- - Enable blast.
2421
-
2422
- - fasta http://www.ncbi.nlm.nih.gov/nuccore/145337669?report=fasta
2423
-
2424
- Input `fasta http://www.ncbi.nlm.nih.gov/nuccore/145337669?report=fasta` not found.
2425
-
2426
- ^^^ this should work.
2427
- fasta_header? http://www.ncbi.nlm.nih.gov/nuccore/145337669?report=fasta
2428
-
2429
- ^^^ and this hmmm.
2430
-
2431
-
2432
- - http://www.biophp.org/minitools/random_seqs/demo.php
2433
-
2434
- ^^^ clone this thingy
2435
-
2436
- composition?
2437
-
2438
- ^^^ also calculate the percentage.... hmm
2439
-
2440
- one part has been cloned finally BUT it has to be
2441
- described... and we may have to show this
2442
- in CGI too hmmmmmmm.
2443
-
2444
- the first one has been done, the second one not quite yet
2445
- and we lack documentation!
2446
-
2447
- Also add this to sinatra yay!
2448
- Well, we have cloned quite a bit so far. Need to finish
2449
- this up eventually.
2450
-
2451
-
2452
- - How do I write Sequences in Fasta format?
2453
-
2454
- FASTA format is a fairly standard bioinformatics output that is convenient and easy to read. BioRuby's Sequence class has a to_fasta method for formatting sequence in FASTA format.
2455
-
2456
- Printing any Bio::Sequence sequence object in FASTA format.
2457
-
2458
- #!/usr/bin/env ruby
2459
-
2460
- require 'bio'
2461
-
2462
- # Generates a sample 100bp sequence.
2463
- seq1 = Bio::Sequence::NA.new("aatgacccgt" * 10)
2464
-
2465
- # Naming this sequence as "testseq" and print in FASTA format
2466
- # (folded by 60 chars per line).
2467
- puts seq1.to_fasta("testseq", 60)
2468
-
2469
- ^^^^ enable this
2470
-
2471
-
2472
-
2473
- -------------------------------------------------------------------------------
2474
- - Identifying amino acid cleavage sites (Sigcleave)
2475
-
2476
- For amino acid sequences we may be interested to know whether
2477
- the amino acid sequence contains a cleavable signal sequence
2478
- for directing the transport of the protein within the cell.
2479
-
2480
- SigCleave is a program (originally part of the EGCG molecular
2481
- biology package) to predict signal sequences, and to identify
2482
- the cleavage site based on the "von Heijne" algorithm.
2483
-
2484
- The threshold setting controls the score reporting. If no
2485
- value for threshold is passed in by the user, the code
2486
- defaults to a reporting value of 3.5.
2487
-
2488
- SigCleave will only return score/position pairs which meet
2489
- the threshold limit.
2490
-
2491
- There are 2 accessor methods for this object.
2492
- signals() will return a perl hash containing the
2493
- sigcleave scores keyed by amino acid position.
2494
- pretty_print() returns a formatted string similar
2495
- to the output of the original sigcleave utility.
2496
-
2497
- The syntax for using Sigcleave is as follows:
2498
-
2499
- # create a Seq object, for example:
2500
- $seqobj = Bio::Seq->new(-seq => "AALLHHHHHHGGGGPPRTTTTTVVVVVVVVVVVVVVV");
2501
-
2502
- use Bio::Tools::Sigcleave;
2503
- $sigcleave_object = new Bio::Tools::Sigcleave
1809
+ - Enable blast.
1810
+ - fasta http://www.ncbi.nlm.nih.gov/nuccore/145337669?report=fasta
1811
+ Input `fasta http://www.ncbi.nlm.nih.gov/nuccore/145337669?report=fasta` not found.
1812
+ ^^^ this should work.
1813
+ fasta_header? http://www.ncbi.nlm.nih.gov/nuccore/145337669?report=fasta
1814
+ ^^^ and this hmmm.
1815
+ - http://www.biophp.org/minitools/random_seqs/demo.php
1816
+ ^^^ clone this thingy
1817
+ composition?
1818
+ ^^^ also calculate the percentage.... hmm
1819
+ one part has been cloned finally BUT it has to be
1820
+ described... and we may have to show this
1821
+ in CGI too hmmmmmmm.
1822
+ the first one has been done, the second one not quite yet
1823
+ and we lack documentation!
1824
+ Also add this to sinatra yay!
1825
+ Well, we have cloned quite a bit so far. Need to finish
1826
+ this up eventually.
1827
+ - How do I write Sequences in Fasta format?
1828
+ FASTA format is a fairly standard bioinformatics output that is convenient and easy to read. BioRuby's Sequence class has a to_fasta method for formatting sequence in FASTA format.
1829
+ Printing any Bio::Sequence sequence object in FASTA format.
1830
+ #!/usr/bin/env ruby
1831
+ require 'bio'
1832
+ # Generates a sample 100bp sequence.
1833
+ seq1 = Bio::Sequence::NA.new("aatgacccgt" * 10)
1834
+ # Naming this sequence as "testseq" and print in FASTA format
1835
+ # (folded by 60 chars per line).
1836
+ puts seq1.to_fasta("testseq", 60)
1837
+ ^^^^ enable this
1838
+ --------------------------------------------------------------------------------
1839
+ (230) Identifying amino acid cleavage sites (Sigcleave)
1840
+ For amino acid sequences we may be interested to know whether
1841
+ the amino acid sequence contains a cleavable signal sequence
1842
+ for directing the transport of the protein within the cell.
1843
+ SigCleave is a program (originally part of the EGCG molecular
1844
+ biology package) to predict signal sequences, and to identify
1845
+ the cleavage site based on the "von Heijne" algorithm.
1846
+ The threshold setting controls the score reporting. If no
1847
+ value for threshold is passed in by the user, the code
1848
+ defaults to a reporting value of 3.5.
1849
+ SigCleave will only return score/position pairs which meet
1850
+ the threshold limit.
1851
+ There are 2 accessor methods for this object.
1852
+ signals() will return a perl hash containing the
1853
+ sigcleave scores keyed by amino acid position.
1854
+ pretty_print() returns a formatted string similar
1855
+ to the output of the original sigcleave utility.
1856
+ The syntax for using Sigcleave is as follows:
1857
+ # create a Seq object, for example:
1858
+ $seqobj = Bio::Seq->new(-seq => "AALLHHHHHHGGGGPPRTTTTTVVVVVVVVVVVVVVV");
1859
+ use Bio::Tools::Sigcleave;
1860
+ $sigcleave_object = new Bio::Tools::Sigcleave
2504
1861
  ( -seq => $seqobj,
2505
- -threshold => 3.5,
2506
- -description => 'test sigcleave protein seq',
1862
+ -threshold => 3.5,
1863
+ -description => 'test sigcleave protein seq',
2507
1864
  );
2508
- %raw_results = $sigcleave_object->signals;
2509
- $formatted_output = $sigcleave_object->pretty_print;
2510
-
2511
- Please see Bio::Tools::Sigcleave for details.
2512
-
2513
- ^^^ add this
2514
- http://doc.bioperl.org/releases/bioperl-current/bioperl-live/Bio/Tools/Sigcleave.html
2515
-
2516
- - enable drawing of images like the following:
2517
-
2518
- http://nar.oxfordjournals.org/content/43/D1/D227/F1.large.jpg
2519
- https://www.researchgate.net/profile/Matt_Oates/publication/268790596/figure/fig1/AS:295477619773440@1447458762108/Summary-of-all-genome-updates-and-additions-at-the-level-of-taxonomic-Class-since-the.png
2520
-
2521
-
2522
- - add reverse showorf
2523
-
2524
- like emboss
2525
- document this
2526
- then upload
2527
- also enable in bioshell
2528
- r1
2529
- r2
2530
- r3
2531
- so wie f1
2532
- f2 f3
2533
- ^^^ da gibt es einen bug. später nochmals probieren.
2534
-
2535
-
2536
- - We will read from NM_001180897.3_Saccharomyces_cerevisiae_S288c_Aga2p_AGA2.fasta
2537
-
2538
- The file NM_001180897.3_Saccharomyces_cerevisiae_S288c_Aga2p_AGA2.fasta has this FASTA header:
2539
-
2540
- >gi|398364826|ref|NM_001180897.3| Saccharomyces cerevisiae S288c Aga2p (AGA2), mRNA
2541
-
2542
- ^^^ this should also (optionally) tell us the organism, via a switch.
2543
- for this we need some way to return the taxonomic ID of an organism
2544
-
2545
- - we have to add expasy...
2546
- functionality to the cmdline too.
2547
- Which one specifically? Let's see...
2548
-
2549
- https://www.expasy.org/
2550
- -------------------------------------------------------------------------------
2551
- - https://biopython.org/wiki/Category%3ACookbook
2552
- ^^^ clone that
2553
- -------------------------------------------------------------------------------
2554
- - include covid genome, and begin to analyse it in bioroebe
2555
- "Das Genom von SARS-CoV-2 sei doppelt so groß wie jenes
2556
- von Influenzaviren, daher scheinen letztere viermal
2557
- so schnell zu mutieren, schrieb Moshiri."
2558
- -------------------------------------------------------------------------------
2559
- - Look at the GUIs that are part of the BioRoebe project.
2560
-
2561
- Polish these part, at the least one widget, then
2562
- make a screenshot, as the first one.
2563
- Then upload the image + new release and docu!
2564
- document also that more images will be added to this
2565
- in the coming weeks and months. Once done, move this to the
2566
- bottom and regularly improve on this part of the bioroebe
2567
- project.
2568
-
2569
- ^^^ also add java gui to it in the long run.
2570
-
2571
- Hmmm. And then, also consider transitioning into gtk3,
2572
- and make mroe screenshots.
2573
- -------------------------------------------------------------------------------
2574
-
2575
- - https://www.ebi.ac.uk/Tools/seqstats/emboss_pepstats/
2576
- http://www.ebi.ac.uk/Tools/services/web/toolresult.ebi?jobId=emboss_pepstats-I20160208-020243-0564-53154194-oy
2577
-
2578
- ^^^^ clone the pepstat functionality
2579
-
2580
- printAA RLAVQYAPLSGCHSTIREDVHNLHFCRARKE*
2581
-
2582
- - Improve on temperature content and how it is calculated
2583
-
2584
- someone googled for it in 2014 so build on it
2585
- -------------------------------------------------------------------------------
2586
- - pfasta /Depot/Temp/bioroebe/NM_000539.3_Homo_sapiens_rhodopsin_RHO.fasta
2587
-
2588
- Will read from the file `/Depot/Temp/bioroebe/NM_000539.3_Homo_sapiens_rhodopsin_RHO.fasta`.
2589
- Bioroebe::ParseFasta: This sequence is assumed to be a protein.
2590
- This sequence has 2768 aminoacids.
2591
- We have identified a total of 1 entries in this fasta dataset.
2592
- Bioroebe::BioShell: We will now assign this data to @_.
2593
- Now assigning aminoacid sequence to:
2594
- AGAGTCATCCAGCTGGAGCCCTGAGTGGCTGAGCTCAGGCCTTCGCAG
2595
- AGAGTCATCCAGCTGGAGCCCTGAGTGGCTGAGCTCAGGCCTTCGCAG
2596
- -------------------------------------------------------------------------------
2597
-
2598
-
2599
- - Formats
2600
-
2601
- BioPerl's SeqIO system understands lot of formats and can interconvert
2602
- all of them. Here is a current listing of formats, as of version 1.6.
2603
-
2604
- ^^^ must implement this too
2605
-
2606
- Name Description File extension
2607
- abi ABI tracefile ab[i1]
2608
- ace Ace database ace
2609
- agave AGAVE XML
2610
- alf ALF tracefile alf
2611
- asciitree write-only, to visualize features
2612
- bsml BSML using bsm,bsml
2613
- bsml_sax BSML, using
2614
- chadoxml CHADO sequence format
2615
- chaos CHAOS sequence format
2616
- chaosxml Chaos XML
2617
- ctf CTF tracefile ctf
2618
- embl EMBL database embl,ebl,emb,dat
2619
- entrezgene Entrez Gene ASN1
2620
- excel Excel
2621
- exp Staden EXP format exp
2622
- fasta FASTA fasta,fast,seq,fa,fsa,nt,aa
2623
- fastq quality score data in FASTA-like format fastq
2624
- flybase_chadoxml variant of Chado XML
2625
- game GAME XML
2626
- gcg GCG gcg
2627
- genbank GenBank gb
2628
- interpro InterProScan XML
2629
- kegg KEGG
2630
- largefasta Large files, fasta format
2631
- lasergene Lasergene format
2632
- locuslink LocusLink
2633
- metafasta
2634
- phd Phred phd,phred
2635
- pir PIR database pir
2636
- pln PLN tracefile pln
2637
- qual Phred
2638
- raw plain text txt
2639
- scf Standard Chromatogram Format scf
2640
- seqxml SeqXML sequence format xml
2641
- strider DNA Strider format
2642
- swiss SwissProt swiss,sp
2643
- tab tab-delimited
2644
- table Table
2645
- tigr TIGR XML
2646
- tigrxml TIGR Coordset XML
2647
- tinyseq NCBI TinySeq XML
2648
- ztr ZTR tracefile ztr
2649
-
2650
- ..........................................................................
2651
- (1) Look at f1 display:
2652
-
2653
-
2654
- 10 20 30 40 50 60 70 80 90 100
2655
- ---------|---------|---------|---------|---------|---------|---------|---------|---------|---------|
2656
- 1 ATGCAGTTACTTCGCTGTTTTTCAATATTTTCTGTTATTGCTTCAGTTTTAGCACAGGAACTGACAACTATATGCGAGCAAATCCCCTCACCAACTTTAG 100
2657
- F1 1 M Q L L R C F S I F S V I A S V L A Q E L T T I C E Q I P S P T L E 34
2658
-
2659
- ^^^ when we do f1
2660
- the aminoacid sequence position is on the next
2661
- line. this is bad.
2662
-
2663
-
2664
- ff1 "ATGCAGTTACTTCGCTGATTTTCTGTTATTGCTTTTTCAATATTTTCTGTTATTGCTTCAGTTTTAGCACAGGAACTGACAACTATATGCGAGCAAATCCCCTCACCAACTTTAG
2665
- ATGCAGTTACTTCGCTTTCTGTTATTGCTTCAGTTTTAGCACAGGAACTGACAACTATATGCGAGCAAATCCCCTCACCAACTTTAG"
2666
-
2667
- this works semi ok...
2668
- we probably have to rewrite the whole thing
2669
-
2670
- BEFORE we add ANY COLOURS.
2671
- OH WELL.
2672
-
2673
- -------------------------------------------------------------------------------
2674
- (100) → Add a primer-design widget
2675
-
2676
- The idea is to be able to manipulate forward and
2677
- reverse primer areas.
2678
-
2679
- AND research how to do this ...
2680
- We now have a ruby-gtk3 widget for this. It's not yet
2681
- perfect but it is a start.
2682
-
2683
-
2684
- https://www.bioinformatics.nl/molbi/SCLResources/sequence_notation.htm
2685
- ^^^ and check what is useful there. perhaps also add
2686
- nicer visual cues to pretty it up a bit.
2687
- -------------------------------------------------------------------------------
2688
- (1) → Compare bioroebe to:
2689
-
2690
- https://www.ncbi.nlm.nih.gov/orffinder
2691
-
2692
- whether both return the same also possibly add a web-gui
2693
- → it must also allow for different tables to be used!
2694
- check this... so that we can search in standard ORF
2695
- but also in different ORFs
2696
- und die länge angeben, zumindest vom längsten ORF start + stop... also so das das ergebnis auch passt
2697
- ...........................................................................
2698
- test reverse complement in bioroebe
2699
- ^^^
2700
- new_WWW/
2701
- ^^^ this should eventually become the new web-related interface.
2702
- Ah well. Perhaps not ... ruby-cgi is soooooo annoying ...
2703
- ...........................................................................
2704
- (154) → the blosum-viewer should be supported in the cgi part
1865
+ %raw_results = $sigcleave_object->signals;
1866
+ $formatted_output = $sigcleave_object->pretty_print;
1867
+ Please see Bio::Tools::Sigcleave for details.
1868
+ ^^^ add this
1869
+ http://doc.bioperl.org/releases/bioperl-current/bioperl-live/Bio/Tools/Sigcleave.html
1870
+ - enable drawing of images like the following:
1871
+ http://nar.oxfordjournals.org/content/43/D1/D227/F1.large.jpg
1872
+ https://www.researchgate.net/profile/Matt_Oates/publication/268790596/figure/fig1/AS:295477619773440@1447458762108/Summary-of-all-genome-updates-and-additions-at-the-level-of-taxonomic-Class-since-the.png
1873
+ - add reverse showorf
1874
+ like emboss
1875
+ document this
1876
+ then upload
1877
+ also enable in bioshell
1878
+ r1
1879
+ r2
1880
+ r3
1881
+ so wie f1
1882
+ f2 f3
1883
+ ^^^ da gibt es einen bug. später nochmals probieren.
1884
+ - we have to add expasy...
1885
+ functionality to the cmdline too.
1886
+ Which one specifically? Let's see...
1887
+ https://www.expasy.org/
1888
+ --------------------------------------------------------------------------------
1889
+ (231) → https://biopython.org/wiki/Category%3ACookbook
1890
+ ^^^ clone that
1891
+ --------------------------------------------------------------------------------
1892
+ (232) → include covid genome, and begin to analyse it in bioroebe
1893
+ "Das Genom von SARS-CoV-2 sei doppelt so groß wie jenes
1894
+ von Influenzaviren, daher scheinen letztere viermal
1895
+ so schnell zu mutieren, schrieb Moshiri."
1896
+ --------------------------------------------------------------------------------
1897
+ (233) Look at the GUIs that are part of the BioRoebe project.
1898
+ Polish these part, at the least one widget, then
1899
+ make a screenshot, as the first one.
1900
+ Then upload the image + new release and docu!
1901
+ document also that more images will be added to this
1902
+ in the coming weeks and months. Once done, move this to the
1903
+ bottom and regularly improve on this part of the bioroebe
1904
+ project.
1905
+ ^^^ also add java gui to it in the long run.
1906
+ Hmmm. And then, also consider transitioning into gtk3,
1907
+ and make mroe screenshots.
1908
+ --------------------------------------------------------------------------------
1909
+ (234) https://www.ebi.ac.uk/Tools/seqstats/emboss_pepstats/
1910
+ http://www.ebi.ac.uk/Tools/services/web/toolresult.ebi?jobId=emboss_pepstats-I20160208-020243-0564-53154194-oy
1911
+ ^^^^ clone the pepstat functionality
1912
+ printAA RLAVQYAPLSGCHSTIREDVHNLHFCRARKE*
1913
+ - Improve on temperature content and how it is calculated
1914
+ someone googled for it in 2014 so build on it
1915
+ --------------------------------------------------------------------------------
1916
+ (235) pfasta /Depot/Temp/bioroebe/NM_000539.3_Homo_sapiens_rhodopsin_RHO.fasta
1917
+ Will read from the file `/Depot/Temp/bioroebe/NM_000539.3_Homo_sapiens_rhodopsin_RHO.fasta`.
1918
+ Bioroebe::ParseFasta: This sequence is assumed to be a protein.
1919
+ This sequence has 2768 aminoacids.
1920
+ We have identified a total of 1 entries in this fasta dataset.
1921
+ Bioroebe::BioShell: We will now assign this data to @_.
1922
+ Now assigning aminoacid sequence to:
1923
+ AGAGTCATCCAGCTGGAGCCCTGAGTGGCTGAGCTCAGGCCTTCGCAG
1924
+ AGAGTCATCCAGCTGGAGCCCTGAGTGGCTGAGCTCAGGCCTTCGCAG
1925
+ --------------------------------------------------------------------------------
1926
+ (236) Formats
1927
+ BioPerl's SeqIO system understands lot of formats and can interconvert
1928
+ all of them. Here is a current listing of formats, as of version 1.6.
1929
+ ^^^ must implement this too
1930
+ Name Description File extension
1931
+ abi ABI tracefile ab[i1]
1932
+ ace Ace database ace
1933
+ agave AGAVE XML
1934
+ alf ALF tracefile alf
1935
+ asciitree write-only, to visualize features
1936
+ bsml BSML using bsm,bsml
1937
+ bsml_sax BSML, using
1938
+ chadoxml CHADO sequence format
1939
+ chaos CHAOS sequence format
1940
+ chaosxml Chaos XML
1941
+ ctf CTF tracefile ctf
1942
+ embl EMBL database embl,ebl,emb,dat
1943
+ entrezgene Entrez Gene ASN1
1944
+ excel Excel
1945
+ exp Staden EXP format exp
1946
+ fasta FASTA fasta,fast,seq,fa,fsa,nt,aa
1947
+ fastq quality score data in FASTA-like format fastq
1948
+ flybase_chadoxml variant of Chado XML
1949
+ game GAME XML
1950
+ gcg GCG gcg
1951
+ genbank GenBank gb
1952
+ interpro InterProScan XML
1953
+ kegg KEGG
1954
+ largefasta Large files, fasta format
1955
+ lasergene Lasergene format
1956
+ locuslink LocusLink
1957
+ metafasta
1958
+ phd Phred phd,phred
1959
+ pir PIR database pir
1960
+ pln PLN tracefile pln
1961
+ qual Phred
1962
+ raw plain text txt
1963
+ scf Standard Chromatogram Format scf
1964
+ seqxml SeqXML sequence format xml
1965
+ strider DNA Strider format
1966
+ swiss SwissProt swiss,sp
1967
+ tab tab-delimited
1968
+ table Table
1969
+ tigr TIGR XML
1970
+ tigrxml TIGR Coordset XML
1971
+ tinyseq NCBI TinySeq XML
1972
+ ztr ZTR tracefile ztr
1973
+ --------------------------------------------------------------------------------
1974
+ (237) → Look at f1 display:
1975
+ 10 20 30 40 50 60 70 80 90 100
1976
+ --------------------------------------------------------------------------------
1977
+ (238) → 1 ATGCAGTTACTTCGCTGTTTTTCAATATTTTCTGTTATTGCTTCAGTTTTAGCACAGGAACTGACAACTATATGCGAGCAAATCCCCTCACCAACTTTAG 100
1978
+ F1 1 M Q L L R C F S I F S V I A S V L A Q E L T T I C E Q I P S P T L E 34
1979
+ ^^^ when we do f1
1980
+ the aminoacid sequence position is on the next
1981
+ line. this is bad.
1982
+ ff1 "ATGCAGTTACTTCGCTGATTTTCTGTTATTGCTTTTTCAATATTTTCTGTTATTGCTTCAGTTTTAGCACAGGAACTGACAACTATATGCGAGCAAATCCCCTCACCAACTTTAG
1983
+ ATGCAGTTACTTCGCTTTCTGTTATTGCTTCAGTTTTAGCACAGGAACTGACAACTATATGCGAGCAAATCCCCTCACCAACTTTAG"
1984
+ this works semi ok...
1985
+ we probably have to rewrite the whole thing
1986
+ BEFORE we add ANY COLOURS.
1987
+ OH WELL.
1988
+ --------------------------------------------------------------------------------
1989
+ (239) → Add a primer-design widget
1990
+ The idea is to be able to manipulate forward and
1991
+ reverse primer areas.
1992
+ AND research how to do this ...
1993
+ We now have a ruby-gtk3 widget for this. It's not yet
1994
+ perfect but it is a start.
1995
+ https://www.bioinformatics.nl/molbi/SCLResources/sequence_notation.htm
1996
+ ^^^ and check what is useful there. perhaps also add
1997
+ nicer visual cues to pretty it up a bit.
1998
+ --------------------------------------------------------------------------------
1999
+ (240) → Compare bioroebe to:
2000
+ https://www.ncbi.nlm.nih.gov/orffinder
2001
+ whether both return the same also possibly add a web-gui
2002
+ it must also allow for different tables to be used!
2003
+ check this... so that we can search in standard ORF
2004
+ but also in different ORFs
2005
+ und die länge angeben, zumindest vom längsten ORF start + stop... also so das das ergebnis auch passt
2006
+ --------------------------------------------------------------------------------
2007
+ (241) → test reverse complement in bioroebe
2008
+ ^^^
2009
+ new_WWW/
2010
+ ^^^ this should eventually become the new web-related interface.
2011
+ Ah well. Perhaps not ... ruby-cgi is soooooo annoying ...
2012
+ --------------------------------------------------------------------------------
2013
+ (242) → the blosum-viewer should be supported in the cgi part
2705
2014
  and sinatra part as well.
2706
2015
  This now works for sinatra. Need to enable this for
2707
2016
  the cgi-part too eventually.
2708
- -------------------------------------------------------------------------------
2709
- (155) → port the sinatra stuff together in bioroebe
2710
- create a dir: web_api
2711
- ^^^ also make params? usable in both sinatra and cgi page
2712
- well ...............
2713
- this is quite tough.
2714
- Hmmmmmm...
2715
- perhaps as a middle-step,
2716
- add tons of HtmlTemplate[]
2717
- and replace the ad-hoc code otherwise...
2718
- ^^^ yeah, finish the HtmlTemplate stuff.
2719
- -------------------------------------------------------------------------------
2720
- (1) https://i.imgur.com/ptcSn12.png
2721
- ^^^ enable such an overview; this shows mass compuation e.g
2722
- peptide mass and such
2723
- -------------------------------------------------------------------------------
2724
- (80) Bioroebe.sanitize_nucleotide_sequence
2725
- ^^^ port this into java. The code has been written for this already,
2726
- but we currently fail to link it.
2727
- -------------------------------------------------------------------------------
2728
- (81) Bioroebe.base_composition
2729
- ^^^^^^^^^ port this into java
2730
- -------------------------------------------------------------------------------
2731
- (82) - work a bit more on tk!!!
2732
- in particular to start it from the bioshell as-is.
2733
- ^^^ this is mostly done for quick
2734
- demonstration purposes
2735
-
2736
- - also add another ruby-tk widget hmm
2737
- a new one...
2738
- ^^^^^ and fix the remaining ones
2739
-
2740
- hamming_distance [PARTIALLY IMPLEMENTED; ~80%]
2741
- protein_to_DNA
2742
- ^^^^ improve both while improving tk_paradise docu as well.
2743
- -------------------------------------------------------------------------------
2744
- (83) Batch-create the .exe files on windows for libui, once
2745
- the first has been added. And then test it too
2746
- AND document it. This should be done with the controller
2747
- eventually. Once this works, we can remove this entry
2748
- here.
2749
- -------------------------------------------------------------------------------
2750
- (84) port more libui stuff in bioroebe. We have two widgets ported so far;
2751
- add more such entries.
2752
- -------------------------------------------------------------------------------
2753
- (85) after libui has been ported, explore how gosu works on windows.
2754
- if possible add things to a gosu-specific UI as well, but
2755
- we may need a common, unified GUI base for that.
2756
- -------------------------------------------------------------------------------
2757
- (86)
2758
-
2759
- add libui bindings AND once done make sure the controller works in
2760
- libui as well. Embed the various things into it.
2761
-
2762
- Tab A set named tabs for placing items in
2763
- ^^^ use this perhaps also in bioroebe hmmm
2764
- yeah.
2765
- -------------------------------------------------------------------------------
2766
- (87) https://github.com/cnjinhao/nana/wiki/User-Works-using-Nana
2767
-
2768
- ^^^ port the "DNA hybrid"
2769
- https://camo.githubusercontent.com/4c27d554ca4d698d288628f21255f917c2c577e35d7e11dd67e21880d56b6b0a/687474703a2f2f6e616e6170726f2e6f72672f696d616765732f73637265656e73686f74732f746864795f7365715f6578706c2e706e67
2770
-
2771
- -------------------------------------------------------------------------------
2772
- (88) Bioroebe::Cell
2773
- ^^^ think about what to do with it. If we don't need it then perhaps
2774
- we should just remove it. Think about this more at 2022, before
2775
- deciding what to do.
2776
- -------------------------------------------------------------------------------
2777
- (89) - Add emboss cgplot functionality.
2778
-
2779
- https://www.bioinformatics.nl/cgi-bin/emboss/cpgplot
2780
- -------------------------------------------------------------------------------
2781
- (90) - integrate calculation of the Instability index (II)
2782
-
2783
- The instability index provides an estimate of the
2784
- stability of your protein in a test tube. Statistical
2785
- analysis of 12 unstable and 32 stable proteins has
2786
- revealed [7] that there are certain dipeptides, the
2787
- occurence of which is significantly different in the
2788
- unstable proteins compared with those in the stable
2789
- ones. The authors of this method have assigned a
2790
- weight value of instability to each of the 400
2791
- different dipeptides (DIWV). Using these weight
2792
- values it is possible to compute an instability
2793
- index (II) which is defined as:
2794
-
2795
-
2796
- i=L-1
2797
-
2798
- II = (10/L) * Sum DIWV(x[i]x[i+1])
2799
-
2800
- i=1
2801
-
2802
-
2803
- where: L is the length of sequence
2804
-
2805
- DIWV(x[i]x[i+1]) is the instability weight value for the dipeptide starting in position i.
2806
-
2807
- A protein whose instability index is smaller than 40
2808
- is predicted as stable, a value above 40 predicts
2809
- that the protein may be unstable.
2810
-
2017
+ --------------------------------------------------------------------------------
2018
+ (243) → port the sinatra stuff together in bioroebe
2019
+ create a dir: web_api
2020
+ ^^^ also make params? usable in both sinatra and cgi page
2021
+ well ...............
2022
+ this is quite tough.
2023
+ Hmmmmmm...
2024
+ perhaps as a middle-step,
2025
+ add tons of HtmlTemplate[]
2026
+ and replace the ad-hoc code otherwise...
2027
+ ^^^ yeah, finish the HtmlTemplate stuff.
2028
+ --------------------------------------------------------------------------------
2029
+ (244) https://i.imgur.com/ptcSn12.png
2030
+ ^^^ enable such an overview; this shows mass compuation e.g
2031
+ peptide mass and such
2032
+ --------------------------------------------------------------------------------
2033
+ (245) Bioroebe.sanitize_nucleotide_sequence
2034
+ ^^^ port this into java. The code has been written for this already,
2035
+ but we currently fail to link it.
2036
+ --------------------------------------------------------------------------------
2037
+ (246) → Batch-create the .exe files on windows for libui, once
2038
+ the first has been added. And then test it too
2039
+ AND document it. This should be done with the controller
2040
+ eventually. Once this works, we can remove this entry
2041
+ here.
2042
+ --------------------------------------------------------------------------------
2043
+ (247) → port more libui stuff in bioroebe. We have two widgets ported so far;
2044
+ add more such entries.
2045
+ --------------------------------------------------------------------------------
2046
+ (248) after libui has been ported, explore how gosu works on windows.
2047
+ if possible add things to a gosu-specific UI as well, but
2048
+ we may need a common, unified GUI base for that.
2049
+ --------------------------------------------------------------------------------
2050
+ (249) → (86)
2051
+ add libui bindings AND once done make sure the controller works in
2052
+ libui as well. Embed the various things into it.
2053
+ Tab A set named tabs for placing items in
2054
+ ^^^ use this perhaps also in bioroebe hmmm
2055
+ yeah.
2056
+ --------------------------------------------------------------------------------
2057
+ (250) → https://github.com/cnjinhao/nana/wiki/User-Works-using-Nana
2058
+ ^^^ port the "DNA hybrid"
2059
+ https://camo.githubusercontent.com/4c27d554ca4d698d288628f21255f917c2c577e35d7e11dd67e21880d56b6b0a/687474703a2f2f6e616e6170726f2e6f72672f696d616765732f73637265656e73686f74732f746864795f7365715f6578706c2e706e67
2060
+ --------------------------------------------------------------------------------
2061
+ (251) → Bioroebe::Cell
2062
+ ^^^ think about what to do with it. If we don't need it then perhaps
2063
+ we should just remove it. Think about this more at 2022, before
2064
+ deciding what to do.
2065
+ --------------------------------------------------------------------------------
2066
+ (252) → Add emboss cgplot functionality.
2067
+ https://www.bioinformatics.nl/cgi-bin/emboss/cpgplot
2068
+ --------------------------------------------------------------------------------
2069
+ (253) integrate calculation of the Instability index (II)
2070
+ The instability index provides an estimate of the
2071
+ stability of your protein in a test tube. Statistical
2072
+ analysis of 12 unstable and 32 stable proteins has
2073
+ revealed [7] that there are certain dipeptides, the
2074
+ occurence of which is significantly different in the
2075
+ unstable proteins compared with those in the stable
2076
+ ones. The authors of this method have assigned a
2077
+ weight value of instability to each of the 400
2078
+ different dipeptides (DIWV). Using these weight
2079
+ values it is possible to compute an instability
2080
+ index (II) which is defined as:
2081
+ i=L-1
2082
+ II = (10/L) * Sum DIWV(x[i]x[i+1])
2083
+ i=1
2084
+ where: L is the length of sequence
2085
+ DIWV(x[i]x[i+1]) is the instability weight value for the dipeptide starting in position i.
2086
+ A protein whose instability index is smaller than 40
2087
+ is predicted as stable, a value above 40 predicts
2088
+ that the protein may be unstable.
2811
2089
  # MEKVQYLTRSAIRRASTIEMPQQARQKLQNLFINFCLILICLLLICIIVMLL
2812
-
2813
- 52
2814
-
2815
- The instability index (II) is computed to be 65.43
2816
- This classifies the protein as unstable.
2817
-
2818
- -------------------------------------------------------------------------------
2819
- (1) → We have now added a method to show all hydrophobic amino acids, via the
2090
+ 52
2091
+ The instability index (II) is computed to be 65.43
2092
+ This classifies the protein as unstable.
2093
+ --------------------------------------------------------------------------------
2094
+ (254) We have now added a method to show all hydrophobic amino acids, via the
2820
2095
  method .hydrophobic_amino_acids?. This works and has been documented
2821
2096
  in May 2022. However had, we also still need a way to PREDICT
2822
2097
  hydrophobic segments in a polypeptide sequence.
2823
- -------------------------------------------------------------------------------
2098
+ --------------------------------------------------------------------------------
2099
+ (255) → <img src="https://i.imgur.com/tkB8MTJ.png" style="margin: 1em">
2100
+ --------------------------------------------------------------------------------
2101
+ (256) → https://www.studocu.com/en-us/document/queens-college-cuny/biochemistry-laboratory/bioinformatics-exercise/13329106
2102
+ ^^^ this enable via a method
2103
+ and add a screenshot
2104
+ we want to colourize an existing string
2105
+ ALSo use javascript for this on the www, otherwise use ruby
2106
+ hmmm. so first the ruby variant + sinatra demo app
2107
+ ^^^ this works now, see the
2108
+ image at:
2109
+ https://i.imgur.com/tkB8MTJ.png
2110
+ However had, we should add a sinatra demo app too,
2111
+ and demonstrate this too and then documen it as-is
2112
+ --------------------------------------------------------------------------------
2113
+ (257) →
2114
+ make sure we have a good fasta-showing widget
2115
+ show how many nucleotides are
2116
+ AND add support to modify this as-is
2117
+ ^^^^
2118
+ The fasta showing widget in ruby-gtk3 is quite nice, but we
2119
+ need to make it more convenient to work with.
2120
+ Perhaps add a context menu that can be customzied?
2121
+ hmm. and perhaps open the sequence in the editor
2122
+ or something ... perhaps also keybindings by default
2123
+ and a help option somewhere to explain all of this.
2124
+ --------------------------------------------------------------------------------
2125
+ (258) → Add a way in bioroebe to store a gene into a yaml file
2126
+ or so, and to also load it up again. Perhaps simplify
2127
+ this automatically. Need some ways to describe that.
2128
+ FastaToYaml
2129
+ ^^^ perhaps?
2130
+ ^^^ describe the why too
2131
+ This class now exists. We have to add more features to it
2132
+ eventually, though.
2133
+ --------------------------------------------------------------------------------
2134
+ (259) → https://pubchem.ncbi.nlm.nih.gov/compound/16131099#section=Top
2135
+ ^^^ this website is quite interesting; try to use components
2136
+ from it.
2137
+ --------------------------------------------------------------------------------
2138
+ (260) → Add some option to show the aminoacid sequence, at the least
2139
+ store it; and optionally show it.
2140
+ possibly always report how many aminoacids are
2141
+ part of that file; and optionally also show
2142
+ the whole sequence.
2143
+ --------------------------------------------------------------------------------
2144
+ (261) → http://insilico.ehu.es/
2145
+ ^^^ check if we have all of this incorporated
2146
+ --------------------------------------------------------------------------------
2147
+ (262) → Integrate these nice GUI parts parts:
2148
+ https://dev.to/kojix2/introduction-to-gr-rb-data-visualization-with-ruby-2c39
2149
+ --------------------------------------------------------------------------------
2150
+ (263) → AND THEN test on windows as well.
2151
+ ^^^^^^^^^^^^^^
2152
+ --------------------------------------------------------------------------------
2153
+ (264) → add mouse chromsoome URL, also in the bioshell
2154
+ and the main README, to be of help for the
2155
+ user. add a mouse subsection.
2156
+ --------------------------------------------------------------------------------
2157
+ (265) → fix the taxonomy stuff...
2158
+ --------------------------------------------------------------------------------
2159
+ (266) → set_dna_sequence alu
2160
+ ^^^ fetch random alu
2161
+ ^^^ alu sequence
2162
+ Ok we started this now adding more details, but we
2163
+ need to become better at searching for this
2164
+ sequence.
2165
+ --------------------------------------------------------------------------------
2166
+ (267) → draw things based on GR
2167
+ --------------------------------------------------------------------------------
2168
+ (268) → https://mycocosm.jgi.doe.gov/help/screenshots/browser_viewer.png
2169
+ ^^^ offer the same functionality
2170
+ --------------------------------------------------------------------------------
2171
+ (269) → https://genome.cshlp.org/content/12/10/1611/F3.expansion.html
2172
+ ^^^ enable this, we must obtain a sequence then store into genbank format
2173
+ so, first fetch; then store as-is.
2174
+ --------------------------------------------------------------------------------
2175
+ (270) → be able to generate nice graphics
2176
+ https://genome.cshlp.org/content/12/10/1611/F1.large.jpg
2177
+ --------------------------------------------------------------------------------
2178
+ (271) → add rmagicks wrappre, perhaps via imageparadise or something
2179
+ the idea is that we can make fancy drawings and generate
2180
+ an image for the end user to see
2181
+ --------------------------------------------------------------------------------
2182
+ (272) → https://bioperl.org/howtos/Beginners_HOWTO.html#item13
2183
+ extend the sequence object and document it
2184
+ also add:
2185
+ class Genome
2186
+ and:
2187
+ def is_circular?
2188
+ @internal_hash[:is_circular]
2189
+ end; alias circular? is_circular? # === circular?
2190
+ def species?
2191
+ @internal_hash[:species] # return the species here
2192
+ end
2193
+ --------------------------------------------------------------------------------
2194
+ (273) → http://lib.ysu.am/open_books/312400.pdf
2195
+ clone:
2196
+ Primer.pl
2197
+ This program was written to support the required informatics for a sequencing
2198
+ lab. The desire was to quickly generate primer pair candidates for use in STS
2199
+ mapping. We use Bioperl modules to fetch the sequences from GenBank.
2200
+ #! /usr/bin/perl
2201
+ #
2202
+ # primers.pl
2203
+ #
2204
+ # Reads a list of
2205
+ % primers.pl AC013798
2206
+ AC013798
2207
+ Left Right Length Penalty
2208
+ CCTCCTGGACAACCTGTGTT TGAAGTCAGGGGACATAGGG 280 0.0823
2209
+ CCTCCTGGACAACCTGTGTT AGGCCAGTAGACTGGGTGTG 298 0.1758
2210
+ CCTCCTGGACAACCTGTGTT GGTGTGAAGTCAGGGGACAT 284 0.1852
2211
+ TTCCCGCATCTCTTAGCAGT AGGCCAGTAGACTGGGTGTG 209 0.1962
2212
+ CTTCCCGCATCTCTTAGCAG GACACTAGTGGCAAGGAGGC 226 0.2362
2213
+ Most of the primers.pl program is extremely simple. The real guts and power
2214
+ of the program lie in the classes and the methods we call. The next section
2215
+ examines the Primer3 module, which is similar to many Bioperl modules
2216
+ --------------------------------------------------------------------------------
2217
+ (274) → Clone all of Emboss. :)
2218
+ → Clone and document the getorf functionality properly.
2219
+ See: http://emboss.sourceforge.net/apps/cvs/emboss/apps/getorf.html
2220
+ http://emboss.sourceforge.net
2221
+ http://emboss.sourceforge.net/apps/cvs/emboss/apps/getorf.html
2222
+ --------------------------------------------------------------------------------
2223
+ (275) → Add useful formulas for bioshell.
2224
+ --------------------------------------------------------------------------------
2225
+ (276) → Polish the GUI sets:
2226
+ https://i.imgur.com/djElIMh.png
2227
+ --------------------------------------------------------------------------------
2228
+ (277) → The taxonomy part should be fully integrated, without it
2229
+ being a standalone part anymore.
2230
+ continue on the taxonomy stuff.
2231
+ ne day this will work again *shake fist*
2232
+ --------------------------------------------------------------------------------
2233
+ (278) → Show the frequency of codons in different tables
2234
+ This works quite ok, but right now the approach is to store
2235
+ this in a .yml file which is not ideal.
2236
+ Thus, we have to add two things:
2237
+ - The ability to store this into a SQL database
2238
+ - The ability to batch-download all of these codons,
2239
+ which first requires that we have a way to obtain all
2240
+ taxonomic ids.
2241
+ Add where this can be found.
2242
+ IMPROVE THIS ALL!!!!!!!
2243
+ --------------------------------------------------------------------------------
2244
+ (279) → improve docu + tests for melting temperature analysis again
2245
+ + usage example + GUI + web-use
2246
+ --------------------------------------------------------------------------------
2247
+ (280) → https://biopython.org/DIST/docs/tutorial/Tutorial.html#sec15
2248
+ ^^^ work through the above, also integrate it + write docs
2249
+ https://raw.githubusercontent.com/biopython/biopython/master/Doc/examples/ls_orchid.fasta
2250
+ --------------------------------------------------------------------------------
2251
+ (281) → work a bit more on tk!!!
2252
+ in particular to start it from the bioshell as-is.
2253
+ ^^^ this is mostly done for quick
2254
+ demonstration purposes
2255
+ - also add another ruby-tk widget hmm
2256
+ a new one...
2257
+ ^^^^^ and fix the remaining ones
2258
+ hamming_distance [PARTIALLY IMPLEMENTED; ~80%]
2259
+ protein_to_DNA
2260
+ ^^^^ improve both while improving tk_paradise docu as well.
2261
+ --------------------------------------------------------------------------------
2262
+ (1) → add 2nd_orf
2263
+ → this shall scan for the 2nd orf
2264
+ → and third ORF as well, then, and document it.
2265
+ --------------------------------------------------------------------------------
2266
+ (299) https://github.com/pjotrp/bigbio
2267
+ ^^^^ include uses cases from that readme