From 317e93cad17b384e3422273f3774dadb3029ad02 Mon Sep 17 00:00:00 2001 From: Zilong-Li Date: Fri, 6 Sep 2024 10:22:21 +0200 Subject: [PATCH] v0.5.0: two api of BcfRecord breaking and copyHeader updated --- readme.org | 21 +- test/index.vcf | 728 ------------------------------------------- test/test.vcf.gz | Bin 0 -> 23728 bytes test/test.vcf.gz.tbi | Bin 0 -> 136 bytes vcfpp.h | 18 +- 5 files changed, 32 insertions(+), 735 deletions(-) delete mode 100644 test/index.vcf create mode 100644 test/test.vcf.gz create mode 100644 test/test.vcf.gz.tbi diff --git a/readme.org b/readme.org index 05d15b8..3eded76 100644 --- a/readme.org +++ b/readme.org @@ -116,21 +116,36 @@ bw.writeLine("chr20\t2006060\trs146931526\tG\tC\t100\tPASS\tAF=0.000998403\tGT\t *** Use an existing VCF as template In this example, we first read VCF file =test/test-vcf-read.vcf.gz=. Second, -we construct an empty variant record and update the record with the input +we construct a variant record and update the record with the input VCF. Third, we construct a BcfWriter object using the meta data in the header of the input VCF, writing out the header and the modified variant record. #+begin_src C++ -BcfReader br("test/test-vcf-read.vcf.gz"); +BcfReader br("test/test.vcf.gz"); +BcfWriter bw("out.vcf.gz", br.header); BcfRecord var(br.header); br.getNextVariant(var); -BcfWriter bw("out.vcf.gz", br.header); bw.writeHeader(); var.setPOS(100001); // update the POS of the variant bw.writeRecord(var); #+end_src +In some cases where we want to modify the samples of a template VCF, we +can associate the =BcfWriter= and =BcfRecord= with another modifiable header +instead of the non-modifiable header as the above. + +#+begin_src C++ +BcfReader br("test/test.vcf.gz"); +BcfWriter bw("out.vcf.gz"); +// copy header of template vcf and restrict to target sample +bw.copyHeader("test/test.vcf.gz", "HG00097"); +BcfRecord var(bw.header); +br.getNextVariant(var); +var.setPOS(100001); +bw.writeRecord(var); // output only sample HG00097 +#+end_src + ** Variants Operation All variants related API can be found [[https://zilongli.org/proj/vcfpp/classvcfpp_1_1_bcf_record][BcfRecord]]. The commonly used are listed below. diff --git a/test/index.vcf b/test/index.vcf deleted file mode 100644 index b32e2d2..0000000 --- a/test/index.vcf +++ /dev/null @@ -1,728 +0,0 @@ -##fileformat=VCFv4.2 -##FILTER= -##bcftoolsVersion=1.8-31-g9ba4024+htslib-1.8-32-g6e87a1e-dirty -##bcftoolsCommand=mpileup --fasta-ref /nfs/srpipe_references/references/Human/1000Genomes_hs37d5/all/fasta/hs37d5.fa test/index.bam -##reference=file:///nfs/srpipe_references/references/Human/1000Genomes_hs37d5/all/fasta/hs37d5.fa -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##contig= -##ALT= -##INFO= -##INFO= -##INFO= -##INFO= -##INFO= -##INFO= -##INFO= -##INFO= -##INFO= -##INFO= -##INFO= -##INFO= -##INFO= -##FORMAT= -#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT ERS220911 -1 9999919 . G <*> 0 . DP=1;I16=1,0,0,0,26,676,0,0,60,3600,0,0,0,0,0,0;QS=1,0;MQ0F=0 PL 0,3,26 -1 9999920 . T <*> 0 . DP=1;I16=1,0,0,0,34,1156,0,0,60,3600,0,0,1,1,0,0;QS=1,0;MQ0F=0 PL 0,3,34 -1 9999921 . A <*> 0 . DP=1;I16=1,0,0,0,33,1089,0,0,60,3600,0,0,2,4,0,0;QS=1,0;MQ0F=0 PL 0,3,33 -1 9999922 . A <*> 0 . DP=1;I16=1,0,0,0,34,1156,0,0,60,3600,0,0,3,9,0,0;QS=1,0;MQ0F=0 PL 0,3,34 -1 9999923 . T <*> 0 . DP=1;I16=1,0,0,0,35,1225,0,0,60,3600,0,0,4,16,0,0;QS=1,0;MQ0F=0 PL 0,3,35 -1 9999924 . C <*> 0 . DP=1;I16=1,0,0,0,35,1225,0,0,60,3600,0,0,5,25,0,0;QS=1,0;MQ0F=0 PL 0,3,35 -1 9999925 . C <*> 0 . DP=1;I16=1,0,0,0,36,1296,0,0,60,3600,0,0,6,36,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -1 9999926 . C <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,7,49,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999927 . A <*> 0 . DP=1;I16=1,0,0,0,36,1296,0,0,60,3600,0,0,8,64,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -1 9999928 . G <*> 0 . DP=1;I16=1,0,0,0,34,1156,0,0,60,3600,0,0,9,81,0,0;QS=1,0;MQ0F=0 PL 0,3,34 -1 9999929 . C <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,10,100,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -1 9999930 . A <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,11,121,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999931 . C <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,12,144,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999932 . T <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,13,169,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999933 . T <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,14,196,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -1 9999934 . T <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,15,225,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999935 . A <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,16,256,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -1 9999936 . G <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,17,289,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -1 9999937 . G <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,18,324,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -1 9999938 . A <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,19,361,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -1 9999939 . G <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,20,400,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -1 9999940 . G <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,21,441,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999941 . C <*> 0 . DP=1;I16=1,0,0,0,41,1681,0,0,60,3600,0,0,22,484,0,0;QS=1,0;MQ0F=0 PL 0,3,41 -1 9999942 . T <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,23,529,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -1 9999943 . A <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,24,576,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -1 9999944 . A <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999945 . G <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -1 9999946 . G <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999947 . C <*> 0 . DP=1;I16=1,0,0,0,35,1225,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,35 -1 9999948 . A <*> 0 . DP=1;I16=1,0,0,0,34,1156,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,34 -1 9999949 . G <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999950 . G <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999951 . C <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -1 9999952 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999953 . G <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999954 . A <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -1 9999955 . T <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -1 9999956 . C <*> 0 . DP=1;I16=1,0,0,0,41,1681,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,41 -1 9999957 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999958 . C <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999959 . T <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999960 . T <*> 0 . DP=1;I16=1,0,0,0,35,1225,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,35 -1 9999961 . G <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999962 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999963 . G <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999964 . A <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999965 . C <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999966 . C <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -1 9999967 . A <*> 0 . DP=1;I16=1,0,0,0,41,1681,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,41 -1 9999968 . G <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999969 . G <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999970 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999971 . G <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999972 . T <*> 0 . DP=1;I16=1,0,0,0,36,1296,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -1 9999973 . T <*> 0 . DP=1;I16=1,0,0,0,36,1296,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -1 9999974 . A <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -1 9999975 . C <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -1 9999976 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999977 . G <*> 0 . DP=1;I16=1,0,0,0,36,1296,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -1 9999978 . A <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -1 9999979 . C <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999980 . C <*> 0 . DP=1;I16=1,0,0,0,33,1089,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,33 -1 9999981 . A <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -1 9999982 . G <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -1 9999983 . C <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999984 . C <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999985 . T <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999986 . G <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999987 . G <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -1 9999988 . C <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -1 9999989 . C <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999990 . G <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999991 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999992 . C <*> 0 . DP=1;I16=1,0,0,0,36,1296,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -1 9999993 . A <*> 0 . DP=1;I16=1,0,0,0,36,1296,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -1 9999994 . C <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,24,576,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 9999995 . G <*> 0 . DP=1;I16=1,0,0,0,33,1089,0,0,60,3600,0,0,23,529,0,0;QS=1,0;MQ0F=0 PL 0,3,33 -1 9999996 . G <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,22,484,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -1 9999997 . C <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,21,441,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 9999998 . G <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,20,400,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -1 9999999 . A <*> 0 . DP=1;I16=1,0,0,0,31,961,0,0,60,3600,0,0,19,361,0,0;QS=1,0;MQ0F=0 PL 0,3,31 -1 10000000 . A <*> 0 . DP=1;I16=1,0,0,0,33,1089,0,0,60,3600,0,0,18,324,0,0;QS=1,0;MQ0F=0 PL 0,3,33 -1 10000001 . A <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,17,289,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -1 10000002 . C <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,16,256,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -1 10000003 . C <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,15,225,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -1 10000004 . C <*> 0 . DP=1;I16=1,0,0,0,29,841,0,0,60,3600,0,0,14,196,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000005 . C <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,13,169,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -1 10000006 . G <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,12,144,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -1 10000007 . T <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,11,121,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -1 10000008 . C <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,10,100,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -1 10000009 . T <*> 0 . DP=1;I16=1,0,0,0,43,1849,0,0,60,3600,0,0,9,81,0,0;QS=1,0;MQ0F=0 PL 0,3,43 -1 10000010 . C <*> 0 . DP=2;I16=1,1,0,0,59,2105,0,0,89,4441,0,0,8,64,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,59 -1 10000011 . T <*> 0 . DP=2;I16=1,1,0,0,76,2888,0,0,89,4441,0,0,8,50,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,67 -1 10000012 . A <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,89,4441,0,0,8,40,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,67 -1 10000013 . C <*> 0 . DP=2;I16=1,1,0,0,66,2250,0,0,89,4441,0,0,8,34,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,66 -1 10000014 . A <*> 0 . DP=2;I16=1,1,0,0,67,2285,0,0,89,4441,0,0,8,32,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,67 -1 10000015 . A <*> 0 . DP=2;I16=1,1,0,0,69,2385,0,0,89,4441,0,0,8,34,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,65 -1 10000016 . T <*> 0 . DP=2;I16=1,1,0,0,75,2817,0,0,89,4441,0,0,8,40,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,68 -1 10000017 . A <*> 0 . DP=2;I16=1,1,0,0,67,2285,0,0,89,4441,0,0,8,50,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,58 -1 10000018 . A <*> 0 . DP=2;I16=1,1,0,0,64,2120,0,0,89,4441,0,0,8,64,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,55 -1 10000019 . A <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,9,81,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000020 . T <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,10,100,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000021 . T <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,29,841,0,0,11,121,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000022 . A <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,29,841,0,0,12,144,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000023 . A <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,13,169,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000024 . A <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,14,196,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000025 . A <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,29,841,0,0,15,225,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000026 . T <*> 0 . DP=1;I16=0,1,0,0,29,841,0,0,29,841,0,0,16,256,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000027 . A <*> 0 . DP=1;I16=0,1,0,0,40,1600,0,0,29,841,0,0,17,289,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000028 . T <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,18,324,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000029 . T <*> 0 . DP=1;I16=0,1,0,0,41,1681,0,0,29,841,0,0,19,361,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000030 . A <*> 0 . DP=1;I16=0,1,0,0,36,1296,0,0,29,841,0,0,20,400,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000031 . G <*> 0 . DP=1;I16=0,1,0,0,40,1600,0,0,29,841,0,0,21,441,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000032 . C <*> 0 . DP=1;I16=0,1,0,0,40,1600,0,0,29,841,0,0,22,484,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000033 . T <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,29,841,0,0,23,529,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000034 . G <*> 0 . DP=1;I16=0,1,0,0,36,1296,0,0,29,841,0,0,24,576,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000035 . G <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000036 . G <*> 0 . DP=1;I16=0,1,0,0,42,1764,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000037 . C <*> 0 . DP=1;I16=0,1,0,0,34,1156,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000038 . A <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000039 . T <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000040 . G <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000041 . G <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000042 . T <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000043 . G <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000044 . G <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000045 . T <*> 0 . DP=1;I16=0,1,0,0,42,1764,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000046 . G <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000047 . T <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000048 . G <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000049 . T <*> 0 . DP=1;I16=0,1,0,0,41,1681,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000050 . G <*> 0 . DP=1;I16=0,1,0,0,31,961,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000051 . C <*> 0 . DP=1;I16=0,1,0,0,16,256,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,16 -1 10000052 . T <*> 0 . DP=1;I16=0,1,0,0,31,961,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000053 . T <*> 0 . DP=1;I16=0,1,0,0,35,1225,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000054 . G <*> 0 . DP=1;I16=0,1,0,0,40,1600,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000055 . T <*> 0 . DP=1;I16=0,1,0,0,33,1089,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000056 . A <*> 0 . DP=1;I16=0,1,0,0,22,484,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,22 -1 10000057 . G <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000058 . T <*> 0 . DP=1;I16=0,1,0,0,40,1600,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000059 . C <*> 0 . DP=1;I16=0,1,0,0,34,1156,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000060 . C <*> 0 . DP=1;I16=0,1,0,0,40,1600,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000061 . C <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000062 . A <*> 0 . DP=1;I16=0,1,0,0,34,1156,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000063 . G <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000064 . C <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000065 . T <*> 0 . DP=1;I16=0,1,0,0,40,1600,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000066 . A <*> 0 . DP=1;I16=0,1,0,0,32,1024,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000067 . C <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000068 . T <*> 0 . DP=1;I16=0,1,0,0,36,1296,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000069 . T <*> 0 . DP=1;I16=0,1,0,0,36,1296,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000070 . G <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000071 . G <*> 0 . DP=1;I16=0,1,0,0,33,1089,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000072 . C <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000073 . G <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000074 . G <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000075 . G <*> 0 . DP=1;I16=0,1,0,0,36,1296,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000076 . C <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000077 . T <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000078 . G <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000079 . A <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000080 . G <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000081 . G <*> 0 . DP=1;I16=0,1,0,0,36,1296,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000082 . T <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000083 . G <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000084 . G <*> 0 . DP=1;I16=0,1,0,0,41,1681,0,0,29,841,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000085 . G <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,24,576,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000086 . A <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,23,529,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000087 . G <*> 0 . DP=1;I16=0,1,0,0,40,1600,0,0,29,841,0,0,22,484,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000088 . A <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,29,841,0,0,21,441,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000089 . A <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,20,400,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000090 . T <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,29,841,0,0,19,361,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000091 . C <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,29,841,0,0,18,324,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000092 . A <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,17,289,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000093 . T <*> 0 . DP=1;I16=0,1,0,0,33,1089,0,0,29,841,0,0,16,256,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000094 . C <*> 0 . DP=1;I16=0,1,0,0,41,1681,0,0,29,841,0,0,15,225,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000095 . C <*> 0 . DP=1;I16=0,1,0,0,36,1296,0,0,29,841,0,0,14,196,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000096 . A <*> 0 . DP=1;I16=0,1,0,0,41,1681,0,0,29,841,0,0,13,169,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000097 . A <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,12,144,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000098 . G <*> 0 . DP=1;I16=0,1,0,0,36,1296,0,0,29,841,0,0,11,121,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000099 . C <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,10,100,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000100 . C <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,29,841,0,0,9,81,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000101 . T <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,8,64,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000102 . T <*> 0 . DP=1;I16=0,1,0,0,36,1296,0,0,29,841,0,0,7,49,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000103 . G <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,29,841,0,0,6,36,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000104 . G <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,29,841,0,0,5,25,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000105 . A <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,29,841,0,0,4,16,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000106 . G <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,29,841,0,0,3,9,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000107 . G <*> 0 . DP=1;I16=0,1,0,0,33,1089,0,0,29,841,0,0,2,4,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000108 . C <*> 0 . DP=1;I16=0,1,0,0,32,1024,0,0,29,841,0,0,1,1,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -1 10000109 . A <*> 0 . DP=1;I16=0,1,0,0,35,1225,0,0,29,841,0,0,0,0,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -2 4999907 . C <*> 0 . DP=1;I16=1,0,0,0,33,1089,0,0,60,3600,0,0,0,0,0,0;QS=1,0;MQ0F=0 PL 0,3,33 -2 4999908 . C <*> 0 . DP=1;I16=1,0,0,0,32,1024,0,0,60,3600,0,0,1,1,0,0;QS=1,0;MQ0F=0 PL 0,3,32 -2 4999909 . A <*> 0 . DP=1;I16=1,0,0,0,33,1089,0,0,60,3600,0,0,2,4,0,0;QS=1,0;MQ0F=0 PL 0,3,33 -2 4999910 . G <*> 0 . DP=1;I16=1,0,0,0,36,1296,0,0,60,3600,0,0,3,9,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -2 4999911 . G <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,4,16,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 4999912 . A <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,5,25,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 4999913 . C <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,6,36,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 4999914 . A <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,7,49,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 4999915 . T <*> 0 . DP=1;I16=1,0,0,0,41,1681,0,0,60,3600,0,0,8,64,0,0;QS=1,0;MQ0F=0 PL 0,3,41 -2 4999916 . C <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,9,81,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 4999917 . C <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,10,100,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 4999918 . A <*> 0 . DP=1;I16=1,0,0,0,35,1225,0,0,60,3600,0,0,11,121,0,0;QS=1,0;MQ0F=0 PL 0,3,35 -2 4999919 . G <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,12,144,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 4999920 . G <*> 0 . DP=1;I16=1,0,0,0,41,1681,0,0,60,3600,0,0,13,169,0,0;QS=1,0;MQ0F=0 PL 0,3,41 -2 4999921 . C <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,14,196,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 4999922 . T <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,15,225,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 4999923 . T <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,16,256,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 4999924 . T <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,17,289,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 4999925 . T <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,18,324,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 4999926 . C <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,19,361,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 4999927 . T <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,20,400,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 4999928 . G <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,21,441,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 4999929 . T <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,22,484,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 4999930 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,23,529,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 4999931 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,24,576,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 4999932 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 4999933 . T <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 4999934 . C <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 4999935 . T <*> 0 . DP=1;I16=1,0,0,0,36,1296,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -2 4999936 . T <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 4999937 . C <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 4999938 . T <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 4999939 . G <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 4999940 . A <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 4999941 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 4999942 . A <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 4999943 . T <*> 0 . DP=1;I16=1,0,0,0,41,1681,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,41 -2 4999944 . C <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 4999945 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 4999946 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 4999947 . G <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 4999948 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 4999949 . C <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 4999950 . A <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 4999951 . G <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 4999952 . A <*> 0 . DP=2;I16=1,1,0,0,64,2120,0,0,120,7200,0,0,25,625,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,64 -2 4999953 . G <*> 0 . DP=2;I16=1,1,0,0,60,1962,0,0,120,7200,0,0,26,626,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,60 -2 4999954 . T <*> 0 . DP=2;I16=1,1,0,0,67,2257,0,0,120,7200,0,0,27,629,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,67 -2 4999955 . C <*> 0 . DP=2;I16=1,1,0,0,57,1805,0,0,120,7200,0,0,28,634,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,57 -2 4999956 . T <*> 0 . DP=2;I16=1,1,0,0,77,2969,0,0,120,7200,0,0,29,641,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -2 4999957 . C <*> 0 . DP=2;I16=1,1,0,0,60,2000,0,0,120,7200,0,0,30,650,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,60 -2 4999958 . A <*> 0 . DP=2;I16=1,1,0,0,71,2561,0,0,120,7200,0,0,31,661,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,71 -2 4999959 . C <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,120,7200,0,0,32,674,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -2 4999960 . A <*> 0 . DP=2;I16=1,1,0,0,71,2561,0,0,120,7200,0,0,33,689,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,71 -2 4999961 . A <*> 0 . DP=2;I16=1,1,0,0,71,2521,0,0,120,7200,0,0,34,706,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,71 -2 4999962 . G <*> 0 . DP=2;I16=1,1,0,0,79,3121,0,0,120,7200,0,0,35,725,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,79 -2 4999963 . C <*> 0 . DP=2;I16=1,1,0,0,70,2452,0,0,120,7200,0,0,36,746,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,70 -2 4999964 . C <*> 0 . DP=2;I16=1,1,0,0,78,3042,0,0,120,7200,0,0,37,769,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,78 -2 4999965 . T <*> 0 . DP=2;I16=1,1,0,0,79,3121,0,0,120,7200,0,0,38,794,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,79 -2 4999966 . T <*> 0 . DP=2;I16=1,1,0,0,76,2888,0,0,120,7200,0,0,39,821,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,76 -2 4999967 . A <*> 0 . DP=2;I16=1,1,0,0,79,3125,0,0,120,7200,0,0,40,850,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,79 -2 4999968 . A <*> 0 . DP=2;I16=1,1,0,0,73,2669,0,0,120,7200,0,0,41,881,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,73 -2 4999969 . T <*> 0 . DP=2;I16=1,1,0,0,73,2669,0,0,120,7200,0,0,42,914,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,73 -2 4999970 . T <*> 0 . DP=2;I16=1,1,0,0,78,3042,0,0,120,7200,0,0,43,949,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,78 -2 4999971 . C <*> 0 . DP=2;I16=1,1,0,0,80,3208,0,0,120,7200,0,0,44,986,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,80 -2 4999972 . T <*> 0 . DP=2;I16=1,1,0,0,75,2817,0,0,120,7200,0,0,45,1025,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,75 -2 4999973 . T <*> 0 . DP=2;I16=1,1,0,0,78,3044,0,0,120,7200,0,0,46,1066,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,78 -2 4999974 . G <*> 0 . DP=2;I16=1,1,0,0,76,2888,0,0,120,7200,0,0,47,1109,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,76 -2 4999975 . C <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,120,7200,0,0,48,1154,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -2 4999976 . A <*> 0 . DP=2;I16=1,1,0,0,72,2594,0,0,120,7200,0,0,49,1201,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,72 -2 4999977 . T <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,120,7200,0,0,50,1250,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -2 4999978 . T <*> 0 . DP=2;I16=1,1,0,0,80,3200,0,0,120,7200,0,0,50,1250,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,80 -2 4999979 . C <*> 0 . DP=2;I16=1,1,0,0,78,3044,0,0,120,7200,0,0,50,1250,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,78 -2 4999980 . T <*> 0 . DP=2;I16=1,1,0,0,81,3281,0,0,120,7200,0,0,50,1250,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,81 -2 4999981 . G <*> 0 . DP=2;I16=1,1,0,0,75,2813,0,0,120,7200,0,0,50,1250,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,75 -2 4999982 . T <*> 0 . DP=2;I16=1,1,0,0,72,2592,0,0,120,7200,0,0,49,1201,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,72 -2 4999983 . G <*> 0 . DP=2;I16=1,1,0,0,75,2813,0,0,120,7200,0,0,48,1154,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,75 -2 4999984 . C <*> 0 . DP=2;I16=1,1,0,0,76,2888,0,0,120,7200,0,0,47,1109,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,76 -2 4999985 . A <*> 0 . DP=2;I16=1,1,0,0,76,2890,0,0,120,7200,0,0,46,1066,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,76 -2 4999986 . C <*> 0 . DP=2;I16=1,1,0,0,76,2890,0,0,120,7200,0,0,45,1025,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,76 -2 4999987 . C <*> 0 . DP=2;I16=1,1,0,0,81,3281,0,0,120,7200,0,0,44,986,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,81 -2 4999988 . T <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,120,7200,0,0,43,949,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -2 4999989 . G <*> 0 . DP=2;I16=1,1,0,0,78,3044,0,0,120,7200,0,0,42,914,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,78 -2 4999990 . C <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,120,7200,0,0,41,881,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -2 4999991 . A <*> 0 . DP=2;I16=1,1,0,0,82,3362,0,0,120,7200,0,0,40,850,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,82 -2 4999992 . G <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,120,7200,0,0,39,821,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -2 4999993 . G <*> 0 . DP=2;I16=1,1,0,0,77,2969,0,0,120,7200,0,0,38,794,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -2 4999994 . T <*> 0 . DP=2;I16=1,1,0,0,76,2890,0,0,120,7200,0,0,37,769,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,76 -2 4999995 . T <*> 0 . DP=2;I16=1,1,0,0,74,2738,0,0,120,7200,0,0,36,746,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,74 -2 4999996 . T <*> 0 . DP=2;I16=1,1,0,0,79,3121,0,0,120,7200,0,0,35,725,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,79 -2 4999997 . A <*> 0 . DP=2;I16=1,1,0,0,76,2888,0,0,120,7200,0,0,34,706,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,76 -2 4999998 . A <*> 0 . DP=2;I16=1,1,0,0,75,2813,0,0,120,7200,0,0,33,689,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,75 -2 4999999 . T C,<*> 0 . DP=2;I16=0,1,1,0,39,1521,38,1444,60,3600,60,3600,25,625,7,49;QS=0.506494,0.493506,0;SGB=-0.379885;RPB=1;MQB=1;MQSB=1;BQB=1;MQ0F=0 PL 32,0,33,35,36,68 -2 5000000 . A <*> 0 . DP=2;I16=1,1,0,0,76,2890,0,0,120,7200,0,0,31,661,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,76 -2 5000001 . T <*> 0 . DP=2;I16=1,1,0,0,76,2896,0,0,120,7200,0,0,30,650,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,76 -2 5000002 . T <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,120,7200,0,0,29,641,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -2 5000003 . A <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,120,7200,0,0,28,634,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -2 5000004 . C <*> 0 . DP=2;I16=1,1,0,0,79,3121,0,0,120,7200,0,0,27,629,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,79 -2 5000005 . A <*> 0 . DP=2;I16=1,1,0,0,71,2525,0,0,120,7200,0,0,26,626,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,71 -2 5000006 . T <*> 0 . DP=2;I16=1,1,0,0,70,2458,0,0,120,7200,0,0,25,625,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,70 -2 5000007 . G <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000008 . G <*> 0 . DP=1;I16=0,1,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 5000009 . A <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000010 . A <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000011 . G <*> 0 . DP=1;I16=0,1,0,0,35,1225,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,35 -2 5000012 . C <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000013 . C <*> 0 . DP=1;I16=0,1,0,0,33,1089,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,33 -2 5000014 . A <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000015 . A <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 5000016 . G <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 5000017 . G <*> 0 . DP=1;I16=0,1,0,0,41,1681,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,41 -2 5000018 . C <*> 0 . DP=1;I16=0,1,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 5000019 . T <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000020 . T <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000021 . A <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000022 . A <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000023 . G <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 5000024 . G <*> 0 . DP=1;I16=0,1,0,0,35,1225,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,35 -2 5000025 . C <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000026 . T <*> 0 . DP=2;I16=1,1,0,0,67,2357,0,0,120,7200,0,0,25,625,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,67 -2 5000027 . T <*> 0 . DP=2;I16=1,1,0,0,72,2600,0,0,120,7200,0,0,25,577,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,72 -2 5000028 . G <*> 0 . DP=2;I16=1,1,0,0,72,2610,0,0,120,7200,0,0,25,533,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,72 -2 5000029 . T <*> 0 . DP=2;I16=1,1,0,0,75,2825,0,0,120,7200,0,0,25,493,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,75 -2 5000030 . G <*> 0 . DP=2;I16=1,1,0,0,72,2594,0,0,120,7200,0,0,25,457,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,72 -2 5000031 . C <*> 0 . DP=2;I16=1,1,0,0,75,2813,0,0,120,7200,0,0,25,425,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,75 -2 5000032 . T <*> 0 . DP=2;I16=1,1,0,0,71,2525,0,0,120,7200,0,0,25,397,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,71 -2 5000033 . C <*> 0 . DP=2;I16=1,1,0,0,74,2738,0,0,120,7200,0,0,25,373,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,74 -2 5000034 . T <*> 0 . DP=2;I16=1,1,0,0,75,2813,0,0,120,7200,0,0,25,353,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,75 -2 5000035 . C <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,120,7200,0,0,25,337,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -2 5000036 . C <*> 0 . DP=2;I16=1,1,0,0,70,2452,0,0,120,7200,0,0,25,325,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,70 -2 5000037 . A <*> 0 . DP=2;I16=1,1,0,0,74,2740,0,0,120,7200,0,0,25,317,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,74 -2 5000038 . A <*> 0 . DP=2;I16=1,1,0,0,84,3530,0,0,120,7200,0,0,25,313,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,84 -2 5000039 . A <*> 0 . DP=2;I16=1,1,0,0,76,2888,0,0,120,7200,0,0,25,313,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,76 -2 5000040 . G <*> 0 . DP=2;I16=1,1,0,0,76,2888,0,0,120,7200,0,0,25,317,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,76 -2 5000041 . C <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,120,7200,0,0,25,325,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -2 5000042 . A <*> 0 . DP=2;I16=1,1,0,0,73,2665,0,0,120,7200,0,0,25,337,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,73 -2 5000043 . G <*> 0 . DP=2;I16=1,1,0,0,75,2813,0,0,120,7200,0,0,25,353,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,75 -2 5000044 . C <*> 0 . DP=2;I16=1,1,0,0,77,2969,0,0,120,7200,0,0,25,373,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -2 5000045 . A <*> 0 . DP=2;I16=1,1,0,0,72,2594,0,0,120,7200,0,0,25,397,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,72 -2 5000046 . G <*> 0 . DP=2;I16=1,1,0,0,75,2813,0,0,120,7200,0,0,25,425,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,75 -2 5000047 . T <*> 0 . DP=2;I16=1,1,0,0,73,2665,0,0,120,7200,0,0,25,457,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,73 -2 5000048 . G <*> 0 . DP=2;I16=1,1,0,0,73,2665,0,0,120,7200,0,0,25,493,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,73 -2 5000049 . A <*> 0 . DP=2;I16=1,1,0,0,73,2689,0,0,120,7200,0,0,25,533,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,73 -2 5000050 . T <*> 0 . DP=2;I16=1,1,0,0,71,2545,0,0,120,7200,0,0,25,577,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,71 -2 5000051 . A <*> 0 . DP=2;I16=1,1,0,0,72,2610,0,0,120,7200,0,0,25,625,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,72 -2 5000052 . T <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 5000053 . G <*> 0 . DP=1;I16=1,0,0,0,34,1156,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,34 -2 5000054 . G <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 5000055 . T <*> 0 . DP=1;I16=1,0,0,0,36,1296,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -2 5000056 . T <*> 0 . DP=1;I16=1,0,0,0,32,1024,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,32 -2 5000057 . T <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 5000058 . G <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000059 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000060 . C <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 5000061 . T <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 5000062 . G <*> 0 . DP=1;I16=1,0,0,0,36,1296,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -2 5000063 . T <*> 0 . DP=1;I16=1,0,0,0,21,441,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,21 -2 5000064 . G <*> 0 . DP=1;I16=1,0,0,0,34,1156,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,34 -2 5000065 . T <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 5000066 . C <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000067 . C <*> 0 . DP=1;I16=1,0,0,0,41,1681,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,41 -2 5000068 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000069 . C <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000070 . A <*> 0 . DP=1;I16=1,0,0,0,34,1156,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,34 -2 5000071 . C <*> 0 . DP=1;I16=1,0,0,0,35,1225,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,35 -2 5000072 . C <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000073 . C <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 5000074 . A <*> 0 . DP=1;I16=1,0,0,0,20,400,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,20 -2 5000075 . A <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000076 . A <*> 0 . DP=1;I16=1,0,0,0,35,1225,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,35 -2 5000077 . T <*> 0 . DP=1;I16=1,0,0,0,35,1225,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,35 -2 5000078 . C <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 5000079 . T <*> 0 . DP=1;I16=1,0,0,0,36,1296,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -2 5000080 . C <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 5000081 . A <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000082 . T <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 5000083 . C <*> 0 . DP=1;I16=1,0,0,0,35,1225,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,35 -2 5000084 . T <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000085 . T <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000086 . G <*> 0 . DP=1;I16=1,0,0,0,32,1024,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,32 -2 5000087 . A <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000088 . A <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000089 . T <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000090 . T <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000091 . T <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 5000092 . T <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000093 . A <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 5000094 . G <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000095 . C <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 5000096 . T <*> 0 . DP=1;I16=1,0,0,0,36,1296,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -2 5000097 . G <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000098 . C <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000099 . C <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000100 . A <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 5000101 . T <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,24,576,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000102 . A <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,23,529,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 5000103 . A <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,22,484,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000104 . C <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,21,441,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000105 . T <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,20,400,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 5000106 . C <*> 0 . DP=1;I16=1,0,0,0,36,1296,0,0,60,3600,0,0,19,361,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -2 5000107 . C <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,18,324,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 5000108 . C <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,17,289,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 5000109 . A <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,16,256,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000110 . T <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,15,225,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 5000111 . G <*> 0 . DP=1;I16=1,0,0,0,33,1089,0,0,60,3600,0,0,14,196,0,0;QS=1,0;MQ0F=0 PL 0,3,33 -2 5000112 . T <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,13,169,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000113 . G <*> 0 . DP=1;I16=1,0,0,0,38,1444,0,0,60,3600,0,0,12,144,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -2 5000114 . T <*> 0 . DP=1;I16=1,0,0,0,26,676,0,0,60,3600,0,0,11,121,0,0;QS=1,0;MQ0F=0 PL 0,3,26 -2 5000115 . T <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,10,100,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000116 . G <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,9,81,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000117 . T <*> 0 . DP=1;I16=1,0,0,0,40,1600,0,0,60,3600,0,0,8,64,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -2 5000118 . A <*> 0 . DP=1;I16=1,0,0,0,41,1681,0,0,60,3600,0,0,7,49,0,0;QS=1,0;MQ0F=0 PL 0,3,41 -2 5000119 . G <*> 0 . DP=1;I16=1,0,0,0,30,900,0,0,60,3600,0,0,6,36,0,0;QS=1,0;MQ0F=0 PL 0,3,30 -2 5000120 . A <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,5,25,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000121 . A <*> 0 . DP=1;I16=1,0,0,0,37,1369,0,0,60,3600,0,0,4,16,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -2 5000122 . G <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,3,9,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000123 . G <*> 0 . DP=1;I16=1,0,0,0,39,1521,0,0,60,3600,0,0,2,4,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -2 5000124 . A <*> 0 . DP=1;I16=1,0,0,0,36,1296,0,0,60,3600,0,0,1,1,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -2 5000125 . A <*> 0 . DP=1;I16=1,0,0,0,33,1089,0,0,60,3600,0,0,0,0,0,0;QS=1,0;MQ0F=0 PL 0,3,33 -10 2999980 . G <*> 0 . DP=1;I16=0,1,0,0,36,1296,0,0,60,3600,0,0,0,0,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -10 2999981 . C <*> 0 . DP=2;I16=1,1,0,0,63,2045,0,0,120,7200,0,0,1,1,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,63 -10 2999982 . A <*> 0 . DP=2;I16=1,1,0,0,73,2677,0,0,120,7200,0,0,3,5,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,73 -10 2999983 . A <*> 0 . DP=2;I16=1,1,0,0,69,2393,0,0,120,7200,0,0,5,13,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,69 -10 2999984 . C <*> 0 . DP=2;I16=1,1,0,0,71,2525,0,0,120,7200,0,0,7,25,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,71 -10 2999985 . A <*> 0 . DP=2;I16=1,1,0,0,68,2314,0,0,120,7200,0,0,9,41,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,68 -10 2999986 . A <*> 0 . DP=2;I16=1,1,0,0,73,2669,0,0,120,7200,0,0,11,61,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,73 -10 2999987 . A <*> 0 . DP=2;I16=1,1,0,0,61,1945,0,0,120,7200,0,0,13,85,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,61 -10 2999988 . G <*> 0 . DP=3;I16=1,2,0,0,108,3896,0,0,180,10800,0,0,15,113,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,102 -10 2999989 . T <*> 0 . DP=3;I16=1,2,0,0,112,4190,0,0,180,10800,0,0,18,146,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,106 -10 2999990 . G <*> 0 . DP=3;I16=1,2,0,0,115,4409,0,0,180,10800,0,0,21,185,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,109 -10 2999991 . A <*> 0 . DP=3;I16=1,2,0,0,110,4082,0,0,180,10800,0,0,24,230,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,104 -10 2999992 . G <*> 0 . DP=3;I16=1,2,0,0,117,4569,0,0,180,10800,0,0,27,281,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,111 -10 2999993 . A <*> 0 . DP=3;I16=1,2,0,0,115,4411,0,0,180,10800,0,0,30,338,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,108 -10 2999994 . C <*> 0 . DP=3;I16=1,2,0,0,111,4145,0,0,180,10800,0,0,33,401,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,105 -10 2999995 . C <*> 0 . DP=3;I16=1,2,0,0,108,3944,0,0,180,10800,0,0,36,470,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,102 -10 2999996 . T <*> 0 . DP=3;I16=1,2,0,0,113,4257,0,0,180,10800,0,0,39,545,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,106 -10 2999997 . C <*> 0 . DP=3;I16=1,2,0,0,111,4109,0,0,180,10800,0,0,42,626,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,105 -10 2999998 . A <*> 0 . DP=3;I16=1,2,0,0,111,4121,0,0,180,10800,0,0,45,713,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,105 -10 2999999 . T <*> 0 . DP=3;I16=1,2,0,0,109,3961,0,0,180,10800,0,0,48,806,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,103 -10 3000000 . C <*> 0 . DP=3;I16=1,2,0,0,109,3979,0,0,180,10800,0,0,51,905,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,103 -10 3000001 . T <*> 0 . DP=3;I16=1,2,0,0,109,3969,0,0,180,10800,0,0,54,1010,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,103 -10 3000002 . C <*> 0 . DP=3;I16=1,2,0,0,109,3961,0,0,180,10800,0,0,57,1121,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,103 -10 3000003 . T <*> 0 . DP=3;I16=1,2,0,0,116,4494,0,0,180,10800,0,0,60,1238,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,110 -10 3000004 . A <*> 0 . DP=3;I16=1,2,0,0,118,4642,0,0,180,10800,0,0,63,1361,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,111 -10 3000005 . C <*> 0 . DP=3;I16=1,2,0,0,117,4569,0,0,180,10800,0,0,66,1490,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,110 -10 3000006 . A <*> 0 . DP=3;I16=1,2,0,0,112,4190,0,0,180,10800,0,0,68,1574,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,106 -10 3000007 . A <*> 0 . DP=3;I16=1,2,0,0,111,4113,0,0,180,10800,0,0,69,1611,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,105 -10 3000008 . A <*> 0 . DP=3;I16=1,2,0,0,115,4417,0,0,180,10800,0,0,70,1650,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,108 -10 3000009 . A <*> 0 . DP=3;I16=1,2,0,0,115,4411,0,0,180,10800,0,0,71,1691,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,109 -10 3000010 . A <*> 0 . DP=3;I16=1,2,0,0,117,4569,0,0,180,10800,0,0,72,1734,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,110 -10 3000011 . A <*> 0 . DP=3;I16=1,2,0,0,116,4494,0,0,180,10800,0,0,73,1779,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,109 -10 3000012 . T <*> 0 . DP=3;I16=1,2,0,0,112,4190,0,0,180,10800,0,0,74,1826,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,106 -10 3000013 . G <*> 0 . DP=3;I16=1,2,0,0,114,4334,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,108 -10 3000014 . A <*> 0 . DP=3;I16=1,2,0,0,117,4577,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,110 -10 3000015 . A <*> 0 . DP=3;I16=1,2,0,0,112,4182,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,105 -10 3000016 . A <*> 0 . DP=3;I16=1,2,0,0,116,4490,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,110 -10 3000017 . A <*> 0 . DP=3;I16=1,2,0,0,113,4259,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,106 -10 3000018 . G <*> 0 . DP=3;I16=1,2,0,0,110,4042,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,104 -10 3000019 . G <*> 0 . DP=3;I16=1,2,0,0,116,4488,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,109 -10 3000020 . T <*> 0 . DP=3;I16=1,2,0,0,115,4409,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,109 -10 3000021 . T <*> 0 . DP=3;I16=1,2,0,0,110,4046,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,104 -10 3000022 . A <*> 0 . DP=3;I16=1,2,0,0,105,3713,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,99 -10 3000023 . G <*> 0 . DP=3;I16=1,2,0,0,114,4334,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,107 -10 3000024 . C <*> 0 . DP=3;I16=1,2,0,0,113,4267,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,107 -10 3000025 . T <*> 0 . DP=3;I16=1,2,0,0,112,4182,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,106 -10 3000026 . T <*> 0 . DP=3;I16=1,2,0,0,119,4725,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,112 -10 3000027 . G <*> 0 . DP=3;I16=1,2,0,0,96,3464,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,93 -10 3000028 . G <*> 0 . DP=3;I16=1,2,0,0,113,4277,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,107 -10 3000029 . T <*> 0 . DP=4;I16=1,2,0,0,116,4490,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,110 -10 3000030 . G <*> 0 . DP=4;I16=1,3,0,0,152,5778,0,0,240,14400,0,0,76,1876,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,134 -10 3000031 . T <*> 0 . DP=4;I16=1,3,0,0,140,5100,0,0,240,14400,0,0,77,1879,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,122 -10 3000032 . A <*> 0 . DP=4;I16=1,3,0,0,137,4965,0,0,240,14400,0,0,78,1884,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,123 -10 3000033 . G <*> 0 . DP=4;I16=1,3,0,0,153,5853,0,0,240,14400,0,0,79,1891,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,135 -10 3000034 . T <*> 0 . DP=4;I16=1,3,0,0,140,4998,0,0,240,14400,0,0,80,1900,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,124 -10 3000035 . G <*> 0 . DP=4;I16=1,3,0,0,152,5794,0,0,240,14400,0,0,81,1911,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,133 -10 3000036 . G <*> 0 . DP=4;I16=1,3,0,0,150,5628,0,0,240,14400,0,0,82,1924,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,132 -10 3000037 . C <*> 0 . DP=4;I16=1,3,0,0,157,6165,0,0,240,14400,0,0,83,1939,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,138 -10 3000038 . A <*> 0 . DP=4;I16=1,2,0,0,110,4034,0,0,180,10800,0,0,75,1875,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,103 -10 3000039 . C <*> 0 . DP=4;I16=1,3,0,0,149,5553,0,0,240,14400,0,0,85,1975,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,131 -10 3000040 . A <*> 0 . DP=4;I16=1,3,0,0,137,4767,0,0,240,14400,0,0,86,1996,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,121 -10 3000041 . T <*> 0 . DP=4;I16=1,3,0,0,148,5514,0,0,240,14400,0,0,87,2019,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,130 -10 3000042 . A <*> 0 . DP=4;I16=1,3,0,0,135,4599,0,0,240,14400,0,0,88,2044,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,119 -10 3000043 . T <*> 0 . DP=4;I16=1,3,0,0,155,6011,0,0,240,14400,0,0,89,2071,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,136 -10 3000044 . C <*> 0 . DP=4;I16=1,3,0,0,150,5646,0,0,240,14400,0,0,90,2100,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,132 -10 3000045 . T <*> 0 . DP=4;I16=1,3,0,0,133,4559,0,0,240,14400,0,0,91,2131,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,118 -10 3000046 . G <*> 0 . DP=4;I16=1,3,0,0,146,5362,0,0,240,14400,0,0,92,2164,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,128 -10 3000047 . T <*> 0 . DP=4;I16=1,3,0,0,150,5630,0,0,240,14400,0,0,93,2199,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,131 -10 3000048 . G <*> 0 . DP=4;I16=1,3,0,0,155,6007,0,0,240,14400,0,0,94,2236,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,137 -10 3000049 . G <*> 0 . DP=4;I16=1,3,0,0,152,5814,0,0,240,14400,0,0,95,2275,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,135 -10 3000050 . T <*> 0 . DP=4;I16=1,3,0,0,147,5465,0,0,240,14400,0,0,96,2316,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,130 -10 3000051 . C <*> 0 . DP=4;I16=1,3,0,0,152,5782,0,0,240,14400,0,0,97,2359,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,133 -10 3000052 . C <*> 0 . DP=4;I16=1,3,0,0,153,5861,0,0,240,14400,0,0,98,2404,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,134 -10 3000053 . C <*> 0 . DP=4;I16=1,3,0,0,149,5595,0,0,240,14400,0,0,99,2451,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,131 -10 3000054 . A <*> 0 . DP=4;I16=1,3,0,0,144,5202,0,0,240,14400,0,0,100,2500,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,127 -10 3000055 . G <*> 0 . DP=4;I16=1,3,0,0,153,5859,0,0,240,14400,0,0,99,2451,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,135 -10 3000056 . C <*> 0 . DP=4;I16=1,3,0,0,157,6163,0,0,240,14400,0,0,97,2355,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,139 -10 3000057 . T <*> 0 . DP=5;I16=2,3,0,0,166,5624,0,0,300,18000,0,0,95,2263,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,145 -10 3000058 . A <*> 0 . DP=5;I16=2,3,0,0,187,7011,0,0,300,18000,0,0,94,2176,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,162 -10 3000059 . T <*> 0 . DP=5;I16=2,3,0,0,177,6285,0,0,300,18000,0,0,93,2095,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,154 -10 3000060 . T <*> 0 . DP=5;I16=2,3,0,0,182,6644,0,0,300,18000,0,0,92,2020,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,159 -10 3000061 . T <*> 0 . DP=5;I16=2,3,0,0,187,7011,0,0,300,18000,0,0,91,1951,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,162 -10 3000062 . G <*> 0 . DP=5;I16=2,3,0,0,188,7084,0,0,300,18000,0,0,90,1888,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,164 -10 3000063 . G <*> 0 . DP=5;I16=2,3,0,0,184,6802,0,0,300,18000,0,0,88,1782,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,160 -10 3000064 . G <*> 0 . DP=5;I16=2,3,0,0,187,6999,0,0,300,18000,0,0,86,1684,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,162 -10 3000065 . A <*> 0 . DP=5;I16=2,3,0,0,185,6865,0,0,300,18000,0,0,84,1594,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,161 -10 3000066 . G A,<*> 0 . DP=5;I16=1,1,1,2,71,2525,112,4186,120,7200,180,10800,39,821,43,691;QS=0.387978,0.612022,0;VDB=0.946712;SGB=-0.511536;RPB=1;MQB=1;MQSB=1;BQB=0.5;MQ0F=0 PL 91,0,56,97,65,152 -10 3000067 . G <*> 0 . DP=5;I16=2,3,0,0,190,7224,0,0,300,18000,0,0,80,1438,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,164 -10 3000068 . C <*> 0 . DP=5;I16=2,3,0,0,185,6877,0,0,300,18000,0,0,78,1372,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,161 -10 3000069 . T <*> 0 . DP=5;I16=2,3,0,0,183,6715,0,0,300,18000,0,0,76,1314,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,158 -10 3000070 . G <*> 0 . DP=5;I16=2,3,0,0,186,6922,0,0,300,18000,0,0,74,1264,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,161 -10 3000071 . A <*> 0 . DP=5;I16=2,3,0,0,183,6707,0,0,300,18000,0,0,72,1222,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,159 -10 3000072 . G <*> 0 . DP=5;I16=2,3,0,0,188,7080,0,0,300,18000,0,0,70,1188,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,164 -10 3000073 . A <*> 0 . DP=5;I16=2,3,0,0,179,6411,0,0,300,18000,0,0,68,1162,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,155 -10 3000074 . T <*> 0 . DP=5;I16=2,3,0,0,175,6187,0,0,300,18000,0,0,66,1144,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,153 -10 3000075 . A <*> 0 . DP=5;I16=2,3,0,0,172,5958,0,0,300,18000,0,0,64,1134,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,149 -10 3000076 . G <*> 0 . DP=5;I16=2,3,0,0,184,6776,0,0,300,18000,0,0,62,1132,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,159 -10 3000077 . G <*> 0 . DP=5;I16=2,3,0,0,182,6654,0,0,300,18000,0,0,60,1138,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,158 -10 3000078 . A <*> 0 . DP=5;I16=2,3,0,0,184,6784,0,0,300,18000,0,0,58,1152,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,159 -10 3000079 . G <*> 0 . DP=5;I16=2,3,0,0,181,6697,0,0,300,18000,0,0,56,1174,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,15,158 -10 3000080 . G <*> 0 . DP=4;I16=2,2,0,0,147,5435,0,0,240,14400,0,0,55,1203,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,12,135 -10 3000081 . A <*> 0 . DP=3;I16=1,2,0,0,114,4334,0,0,180,10800,0,0,55,1237,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,107 -10 3000082 . T <*> 0 . DP=3;I16=1,2,0,0,107,3841,0,0,180,10800,0,0,55,1275,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,101 -10 3000083 . C <*> 0 . DP=3;I16=1,2,0,0,111,4109,0,0,180,10800,0,0,54,1266,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,105 -10 3000084 . A <*> 0 . DP=3;I16=1,2,0,0,94,3054,0,0,180,10800,0,0,53,1259,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,88 -10 3000085 . C <*> 0 . DP=3;I16=1,2,0,0,112,4210,0,0,180,10800,0,0,52,1254,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,106 -10 3000086 . T <*> 0 . DP=3;I16=1,2,0,0,111,4145,0,0,180,10800,0,0,51,1251,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,105 -10 3000087 . T <*> 0 . DP=3;I16=1,2,0,0,112,4214,0,0,180,10800,0,0,50,1250,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,106 -10 3000088 . G <*> 0 . DP=2;I16=1,1,0,0,70,2450,0,0,120,7200,0,0,50,1250,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,70 -10 3000089 . A <*> 0 . DP=2;I16=1,1,0,0,75,2813,0,0,120,7200,0,0,50,1250,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,75 -10 3000090 . G <*> 0 . DP=2;I16=1,1,0,0,69,2385,0,0,120,7200,0,0,50,1250,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,69 -10 3000091 . C <*> 0 . DP=3;I16=1,2,0,0,109,4001,0,0,180,10800,0,0,50,1250,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,103 -10 3000092 . C <*> 0 . DP=3;I16=1,2,0,0,110,4058,0,0,180,10800,0,0,51,1251,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,104 -10 3000093 . C <*> 0 . DP=3;I16=1,2,0,0,107,3817,0,0,180,10800,0,0,52,1254,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,101 -10 3000094 . A <*> 0 . DP=3;I16=1,2,0,0,92,3026,0,0,180,10800,0,0,53,1259,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,88 -10 3000095 . A <*> 0 . DP=3;I16=1,2,0,0,102,3518,0,0,180,10800,0,0,54,1266,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,96 -10 3000096 . G <*> 0 . DP=3;I16=1,2,0,0,114,4332,0,0,180,10800,0,0,55,1275,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,108 -10 3000097 . A <*> 0 . DP=3;I16=1,2,0,0,110,4038,0,0,180,10800,0,0,56,1286,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,104 -10 3000098 . G <*> 0 . DP=3;I16=1,2,0,0,111,4109,0,0,180,10800,0,0,57,1299,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,105 -10 3000099 . G <*> 0 . DP=3;I16=1,2,0,0,113,4261,0,0,180,10800,0,0,58,1314,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,106 -10 3000100 . T <*> 0 . DP=3;I16=1,2,0,0,112,4230,0,0,180,10800,0,0,59,1331,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,106 -10 3000101 . C <*> 0 . DP=3;I16=1,2,0,0,111,4109,0,0,180,10800,0,0,60,1350,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,104 -10 3000102 . A <*> 0 . DP=3;I16=1,2,0,0,116,4488,0,0,180,10800,0,0,61,1371,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,110 -10 3000103 . A <*> 0 . DP=3;I16=1,2,0,0,114,4332,0,0,180,10800,0,0,62,1394,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,108 -10 3000104 . G <*> 0 . DP=3;I16=1,2,0,0,117,4565,0,0,180,10800,0,0,62,1370,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,111 -10 3000105 . G <*> 0 . DP=3;I16=1,2,0,0,118,4650,0,0,180,10800,0,0,62,1350,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,111 -10 3000106 . C <*> 0 . DP=3;I16=1,2,0,0,114,4370,0,0,180,10800,0,0,62,1334,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,107 -10 3000107 . T <*> 0 . DP=3;I16=1,2,0,0,111,4109,0,0,180,10800,0,0,62,1322,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,105 -10 3000108 . G <*> 0 . DP=3;I16=1,2,0,0,119,4721,0,0,180,10800,0,0,62,1314,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,112 -10 3000109 . C <*> 0 . DP=3;I16=1,2,0,0,110,4054,0,0,180,10800,0,0,62,1310,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,104 -10 3000110 . A <*> 0 . DP=3;I16=1,2,0,0,107,3821,0,0,180,10800,0,0,62,1310,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,101 -10 3000111 . G <*> 0 . DP=3;I16=1,2,0,0,121,4893,0,0,180,10800,0,0,62,1314,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,114 -10 3000112 . T <*> 0 . DP=3;I16=1,2,0,0,110,4034,0,0,180,10800,0,0,62,1322,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,104 -10 3000113 . G <*> 0 . DP=3;I16=1,2,0,0,112,4186,0,0,180,10800,0,0,62,1334,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,105 -10 3000114 . A <*> 0 . DP=3;I16=1,2,0,0,115,4409,0,0,180,10800,0,0,62,1350,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,108 -10 3000115 . G <*> 0 . DP=3;I16=1,2,0,0,113,4261,0,0,180,10800,0,0,62,1370,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,107 -10 3000116 . C <*> 0 . DP=3;I16=1,2,0,0,115,4417,0,0,180,10800,0,0,62,1394,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,108 -10 3000117 . C <*> 0 . DP=3;I16=1,2,0,0,115,4409,0,0,180,10800,0,0,61,1371,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,109 -10 3000118 . A <*> 0 . DP=3;I16=1,2,0,0,118,4642,0,0,180,10800,0,0,60,1350,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,112 -10 3000119 . T <*> 0 . DP=3;I16=1,2,0,0,108,3890,0,0,180,10800,0,0,59,1331,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,102 -10 3000120 . G <*> 0 . DP=3;I16=1,2,0,0,113,4259,0,0,180,10800,0,0,58,1314,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,107 -10 3000121 . A <*> 0 . DP=3;I16=1,2,0,0,106,3750,0,0,180,10800,0,0,57,1299,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,100 -10 3000122 . T <*> 0 . DP=3;I16=1,2,0,0,111,4109,0,0,180,10800,0,0,56,1286,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,105 -10 3000123 . T <*> 0 . DP=3;I16=1,2,0,0,115,4411,0,0,180,10800,0,0,55,1275,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,108 -10 3000124 . G <*> 0 . DP=3;I16=1,2,0,0,108,3890,0,0,180,10800,0,0,54,1266,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,101 -10 3000125 . C <*> 0 . DP=3;I16=1,2,0,0,113,4275,0,0,180,10800,0,0,53,1259,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,106 -10 3000126 . A <*> 0 . DP=3;I16=1,2,0,0,109,3977,0,0,180,10800,0,0,52,1254,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,103 -10 3000127 . T <*> 0 . DP=3;I16=1,2,0,0,111,4121,0,0,180,10800,0,0,51,1251,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,105 -10 3000128 . C <*> 0 . DP=3;I16=1,2,0,0,101,3489,0,0,180,10800,0,0,50,1250,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,9,96 -10 3000129 . C <*> 0 . DP=2;I16=1,1,0,0,78,3042,0,0,120,7200,0,0,50,1250,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,78 -10 3000130 . C <*> 0 . DP=2;I16=1,1,0,0,76,2888,0,0,120,7200,0,0,50,1250,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,76 -10 3000131 . T <*> 0 . DP=2;I16=1,1,0,0,76,2890,0,0,120,7200,0,0,50,1250,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,76 -10 3000132 . G <*> 0 . DP=2;I16=1,1,0,0,80,3202,0,0,120,7200,0,0,49,1201,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,80 -10 3000133 . C <*> 0 . DP=2;I16=1,1,0,0,74,2738,0,0,120,7200,0,0,48,1154,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,74 -10 3000134 . A <*> 0 . DP=2;I16=1,1,0,0,75,2817,0,0,120,7200,0,0,47,1109,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,75 -10 3000135 . C <*> 0 . DP=2;I16=1,1,0,0,74,2738,0,0,120,7200,0,0,46,1066,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,74 -10 3000136 . T <*> 0 . DP=2;I16=1,1,0,0,74,2756,0,0,120,7200,0,0,45,1025,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,74 -10 3000137 . C <*> 0 . DP=2;I16=1,1,0,0,75,2817,0,0,120,7200,0,0,44,986,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,75 -10 3000138 . C <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,120,7200,0,0,43,949,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -10 3000139 . A <*> 0 . DP=2;I16=1,1,0,0,78,3044,0,0,120,7200,0,0,42,914,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,78 -10 3000140 . A <*> 0 . DP=2;I16=1,1,0,0,73,2677,0,0,120,7200,0,0,41,881,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,73 -10 3000141 . C <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,120,7200,0,0,40,850,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -10 3000142 . C <*> 0 . DP=2;I16=1,1,0,0,81,3281,0,0,120,7200,0,0,39,821,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,81 -10 3000143 . T <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,120,7200,0,0,38,794,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -10 3000144 . G <*> 0 . DP=2;I16=1,1,0,0,80,3208,0,0,120,7200,0,0,37,769,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,80 -10 3000145 . G <*> 0 . DP=2;I16=1,1,0,0,78,3044,0,0,120,7200,0,0,36,746,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,78 -10 3000146 . G <*> 0 . DP=2;I16=1,1,0,0,78,3044,0,0,120,7200,0,0,35,725,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,78 -10 3000147 . T <*> 0 . DP=2;I16=1,1,0,0,73,2677,0,0,120,7200,0,0,34,706,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,73 -10 3000148 . G <*> 0 . DP=2;I16=1,1,0,0,79,3121,0,0,120,7200,0,0,33,689,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,79 -10 3000149 . A <*> 0 . DP=2;I16=1,1,0,0,71,2545,0,0,120,7200,0,0,32,674,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,71 -10 3000150 . C <*> 0 . DP=2;I16=1,1,0,0,74,2746,0,0,120,7200,0,0,31,661,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,74 -10 3000151 . A <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,120,7200,0,0,30,650,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -10 3000152 . G <*> 0 . DP=2;I16=1,1,0,0,77,2965,0,0,120,7200,0,0,29,641,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,77 -10 3000153 . A <*> 0 . DP=2;I16=1,1,0,0,72,2594,0,0,120,7200,0,0,28,634,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,72 -10 3000154 . G <*> 0 . DP=2;I16=1,1,0,0,75,2817,0,0,120,7200,0,0,27,629,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,75 -10 3000155 . C <*> 0 . DP=2;I16=1,1,0,0,56,1730,0,0,120,7200,0,0,26,626,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,56 -10 3000156 . A <*> 0 . DP=2;I16=1,1,0,0,67,2285,0,0,120,7200,0,0,25,625,0,0;QS=1,0;MQSB=1;MQ0F=0 PL 0,6,67 -10 3000157 . A <*> 0 . DP=1;I16=0,1,0,0,31,961,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,31 -10 3000158 . G <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -10 3000159 . A <*> 0 . DP=1;I16=0,1,0,0,29,841,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -10 3000160 . C <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -10 3000161 . C <*> 0 . DP=1;I16=0,1,0,0,36,1296,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -10 3000162 . C <*> 0 . DP=1;I16=0,1,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -10 3000163 . T <*> 0 . DP=1;I16=0,1,0,0,36,1296,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -10 3000164 . G <*> 0 . DP=1;I16=0,1,0,0,40,1600,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -10 3000165 . T <*> 0 . DP=1;I16=0,1,0,0,36,1296,0,0,60,3600,0,0,25,625,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -10 3000166 . T <*> 0 . DP=1;I16=0,1,0,0,40,1600,0,0,60,3600,0,0,24,576,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -10 3000167 . A <*> 0 . DP=1;I16=0,1,0,0,36,1296,0,0,60,3600,0,0,23,529,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -10 3000168 . A <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,60,3600,0,0,22,484,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -10 3000169 . T <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,60,3600,0,0,21,441,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -10 3000170 . A <*> 0 . DP=1;I16=0,1,0,0,34,1156,0,0,60,3600,0,0,20,400,0,0;QS=1,0;MQ0F=0 PL 0,3,34 -10 3000171 . T <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,60,3600,0,0,19,361,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -10 3000172 . T <*> 0 . DP=1;I16=0,1,0,0,40,1600,0,0,60,3600,0,0,18,324,0,0;QS=1,0;MQ0F=0 PL 0,3,40 -10 3000173 . T <*> 0 . DP=1;I16=0,1,0,0,35,1225,0,0,60,3600,0,0,17,289,0,0;QS=1,0;MQ0F=0 PL 0,3,35 -10 3000174 . A <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,60,3600,0,0,16,256,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -10 3000175 . A <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,60,3600,0,0,15,225,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -10 3000176 . C <*> 0 . DP=1;I16=0,1,0,0,36,1296,0,0,60,3600,0,0,14,196,0,0;QS=1,0;MQ0F=0 PL 0,3,36 -10 3000177 . A <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,60,3600,0,0,13,169,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -10 3000178 . T <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,60,3600,0,0,12,144,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -10 3000179 . T <*> 0 . DP=1;I16=0,1,0,0,38,1444,0,0,60,3600,0,0,11,121,0,0;QS=1,0;MQ0F=0 PL 0,3,38 -10 3000180 . A <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,60,3600,0,0,10,100,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -10 3000181 . T <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,60,3600,0,0,9,81,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -10 3000182 . C <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,60,3600,0,0,8,64,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -10 3000183 . T <*> 0 . DP=1;I16=0,1,0,0,39,1521,0,0,60,3600,0,0,7,49,0,0;QS=1,0;MQ0F=0 PL 0,3,39 -10 3000184 . A <*> 0 . DP=1;I16=0,1,0,0,35,1225,0,0,60,3600,0,0,6,36,0,0;QS=1,0;MQ0F=0 PL 0,3,35 -10 3000185 . C <*> 0 . DP=1;I16=0,1,0,0,37,1369,0,0,60,3600,0,0,5,25,0,0;QS=1,0;MQ0F=0 PL 0,3,37 -10 3000186 . C <*> 0 . DP=1;I16=0,1,0,0,35,1225,0,0,60,3600,0,0,4,16,0,0;QS=1,0;MQ0F=0 PL 0,3,35 -10 3000187 . A <*> 0 . DP=1;I16=0,1,0,0,35,1225,0,0,60,3600,0,0,3,9,0,0;QS=1,0;MQ0F=0 PL 0,3,35 -10 3000188 . G <*> 0 . DP=1;I16=0,1,0,0,33,1089,0,0,60,3600,0,0,2,4,0,0;QS=1,0;MQ0F=0 PL 0,3,33 -10 3000189 . A <*> 0 . DP=1;I16=0,1,0,0,29,841,0,0,60,3600,0,0,1,1,0,0;QS=1,0;MQ0F=0 PL 0,3,29 -10 3000190 . A <*> 0 . DP=1;I16=0,1,0,0,26,676,0,0,60,3600,0,0,0,0,0,0;QS=1,0;MQ0F=0 PL 0,3,26 diff --git a/test/test.vcf.gz b/test/test.vcf.gz new file mode 100644 index 0000000000000000000000000000000000000000..9438b0c3de77759dde856c75801735cdc18f38cd GIT binary patch literal 23728 zcmZsBWmHt}7w!xU-Cfe%k}3!cB}kW~fRaOtbmK@VBB+!gFm(3-N((3=0z-ExNQX!W z2;6u4{n!0)*Y$&I=A3=@yW`o<-up1@@q{o4_TLpgL`4+`@{=2FZXNLbCUsbjt)cTl z&3wyU)5}I@`IgwbDlzp@Je@@WKY2B#@FrZiKkzbj{XM%?S9|)HCOG(P0Df=2ynICT zZGF99|0AK5cbnyFTf1Ft=!l#0TgPWVW|zlf7Y4@dou%bEl0tWsiUx&T${qqVEf%Lh(KrJlB53Vj3;K__J0#{76JF(RiE8)AQWW_ih|LJ!e|f zU;CzS=JaxS(f1_==yj^$W^a?0air3Ti5#iY<6%d@xFe?1_UY9qc_?@6 z>iy<$MiXDQn3JqB7kQ-|isz@W)f7KPo7KSa0oXV}nAe$)iquzokteS{3l@|-jJdyl z?7(pwi}J3{2} zENggLhNdy$hZ6N0R{2}&tsT^Ldy*Bm+@<*fOlLGcvv%&-ncESY-N##>EuZw5fAJGe zJC)hgdCsL|uH0eMIbXXYeP8NWo^NeGT}wyIL`O=8!#?M_eJSN#zfd1SmNn|2H;BjoI zll#Uc&@}COkDyLTmmfu4*#6dw6wlM~Kl4q$7~K9Q@Cl|6_25&#c|F0F>BUrXg~R)E zoPqeKFP@+^hYn?sl=4aE9T~<_GwqYm!+h_%J_+=r&1Nr7l(|w{gY_=!tq%Pk`Fy3m zW4ZFv(*r=O5ccq)d$k9-DeZ?B!juZ1!kr3Kd{eDmxW zSh9TN_>w&FhMh#hcvy^UI?2-(TGS&x?>k3>hjX}lzn-b-Fy>g3ywh8}v&aK?mo)$O z0@nzy#CKZYuHOUF>d!$<8h6{hF^mShE`gSba#g6fRl^alL6io&!947fD@0;3qwI6i^iEirhBVCYiAct z%sCdlMB(SorV^C{r%}?)&EJDxXrE*>v(FC?Re%;M9)Z89>WjqtB0Bo_(J`9{SRa1Y<-Kt zZMJgvZ-2uxv1pSSom;=D5fc_20Zn3Ctp^8R<9UtlO_VEovNxXD;pJg?cOG*))BQ5B zTw^+Tt5@xQ!yw7B`v#MGPij}aTp|Q-S znep)Z%g4vT`6~~_%ZIRQ{o>sw@QB>Li|_N12>&7L)pb5syLXspbDhMN(ma_Q{BcNA zDp5C@l~n>L;FEV*bGW z{bufuU{Q&Lxm^2@CtH@H5)Wxq-*wh@H_9c<*{oYV@oLzPn=jdEsoJs($54;P?VLfM zF1c~5BT0yxTwG%HxqiiTBT1|>ru}8 z$tREA&ND?t^L@>Ukvlb>2mOyqL(5>^x(nH~xXnVCC^5~V;%o1?<)|K?Lh{kRNm;9iZS zQu;I_@qk5{bC+df!2Dz3B@0xFPufwUNS~?c`q@<`q=@a7{O#kQZ%m^k@o%nAV@f;T zl9?vZ%HeQK6(x^2;q>5WYk!WMvY^l8oG0oc=GD47iho2osls}=s&uty8z*|P7BEu1 ztTV=*_gI#u!OKw4jIWtJk4?uUr}=sQ-(@TRu>NKenv~d*b!xug?61vI!x#1sN^ugy z@A1A#E-s5PLV8XANo}~P-rAYN9H%BR97`Ph_xf;gci@O^Wh{<%DdKV;sjHTwhbrn) zyXdVR`S?6mK<9c##Jz|=B;!XCjzXnnLriZrbPjJ-_)k@zsooQ|Hu>|8JshEPSn2^A zJ?_z`iz`tGZ=l^>Z+HB?G1tr>)g44AKj+q_b|U2|Ft^dBZl=I$`1)d#;$*Np{Ajv% zrYro&<+@e;vrBU-^32BbUUUPk}mMUvG~|tB5r`T zq_{-1s8h=6>$Hhk78Bn((5pciae;^*g4B`agOtU)Zk>In;|CQy5o2a_kyPu!^eR z-HAPVNnh^pDV#4?yF}r*+uvU$p)Y&w#ezBk7h7E(`}d(Vma8wg+)=?lrC>jwQ7B#= z^{?feDd-L@*9%6#Al|-aOWN`ISILW)vsY&wiyz)}Z+^GB{F^)(Mq6p?b$VoMW9kv8 zzr@=iemeIj>awZhvU?<4-12X-S!^tw_-XTs<3``mLGxv)(&x2zS@G+8N)Zv_^A#<} z#;q`Za`Dr!(!J&(yKZr_fFV6qXA#>2@63Mkitin(;urRg_*5NtG{#@<4jrsDriz6! z{BW#nY4J4G32SSKZ_4mG4eMWf&2`YESXW^lQQucBG`E;#x%T&PP|i*zW^K31QfO=U z!%Fw+wY>1T1wIzXhsU#vSN_`-Z4qVS50N+(VHe?+Bz(xZ9nlWY03Qp3oa$d6&NH`Q z!~DCKF|6WT5pT^D=*Bq~&TgF4Mmz1DpKP67H0Vz4EnYqvOIgnePfV0AJ832FEl;>I zGWl~b`N#H$y*Db%wj<}{Pdx3+Kzr5A?LQk{@)WoGL#{3rSOe(g!a}ze|1LJwE-p4* z-uq(Xb?RxG4sToC`!o}l8}>PTdMn#Dc97t@|d0H7c$xY7vZj` zxrbNBlM%_2draf{0rJm^CaqZae$KFmFILQL5u0hC6$KAlE=5nN)cgCba_PNf=+7@g zZg*evfG%__ri$N1=O|R4KIT;lS0THX-ceP$(ky>X-jHBWJ2xb3QHXrCLvHTc81Aig zs~opT$<+Mq9npv&gY9T7+_k2w{j|czbKeSD zy1Lv!i3#e2&l5fA{;}KimsS5;BUKPM;u`;a%K9nfV6Exzhgcu}=T}~*^f~v2&upqs z{Q^z|e23CnZRg3Jq#Z7N_Xwg=UtSiC2i4`CO0={VLwk;I`Qqk&TPF*qk@yet zL~B(>%Y*rSb!6w_6<G{2Pd#ZMbY<#D_?w!9)PH10Ue zp5|Paz`rm0N9(@8XxjXDe0m}4xb_N+8}~C=`$e+3__u}TujqFV_4V%Z;9B0|rD*lO zz~6R#=!U=3t&m^XTA~#|s(nLnqRV<%ZeL}l>96BjSdRh4w)W-Ak4BGQLq)%=NupyS znlNQNsr#IREC}_k?_o6|IqiBM`}RGTt0ZId!Y;N@Fa7-`o5;l&%=PzC?fx%jXTR6Tn3J4*U$UXDh%PU3 zeoN(>Tt&3pxTNrO(kia?oZ~0X&OW|aKdbOfO^^-&nSkq^o-I4iqtEwcw3?(+&J6MI zlV|(k^mR1a^j^9sInOC>FKaYDMNvP^3dgNX_aWC(T&@=@Tr9HE5#qgA2Mk#-RKYBtO{E-#!pP~(F5kR zp@b)<8-uC-9u%?_6IPL1ciSS%DKv1Zw!BXJ&2F?`2pw&Hv=;jk_TXC}?oOw(U}(f8 zzGZf{b{4b}*7zHvxERz-)coO4G@_-<;+?+rpEjYep~e@TdwsUnac)V{q7wUsM)IS2 zf996S=H5J9Q@#c8tS~!r&BOAZjrIGk6^V4|l2&F!${!j}()Yzkt6aPjtNe2RYV86f zB*-JYYskr3?}J+wLHmP0yDIW?H`TLw8(B+sg}m|6Crbj()SAT zslyyOZR;IL#UL9m__Ry4Dnp-ot%u=y<2I->VaP2{WQk2!3URj};wVc(pW^PBK_PqA z$G`Z?a|(qw>)TiMveIvcQvOQyiti}tY$3HhXxYHsOFz@T6brl40xgb{W|(&`3IeWg6kG%to+XPTKaeWwi_muvt9nKQm5xN@BR}hO6i|i`VRnjSiyu8!DWayQ}rfc|U`=|HZnM(a_%45o7yinEhwT!6B zxy*g7k2fAoLX7WaL7TQ5lZWZI@xOi%MUP~^^YYc&6hD=B*d|lrJS0Cf8~o{fOK0X_ zyfZ5i;qschwpl;c+LZ;t+Mqqdj+bMF3hI!lWy0leL(aBC6<(?$5ugQe9B zFI1N9E^4$lxQduP`$Dm*BjVN1c4fMUf;rYh8C!4nc&Z$Bx{R*faN4lSzlSrlz9=+r zSkYzBw_o|%R%mtao=%>Tyua02r`0#Fl--*__c3S3rY=oH0U~cod78AeDVk6BO-tva z=A3IC_U;R%q=r{0TvUy$R-|xl{T{D+-(Ti!CwzX1=hj&$gWOJB!%e?*Xb z+rg-I%GICqt!KaWPj@2z4qaVFTrLLGoFDpMFi_rjz`mGxg;NLR`Wo2`Jv|)RVb~7| zKcT6y6S~M|xqvws!dYq*6*$lF%=exebe;grOo*TXs1Y4=H)1>REV#V+#bh#92Xd}^wxT9W2S_HVYDk@d*6 z0^2o4wDObXK3m)J@Z1fIXxME=)V$R^;{&qDwcm@78nh_>HaQG0ve4!J!TN8`6zIT6 z9fi#OCSRDbOMteQbJoJZNOKP(W)0o*5S}w?VN%$x#~y$$t?kCcTIjjK`!RC8^Dh2M z2VAY&T>FZkWp)^AWKlJaLg- z>FIBhFN#`Y*7X*_zsuk<+E$qDmgYh!K@=-x8TFUGGDeJ#!gvo@npzlPOM#3Qo)3q!e;Wx_EuYL5~G;#^`TcGUTuQA$SP zy!l6g;1nH?Vfu_h-KR`dW_|u5*r;_KKZlk%Y=YPsvH0v!Qnr=;zMC`J@4_B@I9MQb zbcEY=fWfRWc37$$%he0Z$;Bs*+R^&BHAYD2m{ zeInr+yB&5@%g?ol{D!OsnT~7YG8ryOL!BOD`WQ1QqUfcrS+WgZc_O>s1zzNK%eugm zgo)g+^*0&%x3=EQl$;z%z7@_=M(H$4r%;9Ihao!m{3H{zNFl}7d`(Z}-injEKZxs& zv^lqNpJZEbK6*Y78J>OgtL0lqq+)5>b>~e#$*VtwO}3G~DSxj+db~1jAIiNYw$=#A z(g_m7Y&p!{#Q4;`cmR?6P!I;9uN6x-)}Au(ew4ofb00O?Hx`96F!D1*0%dg;C519m zi_hA-Q~f5Z(4=rwUe}D)fPXf#Y;5gtIu^RL@>)!}M$eRV-6tcpE`PGGZYo1M z?9a4E<=6yK+2S4fc{A7kfk?IveI8u+r7pLzDZx6k35TF?kAwD$fm7 zpLrX8l0t3bk~iLKxPFqwgOghhRdeL?Ben7Vd1-14!pV-}$)%RpQTm&98v(~EZ`n)* z`fUVeooLrLnrc^WhUB`pjM8V6Cf23gP|n>3C=(yo-l#b~Aj+bgCUI;}lz=4bo2$HO zGjQ6-grazuzk*7C5PK@vIidmzhqX=+LhKlOoU89@J@!Z;l*I~J$MAIoK_b*i9t#cA z+nXX?pX?K&$L;i!l>0%OqUJ)f4v)uCxDb=TkNkc^ep1gj>)TX~WzOhbQIF0M>EL|1 z>0F^F_M56!khB-;?__<|O_(qYu7xqcxrc%l$SsFYjblVrZ$#}yK~7x^Q5JGamju1y zD0T$vf{K2h;^pvwILb;6!M7Bq=jG}`FxfsmDDu%yQpOqbd?6FMzWzF41V8P`H`edztsN^yXq{KCT<#_j>twhF`mq#5XH)kc4vCC4QZHxhdmn_d z-I@ub8z!3&)@DF|IN}QeO3cdlC|9S438S4?had?v9nCZf!i;|$h;-U4(AXg&cGMQY zIQtE>w_m9#Hcc0Z`)YPa%GF>iO{MHy?)ExSI7EUqd;< zTO_d0WbHai`ZeQ&lE4<2k7V$J%q6+$QagC${Z1x`P@Wd_kvC>qDE6zvTT+O;tLp~N zxrh~TX>j!uq-H4Zmc^4>S%Z<;VVUbru=f=1W0Y-yuYH$y93`#pT9!oHAc$}aPBKK5 z@xD3~qNANQ>)+rksm1GVZaPX9Hno+d_5IN;-E^!&J06%ZSnSa33;kGx!;wXq33Pn+XMmpX9P5vOYl$cO@-~ zSwcIN&RdzA@S+XG0a}|Fc3@L3(Nd^j{^!somA512pg?|(Su5wd zkGg}}e@PD*CU;vPK?f8|A9+_f!&k+tj|buq+OOM~CmxMetZ15XKWf>UHlv?)@1RTe zjH~Ep;oN^PI<;wE-Ffjk1N2|`+XbvqK7o%nB)8AESHcC{NWb>I`>0XnZ||PwAHO z($V){3d&1}hDG1!PqmRGV&dl z=97_>l*A-o#N7L1;?yu9uG<{AN&1w*1oDlrcCm)XNO+b#sYqLy<1M`t$cg&#RPiPz zoqA<6X@`g!q8p~Q|HtXY9Ty1mxDTWc_Qdk|{TS-&B=$rqnZ{(+x^&bgriASe{xs8@ zp#AoP7smN$)gZnva?(Fn@!HvfLy`WmpB>ol^T|t+-wJ#vFXp%V1i?6XrzkOx5!*WzGFx&}~$WYv= zNfF(_f^3)A(5MPkdK%bmyF&$Bw%R)-Aw6Za`b*l{hBPp?Y5F`vP zr@UmjRYtVQzxBiAX=LW32nY`UwhN7~oHR2-sH2C7Xxv)P6Jq$Aw5P=d4PpO1aD+?{ zg(jwNuHv+JrZ(^+0tD=loi~g%5~rZ3T74pv_m^^2XV(o8NES%R4v`V$6qXwyP)hU| zf_;vWrVRJCN-rk`Bw-VpnB>ht3^$%5S)#HRp@1goQJ(eUzUA05co_ze{N}oBq#cqx zast4A@a75;%CHRwBOSbCYV#+X11r2wC=dKpaYJ5%+iQ1nJGgIvs30n zhzW;cOUP?rMkNsrI|&QJrpn^+HO@WE`_SWzRrV+3k?zP&vkxcP9QzJ~VMt8n{6}tI zXkrc3Yb9LZ;lzDwe2z>G0)UwPIKeyXIBd1YVEq29$Q$W^Y}VM$h`CP>41~d*k+g6{ z<%!&DR*!+hGmttWIo9XiUp}5h7T;F+SU;l#fhJD9VG-OX%sGq7o*GQTS=UK|BJ?Tu z<6Lse*Ojh)<0B00VN8&|Bn9q#ej)&nED7A;(}b*GLFNW3i;vd{_f!P@xOH`8Nu>Ck06(~n03WmA4GSBV) zE(+!3(!h^Z4O@tD)JuYD@_&Ky(4e7cfI+&w7h%MZg#0w%1P+@-)h5C+VODdD2WRyF(KUjUl0yKtB$bya%!|H_blUouybyh&WcPN&a7#wJRO+P- z1Ysf$3@o{Xp`U2z$h*$saAW~n*k~O7ao{w}1Mnb#bge$NI6<%rofrTu4W3CQLB8;5 zg1iHwxIgw17om^2l)^wUH-n)ZSpX|o@chVNA%Pp`;DlNJwb#C4j}GT~0=a@f*#DoO z`8=b^nr|x;?WKxYo@>e-Qx5PuV*hhf@{YfIFY(K0hwB;-qCN)f^`KC;fhScYf$19= zUX?PVZp~-s6IO2ay%A>=SATwOiCN@P1`6LEbU~d^&lxwYAj1=Q2U+|!VJxiZ7fLO-d!Q(|qOG9PA-z3LWwu$8;8#y(5 zlDHk+B1b5BxSBaN)pBC+ITWaQ@>o^tk@P4Jgc$}(8v~TsctvL9-k>_yQ-unUk-;Ad z*!sTAB|VhdT-q2pl#_rEA+e^f4&WnLHtA@AtA@%%2*wO?Ik_S{1k+4p2B&7! zW%G~~yc{YAj*_OOlZN;KrwgOgfiNC|5vlggV;wTWZ4EyzuVq0w5mM&{Jd_tgKt@KN zs=csB;8@1cbBLbekCCy$EeSQF7DC_Q3y>*uHM0pqL{)Jy(#es{xfaSvaYgnsk-j?E zILA|cv(Rfhjl9)Qfr6kxjQS05GGh2p*U>huE`g^U2tyf|kAN1;F;IklMB`7F@ zNee!TfhEV#12HJDkg&!6?G8>I_mQbGtAP|PHgPTj&}fhtHdPXgBpT|{$4RA{2%`8R z4@#`shQu+Df<3DV$!7w&$cP-l;1uvK()mJM`n;fInra;}avUt9k!omrgvmc^I9NwB z(JyWWrpVgw`{}JEDwNA2lPa# ziLi`rN50Jc=i8_Oj+QD7<~kai!-f^C|1Emu_Y4S{u*c2gPs?Q^ldBCIvt-HFto#)&FovpX8@VtS*gsdi}Ly(|PQr z2={+q0-GKTKAws?1eYjKWSZ5ZP&An}-1J}pfTIn`L3Pf|9NCx)Lz4~=W9#Ri1Qz`* zdczQze-ds?%#ksyYrA;=w9@!MJ=Rnvx+2&tDeajdhO})bf8`TPTxY zWU;4Y07JBUT*2|=0;KAAWbkdE38^Nc(1xd}-45N@0M#o4 zfcDj^9A%>npvnwF0-)-ICkpSz;E#doJ^%-<$n8uBmEg|A;21D6)BRT#Fi7TE2 zNjAkKU#g6~K<$$wn2F&o0EVU*m^7&<{2{^&SkIEDCany87H5yukcY%aISN1vL1B?3 zrG}bO@^m>V`rii$8U{tjijk%T053sRpg^r#Ms`v1?^8e2*a1`lvca5rPV ztFHkBq~?MXG)Zb9WU4Jla@2KS=(OJq7M1D%m&Y>6pa%YLQbsoCapWqI|6{Y-M(`;1 z-v5P#i+wX0MsiG808)tMs1lY8Iy8V%3kXmt!$~+>)&$lw1Lo(%|8PYM_`(w)`d@V2 z2M^FBU$hpi@caPZ0*Hd5!j@5-2Lft|Z2+gx0^SZhS~7zhfC7lnG=|Fg$T}n+${u?k zOUk1dQPi+d9QCUh0jh3H3ARk6h5ywO**YU2|6D8hv`~o|wlm4`k$4m2a2=RG0D6k&znVYK&JwQDR4IGEb|^X-^ywQYG{9eu8e;g;2aLCLk9n2D!< z2M@C9uk{WGVJ3Q80tO#%qqp09QtjPdEZpj^wPSuxFz`O+bNnI8F_5y~eDyQ_F!l_5 zY#!O%z2H(B2{N4AI^@fsh}fT_5g&0rk+Lc1(B}sg#!uxfb7$|J<6N%B z3}BSsC>qOybB=l6lXgHHEIA9n4qy3^^FfQx6DqCX(mBn)OEirkYw?HQ7j96M<`lQs z`l46X)U`>xZF%-;4<{DX%-y~&SVp|kSva2ft}fX-V#Ol*Oe8egfP~4Wf zb}!nNM-)HUZB1QwI5Vz&NABo2&&7TdTi)VY3L}Ar0$&;YM0)^9v^}diF%RB4SJ&3x zGr@$tvYA7OzrK6CkcaCEN>z0@(6Ask^p%LpHrY>+>c5afe4o^l^;gI7jW)%<^GYW- zaO#q~->nWvkJ(D6<^oDxIF0k6U5!7Mz%uWa`Y_oJ!$uS6+S&e~3;27xEbh0E`N8KY zw^_nzKejKQIo>^w5f_DRvt)OJzG|Nr)?E!o7!q&j^u}=ea<5+EF`2%+!px2Po0u-VkV-m_-hYNAX zOT)GSGf(;5blQHjG->$k!A3W3*#t8^X$f$10o1)J!tosZllKiF;$Q0}~^ zjzRR&0q?vWne$l9nmO0w+ZuBL!46+CTN2LLs$v7;m!}ZuB&BxlHkFQ`$VN*Hv)wK& z);M<)$eX9F`X}s3Y^kuB)j<~49uU5VCP?Jo%u*sX-23Yp7rb^u6;rvgIJ$Z}6Q+q$FU5K2>sAY|K=!^{+b@Z3!6moASx04zcA7Ysn+s$(B zZ#uu7q&)iM8B;=1$M&q5LtoWGIcEY+0JrP)S>1WbNGI_~nI2Loj%zTw|VO=ddJa)gq=OThM8S;zsYSot&P|;nhBUc1SV7>d zsvoH#dzn*QGkIZd1;?h`7}>|*gz{)qm1;kwAM1#O;Up)MF=pn`xfBn=oDGdrMx9y&?iPS@^Vk=I_+{CiPaB z&YH5E4EtVzU zc75NkWDhXE*6!gN-uZpsk)#xBm$6~%m+Y6=Pu7XB9h@2L7R08|cEojVZx5tg_O{4@ z%{zSl;0rpf1Yd-FC=4{-ZiCh@e4^m7AyU=TykXJ-20sM_1jjHs*ozCT+6~+|qu&Al z+fSu~KoH7Pi;z#Az=1rw2X?YF)2Cq8oT4}i2nAyn2p$(sMIrH@*vl z5Z@d%h%hvUk@0xmz_vLVRSdlA(uBpn4I!$^l4QCErell|&(7h(U2U%piG{oU>u1{O zNrV?M15izUaqP&^Ksz0SDvY~q?s$`$dzF|_kjjrSYwIYI!@}c6E;gHq%?S^RZ zX`#>_$tT5Z9yn~|TuT7%w^Cw^9IPHJtM|Y058lsPB`|13kYaNI#AKzEmel+6*#m$! z+_=lBbzs_;;n%isQ^LUE$6KTPt%wgzY(yE~*4wDS$Kq@ooEb2*Nuwbf>9YT3Y4$GM zXBdKHl4qFenprSoCdE}DKoG#}hg6%2sQ)1bmKu!1PuITOz8B$hzzD}7{ppe&Fxlxk zZO8M13lv~6WYL26V+Rt{PEy-tx(Duv$~DkW@EjPLnuw8VIxNN%;C=}TNN`YnAQJ_H z7}!;nxnffP;?yzX@pOelBsr-3Ra2<`&B?l|V4#I2lx(-$2Ks|&J>b2dH$}gX+?M+o zpUMsgyCy85R)MLLL8`1qXv*X@O43 zqtcP;1Dyd&xy$;$RsErT2AKh4o|B(TT*kPy>R5#*aAbf|rhIhNicW1{(^m&vVgmP5 z3DfJp*=vAJY>Ztdy?#t-i{__@kF6?^z`x}SEbu)eui4?22Veqj=cGEDs*y)|V2>_< zY9=OwX%x%{fC<#r_h@9$*B*j-tJ6qTTVt$!A3IBfJtXzt7GX+3J$$73Z?m+^jY)lp zeHQ@8)fF)f)|U*hp#pSP?6Rq{U}Swk21{}aU$_80EvM;In8#ov#(H`~jZCyU7rV`Y zA>LY<0CU}!B#zl~_xjks1un!tXACpqS;8lSzV8xe`OaUp99in6sOUt>3U^? zadYg@|KB95GCKRV5=pG<#sAHMyu~l`x^rpPV`n~;Q(VRZn<`dnCbEdrHMo58S;oSQ zUmoxjD=dSa>%~CE$Q0nC)XZU0=+T(lOiXr`Nym|CeiYFbi12 zz%uZ93vG!D_R4tUbDF@9?;M*lPw{+GJ94j}KIBdO|GzhJpXlFZStqjIcpS5L;`K+b zhI@jb?>h%E$L#dflU(@VZUgEI_jIHP6Z5Ues|y_F4e6Hh&DMALp^mK-79#Y0b}~cj zKfJuM|DGoA@ zTP8SPDAlm;z2qQ5J?Bus1(uU=sQ*($g#1ht*Uu#IobZkTEmaioACJlqAL$BVLe4Oi z5Ey7AWj*;!q^Sjwff<6pctePQ)r881YP7J^I7XG4a+KQm@^j+yR5Bwmfjb%unBD%= z{Jeh46UTneGD0qbk3+8WX`%!lOAj2ix3VY;$X#)dO;}Ba8*0w$3O@>)0J()Nc+#yB zOLfenB$d6>I>$+T&c{gt+;OZu*H?c#_SpM$?z>%3 zQ@o5C6nGFeT_^LD@>YqdGNruc9o3{jLTkT}-&7xQRmN4xGS;|}a=vO594bjcgf=&3 z_9<^KwIs$#;v?^?nXn)^jw!pddz4vE3fMNul*d)+GA3`U)fDe?roTF?bFAFwIdtrQ zuCHR-SV5N^^^;qLNJWvWzzp^*_QAnJMXvW<*`aEbbN!V;ny4T#E??huh2NC2UV)la zEB&d@6ZLHyP24E$C_dtLrinHv1K$c|d13bm#y=9evZ)kKDK^K#jvd&9YWUfyp1;%= z{-nBI^)ULj8h9jEFv+fa+r#_ZJP=H|`2w3Ym4n-;_dwR!nRvlL6f9 zRQn-qR6%h<$SKuF*`d=$_`V;+t+X`aG2pWtskv=qvJ;u|yeg{-fFBiiWD+;Ho`C|A zkMZm%sG@pQ1*~L`sh|{)ADf&vKg{81QjICV0=Q)zgM#jW0*;Mg@LR93M;^<*3!&CmwxkG3K zG`Rx6%23eVyO{jW6u<^F2;Whk6l%)AK3MRSS8L_aJr6DEo-=S)U!&K6@VZRlG*^sB% zuVRKl6<#20bF<)6fX+)7vy}>t7 zH8`SdR!SJTJ+_?Q9E)L-cpg;J1-vy-4mS8j*<4DN2T}P+6%zMQUs(+UPr^^)D8a`P zCZX^7wH0m^{A+LU3{)|M94mxkY?5O%#t=c|%}7S&$M{ARxmE*9%pl#OfLjzmwJNKb z-29-dCh%NA5b8jf2M|TfrGg5g2F%0^M*@a=q@%3X8K#en<5+*pPDPs5@&cj)E&rq% z^SNLpw~7$(XCR?iz&LoYjIFGu_$Ozyf*-}l!d!azqsJN$&fl(B+W7dTm01lY1G@_> z!O|5nP5{XP5-x`fhA3-HF(Cp=@~22#Q=>NTaFGgfhE)zqR5 zC^(dp@F+n$AXNW(L=C4?rUg=pKs)qzT-ATP%7U~6a5=-0Tj0i~vKlQdmGVFxsS0J# zmk0XSfKSwLtkMcQ2H% zx>8!6X^SRYZybZ{M-r-hTcXWbfAuL5xhqi9ztuzxz^j_XtflQeu>#tdye+Ta?#V z>0Vj8lnPn?U7>Um`CT)HwEA^U4 z&ntEP+WXflL`*V!&|52Og#TPaGR*!HtM!SH_B~IxSj*N=mcKNLu6~`YXNCtmH3 zj(W9+nVdaT&@>lsVVDRq5&f368||U4?)T=D9eLk@t3{&!OO#v)-_E__SuMtr91~;e zH*!7XqnzJLAGJT!9B|;deQWo*yPWElrlWCq*e6*r;j+>ZAN7Wm-jQeg-kzf-UnnSv z-_;Iv3I57`MeI?yD9GJDvf=Ny=5#l)_lI?yF=sXLIXaPVE2nJHOpO`iS$fl3f+O4Q zS5Wu4@?FoJ-EFN43yQZjS5h;XHR9yiMW}tHXIT#M^dn`H8`qt*qdo3TQTXh3M6orr z7Bqa~{MEPv>o}cp@RJaEFD_>tIIOaIH$~6CA*FG~cO*rca!tr*A+@hyv-T@eJ!?l_ z`7%^jKFeUjD5V0eK6}}dpQ!WMqp&vD$NtyH!HXJU!bubJ6Sfb01*bJAld05HEV_s7qOt~v{g(jB<&X*)j~v;QftC4~M&uaxyeSoR=! zf__z}hy-^x{_Ujm^1B8e(HXl5uEz9+INfTh_wk>~P~n<)k@YUccC8W3fwUyDe%GQV zMNy3c4wjo|j~M!`+*=e{-1z1ur?u(Tei+NvMT@OHed}CTC^TGJIVM-T_6QwZ(oqzU zG^7%`LR-`RjZlH_2447eMjB_-Uzc*WD)ycqVHDr}M zF44D;B}%8KZF{GTc{FB9icVKv?cmjpjamP5m)ji{eX-4}MefrXqw`Jqt)lavEgpM~ zZ7|umM)=(b){d4deQh&q_hhNQ##51rmZU|gGoDM-zOu5a)t>6MU*UjlMO7(jfZWRo zm%b>UvaAwU3-_ujzLoKF3nR|z#wmtqxv&<*wer9V|5}$*sHMfYn6yxN>Q^$|V#Mpb zdw%6fQp^6Fj`*>N>s@}Xi8G(hjCg@T;pR%S_MEZXteKTV2(Q6}Eh#J3Ad_iqC zH~pXl zozL%HcwG2FaGcmSdZFCCuyVsGxKOl+Zv6J=*4_<8nT$_m=Pn95<^+~D!+QEU7acsc zr9UB}rDZ#kFt)cfy6iJNJ7b$)Ar+#fMa8_j_T$UtT(j?mZhwn`it6vkqEh}(0A?4N z=%sdL)c}_sKrOFor&0BOO{i0l=fkf{&TlI!D%yPO8+HPp%jw9^lR!sC6`)@k#hE7~?YG&ew2 zJb!S{;7OLIBmDr zd*|taVV1Ie5^7y)ksmPj6B+J!C+hU7NmxQ|NQ{Atj+{PsX zX)du2OJ%m=>nSFg)8L&RFO)Vim6Vl6Mpk=2hEZx^0)nw3S}{%HZ8si((dbL zX{nj>F9B-W7L66_=J!6yyb|Axg|>9D5=Od`iIEnaZSFhBnm5*7hyfjNF9RZlNqC9`DTn|Yzy^Mn8yegMt>q!Xinsz@sg|BI4ua#@svC!_%-(FIKiJC^QPB006*yuRe7qcxwv~{8pXTzUeA-BbG_|fJ(ARx`+Q87 zUuyAVzfVdkM+eqIsn~m)k$5sV{!DTa>%L{mZRVpH@!xv-BHk_yMEx36zM+`fBt`{ea+Z*OlNKYsk?kJopPU%mhR_xJDa9{=+4 z8MaWz?0j6w^wL z=wiCPj6yIFO^DgpK-;HZIfX7ICfE+SY0}0z-_QA))&~(T?=TJ1Xvrdy+nXtFz$Ip7t!*72cBDtsvBp2n7T(1gEu9w$@D$#^QnqVcG z;Gw4Oej?lgBsZ%JLb6I2FLk2--TnUS*UdT8-S4(`-F$Dm`@Pk!i$r|k8TbnNc25>M zz02y&$%4GU&B%(L?qi0d?oBa{z$cEQ(c3~CQJ*;uI

hOb9&rxNhLJChnSRJH z!>N4CIFfE#jw6X0gPDHNjv@)8hM9iUjv}FBn&8%jd>4=_yG)_L&#$$Pl# z*v;$q&6U;q@eA0uGeRbpz`7j=69I^f{x2XdQU;Og*glBssVhOMN|25csHy}iD^UPY zsoZ-F#eYNu03%);zlNC7z54}?*zPx63Q49RCy-nugCv!&VkFm-(uB%uLXw(bc}-A8 zQyG-P@}5MJJ)BH%90du$X-ZO%1`KoRER6^Sn~2RpLX!F-qBKfopg+oXQ#8NFeIEL& zB*k$xdus6AX<7}CAp{Ih01Xz)7}Fo3&31J}_zg>zi0f77NOpUcC4wq+1SvTJD|7@( zJJQXEwaj*30rs#FWu^i0u(=?X85YRHqO++*N^0TR)T%^kJA=@JJ~fF@q@ip4sX>Jr z4ei*Snq#QZK$Wf|MPB z6*(dwR(Le1@-q+xtRv3Lk4WUPj>S!Wf+CNF%IQK%x?nk7(%IM!W@N!~ z4h5#_jCFK46op7<>Z8K}oz(}Z=L0O~!#JEWLSpK#BfR+sr<5buhJF2W&S)Hs>-_w1 zU6(8tK{C1J=b%1#2Xz9{U}Mafu@0xk>>bP*ThfkcPTA7fd!^gJTfW?KN{ zWRD!SP6I^=WhQ=Q+*|Vt%C;&BSXR|;0Nd*DC$8k_XJ!Tr%cNqrcE@7(J*Yo`pBLTx z;lS%$JPlCiV%US73-@t!Xodo@gQ&LmNVU)7ztyx`TSB81240W?rv*OVw$#W7HHv{K zxbIt*b#0;KSpnH>~knP7Mx*o;#{-8E%i)+vCk2{{iZT#$vh1001A0 z2m}BC000301^_}s0sveN?c7_B+r||D;Ai!>j3(zg7ZQ3<{-1~ao3}CA06%|301L3%?lw?wZ#QWL96BYgw{JF5zj(VLiJ`vT zZ=int>Snv!VDsAT%MiGJ{ceNag0Ek{gALU0UZho}djIwPhAnBRzj?QR_4SL*n#sEQ z{VzZK^y&MbH?QB-+c)3d{P}+4G&GCb{YIr(+`hiwC_!4>Z}%HV(v05TZ$gxMd)pia z&BAuyE~G=!&%{9-Y8u%On{Jz(cl}Yj{-|AZ6s3jzFw*b*)#2yyfM%uNC9O1X*RSj+ zTe1^$3m<>X&XiV${XhOV{s?JhKiP6S9KkxgW4%^}FJ{xXU(s2Ko zU?k^5zden%+xPil`|+^b>`mI0ul{l~UU7BZ#&k2>Y@oZ2-~RI9)8U7WLVq>80U54o zH^AIq*?y1EU*B%PIoDm~ZlJ%``RaH==q~qU>ut7mH$ByEdaB*@RJ-Y^cJotpAN^uz zdQkUC&xiJfy=JFx@0P{T?7;TUnGfv?dyS`G-91#>7uAcm)!n=5?W>jj!d76Qg34BQ z1z$jg1}dz$0>g&286pD}RjvZ#jE6tH&%_B6aYD7L9#*^6!-~A;3(TR;905RbtvFwP z_wU2I`|m#dYyZoS8?c^Tz5C(&pMLuO!#}2@i6*bm()9AfUq2lVfBXK!KXxAvfB$=X z(`zpl`q3=*jaSif$JOic-4a~{mScT?`0j6SKmL3Aj;&{ogj_nRm*h-MMgu?(N@a`n<>9@)eEzTMCSKI!2@9O-iwFa|oTDIP)4atKZj zYdXRq+w|-0utq3z$Tq67Kde#9FtX*Sgoib_0+CY~Er@qCVkRgXur;rj(C6%S5 zekmy~CB>zre37JFpp5f#V@c*vQh)5t!(AQI0Fs1 zwnNsoe;Q-$psZ~UoJotInY0L6?W*l+x7w~aL~nh#R9aRF>C1n)wB(btEZ&u+T~Dbc z$*7f%TIZM6Ht&~0#b^-=-jXlJkA$TFvi%CBZ+atN$E^=Nw|)#D5t7Zlw-8+*&A}Cl zkzEchNkq=UMcp_POOlH#VInR~#GOQ%llwXo7i0`J zElLo7Sd)~+O_#)z9@glNi5oPyBr283CCeh3&LCxS$vOd`3(4g&paD`n<=mcER^k)e!6M0>L%PvTJk{P%t$Z;yJ;7nY3sR1g<_ z7vhpYCN2q7+eGd*nd^clE1aAA^398zS20>B=G4ee(#3qUaZ`2IU{#muCUG%<^GAbQ z=cE(YyAfP$a&d~*X#rezSiiZs+djPf=H}an`|Zm&-|k;N+%yQwQef852#l7127$@y z81TM$+XIZ~`Y&0E*IZ9<1lJRcC_~7L&ciZ`w}8G0mu~y)dLUPYBKYIEJPpZwpTBL= ziOyy7=yGT$={&olv{^JePv|9gb{g2aoY^6xXE>MAr;h|Y1kGgy8)I}R1F(1`v9~dL zae1hFHpM3^K*WbN#yP8qfv2!Wn!y^~3|7Z8Si3{U?=%vh$@DE?(5R;q{uF2{cK6PO ziGNGX(ac|@enjez1wRY9Ke>wlEkvxcNtZF1hpGTU z(TBJ!dmgWL4zcJG$a1w4b04p*cH&VS;V6#fY9~Ch+KEpsa^k7dX6$0i(K8R{#ABF? zC$fko6Y~_^^A$N~q2g1EoVb+kw_W7?av_~9T4!h45yWnM;hj#r^Lg!)e5!rkt#+%s zRdB&Zx!5}l=$q+}0(YLV#9d$gB#beinaQf0fT&HOlY$NLc=k=nt)*%y^h~OWMba_* z)?0%yU~=bX&306seOp?^OMh?vgq8TcV-RNLqwz}PKf|?WViT z&bnKkbtbMUm^-~GmwPVWWEN{?H{+j6gJWK8O8%OXz@}g@nB#G$GceDn(~BKa>jvs6 z)fCf|)SdWzJJ;iP8t5u8>-%8l!X+#{$I>V_$)!oMNx#)HWgC0$yCpiQtaamMg5kIJfWe8nTXs~L1}YGXozVRWD%5#g@C~y z*5C_4DW45W!}VFXgR6ql^p8fCa5gwa4W8>OVpFP5#imq$*3k6P$vD-Qgr;;9#T04k z5Q3tYNyyqI={y9cDPl!eNg`B&QE|C=>T>baRd-cN#dj<6QPNb+f;@}?6|bu@hS-wl zJRPqxrkuQ=3(v#T7^n+M`FQ71slu3wn_WaoK4%c!RSX`Ph`Bx$_YgR6eN8HUA%DkH z@!<^F1orQ-#t1e#OuYh=QX|+HB1*B~cs+Hr2oNLq7RspmPKQjgB5hlNxeg**8Q6VU2=_rfInWBkLtKivDnz5iSHBIfQztopWd_R$7k z$ZwemK+Eo>wwQu@;JXAD`UJP=Ug+J+*-P)2*kiEbN)<@-P1>s#7n3jDPt?liAdd2- zo};pSSq(CuQ-)DKmZ!*<%A8KFk?+ay4etHk*^5yiF4&7PaM3*q^|ces^x5ZIP$RK> zKKMR!PD{Twe6NaA$zye|@K8+#oh0PMILK34kvydpal`yFGpC18mdZKu(HB~~*J+N| z6(^FZHI75XqpgOi_qii3X2waR-nby8+Ft79Yg!Afl#m18hj&bI65dG>IT;&gAokyh zO=zbh_xZ$ggPi7#GmPSb7OQ+|UfRd{X#Y%i+vnP0eY$~$j&gTNd2ce;ob`1-0z4Wk ze8GwnPLI9#U$>pcSw-=DSd%Uadye8vKZKcnGM{b6xqi58C|tH`lNV2Vw!nuDL%MvI zx7oChPaS0t>nr!)4yKlxwhMdusG+SW6#+>rpC&Xsnc(7tz%U^& z4D%SNjmFb8Erl^rS9t1WsUkMbFU@n$K;cw3VJLSKc9YNE`clTLzEzQ|E>uK2l=_0p znh7-tdd|cUMib?%i5CWb;8M~1#7sva=Y4X@J5cz0j{GiZCR}PakM#)7YA5i#u^yyQ zU#Wh+=9|y<*YUD! z){hTQXB{;sI%yzyY0~78If@a4i%m6fq^AaZ@APC(4P4Vx!(&ekpX{lDGd(pE1!?T5 zNsAZs)a0A;CYW8I(Sa_HZU3#AHY$pv|&`2|TrjIUeq=73MX>iFbUED|mpRbWdn+_Nq zyJSN1V$;I-P~YUgt1!-3O*3-w zV&3#Pc6)n-`gEra1)_4N4KYUy^%b2qQePB*a(yiR>`&0;V!xUAGhyj`_t@Vx1<2Gl zyQn`pb)}G07SFfLJ(YO^E{pP`zb+(AcK`$@bT&3+l<2U2w%!`49~qibeMx94_5Ua| zjRvY%kx!E@-Z1c)o*UwXB3GPBku1g&{OP}9uv(H?OaK5MiwFb&00000{{{d;LjnLB O00RI3000000000M9NZHC literal 0 HcmV?d00001 diff --git a/test/test.vcf.gz.tbi b/test/test.vcf.gz.tbi new file mode 100644 index 0000000000000000000000000000000000000000..63156e73682feb684e35c3c613d0e05c358bf0ef GIT binary patch literal 136 zcmb2|=3rp}f&Xj_PR>jW*$mu;pHfm%5)u-U5)v9N@&Li9fLXnbjf$3wPn;0Y30CF_ zJbEDE9s9)jLWd6=IPk!fbD41HocWBRGYphgOD=VP8CW$%4o0!F^D;6#xDt3xii?3k O9?dRk24=9;AOZjpSSa`a literal 0 HcmV?d00001 diff --git a/vcfpp.h b/vcfpp.h index c8d4819..b6c263e 100644 --- a/vcfpp.h +++ b/vcfpp.h @@ -385,7 +385,7 @@ class BcfHeader + "-th sample are not in the VCF.\nparameter samples:" + samples); } } - + /** @brief set the VCF version */ inline void setVersion(const std::string & version) const { @@ -1769,13 +1769,23 @@ class BcfWriter hp = &h; } - /// copy header of given VCF and restrict on samples + /// copy header of given VCF and restrict on samples. if samples=="", FORMAT removed and only sites left void copyHeader(const std::string & vcffile, std::string samples = "-") { htsFile * fp2 = hts_open(vcffile.c_str(), "r"); if(!fp2) throw std::invalid_argument("I/O error: input file is invalid"); - header.hdr = bcf_hdr_read(fp2); - header.setSamples(samples); + if(samples == "") + { // site-only + bcf_hdr_t * hfull = bcf_hdr_read(fp2); + header.hdr = bcf_hdr_subset(hfull, 0, 0, 0); + bcf_hdr_remove(header.hdr, BCF_HL_FMT, NULL); + bcf_hdr_destroy(hfull); + } + else + { + header.hdr = bcf_hdr_read(fp2); + header.setSamples(samples); + } hts_close(fp2); initalHeader(header); }