| Previous changeset 28:6322be79bd8e (2024-09-28) Next changeset 30:c7275bd8b4d6 (2024-11-18) |
|
Commit message:
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tool_collections/snpeff commit 4cc10cfe702828a91ecf8bb58d3f84a36b0578f7 |
|
modified:
gbk2fa.py snpEff.xml snpEff_create_db.xml snpEff_databases.xml snpEff_download.xml snpEff_macros.xml snpeff_get_chr_names.xml test-data/input.vcf |
| b |
| diff -r 6322be79bd8e -r ca2b512e8d7c gbk2fa.py --- a/gbk2fa.py Sat Sep 28 16:27:56 2024 +0000 +++ b/gbk2fa.py Mon Oct 21 13:56:15 2024 +0000 |
| [ |
| @@ -6,6 +6,9 @@ def get_opener(gbk_filename): + """Determines the appropriate opener for a given file, supporting + bzip2, gzip, or standard open. + """ try: bz2.open(gbk_filename).read(1) return bz2.open @@ -18,30 +21,38 @@ return open -parser = argparse.ArgumentParser() -parser.add_argument( - "genbank_file", - help="GenBank input file. Can be compressed with gzip or bzip2" -) -parser.add_argument( - "fasta_file", help="FASTA output datset" -) -parser.add_argument( - "--remove_version", action="store_true", - help="Remove version number from NCBI form formatted accession numbers. " - "For example, this would convert 'B000657.2' to 'B000657'" -) -args = parser.parse_args() +def main(): + parser = argparse.ArgumentParser( + description="Convert GenBank files to FASTA format. " + "Supports gzip and bzip2 compressed files." + ) + parser.add_argument( + "genbank_file", + help="GenBank input file. Can be compressed with gzip or bzip2" + ) + parser.add_argument( + "fasta_file", + help="FASTA output dataset" + ) + parser.add_argument( + "--remove_version", action="store_true", + help="Remove version number from NCBI formatted accession numbers. " + "For example, this converts 'B000657.2' to 'B000657'." + ) + args = parser.parse_args() + + gbk_open = get_opener(args.genbank_file) + with gbk_open(args.genbank_file, 'rt') as input_handle, \ + open(args.fasta_file, 'w') as output_handle: + for seq_record in SeqIO.parse(input_handle, 'genbank'): + if args.remove_version: + seq_id = seq_record.id.split('.')[0] + else: + seq_id = seq_record.id + print(f'Writing FASTA record: {seq_id}') + output_handle.write(f'>{seq_id}\n') + output_handle.write(f'{seq_record.seq}\n') -gbk_open = get_opener(args.genbank_file) -with gbk_open(args.genbank_file, 'rt') as input_handle, \ - open(args.fasta_file, 'w') as output_handle: - for seq_record in SeqIO.parse(input_handle, 'genbank'): - if args.remove_version: - seq_id = seq_record.id.split('.')[0] - else: - seq_id = seq_record.id - print('Writing FASTA record: {}'.format(seq_id)) - print('>' + seq_id, file=output_handle) - print(seq_record.seq, file=output_handle) +if __name__ == "__main__": + main() |
| b |
| diff -r 6322be79bd8e -r ca2b512e8d7c snpEff.xml --- a/snpEff.xml Sat Sep 28 16:27:56 2024 +0000 +++ b/snpEff.xml Mon Oct 21 13:56:15 2024 +0000 |
| [ |
| b'@@ -1,4 +1,4 @@\n-<tool id="snpEff" name="SnpEff eff:" version="@WRAPPER_VERSION@.galaxy2">\n+<tool id="snpEff" name="SnpEff eff:" version="@SNPEFF_VERSION@+galaxy@WRAPPER_VERSION@" profile="23.0">\n <description> annotate variants</description>\n <macros>\n <import>snpEff_macros.xml</import>\n@@ -13,7 +13,7 @@\n ln -s \'${intervals}\' intervals.bed &&\n #end if \n snpEff @JAVA_OPTIONS@ eff\n- -i $inputFormat -o ${outputConditional.outputFormat} -upDownStreamLen $udLength\n+ -i $inputFormat -o ${outputFormat} -upDownStreamLen $udLength\n #if $spliceSiteSize and str($spliceSiteSize) != \'\':\n -spliceSiteSize "$spliceSiteSize"\n #end if\n@@ -53,9 +53,6 @@\n #if $csvStats:\n -csvStats \'$csvFile\'\n #end if\n- #if str($offset) != \'default\':\n- ${offset}\n- #end if\n #if str($chr).strip() != \'\':\n -chr \'$chr\'\n #end if\n@@ -103,35 +100,21 @@\n mkdir \'$statsFile.files_path\' &&\n mv \'$genes_file\' \'#echo os.path.join($statsFile.files_path, $genes_file_name)#\'\n #end if\n- #if $outputConditional.outputFormat == \'gatk\' and $outputConditional.gatk_v1\n- &&\n- ## Replace real SnpEff version with 2.0.5 to prevent this GATK 1.x error: "The version of SnpEff used to generate the SnpEff input file (x.x) is not currently supported by the GATK. Supported versions are: [2.0.5]"\n- sed -i.bak -e \'s/^\\#\\#SnpEffVersion="\\(\\S*\\s\\)/\\#\\#SnpEffVersion="2.0.5 - real is \\1/\' \'$snpeff_output\'\n- #end if\n ]]></command>\n <inputs>\n <param name="input" type="data" format="vcf,bed" label="Sequence changes (SNPs, MNPs, InDels)"/>\n-\n- <param name="inputFormat" type="select" label="Input format">\n+ <param argument="-i" name="inputFormat" type="select" label="Input format" help="Specify the format of input dataset(s)">\n <option value="vcf" selected="true">VCF</option>\n- <option value="bed">BED (Deprecated)</option>\n+ <option value="bed">BED</option>\n </param>\n-\n- <conditional name="outputConditional">\n- <param name="outputFormat" type="select" label="Output format">\n- <option value="vcf" selected="true">VCF (only if input is VCF)</option>\n- <option value="gatk">GATK-compatible VCF (only if input is VCF)</option>\n- <option value="bed">BED</option>\n- <option value="bedAnn">BED annotations</option>\n- </param>\n- <when value="vcf" />\n- <when value="gatk">\n- <param name="gatk_v1" type="boolean" checked="true" label="Compatible with GATK 1.x" />\n- </when>\n- <when value="bed" />\n- <when value="bedAnn" />\n- </conditional>\n- <param name="csvStats" type="boolean" truevalue="-csvStats" falsevalue="" checked="false" label="Create CSV report, useful for downstream analysis (-csvStats)" />\n+ <param argument="-o" name="outputFormat" type="select" label="Output format" help="Specify output format">\n+ <option value="vcf" selected="true">VCF (only if input is VCF)</option>\n+ <option value="gatk">GATK-compatible VCF (only if input is VCF)</option>\n+ <option value="bed">BED</option>\n+ <option value="bedAnn">BED annotations</option>\n+ </param>\n+ <param argument="-csvStats" type="boolean" truevalue="-csvStats" falsevalue="" checked="false" label="Create CSV report?" help="Useful for downstream analyses and report generation" />\n+ <param argument="-noStats" name="generate_stats" type="boolean" truevalue="" falsevalue="-noStats" checked="true" label="Produce Summary Stats?" help="Generates an HTML summary of results"/>\n <conditional name="snpDb">\n <param name="genomeSrc" type="select" label="Genome source">\n <!-- These options are referenced in the'..b' <option value="20000">20000 bases</option>\n </param>\n-\n <param name="spliceSiteSize" argument="-ss" type="select" optional="true" label="Set size for splice sites (donor and acceptor) in bases">\n <option value="1">1 base</option>\n <option value="2" selected="true">2 bases</option>\n@@ -232,7 +212,6 @@\n <option value="8">8 bases</option>\n <option value="9">9 bases</option>\n </param>\n-\n <conditional name="spliceRegion">\n <param name="setSpliceRegions" type="select" label="spliceRegion Settings">\n <option value="no">Use Defaults</option>\n@@ -245,7 +224,6 @@\n <param argument="-spliceRegionIntronMax" type="integer" value="" min="1" max="10" optional="true" label="Set maximum number of bases for splice site region within intron. Default: 8 bases" />\n </when>\n </conditional>\n-\n <param name="annotations" type="select" display="checkboxes" multiple="true" label="Annotation options">\n <option value="-formatEff">Use \'EFF\' field compatible with older versions (instead of \'ANN\')</option>\n <option value="-classic">Use Classic Effect names and amino acid variant annotations (NON_SYNONYMOUS_CODING vs missense_variant and G180R vs p.Gly180Arg/c.538G>C)</option>\n@@ -334,12 +312,6 @@\n </param>\n </when>\n </conditional>\n-\n- <param name="offset" type="select" display="radio" label="Chromosomal position">\n- <option value="default" selected="true">Use default (based on input type)</option>\n- <option value="-0">Force zero-based positions (both input and output)</option>\n- <option value="-1">Force one-based positions (both input and output)</option>\n- </param>\n <param argument="-chr" type="text" label="Text to prepend to chromosome name">\n <help>\n By default SnpEff simplifies all chromosome names. For instance \'chr1\' is just \'1\'.\n@@ -347,7 +319,6 @@\n </help>\n <validator type="regex" message="No whitespace allowed">^\\S*$</validator>\n </param>\n- <param name="generate_stats" argument="-noStats" type="boolean" truevalue="" falsevalue="-noStats" checked="true" label="Produce Summary Stats" />\n <param argument="-noLog" type="boolean" truevalue="-noLog" falsevalue="" checked="true" label="Suppress reporting usage statistics to server" />\n </inputs>\n <outputs>\n@@ -375,8 +346,8 @@\n <param name="generate_stats" value="true"/>\n <output name="snpeff_output">\n <assert_contents>\n- <has_text_matching expression="KJ660346\\t572\\t.*missense_variant" />\n- <has_text_matching expression="KJ660346\\t1024\\t.*synonymous_variant" />\n+ <has_text_matching expression="KJ660346.1\\t572\\t.*missense_variant" />\n+ <has_text_matching expression="KJ660346.1\\t1024\\t.*synonymous_variant" />\n </assert_contents>\n </output>\n <output name="statsFile">\n@@ -398,13 +369,13 @@\n <param name="csvStats" value="true"/>\n <output name="snpeff_output">\n <assert_contents>\n- <has_text_matching expression="KJ660346\\t572\\t.*missense_variant" />\n- <has_text_matching expression="KJ660346\\t1024\\t.*synonymous_variant" />\n+ <has_text_matching expression="KJ660346.1\\t572\\t.*missense_variant" />\n+ <has_text_matching expression="KJ660346.1\\t1024\\t.*synonymous_variant" />\n </assert_contents>\n </output>\n <output name="csvFile">\n <assert_contents>\n- <has_n_lines n="185"/>\n+ <has_n_lines n="134"/>\n <has_n_columns n="1" sep=","/>\n </assert_contents>\n </output>\n' |
| b |
| diff -r 6322be79bd8e -r ca2b512e8d7c snpEff_create_db.xml --- a/snpEff_create_db.xml Sat Sep 28 16:27:56 2024 +0000 +++ b/snpEff_create_db.xml Mon Oct 21 13:56:15 2024 +0000 |
| [ |
| @@ -1,11 +1,11 @@ -<tool id="snpEff_build_gb" name="SnpEff build:" version="@WRAPPER_VERSION@.galaxy6" profile="22.01"> +<tool id="snpEff_build_gb" name="SnpEff build:" version="@SNPEFF_VERSION@+galaxy@WRAPPER_VERSION@" profile="23.0"> <description> database from Genbank or GFF record</description> <macros> <import>snpEff_macros.xml</import> </macros> <requirements> <expand macro="requirement" /> - <requirement type="package" version="1.79">biopython</requirement> + <requirement type="package" version="1.84">biopython</requirement> </requirements> <expand macro="stdio" /> <expand macro="version_command" /> @@ -36,7 +36,7 @@ ln -s '${input_type.input}' 'snpeff_output/${genome_version}/genes.${input_type.input_type_selector}' && #end if - snpEff @JAVA_OPTIONS@ build -v + snpEff @JAVA_OPTIONS@ build -noCheckCds -noCheckProtein -v -configOption '${genome_version}'.genome='${genome_version}' -configOption '${genome_version}'.codonTable='${codon_table}' #if str($input_type.input_type_selector) == "gb": @@ -186,7 +186,7 @@ <help><![CDATA[ **What it does** -This tool uses `"snpEff build -genbank"` or `"snpEff build -gff3"` commands to create a snpEff database. +This tool uses `snpEff build` to create a snpEff database. ------ @@ -201,7 +201,7 @@ .. class:: warningmark - SnpEff errors out on highly fragmented genomes containing multiple scaffolds. This is because a single gene may be split between multiple scaffolds causing SnpEff to crash. If this is happening use GFF route described below. + SnpEff errors out on highly fragmented genomes containing multiple scaffolds. This is because a single gene may be split between multiple scaffolds causing SnpEff to crash. If this is happening use the GFF route described below. ------- |
| b |
| diff -r 6322be79bd8e -r ca2b512e8d7c snpEff_databases.xml --- a/snpEff_databases.xml Sat Sep 28 16:27:56 2024 +0000 +++ b/snpEff_databases.xml Mon Oct 21 13:56:15 2024 +0000 |
| [ |
| @@ -1,4 +1,4 @@ -<tool id="snpEff_databases" name="SnpEff databases:" version="@WRAPPER_VERSION@.galaxy2"> +<tool id="snpEff_databases" name="SnpEff databases:" version="@SNPEFF_VERSION@+galaxy@WRAPPER_VERSION@" profile="23.0"> <description> list available databases</description> <macros> <import>snpEff_macros.xml</import> @@ -19,6 +19,10 @@ | grep -v '${exclude_pattern}' #end if + #if str($include_download_path) == "no": + | cut -f 1,2,3,4 + #end if + > '${snpeff_dbs}' ]]></command> <inputs> @@ -38,7 +42,10 @@ </valid> </sanitizer> </param> - + <param name="include_download_path" type="select" display="radio" label="Include download paths?" help="When snpEff dumps the list of available databases, it includes their download paths. These are not needed in the Galaxy context."> + <option value="yes">Yes</option> + <option value="no" selected="true">No</option> + </param> </inputs> <outputs> <data name="snpeff_dbs" format="tabular" label="${tool.name} @SNPEFF_VERSION@ available databases" /> @@ -63,12 +70,16 @@ <help><![CDATA[ **What it does** -This tool downloads the master list of snpEff databases from @SNPEFF_DATABASE_URL@. You can then look at this list and decide which database to use for your analysis. For example, if **List entries matching the following expression** parameter of this tool is set to *Mouse* the it will produce a tabular dataset with the following content:: +This tool downloads the master list of snpEff databases from a remote SnpEff repository. You can then look at this list and decide which database to use for your analysis. For example, if **List entries matching the following expression** parameter of this tool is set to *Mouse*, it will produce a tabular dataset with the following content:: - mm10 Mouse http://downloads.sourceforge.net/project/snpeff/databases/v4_3/snpEff_v4_3_mm10.zip - mm9 Mouse http://downloads.sourceforge.net/project/snpeff/databases/v4_3/snpEff_v4_3_mm9.zip + mm10 Mouse + mm39 Mouse + mm9 Mouse -This means that there two available snpEff databases for mouse genome versions mm9 and mm10. In order to download these databases you should use identifier from the first column (e.g., mm9 or mm10 in this case). +This means that there are three available snpEff databases for mouse genome. If you want to use mm39 in you analysis: + + - set **Genome source** option of **SnpEff eff** Galaxy tool to *Download on demand* + - enter 'mm39' into **Snpff Genome Version Name** text box ------- @@ -83,6 +94,7 @@ @SNPEFF_IN_GALAXY_INFO@ @EXTERNAL_DOCUMENTATION@ + ]]></help> <expand macro="citations" /> </tool> |
| b |
| diff -r 6322be79bd8e -r ca2b512e8d7c snpEff_download.xml --- a/snpEff_download.xml Sat Sep 28 16:27:56 2024 +0000 +++ b/snpEff_download.xml Mon Oct 21 13:56:15 2024 +0000 |
| [ |
| @@ -1,4 +1,4 @@ -<tool id="snpEff_download" name="SnpEff download:" version="@WRAPPER_VERSION@.galaxy2"> +<tool id="snpEff_download" name="SnpEff download:" version="@SNPEFF_VERSION@+galaxy@WRAPPER_VERSION@" profile="23.0"> <description> download a pre-built database</description> <macros> <import>snpEff_macros.xml</import> @@ -42,7 +42,7 @@ <help><![CDATA[ **What it does** -This tool downloads a specified database from @SNPEFF_DATABASE_URL@. It deposits it into the history. +This tool downloads a specified database from a remote SnpEff repository. It deposits it into the history. ------- |
| b |
| diff -r 6322be79bd8e -r ca2b512e8d7c snpEff_macros.xml --- a/snpEff_macros.xml Sat Sep 28 16:27:56 2024 +0000 +++ b/snpEff_macros.xml Mon Oct 21 13:56:15 2024 +0000 |
| b |
| @@ -1,6 +1,6 @@ <macros> <xml name="requirement"> - <requirement type="package" version="4.3.1t">snpeff</requirement> + <requirement type="package" version="5.2">snpeff</requirement> <yield/> </xml> <xml name="stdio"> @@ -14,9 +14,8 @@ snpEff -version ]]></version_command> </xml> - <token name="@WRAPPER_VERSION@">4.3+T</token> - <token name="@SNPEFF_VERSION@">SnpEff4.3</token> - <token name="@SNPEFF_DATABASE_URL@">https://sourceforge.net/projects/snpeff/files/databases/v4_3/</token> + <token name="@WRAPPER_VERSION@">0</token> + <token name="@SNPEFF_VERSION@">5.2</token> <token name="@JAVA_OPTIONS@">-Xmx\${GALAXY_MEMORY_MB:-8192}m</token> <xml name="ref_select"> <conditional name="reference_source"> @@ -59,7 +58,7 @@ **Download pre-built databases** -SnpEff project generates large numbers of pre-build databases. These are available at @SNPEFF_DATABASE_URL@ and can downloaded. Follow these steps: +SnpEff project generates large numbers of pre-build databases. To obtain and use them follow these steps: #. Use **SnpEff databases** tool to generate a list of existing databases. Note the name of the database you need. #. Use **SnpEff download** tool to download the database. |
| b |
| diff -r 6322be79bd8e -r ca2b512e8d7c snpeff_get_chr_names.xml --- a/snpeff_get_chr_names.xml Sat Sep 28 16:27:56 2024 +0000 +++ b/snpeff_get_chr_names.xml Mon Oct 21 13:56:15 2024 +0000 |
| b |
| @@ -1,4 +1,4 @@ -<tool id="snpEff_get_chr_names" name="SnpEff chromosome-info:" version="@WRAPPER_VERSION@.galaxy2"> +<tool id="snpEff_get_chr_names" name="SnpEff chromosome-info:" version="@SNPEFF_VERSION@+galaxy@WRAPPER_VERSION@" profile="23.0"> <description>list chromosome names/lengths</description> <macros> <import>snpEff_macros.xml</import> |
| b |
| diff -r 6322be79bd8e -r ca2b512e8d7c test-data/input.vcf --- a/test-data/input.vcf Sat Sep 28 16:27:56 2024 +0000 +++ b/test-data/input.vcf Mon Oct 21 13:56:15 2024 +0000 |
| b |
| b'@@ -2,59 +2,9 @@\n ##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">\n ##contig=<ID="KJ660346",length=18959>\n ##reference=http://www.ncbi.nlm.nih.gov/nuccore/KJ660346.2\n-#CHROM\tPOS\tID\tREF\tALT\tQUAL\tFILTER\tINFO\tFORMAT\tEBOV_2014_EM095\tEBOV_2014_G3676\tEBOV_2014_G3680\tEBOV_2014_G3686\tEBOV_2014_G3687\tEBOV_2014_G3683\tEBOV_2014_EM096\tEBOV_2014_EM098\tEBOV_2014_EM104\tEBOV_2014_EM106\tEBOV_2014_EM112\tEBOV_2014_EM113\tEBOV_2014_EM124\tEBOV_2014_G3677\tEBOV_2014_G3682\tEBOV_2014_G3707\tEBOV_2014_G3713\tEBOV_2014_G3724\tEBOV_2014_G3734\tEBOV_2014_G3735\tEBOV_2014_G3750\tEBOV_2014_G3758\tEBOV_2014_G3764\tEBOV_2014_G3769\tEBOV_2014_G3782\tEBOV_2014_G3786\tEBOV_2014_G3788\tEBOV_2014_G3796\tEBOV_2014_G3798\tEBOV_2014_G3799\tEBOV_2014_G3800\tEBOV_2014_G3805\tEBOV_2014_G3807\tEBOV_2014_G3810\tEBOV_2014_G3820\tEBOV_2014_G3838\tEBOV_2014_G3840\tEBOV_2014_G3841\tEBOV_2014_G3848\tEBOV_2014_NM042\tEBOV_2014_G3850\tEBOV_2014_EM110\tEBOV_2014_EM111\tEBOV_2014_EM119\tEBOV_2014_G3729\tEBOV_2014_G3765\tEBOV_2014_G3770\tEBOV_2014_G3789\tEBOV_2014_G3825\tEBOV_2014_G3845\tEBOV_2014_G3851\tEBOV_2014_G3857\tEBOV_2014_EM115\tEBOV_2014_EM120\tEBOV_2014_G3752\tEBOV_2014_G3795\tEBOV_2014_G3808\tEBOV_2014_G3823\tEBOV_2014_EM121\tEBOV_2014_G3771\tEBOV_2014_G3816\tEBOV_2014_G3829\tEBOV_2014_G3846\tEBOV_2014_G3856\tEBOV_2014_G3826\tEBOV_2014_G3827\tEBOV_2014_G3809\tEBOV_2014_G3814\tEBOV_2014_G3821\tEBOV_2014_G3822\tEBOV_2014_G3679\tEBOV_2014_G3819\tEBOV_2014_G3817\tEBOV_2014_G3834\tEBOV_2014_G3818\tEBOV_2014_G3787\tEBOV_2014_G3831\tEBOV_2014_G3670\tEBOV_2014_KJ660346\tEBOV_2014_KJ660347\tEBOV_2014_KJ660348\n-KJ660346\t572\t.\tA\tG\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t800\t.\tC\tT\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t0\t0\t0\n-KJ660346\t1024\t.\tA\tC\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t1288\t.\tA\tT\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\n-KJ660346\t1492\t.\tA\tG\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t1849\t.\tC\tT\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t1\t1\n-KJ660346\t2124\t.\tA\tG\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\n-KJ660346\t2185\t.\tG\tA\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\n-KJ660346\t2341\t.\tA\tG\t.\t.\t.\tGT\t0\t0\t0\t0\t.\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t2364\t.\tA\tG\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t2497\t.\tA\tG\t.\t.\t.\tGT\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t2931\t.\tG\tA\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t1\n-KJ660346\t3116\t.\tC\tG\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0'..b'\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t11142\t.\tG\tA\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t11811\t.\tT\tC\t.\t.\t.\tGT\t1\t1\t1\t0\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\n-KJ660346\t11943\t.\tG\tA\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t1\t0\t0\t0\t0\n-KJ660346\t12878\t.\tG\tA\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t12885\t.\tA\tC\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t.\t.\t.\n-KJ660346\t13856\t.\tG\tA\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t1\t1\n-KJ660346\t13923\t.\tT\tC\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t14019\t.\tT\tC\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t1\t1\t0\t0\t1\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t14232\t.\tC\tT\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t15599\t.\tG\tA\t.\t.\t.\tGT\t1\t1\t1\t0\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\n-KJ660346\t15660\t.\tC\tT\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t1\t1\n-KJ660346\t15963\t.\tG\tA\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t0\t0\t0\n-KJ660346\t16054\t.\tT\tA\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t1\t0\t0\t0\t0\n-KJ660346\t16455\t.\tT\tC\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t16750\t.\tT\tC\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t17142\t.\tT\tC\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t1\t0\t0\t0\n-KJ660346\t17985\t.\tT\tC\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t18412\t.\tT\tC\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\n-KJ660346\t18895\t.\tC\tT\t.\t.\t.\tGT\t0\t0\t0\t0\t0\t0\t0\t.\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t.\t0\t0\t0\t0\t0\t0\t0\t0\t.\t0\t0\t0\t0\t.\t0\t.\t0\t0\t0\t.\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t.\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t0\t1\t1\t1\n+#CHROM\tPOS\tID\tREF\tALT\tQUAL\tFILTER\tINFO\tFORMAT\tEBOV_2014_EM095\tEBOV_2014_G3676\n+KJ660346.1\t572\t.\tA\tG\t.\t.\t.\tGT\t0\t0\n+KJ660346.1\t800\t.\tC\tT\t.\t.\t.\tGT\t0\t0\n+KJ660346.1\t1024\t.\tA\tC\t.\t.\t.\tGT\t0\t0\n+KJ660346.1\t1288\t.\tA\tT\t.\t.\t.\tGT\t0\t0\n+KJ660346.1\t1492\t.\tA\tG\t.\t.\t.\tGT\t0\t0\n' |