bio-polymarker 1.3.2

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (177) hide show
  1. checksums.yaml +7 -0
  2. data/.travis.yml +24 -0
  3. data/Gemfile +23 -0
  4. data/README.md +205 -0
  5. data/Rakefile +61 -0
  6. data/SECURITY.md +16 -0
  7. data/VERSION +1 -0
  8. data/bin/bfr.rb +128 -0
  9. data/bin/blast_triads.rb +166 -0
  10. data/bin/blast_triads_promoters.rb +192 -0
  11. data/bin/count_variations.rb +36 -0
  12. data/bin/filter_blat_by_target_coverage.rb +69 -0
  13. data/bin/filter_exonerate_by_identity.rb +38 -0
  14. data/bin/find_best_blat_hit.rb +33 -0
  15. data/bin/find_best_exonerate.rb +17 -0
  16. data/bin/get_longest_hsp_blastx_triads.rb +66 -0
  17. data/bin/hexaploid_primers.rb +168 -0
  18. data/bin/homokaryot_primers.rb +183 -0
  19. data/bin/mafft_triads.rb +120 -0
  20. data/bin/mafft_triads_promoters.rb +403 -0
  21. data/bin/map_markers_to_contigs.rb +66 -0
  22. data/bin/marker_to_vcf.rb +241 -0
  23. data/bin/markers_in_region.rb +42 -0
  24. data/bin/mask_triads.rb +169 -0
  25. data/bin/polymarker.rb +410 -0
  26. data/bin/polymarker_capillary.rb +443 -0
  27. data/bin/polymarker_deletions.rb +350 -0
  28. data/bin/snp_position_to_polymarker.rb +101 -0
  29. data/bin/snps_between_bams.rb +107 -0
  30. data/bin/tag_stats.rb +75 -0
  31. data/bin/vcfLineToTable.rb +56 -0
  32. data/bin/vcfToPolyMarker.rb +82 -0
  33. data/bio-polymarker.gemspec +227 -0
  34. data/conf/defaults.rb +1 -0
  35. data/conf/primer3_config/dangle.dh +128 -0
  36. data/conf/primer3_config/dangle.ds +128 -0
  37. data/conf/primer3_config/interpretations/dangle_i.dh +131 -0
  38. data/conf/primer3_config/interpretations/dangle_i.ds +131 -0
  39. data/conf/primer3_config/interpretations/loops_i.dh +34 -0
  40. data/conf/primer3_config/interpretations/loops_i.ds +31 -0
  41. data/conf/primer3_config/interpretations/stack_i.dh +257 -0
  42. data/conf/primer3_config/interpretations/stack_i.ds +256 -0
  43. data/conf/primer3_config/interpretations/stackmm_i_mm.dh +257 -0
  44. data/conf/primer3_config/interpretations/stackmm_i_mm.ds +256 -0
  45. data/conf/primer3_config/interpretations/tetraloop_i.dh +79 -0
  46. data/conf/primer3_config/interpretations/tetraloop_i.ds +81 -0
  47. data/conf/primer3_config/interpretations/triloop_i.dh +21 -0
  48. data/conf/primer3_config/interpretations/triloop_i.ds +18 -0
  49. data/conf/primer3_config/interpretations/tstack2_i.dh +256 -0
  50. data/conf/primer3_config/interpretations/tstack2_i.ds +256 -0
  51. data/conf/primer3_config/interpretations/tstack_i.dh +256 -0
  52. data/conf/primer3_config/interpretations/tstack_i.ds +256 -0
  53. data/conf/primer3_config/interpretations/tstack_tm_inf_i.dh +256 -0
  54. data/conf/primer3_config/interpretations/tstack_tm_inf_i.ds +256 -0
  55. data/conf/primer3_config/loops.dh +30 -0
  56. data/conf/primer3_config/loops.ds +30 -0
  57. data/conf/primer3_config/stack.dh +256 -0
  58. data/conf/primer3_config/stack.ds +256 -0
  59. data/conf/primer3_config/stackmm.dh +256 -0
  60. data/conf/primer3_config/stackmm.ds +256 -0
  61. data/conf/primer3_config/tetraloop.dh +77 -0
  62. data/conf/primer3_config/tetraloop.ds +77 -0
  63. data/conf/primer3_config/triloop.dh +16 -0
  64. data/conf/primer3_config/triloop.ds +16 -0
  65. data/conf/primer3_config/tstack.dh +256 -0
  66. data/conf/primer3_config/tstack2.dh +256 -0
  67. data/conf/primer3_config/tstack2.ds +256 -0
  68. data/conf/primer3_config/tstack_tm_inf.ds +256 -0
  69. data/lib/bio/BFRTools.rb +465 -0
  70. data/lib/bio/BIOExtensions.rb +153 -0
  71. data/lib/bio/PolyploidTools/ChromosomeArm.rb +63 -0
  72. data/lib/bio/PolyploidTools/ExonContainer.rb +245 -0
  73. data/lib/bio/PolyploidTools/Marker.rb +175 -0
  74. data/lib/bio/PolyploidTools/Mask.rb +116 -0
  75. data/lib/bio/PolyploidTools/NoSNPSequence.rb +292 -0
  76. data/lib/bio/PolyploidTools/PrimerRegion.rb +30 -0
  77. data/lib/bio/PolyploidTools/SNP.rb +804 -0
  78. data/lib/bio/PolyploidTools/SNPMutant.rb +86 -0
  79. data/lib/bio/PolyploidTools/SNPSequence.rb +55 -0
  80. data/lib/bio/db/blast.rb +114 -0
  81. data/lib/bio/db/exonerate.rb +333 -0
  82. data/lib/bio/db/primer3.rb +820 -0
  83. data/lib/bio-polymarker.rb +28 -0
  84. data/test/data/7B_amplicon_test.fa +12 -0
  85. data/test/data/7B_amplicon_test.fa.fai +1 -0
  86. data/test/data/7B_amplicon_test_reference.fa +110 -0
  87. data/test/data/7B_amplicon_test_reference.fa.fai +3 -0
  88. data/test/data/7B_marker_test.txt +1 -0
  89. data/test/data/BS00068396_51.fa +2 -0
  90. data/test/data/BS00068396_51_blast.tab +4 -0
  91. data/test/data/BS00068396_51_contigs.aln +1412 -0
  92. data/test/data/BS00068396_51_contigs.dnd +7 -0
  93. data/test/data/BS00068396_51_contigs.fa +8 -0
  94. data/test/data/BS00068396_51_contigs.fa.fai +4 -0
  95. data/test/data/BS00068396_51_contigs.fa.nhr +0 -0
  96. data/test/data/BS00068396_51_contigs.fa.nin +0 -0
  97. data/test/data/BS00068396_51_contigs.fa.nsq +0 -0
  98. data/test/data/BS00068396_51_contigs.nhr +0 -0
  99. data/test/data/BS00068396_51_contigs.nin +0 -0
  100. data/test/data/BS00068396_51_contigs.nsq +0 -0
  101. data/test/data/BS00068396_51_exonerate.tab +6 -0
  102. data/test/data/BS00068396_51_for_polymarker.txt +1 -0
  103. data/test/data/BS00068396_51_genes.txt +14 -0
  104. data/test/data/IWGSC_CSS_1AL_scaff_1455974.fa +112 -0
  105. data/test/data/IWGSC_CSS_1AL_scaff_1455974_aln_contigs.fa +2304 -0
  106. data/test/data/IWGSC_CSS_1AL_scaff_1455974_aln_contigs.fa.fai +11 -0
  107. data/test/data/LIB1716.bam +0 -0
  108. data/test/data/LIB1716.bam.bai +0 -0
  109. data/test/data/LIB1719.bam +0 -0
  110. data/test/data/LIB1719.bam.bai +0 -0
  111. data/test/data/LIB1721.bam +0 -0
  112. data/test/data/LIB1721.bam.bai +0 -0
  113. data/test/data/LIB1722.bam +0 -0
  114. data/test/data/LIB1722.bam.bai +0 -0
  115. data/test/data/PST130_7067.csv +1 -0
  116. data/test/data/PST130_7067.fa +2 -0
  117. data/test/data/PST130_7067.fa.fai +1 -0
  118. data/test/data/PST130_7067.fa.ndb +0 -0
  119. data/test/data/PST130_7067.fa.nhr +0 -0
  120. data/test/data/PST130_7067.fa.nin +0 -0
  121. data/test/data/PST130_7067.fa.not +0 -0
  122. data/test/data/PST130_7067.fa.nsq +0 -0
  123. data/test/data/PST130_7067.fa.ntf +0 -0
  124. data/test/data/PST130_7067.fa.nto +0 -0
  125. data/test/data/PST130_reverse_primer.csv +1 -0
  126. data/test/data/S22380157.fa +16 -0
  127. data/test/data/S22380157.fa.fai +1 -0
  128. data/test/data/S22380157.vcf +67 -0
  129. data/test/data/S58861868/LIB1716.bam +0 -0
  130. data/test/data/S58861868/LIB1716.sam +651 -0
  131. data/test/data/S58861868/LIB1719.bam +0 -0
  132. data/test/data/S58861868/LIB1719.sam +805 -0
  133. data/test/data/S58861868/LIB1721.bam +0 -0
  134. data/test/data/S58861868/LIB1721.sam +1790 -0
  135. data/test/data/S58861868/LIB1722.bam +0 -0
  136. data/test/data/S58861868/LIB1722.sam +1271 -0
  137. data/test/data/S58861868/S58861868.fa +16 -0
  138. data/test/data/S58861868/S58861868.fa.fai +1 -0
  139. data/test/data/S58861868/S58861868.vcf +76 -0
  140. data/test/data/S58861868/header.txt +9 -0
  141. data/test/data/S58861868/merged.bam +0 -0
  142. data/test/data/S58861868/merged_reheader.bam +0 -0
  143. data/test/data/S58861868/merged_reheader.bam.bai +0 -0
  144. data/test/data/Test3Aspecific.csv +2 -0
  145. data/test/data/Test3Aspecific_contigs.fa +6 -0
  146. data/test/data/bfr_out_test.csv +5 -0
  147. data/test/data/chr1A_C1145499T/chr1A_C1145499T.csv +1 -0
  148. data/test/data/chr1A_G540414846C/chr1A_G540414846C.csv +1 -0
  149. data/test/data/chr1A_G540414846C/chr1A_G540414846C.fa +2 -0
  150. data/test/data/chr1A_T517634750C/chr1A_T517634750C.csv +1 -0
  151. data/test/data/chr2D_C112180134A/chr2D_C112180134A.csv +1 -0
  152. data/test/data/chr4D_C14473543T/chr4D_C14473543T.csv +1 -0
  153. data/test/data/chr4D_C14473543T/chr4D_C14473543T.fa +2 -0
  154. data/test/data/headerMergeed.txt +9 -0
  155. data/test/data/headerS2238015 +1 -0
  156. data/test/data/mergedLibs.bam +0 -0
  157. data/test/data/mergedLibsReheader.bam +0 -0
  158. data/test/data/mergedLibsSorted.bam +0 -0
  159. data/test/data/mergedLibsSorted.bam.bai +0 -0
  160. data/test/data/patological_cases5D.csv +1 -0
  161. data/test/data/primer_3_input_header_test +5 -0
  162. data/test/data/short_primer_design_test.csv +10 -0
  163. data/test/data/some_tests/some_tests.csv +201 -0
  164. data/test/data/test_from_mutant.csv +3 -0
  165. data/test/data/test_iselect.csv +196 -0
  166. data/test/data/test_iselect_reference.fa +1868 -0
  167. data/test/data/test_iselect_reference.fa.fai +934 -0
  168. data/test/data/test_primer3_error.csv +4 -0
  169. data/test/data/test_primer3_error_contigs.fa +10 -0
  170. data/test/test_bfr.rb +135 -0
  171. data/test/test_blast.rb +47 -0
  172. data/test/test_exon_container.rb +17 -0
  173. data/test/test_exonearate.rb +48 -0
  174. data/test/test_integration.rb +76 -0
  175. data/test/test_snp_parsing.rb +121 -0
  176. data/test/test_wrong_selection.sh +5 -0
  177. metadata +356 -0
data/conf/primer3_config/loops.dh ADDED
@@ -0,0 +1,30 @@
1
+ 1 inf 0.0 inf
2
+ 2 inf 0.0 inf
3
+ 3 0.0 0.0 0.0
4
+ 4 0.0 0.0 0.0
5
+ 5 0.0 0.0 0.0
6
+ 6 0.0 0.0 0.0
7
+ 7 0.0 0.0 0.0
8
+ 8 0.0 0.0 0.0
9
+ 9 0.0 0.0 0.0
10
+ 10 0.0 0.0 0.0
11
+ 11 0.0 0.0 0.0
12
+ 12 0.0 0.0 0.0
13
+ 13 0.0 0.0 0.0
14
+ 14 0.0 0.0 0.0
15
+ 15 0.0 0.0 0.0
16
+ 16 0.0 0.0 0.0
17
+ 17 0.0 0.0 0.0
18
+ 18 0.0 0.0 0.0
19
+ 19 0.0 0.0 0.0
20
+ 20 0.0 0.0 0.0
21
+ 21 0.0 0.0 0.0
22
+ 22 0.0 0.0 0.0
23
+ 23 0.0 0.0 0.0
24
+ 24 0.0 0.0 0.0
25
+ 25 0.0 0.0 0.0
26
+ 26 0.0 0.0 0.0
27
+ 27 0.0 0.0 0.0
28
+ 28 0.0 0.0 0.0
29
+ 29 0.0 0.0 0.0
30
+ 30 0.0 0.0 0.0
data/conf/primer3_config/loops.ds ADDED
@@ -0,0 +1,30 @@
1
+ 1 -1.0 -12.89 -1.0
2
+ 2 -1.0 -9.35 -1.0
3
+ 3 -10.31 -9.99 -11.28
4
+ 4 -11.6 -10.31 -11.28
5
+ 5 -12.89 -10.64 -10.64
6
+ 6 -14.18 -11.28 -12.89
7
+ 7 -14.83 -11.92 -13.54
8
+ 8 -15.47 -12.57 -13.86
9
+ 9 -15.79 -13.21 -14.5
10
+ 10 -15.79 -13.86 -14.83
11
+ 11 -16.26 -14.32 -15.29
12
+ 12 -16.76 -14.5 -16.12
13
+ 13 -17.15 -14.89 -16.5
14
+ 14 -17.41 -15.47 -16.44
15
+ 15 -17.74 -15.81 -16.77
16
+ 16 -18.05 -16.12 -17.08
17
+ 17 -18.34 -16.41 -17.38
18
+ 18 -18.7 -16.76 -17.73
19
+ 19 -18.96 -17.02 -17.99
20
+ 20 -19.02 -17.08 -18.37
21
+ 21 -19.25 -17.32 -18.61
22
+ 22 -19.48 -17.55 -18.84
23
+ 23 -19.7 -17.76 -19.05
24
+ 24 -19.9 -17.97 -19.26
25
+ 25 -20.31 -18.05 -19.66
26
+ 26 -20.5 -18.24 -19.85
27
+ 27 -20.68 -18.42 -20.04
28
+ 28 -20.86 -18.6 -20.21
29
+ 29 -21.03 -18.77 -20.38
30
+ 30 -21.28 -19.02 -20.31
data/conf/primer3_config/stack.dh ADDED
@@ -0,0 +1,256 @@
1
+ inf
2
+ inf
3
+ inf
4
+ inf
5
+ inf
6
+ inf
7
+ inf
8
+ inf
9
+ inf
10
+ inf
11
+ inf
12
+ inf
13
+ inf
14
+ inf
15
+ inf
16
+ -7900
17
+ inf
18
+ inf
19
+ inf
20
+ inf
21
+ inf
22
+ inf
23
+ inf
24
+ inf
25
+ inf
26
+ inf
27
+ inf
28
+ inf
29
+ inf
30
+ inf
31
+ -8400
32
+ inf
33
+ inf
34
+ inf
35
+ inf
36
+ inf
37
+ inf
38
+ inf
39
+ inf
40
+ inf
41
+ inf
42
+ inf
43
+ inf
44
+ inf
45
+ inf
46
+ -7800
47
+ inf
48
+ inf
49
+ inf
50
+ inf
51
+ inf
52
+ inf
53
+ inf
54
+ inf
55
+ inf
56
+ inf
57
+ inf
58
+ inf
59
+ inf
60
+ inf
61
+ -7200
62
+ inf
63
+ inf
64
+ inf
65
+ inf
66
+ inf
67
+ inf
68
+ inf
69
+ inf
70
+ inf
71
+ inf
72
+ inf
73
+ inf
74
+ inf
75
+ inf
76
+ -8500
77
+ inf
78
+ inf
79
+ inf
80
+ inf
81
+ inf
82
+ inf
83
+ inf
84
+ inf
85
+ inf
86
+ inf
87
+ inf
88
+ inf
89
+ inf
90
+ inf
91
+ -8000
92
+ inf
93
+ inf
94
+ inf
95
+ inf
96
+ inf
97
+ inf
98
+ inf
99
+ inf
100
+ inf
101
+ inf
102
+ inf
103
+ inf
104
+ inf
105
+ inf
106
+ -10600
107
+ inf
108
+ inf
109
+ inf
110
+ inf
111
+ inf
112
+ inf
113
+ inf
114
+ inf
115
+ inf
116
+ inf
117
+ inf
118
+ inf
119
+ inf
120
+ inf
121
+ -7800
122
+ inf
123
+ inf
124
+ inf
125
+ inf
126
+ inf
127
+ inf
128
+ inf
129
+ inf
130
+ inf
131
+ inf
132
+ inf
133
+ inf
134
+ inf
135
+ inf
136
+ -8200
137
+ inf
138
+ inf
139
+ inf
140
+ inf
141
+ inf
142
+ inf
143
+ inf
144
+ inf
145
+ inf
146
+ inf
147
+ inf
148
+ inf
149
+ inf
150
+ inf
151
+ -9800
152
+ inf
153
+ inf
154
+ inf
155
+ inf
156
+ inf
157
+ inf
158
+ inf
159
+ inf
160
+ inf
161
+ inf
162
+ inf
163
+ inf
164
+ inf
165
+ inf
166
+ -8000
167
+ inf
168
+ inf
169
+ inf
170
+ inf
171
+ inf
172
+ inf
173
+ inf
174
+ inf
175
+ inf
176
+ inf
177
+ inf
178
+ inf
179
+ inf
180
+ inf
181
+ -8400
182
+ inf
183
+ inf
184
+ inf
185
+ inf
186
+ inf
187
+ inf
188
+ inf
189
+ inf
190
+ inf
191
+ inf
192
+ inf
193
+ inf
194
+ inf
195
+ inf
196
+ -7200
197
+ inf
198
+ inf
199
+ inf
200
+ inf
201
+ inf
202
+ inf
203
+ inf
204
+ inf
205
+ inf
206
+ inf
207
+ inf
208
+ inf
209
+ inf
210
+ inf
211
+ -8200
212
+ inf
213
+ inf
214
+ inf
215
+ inf
216
+ inf
217
+ inf
218
+ inf
219
+ inf
220
+ inf
221
+ inf
222
+ inf
223
+ inf
224
+ inf
225
+ inf
226
+ -8500
227
+ inf
228
+ inf
229
+ inf
230
+ inf
231
+ inf
232
+ inf
233
+ inf
234
+ inf
235
+ inf
236
+ inf
237
+ inf
238
+ inf
239
+ inf
240
+ inf
241
+ -7900
242
+ inf
243
+ inf
244
+ inf
245
+ inf
246
+ inf
247
+ inf
248
+ inf
249
+ inf
250
+ inf
251
+ inf
252
+ inf
253
+ inf
254
+ inf
255
+ inf
256
+ inf
data/conf/primer3_config/stack.ds ADDED
@@ -0,0 +1,256 @@
1
+ inf
2
+ inf
3
+ inf
4
+ inf
5
+ inf
6
+ inf
7
+ inf
8
+ inf
9
+ inf
10
+ inf
11
+ inf
12
+ inf
13
+ inf
14
+ inf
15
+ inf
16
+ -22.2
17
+ inf
18
+ inf
19
+ inf
20
+ inf
21
+ inf
22
+ inf
23
+ inf
24
+ inf
25
+ inf
26
+ inf
27
+ inf
28
+ inf
29
+ inf
30
+ inf
31
+ -22.4
32
+ inf
33
+ inf
34
+ inf
35
+ inf
36
+ inf
37
+ inf
38
+ inf
39
+ inf
40
+ inf
41
+ inf
42
+ inf
43
+ inf
44
+ inf
45
+ inf
46
+ -21.0
47
+ inf
48
+ inf
49
+ inf
50
+ inf
51
+ inf
52
+ inf
53
+ inf
54
+ inf
55
+ inf
56
+ inf
57
+ inf
58
+ inf
59
+ inf
60
+ inf
61
+ -20.4
62
+ inf
63
+ inf
64
+ inf
65
+ inf
66
+ inf
67
+ inf
68
+ inf
69
+ inf
70
+ inf
71
+ inf
72
+ inf
73
+ inf
74
+ inf
75
+ inf
76
+ -22.7
77
+ inf
78
+ inf
79
+ inf
80
+ inf
81
+ inf
82
+ inf
83
+ inf
84
+ inf
85
+ inf
86
+ inf
87
+ inf
88
+ inf
89
+ inf
90
+ inf
91
+ -19.9
92
+ inf
93
+ inf
94
+ inf
95
+ inf
96
+ inf
97
+ inf
98
+ inf
99
+ inf
100
+ inf
101
+ inf
102
+ inf
103
+ inf
104
+ inf
105
+ inf
106
+ -27.2
107
+ inf
108
+ inf
109
+ inf
110
+ inf
111
+ inf
112
+ inf
113
+ inf
114
+ inf
115
+ inf
116
+ inf
117
+ inf
118
+ inf
119
+ inf
120
+ inf
121
+ -21.0
122
+ inf
123
+ inf
124
+ inf
125
+ inf
126
+ inf
127
+ inf
128
+ inf
129
+ inf
130
+ inf
131
+ inf
132
+ inf
133
+ inf
134
+ inf
135
+ inf
136
+ -22.2
137
+ inf
138
+ inf
139
+ inf
140
+ inf
141
+ inf
142
+ inf
143
+ inf
144
+ inf
145
+ inf
146
+ inf
147
+ inf
148
+ inf
149
+ inf
150
+ inf
151
+ -24.4
152
+ inf
153
+ inf
154
+ inf
155
+ inf
156
+ inf
157
+ inf
158
+ inf
159
+ inf
160
+ inf
161
+ inf
162
+ inf
163
+ inf
164
+ inf
165
+ inf
166
+ -19.9
167
+ inf
168
+ inf
169
+ inf
170
+ inf
171
+ inf
172
+ inf
173
+ inf
174
+ inf
175
+ inf
176
+ inf
177
+ inf
178
+ inf
179
+ inf
180
+ inf
181
+ -22.4
182
+ inf
183
+ inf
184
+ inf
185
+ inf
186
+ inf
187
+ inf
188
+ inf
189
+ inf
190
+ inf
191
+ inf
192
+ inf
193
+ inf
194
+ inf
195
+ inf
196
+ -21.3
197
+ inf
198
+ inf
199
+ inf
200
+ inf
201
+ inf
202
+ inf
203
+ inf
204
+ inf
205
+ inf
206
+ inf
207
+ inf
208
+ inf
209
+ inf
210
+ inf
211
+ -22.2
212
+ inf
213
+ inf
214
+ inf
215
+ inf
216
+ inf
217
+ inf
218
+ inf
219
+ inf
220
+ inf
221
+ inf
222
+ inf
223
+ inf
224
+ inf
225
+ inf
226
+ -22.7
227
+ inf
228
+ inf
229
+ inf
230
+ inf
231
+ inf
232
+ inf
233
+ inf
234
+ inf
235
+ inf
236
+ inf
237
+ inf
238
+ inf
239
+ inf
240
+ inf
241
+ -22.2
242
+ inf
243
+ inf
244
+ inf
245
+ inf
246
+ inf
247
+ inf
248
+ inf
249
+ inf
250
+ inf
251
+ inf
252
+ inf
253
+ inf
254
+ inf
255
+ inf
256
+ inf