Export bed file.

[lightning.git] / example-su92l-1kg.sh

diff --git a/example-su92l-1kg.sh b/example-su92l-1kg.sh
index b73fcfed50c0368b76f7bb8299fb187714783ce3..1e5ecbc1b348535695714d096c7e88c173202c24 100755 (executable)
--- a/example-su92l-1kg.sh
+++ b/example-su92l-1kg.sh
@@ -16,10 +16,12 @@ tagset=su92l-4zz18-92bx4zjg5hgs3yc/tagset.fa.gz
 
 genome=$(lightning     ref2genome   -project ${project} -priority ${priority} -ref ${ref_fa})
 fasta=$(lightning      vcf2fasta    -project ${project} -priority ${priority} -ref ${ref_fa} -genome ${genome} -mask=true ${gvcf})
-unfiltered=$(lightning import       -project ${project} -priority ${priority} -tag-library ${tagset} ${fasta})
+unfiltered=$(lightning import       -project ${project} -priority ${priority} -tag-library ${tagset} -skip-ooo=true ${fasta})
+stats=$(lightning      stats        -project ${project} -priority ${priority} -i ${unfiltered})
 filtered=$(lightning   filter       -project ${project} -priority ${priority} -i ${unfiltered} -min-coverage "0.9" -max-variants "30")
-numpy=$(lightning      export-numpy -project ${project} -priority ${priority} -i ${filtered})
-pca=$(lightning        pca          -project ${project} -priority ${priority} -i ${numpy})
+#numpy=$(lightning     export-numpy -project ${project} -priority ${priority} -i ${filtered} -one-hot)
+#pca=$(lightning       pca-py       -project ${project} -priority ${priority} -i ${numpy})
+pca=$(lightning        pca-go       -project ${project} -priority ${priority} -i ${filtered} -one-hot)
 plot=$(lightning       plot         -project ${project} -priority ${priority} -i ${pca} -labels-csv ${info}/sample_info.csv -sample-fasta-dir ${fasta})
 echo >&2 "https://workbench2.${plot%%-*}.arvadosapi.com/collections/${plot}"
 echo ${plot%%/*}