@@ -123,11 +123,11 @@ Here our analysis results are stored in potato directory. File descriptions are
123
123
124
124
$ cd /home/zhangtaolab/data/chorus2_project/analysis/potato
125
125
$ head -5 potato_dm_v404_all_pm_un.fasta_all.bed
126
- chr01 97858 97902 ATTTTCCATGGACCTCATTAAGATTAGCTATTGAACCAGTTACCC
127
- chr01 103008 103052 ACAGCCAAATCGTCCCATATTCAAGGATAAACGACCCACGAATCA
128
- chr01 127095 127139 ATCTATATCTACTACACCAGAATATTCATACACAAATAAATTACT
129
- chr01 127101 127145 ATCTACTACACCAGAATATTCATACACAAATAAATTACTACTATT
130
- chr01 127815 127859 CTCAGATTTACCGAATTATTCCTGTGACAAAATATTACTTCCAGT
126
+ chr01 97858 97902 ATTTTCCATGGACCTCATTAAGATTAGCTATTGAACCAGTTACCC
127
+ chr01 103008 103052 ACAGCCAAATCGTCCCATATTCAAGGATAAACGACCCACGAATCA
128
+ chr01 127095 127139 ATCTATATCTACTACACCAGAATATTCATACACAAATAAATTACT
129
+ chr01 127101 127145 ATCTACTACACCAGAATATTCATACACAAATAAATTACTACTATT
130
+ chr01 127815 127859 CTCAGATTTACCGAATTATTCCTGTGACAAAATATTACTTCCAGT
131
131
132
132
There are four columns in each row, first column is chromosome name, second is oligo start site, third is oligo end site, the last one is oligo probe sequence. Users can use excel or text editor to open this file.
133
133
@@ -154,13 +154,13 @@ After running NGS filtering, three files (*\*.jf, \*.bw, \*.bed*) will output to
154
154
.. code-block :: bash
155
155
156
156
$ head -5 potato_dm_v404_all_pm_un.fasta_kmer.bed
157
- chr01 97858 97902 ATTTTCCATGGACCTCATTAAGATTAGCTATTGAACCAGTTACCC 324 +
158
- chr01 103008 103052 ACAGCCAAATCGTCCCATATTCAAGGATAAACGACCCACGAATCA 146 +
159
- chr01 127095 127139 ATCTATATCTACTACACCAGAATATTCATACACAAATAAATTACT 318 +
160
- chr01 127101 127145 ATCTACTACACCAGAATATTCATACACAAATAAATTACTACTATT 373 +
161
- chr01 127815 127859 CTCAGATTTACCGAATTATTCCTGTGACAAAATATTACTTCCAGT 281 +
157
+ chr01 97858 97902 ATTTTCCATGGACCTCATTAAGATTAGCTATTGAACCAGTTACCC 324 +
158
+ chr01 103008 103052 ACAGCCAAATCGTCCCATATTCAAGGATAAACGACCCACGAATCA 146 +
159
+ chr01 127095 127139 ATCTATATCTACTACACCAGAATATTCATACACAAATAAATTACT 318 +
160
+ chr01 127101 127145 ATCTACTACACCAGAATATTCATACACAAATAAATTACTACTATT 373 +
161
+ chr01 127815 127859 CTCAGATTTACCGAATTATTCCTGTGACAAAATATTACTTCCAGT 281 +
162
162
163
- There are six columns in each row, first four columns are the same as *Tpotato_dm_v404_all_pm_un .fasta_all.bed *.
163
+ There are six columns in each row, first four columns are the same as *potato_dm_v404_all_pm_un .fasta_all.bed *.
164
164
The fifth column is the k-mer score, last column is target strand of probes.
165
165
166
166
@@ -172,11 +172,11 @@ Confident probes need to be further selected based on the k-mer of each oligo. C
172
172
173
173
$ ChorusNGSselect -i potato_DM_v404_filtered.bed -o potato_DM_v404_selected.bed
174
174
$ head -5 potato_DM_v404_selected.bed
175
- chr01 97858 97902 ATTTTCCATGGACCTCATTAAGATTAGCTATTGAACCAGTTACCC 324 +
176
- chr01 103008 103052 TGATTCGTGGGTCGTTTATCCTTGAATATGGGACGATTTGGCTGT 146 -
177
- chr01 127095 127139 ATCTATATCTACTACACCAGAATATTCATACACAAATAAATTACT 318 +
178
- chr01 127815 127859 ACTGGAAGTAATATTTTGTCACAGGAATAATTCGGTAAATCTGAG 281 -
179
- chr01 133524 133568 ATTAATCAATAAAGGAAAAGCAAGGTTGGAATACGGTTTCATCCT 389 +
175
+ chr01 97858 97902 ATTTTCCATGGACCTCATTAAGATTAGCTATTGAACCAGTTACCC 324 +
176
+ chr01 103008 103052 TGATTCGTGGGTCGTTTATCCTTGAATATGGGACGATTTGGCTGT 146 -
177
+ chr01 127095 127139 ATCTATATCTACTACACCAGAATATTCATACACAAATAAATTACT 318 +
178
+ chr01 127815 127859 ACTGGAAGTAATATTTTGTCACAGGAATAATTCGGTAAATCTGAG 281 -
179
+ chr01 133524 133568 ATTAATCAATAAAGGAAAAGCAAGGTTGGAATACGGTTTCATCCT 389 +
180
180
181
181
There are six columns in each row, which are the same as ChorusNGSfilter output.
182
182
The final probes can be synthesized directly for oligo-FISH or imported into ChorusPBGUI for further selection.
0 commit comments