You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
There are actually identical to the equivalent module at [scripts folder](https://github.com/jsarriaa/PHGv2Tools/tree/main/Modules) but it is coded to work different just if the file is ending in .py or .ipynb.
When defining the main function, note the block of code where input variables are defined. When working in a notebook declare there whatever you need. Whole script will run with them.
Just on top of the main function, see the "AmIaNotebook" function. You can declare there whatever is provide as argument on the command line function. It will work on the same way.
Find examples of usage at the [Example database](https://github.com/jsarriaa/PHGv2Tools/tree/main/Example_database)
51
+
52
+
>>>>>>> origin/main
48
53
## Installation
49
54
To use this package, a conda environment is used. It is a modified version of the original phg one, adding the pygenometracks package for python plotting. For everything further needed, it will be updated.
50
55
@@ -112,7 +117,7 @@ Extrapolating this into the PHG ranges, it is useful to have the information of
112
117
- sector diagram plotting the % of core, accessory and unique haplotypes.
113
118
It takes as arguments:
114
119
```
115
-
--pangenome-hvcf // -hvcf <Merged of all hvcf of the pangenome database> #Built with phg merge-hvcfs
120
+
--pangenome-hvcf // -phvcf <Merged of all hvcf of the pangenome database> #Built with phg merge-hvcfs
@@ -150,12 +155,12 @@ For those who are looking for genes alleles; Ranges are built with genes annotat
150
155
151
156
These are the mandatory arguments:
152
157
```
153
-
--hvcf // -hf <Merged of all hvcf of the pangenome database> #Built with phg merge-hvcfs
154
-
--reference-fasta // -ref <Reference.fasta> #Built with phg prepare-assemblies
158
+
--pangenome-hvcf // -phvcf <Merged of all hvcf of the pangenome database> #Built with phg merge-hvcfs
159
+
--reference-fasta // -ref-fa <Reference.fasta> #Built with phg prepare-assemblies
155
160
```
156
161
And then, provide the coordinates. These are optional, but will be asked ahead if there are not as a command line argument.
157
162
```
158
-
--chromosome // -c <INT> i.e. <2>
163
+
--chromosome // -chr <INT> i.e. <2>
159
164
--start // -s <START> i.e. <1000>
160
165
--end // -e <END> i.e. <5000>
161
166
```
@@ -189,18 +194,18 @@ Taking an imputed h.vcf from ```phg map-kmers``` and ```find-paths``` and plots
189
194
190
195
```plot-imputed-hvcf``` takes as argument:
191
196
```
192
-
--input-hvcf // -hvcf <Imputed h.vcf> #Built with phg find-paths
193
-
--reference-hvcf // -ref <Reference.h.vcf> #Built with phg create-ref-vcf
194
-
--pangenome-hvcf-folder // -folder <Folder path to the h.vcf files of the pangenome database> #Built with phg create-maf-vcf
197
+
--imputed-hvcf // -ihvcf <Imputed h.vcf> #Built with phg find-paths
198
+
--reference-hvcf // -refhvcf <Reference.h.vcf> #Built with phg create-ref-vcf
199
+
--pangenome-hvcf-folder // -pfolder <Folder path to the h.vcf files of the pangenome database> #Built with phg create-maf-vcf
195
200
```
196
201
197
202
198
203
199
204
-#### Check identity against pangenome
200
205
Reads all pangenome haplotypes and compare with an imputated h.vcf, to extract the percentage of the genome is used to build up the imputed haplotype. It is now only based in nº ranges itself, not based in base pairs absolute amount. Call the script with ```check-imputated-haplotype``` and provide as arguments:
201
206
```
202
-
--hvcf-folder // -folder <Folder path to the h.vcf files of the pangenome database> #Built with phg create-maf-vcf
203
-
--hvcf-file // -file <Imputed h.vcf> #Built with phg find-paths
207
+
--pangenome-hvcf-folder // -pf <Folder path to the h.vcf files of the pangenome database> #Built with phg create-maf-vcf
208
+
--imputed-hvcf // -ihvcf <Imputed h.vcf> #Built with phg find-paths
204
209
```
205
210
Working on plotting the results, for now we get a list for the results:
0 commit comments