bio-polymarker 1.3.2

Sign up to get free protection for your applications and to get access to all the features.
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