# HG changeset patch
# User sanbi-uwc
# Date 1516881164 18000
# Node ID d1116f934d3d3c45f0dbb1b54fa2287cd89bf4cc
planemo upload for repository https://github.com/SANBI-SA/tools-sanbi-uwc/tree/master/tools/freebayes commit fa47457b312f36cead164f993d46bcc31d13fe57
diff -r 000000000000 -r d1116f934d3d freebayes.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/freebayes.xml Thu Jan 25 06:52:44 2018 -0500
@@ -0,0 +1,705 @@
+
+ bayesian genetic variant detector
+
+ macros.xml
+
+
+ gawk
+ parallel
+
+ &1 || echo "Error running samtools faidx for FreeBayes" >&2 &&
+ #else:
+ #set $reference_fasta_filename = str( $reference_source.ref_file.fields.path )
+ #end if
+
+ #if $reference_source.batchmode.processmode == 'merge':
+ #set $input_bamfiles = $reference_source.batchmode.input_bams
+ #else:
+ #set $input_bamfiles = [ $reference_source.batchmode.input_bams ]
+ #end if
+
+ #for $bam_count, $input_bam in enumerate( $input_bamfiles ):
+ ln -s -f '${input_bam}' 'b_${bam_count}.bam' &&
+ ln -s -f '${input_bam.metadata.bam_index}' 'b_${bam_count}.bam.bai' &&
+ #end for
+
+ ## Tabixize optional input_variant_vcf file (for --variant-input option)
+ #if ( str( $options_type.options_type_selector ) == 'cline' or str( $options_type.options_type_selector ) == 'full' ) and str( $options_type.optional_inputs.optional_inputs_selector ) == 'set' and str( $options_type.optional_inputs.input_variant_type.input_variant_type_selector ) == "provide_vcf":
+ ln -s -f '${options_type.optional_inputs.input_variant_type.input_variant_vcf}' input_variant_vcf.vcf.gz &&
+ ln -s -f '${Tabixized_input}' input_variant_vcf.vcf.gz.tbi &&
+ #end if
+
+ ##if the user has specified a region or target file, just use that instead of calculating a set of unique regions
+ #if str( $target_limit_type.target_limit_type_selector ) == "limit_by_target_file":
+ ln -s '${target_limit_type.input_target_bed}' regions_all.bed &&
+ #elif str( $target_limit_type.target_limit_type_selector ) == "limit_by_region":
+ printf '${target_limit_type.region_chromosome}\t${target_limit_type.region_start}\t${target_limit_type.region_end}' > regions_all.bed &&
+ #else
+ ##divide up the regions in the bam file for efficient processing
+ #for $bam_count, $input_bam in enumerate( $input_bamfiles ):
+ samtools view -H b_${bam_count}.bam |
+ grep '^@SQ' |
+ cut -f 2- |
+ awk '{ gsub("^SN:","",$1); gsub("^LN:","",$2); print $1"\t0\t"$2; }' >> regions_all.bed &&
+ #end for
+ #end if
+
+ sort -u regions_all.bed > regions_uniq.bed &&
+ ## split into even small chunks, this has some disatvantages and will not be used for the moment
+ ## bedtools makewindows -b regions_uniq.bed -w 10000000 -s 9990000 > regions.bed &&
+
+ mkdir vcf_output failed_alleles trace &&
+
+ ## Finished setting up inputs
+
+ for i in `cat regions_uniq.bed | awk '{print $1":"$2".."$3}'`;
+ do
+ echo "
+
+ ## COMMAND LINE STARTS HERE
+
+ freebayes
+
+ --region '\$i'
+
+ #for $bam_count, $input_bam in enumerate( $input_bamfiles ):
+ --bam 'b_${bam_count}.bam'
+ #end for
+ --fasta-reference '${reference_fasta_filename}'
+
+ ## Outputs
+ --vcf './vcf_output/part_\$i.vcf'
+
+ ##advanced options
+ #if str( $options_type.options_type_selector ) == "simple":
+ #pass
+ #elif str( $options_type.options_type_selector ) == "simple_w_filters":
+ --standard-filters
+ --min-coverage ${options_type.min_coverage}
+ #elif str( $options_type.options_type_selector ) == "naive":
+ --haplotype-length 0
+ --min-alternate-count 1
+ --min-alternate-fraction 0
+ --pooled-continuous
+ --report-monomorphic
+ #elif str( $options_type.options_type_selector ) == "naive_w_filters":
+ --haplotype-length 0
+ --min-alternate-count 1
+ --min-alternate-fraction 0
+ --pooled-continuous
+ --report-monomorphic
+ --standard-filters
+ --min-coverage ${options_type.min_coverage}
+ #elif str( $options_type.options_type_selector ) == "full":
+ #if str( $options_type.optional_inputs.optional_inputs_selector ) == 'set':
+ ${options_type.optional_inputs.report_monomorphic}
+
+ #if $options_type.optional_inputs.output_trace_option:
+ --trace ./trace/part_'\$i'.txt
+ #end if
+ #if $options_type.optional_inputs.output_failed_alleles_option:
+ --failed-alleles ./failed_alleles/part_'\$i'.bed
+ #end if
+ #if $options_type.optional_inputs.samples:
+ --samples '${options_type.optional_inputs.samples}'
+ #end if
+ #if $options_type.optional_inputs.populations:
+ --populations '${options_type.optional_inputs.populations}'
+ #end if
+ #if $options_type.optional_inputs.A:
+ --cnv-map '${options_type.optional_inputs.A}'
+ #end if
+ #if str( $options_type.optional_inputs.input_variant_type.input_variant_type_selector ) == "provide_vcf":
+ --variant-input 'input_variant_vcf.vcf.gz' ## input_variant_vcf.vcf.gz is symlinked to a galaxy-generated dataset in "Tabixize optional input_variant_vcf file" section of the command line above
+ ${options_type.optional_inputs.input_variant_type.only_use_input_alleles}
+ #end if
+ #if $options_type.optional_inputs.haplotype_basis_alleles:
+ --haplotype-basis-alleles '${options_type.optional_inputs.haplotype_basis_alleles}'
+ #end if
+ #if $options_type.optional_inputs.observation_bias:
+ --observation-bias '${options_type.optional_inputs.observation_bias}'
+ #end if
+ #if $options_type.optional_inputs.contamination_estimates:
+ --contamination-estimates '${options_type.optional_inputs.contamination_estimates}'
+ #end if
+ #end if
+
+ ## REPORTING
+ #if str( $options_type.reporting.reporting_selector ) == "set":
+ --pvar ${options_type.reporting.pvar}
+ #end if
+ ## POPULATION MODEL
+ #if str( $options_type.population_model.population_model_selector ) == "set":
+ --theta ${options_type.population_model.T}
+ --ploidy ${options_type.population_model.P}
+ ${options_type.population_model.J}
+ ${options_type.population_model.K}
+ #end if
+
+ ## REFERENCE ALLELE
+ #if str( $options_type.reference_allele.reference_allele_selector ) == "set":
+ ${options_type.reference_allele.Z}
+ --reference-quality '${options_type.reference_allele.reference_quality}'
+ #end if
+
+ ## ALLELE SCOPE
+ #if str( $options_type.allele_scope.allele_scope_selector ) == "set":
+ ${options_type.allele_scope.I}
+ ${options_type.allele_scope.i}
+ ${options_type.allele_scope.X}
+ ${options_type.allele_scope.u}
+ ${options_type.allele_scope.no_partial_observations}
+
+ -n ${options_type.allele_scope.n}
+
+ --haplotype-length ${options_type.allele_scope.haplotype_length}
+ --min-repeat-size ${options_type.allele_scope.min_repeat_length}
+ --min-repeat-entropy ${options_type.allele_scope.min_repeat_entropy}
+ #end if
+
+ ## REALIGNMENT
+ ${options_type.O}
+
+ ##INPUT FILTERS
+ #if str( $options_type.input_filters.input_filters_selector ) == "set":
+ ${options_type.input_filters.use_duplicate_reads}
+ -m ${options_type.input_filters.m}
+ -q ${options_type.input_filters.q}
+ -R ${options_type.input_filters.R}
+ -Y ${options_type.input_filters.Y}
+ -e ${options_type.input_filters.e}
+ -F ${options_type.input_filters.F}
+ -C ${options_type.input_filters.C}
+ -G ${options_type.input_filters.G}
+
+ #if str( $options_type.input_filters.mismatch_filters.mismatch_filters_selector ) == "set":
+ -Q ${options_type.input_filters.mismatch_filters.Q}
+ #if str($options_type.input_filters.mismatch_filters.U)
+ -U ${options_type.input_filters.mismatch_filters.U}
+ #end if
+ -z ${options_type.input_filters.mismatch_filters.z}
+
+ --read-snp-limit ${options_type.input_filters.mismatch_filters.read_snp_limit}
+ #end if
+
+ --min-coverage ${options_type.input_filters.min_coverage}
+ --min-alternate-qsum ${options_type.input_filters.min_alternate_qsum}
+ #end if
+
+ ## POPULATION AND MAPPABILITY PRIORS
+ #if str( $options_type.population_mappability_priors.population_mappability_priors_selector ) == "set":
+ ${options_type.population_mappability_priors.k}
+ ${options_type.population_mappability_priors.w}
+ ${options_type.population_mappability_priors.V}
+ ${options_type.population_mappability_priors.a}
+ #end if
+
+ ## GENOTYPE LIKELIHOODS
+ #if str( $options_type.genotype_likelihoods.genotype_likelihoods_selector ) == "set":
+ ${$options_type.genotype_likelihoods.experimental_gls}
+
+ --base-quality-cap ${$options_type.genotype_likelihoods.base_quality_cap}
+ --prob-contamination ${$options_type.genotype_likelihoods.prob_contamination}
+ #end if
+
+ ## ALGORITHMIC FEATURES
+ #if str( $options_type.algorithmic_features.algorithmic_features_selector ) == "set":
+ -B '${options_type.algorithmic_features.B}'
+ -W '${options_type.algorithmic_features.W}'
+ -D '${options_type.algorithmic_features.D}'
+
+ #if str($options_type.algorithmic_features.genotype_variant_threshold)
+ -S ${options_type.algorithmic_features.genotype_variant_threshold}
+ #end if
+
+ ${options_type.algorithmic_features.N}
+ ${options_type.algorithmic_features.j}
+ ${options_type.algorithmic_features.H}
+ ${options_type.algorithmic_features.genotype_qualities}
+ ${options_type.algorithmic_features.report_genotype_likelihood_max}
+
+ --genotyping-max-banddepth ${options_type.algorithmic_features.genotyping_max_banddepth}
+ #end if
+ #end if
+
+ ";
+ done > freebayes_commands.sh &&
+
+ cat freebayes_commands.sh |
+ parallel --will-cite -j \${GALAXY_SLOTS:-1} &&
+
+ ## make VCF header
+ grep "^#" "./vcf_output/part_\$i.vcf" > header.txt &&
+
+ for i in `cat regions_uniq.bed | awk '{print $1":"$2".."$3}'`;
+ do
+ ## if this fails then it bails out the script
+ cat "./vcf_output/part_\$i.vcf" | grep -v "^#" || true
+ ;
+ done | sort -k1,1 -k2,2n -k5,5 -u | cat header.txt - > '${output_vcf}'
+
+ #if str( $options_type.options_type_selector ) == "full":
+ #if str( $options_type.optional_inputs.optional_inputs_selector ) == 'set':
+ #if $options_type.optional_inputs.output_failed_alleles_option:
+ &&
+ for i in `cat regions.bed | awk '{print $1":"$2".."$3}'`;
+ do
+ cat "./failed_alleles/part_\$i.bed"
+ ;
+ done > '${output_failed_alleles_bed}'
+ #end if
+
+ #if $options_type.optional_inputs.output_trace_option:
+ &&
+ for i in `cat regions.bed | awk '{print $1":"$2".."$3}'`;
+ do
+ cat './trace/part_\$i.txt'
+ ;
+ done > '${output_trace}'
+ #end if
+ #end if
+ #end if
+ ]]>
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ ( options_type['options_type_selector'] == 'cline' or options_type['options_type_selector'] == 'full' ) and options_type['optional_inputs']['optional_inputs_selector'] == 'set' and options_type['optional_inputs']['output_failed_alleles_option'] is True
+
+
+ ( options_type['options_type_selector'] == 'cline' or options_type['options_type_selector'] == 'full' ) and options_type['optional_inputs']['optional_inputs_selector'] == 'set' and options_type['optional_inputs']['output_trace_option'] is True
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+**What it does**
+
+FreeBayes is a Bayesian genetic variant detector designed to find small polymorphisms, specifically SNPs (single-nucleotide polymorphisms), indels (insertions and deletions), MNPs (multi-nucleotide polymorphisms), and complex events (composite insertion and substitution events) smaller than the length of a short-read sequencing alignment.
+
+See https://github.com/ekg/freebayes for details on FreeBayes.
+
+------
+
+**Description**
+
+Provided some BAM dataset(s) and a reference sequence, FreeBayes will produce a VCF dataset describing SNPs, indels, and complex variants in samples in the input alignments.
+
+By default, FreeBayes will consider variants supported by at least 2 observations in a single sample (-C) and also by at least 20% of the reads from a single sample (-F). These settings are suitable to low to high depth sequencing in haploid and diploid samples, but users working with polyploid or pooled samples may wish to adjust them depending on the characteristics of their sequencing data.
+
+FreeBayes is capable of calling variant haplotypes shorter than a read length where multiple polymorphisms segregate on the same read. The maximum distance between polymorphisms phased in this way is determined by the --max-complex-gap, which defaults to 3bp. In practice, this can comfortably be set to half the read length.
+
+Ploidy may be set to any level (-p), but by default all samples are assumed to be diploid. FreeBayes can model per-sample and per-region variation in copy-number (-A) using a copy-number variation map.
+
+FreeBayes can act as a frequency-based pooled caller and describe variants and haplotypes in terms of observation frequency rather than called genotypes. To do so, use --pooled-continuous and set input filters to a suitable level. Allele observation counts will be described by AO and RO fields in the VCF output.
+
+-------
+
+**Galaxy-specific options**
+
+Galaxy allows five levels of control over FreeBayes options, provided by the **Choose parameter selection level** menu option. These are:
+
+ 1. *Simple diploid calling*: The simplest possible FreeBayes application. Equivalent to using FreeBayes with only a BAM input and no other parameter options.
+ 2. *Simple diploid calling with filtering and coverage*: Same as #1 plus two additional options: -0 (standard filters: --min-mapping-quality 30 --min-base-quality 20 --min-supporting-allele-qsum 0 --genotype-variant-threshold 0) and --min-coverage.
+ 3. *Frequency-based pooled calling*: This is equivalent to using FreeBayes with the following options: --haplotype-length 0 --min-alternate-count 1 --min-alternate-fraction 0 --pooled-continuous --report-monomorphic. This is the best choice for calling variants in mixtures such as viral, bacterial, or organellar genomes.
+ 4. *Frequency-based pooled calling with filtering and coverage*: Same as #3 but adds -0 and --min-coverage like in #2.
+ 5. *Complete list of all options*: Gives you full control by exposing all FreeBayes options as Galaxy parameters.
+
+------
+
+**Acknowledgments**
+
+The initial version of the wrapper was produced by Dan Blankenberg and upgraded by Anton Nekrutenko.
+TNG was developed by Bjoern Gruening.
+
+
+
+ @article{Tange2011a,
+ title = {GNU Parallel - The Command-Line Power Tool},
+ author = {O. Tange},
+ address = {Frederiksberg, Denmark},
+ journal = {;login: The USENIX Magazine},
+ month = {Feb},
+ number = {1},
+ volume = {36},
+ url = {http://www.gnu.org/s/parallel},
+ year = {2011},
+ pages = {42-47}
+ }
+
+
+
diff -r 000000000000 -r d1116f934d3d leftalign.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/leftalign.xml Thu Jan 25 06:52:44 2018 -0500
@@ -0,0 +1,70 @@
+
+
+ indels in BAM datasets
+
+ macros.xml
+
+
+ &1 || echo "Error running samtools faidx for leftalign" >&2 &&
+ #else:
+ #set $reference_fasta_filename = str( $reference_source.ref_file.fields.path )
+ #end if
+
+ ##start leftalign commandline
+ cat '${input_bam}' |
+ bamleftalign
+ --fasta-reference '${reference_fasta_filename}'
+ -c
+ --max-iterations "${iterations}"
+ > '${output_bam}'
+ ]]>
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+When calling indels, it is important to homogenize the positional distribution of insertions and deletions in the input by using left realignment. Left realignment will place all indels in homopolymer and microsatellite repeats at the same position, provided that doing so does not introduce mismatches between the read and reference other than the indel. This method is computationally inexpensive and handles the most common classes of alignment inconsistency.
+
+This is leftalign utility from FreeBayes package.
+
+
+
diff -r 000000000000 -r d1116f934d3d macros.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/macros.xml Thu Jan 25 06:52:44 2018 -0500
@@ -0,0 +1,49 @@
+
+ 1.1.0.46
+
+
+ freebayes
+ samtools
+
+
+
+
+
+
+ @misc{1207.3907,
+ Author = {Erik Garrison},
+ Title = {Haplotype-based variant detection from short-read sequencing},
+ Year = {2012},
+ Eprint = {arXiv:1207.3907},
+ url = {http://arxiv.org/abs/1207.3907}
+ }
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff -r 000000000000 -r d1116f934d3d test-data/freebayes-phix174-test1.vcf
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/freebayes-phix174-test1.vcf Thu Jan 25 06:52:44 2018 -0500
@@ -0,0 +1,85 @@
+##fileformat=VCFv4.2
+##reference=localref.fa
+##contig=
+##phasing=none
+##commandline="freebayes --region phiX174:0..5386 --bam b_0.bam --fasta-reference localref.fa --vcf ./vcf_output/part_phiX174:0..5386.vcf"
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT unknown
+phiX174 311 . A G 0.00392026 . AB=0.222222;ABP=9.04217;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=9;DPB=9;DPRA=0;EPP=7.35324;EPPR=5.80219;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=33.5714;NS=1;NUMALT=1;ODDS=7.0097;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=70;QR=478;RO=7;RPL=0;RPP=7.35324;RPPR=3.32051;RPR=2;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=5;SRP=5.80219;SRR=2;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:9:7,2:7:478:2:70:-1.91487,0,-18.7749
+phiX174 374 . T G 23.3387 . AB=0.5;ABP=3.0103;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=4;DPB=4;DPRA=0;EPP=7.35324;EPPR=7.35324;GTI=0;LEN=1;MEANALT=1;MQM=37;MQMR=25;NS=1;NUMALT=1;ODDS=1.58025;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=75;QR=141;RO=2;RPL=1;RPP=3.0103;RPPR=7.35324;RPR=1;RUN=1;SAF=1;SAP=3.0103;SAR=1;SRF=0;SRP=7.35324;SRR=2;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:4:2,2:2:141:2:75:-5.19433,0,-3.54586
+phiX174 913 . A C 2.27318 . AB=0.4;ABP=3.44459;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=5;DPB=5;DPRA=0;EPP=3.0103;EPPR=3.73412;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=29;NS=1;NUMALT=1;ODDS=0.380296;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=67;QR=171;RO=3;RPL=1;RPP=3.0103;RPPR=3.73412;RPR=1;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=3;SRP=9.52472;SRR=0;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:5:3,2:3:171:2:67:-3.54411,0,-6.47921
+phiX174 1205 . A C 0.00388031 . AB=0.2;ABP=10.8276;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=10;DPB=10;DPRA=0;EPP=7.35324;EPPR=7.35324;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=34;NS=1;NUMALT=1;ODDS=7.01995;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=67;QR=456;RO=8;RPL=2;RPP=7.35324;RPPR=7.35324;RPR=0;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=2;SRP=7.35324;SRR=6;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:10:8,2:8:456:2:67:-2.11454,0,-21.3383
+phiX174 1245 . G T 0.0324554 . AB=0.222222;ABP=9.04217;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=9;DPB=9;DPRA=0;EPP=7.35324;EPPR=3.32051;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=28.4286;NS=1;NUMALT=1;ODDS=4.8927;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=70;QR=389;RO=7;RPL=2;RPP=7.35324;RPPR=3.32051;RPR=0;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=0;SRP=18.2106;SRR=7;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:9:7,2:7:389:2:70:-2.83427,0,-14.5618
+phiX174 1249 . T G 0.0166698 . AB=0.222222;ABP=9.04217;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=9;DPB=9;DPRA=0;EPP=3.0103;EPPR=3.32051;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=30.1429;NS=1;NUMALT=1;ODDS=5.56079;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=74;QR=464;RO=7;RPL=1;RPP=3.0103;RPPR=5.80219;RPR=1;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=2;SRP=5.80219;SRR=5;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:9:7,2:7:464:2:74:-1.94207,0,-16.5662
+phiX174 1445 . C A 0.147157 . AB=0.285714;ABP=5.80219;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=7;DPB=7;DPRA=0;EPP=7.35324;EPPR=6.91895;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=32.2;NS=1;NUMALT=1;ODDS=3.36782;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=76;QR=273;RO=5;RPL=1;RPP=3.0103;RPPR=6.91895;RPR=1;RUN=1;SAF=1;SAP=3.0103;SAR=1;SRF=2;SRP=3.44459;SRR=3;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:7:5,2:5:273:2:76:-2.52649,0,-12.4911
+phiX174 1577 . A C 0.0123232 . AB=0.222222;ABP=9.04217;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=9;DPB=9;DPRA=0;EPP=7.35324;EPPR=3.32051;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=35.2857;NS=1;NUMALT=1;ODDS=5.8634;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=60;QR=460;RO=7;RPL=1;RPP=3.0103;RPPR=10.7656;RPR=1;RUN=1;SAF=1;SAP=3.0103;SAR=1;SRF=4;SRP=3.32051;SRR=3;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:9:7,2:7:460:2:60:-1.81064,0,-19.8257
+phiX174 1631 . T G 0.00100612 . AB=0.2;ABP=10.8276;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=10;DPB=10;DPRA=0;EPP=7.35324;EPPR=4.09604;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=28;NS=1;NUMALT=1;ODDS=8.3701;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=68;QR=500;RO=8;RPL=0;RPP=7.35324;RPPR=3.0103;RPR=2;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=3;SRP=4.09604;SRR=5;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:10:8,2:8:500:2:68:-1.52818,0,-17.3788
+phiX174 1772 . T G 0.0180574 . AB=0.2;ABP=10.8276;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=10;DPB=10;DPRA=0;EPP=3.0103;EPPR=3.32051;GTI=0;LEN=1;MEANALT=2;MQM=31;MQMR=31.8571;NS=1;NUMALT=1;ODDS=5.48067;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=59;QR=425;RO=7;RPL=1;RPP=3.0103;RPPR=3.32051;RPR=1;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=1;SRP=10.7656;SRR=6;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:10:7,2:7:425:2:59:-1.97686,0,-17.3816
+phiX174 1945 . T G 1.01422 . AB=0.333333;ABP=4.45795;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=6;DPB=6;DPRA=0;EPP=3.0103;EPPR=5.18177;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=37;NS=1;NUMALT=1;ODDS=1.3354;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=59;QR=263;RO=4;RPL=1;RPP=3.0103;RPPR=5.18177;RPR=1;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=2;SRP=3.0103;SRR=2;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:6:4,2:4:263:2:59:-3.25425,0,-11.8637
+phiX174 2230 . T G 0.00580188 . AB=0.2;ABP=10.8276;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=10;DPB=10;DPRA=0;EPP=7.35324;EPPR=3.0103;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=35.5;NS=1;NUMALT=1;ODDS=6.61746;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=80;QR=491;RO=8;RPL=0;RPP=7.35324;RPPR=4.09604;RPR=2;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=3;SRP=4.09604;SRR=5;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:10:8,2:8:491:2:80:-2.28934,0,-22.3591
+phiX174 2699 . C A 2.04855 . AB=0.5;ABP=3.0103;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=4;DPB=4;DPRA=0;EPP=7.35324;EPPR=3.0103;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=37;NS=1;NUMALT=1;ODDS=0.517362;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=76;QR=109;RO=2;RPL=1;RPP=3.0103;RPPR=3.0103;RPR=1;RUN=1;SAF=1;SAP=3.0103;SAR=1;SRF=0;SRP=7.35324;SRR=2;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:4:2,2:2:109:2:76:-3.4054,0,-5.58455
+phiX174 2722 . T G 1.01772 . AB=0.4;ABP=3.44459;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=5;DPB=5;DPRA=0;EPP=7.35324;EPPR=3.73412;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=33;NS=1;NUMALT=1;ODDS=1.33197;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=67;QR=150;RO=3;RPL=1;RPP=3.0103;RPPR=3.73412;RPR=1;RUN=1;SAF=1;SAP=3.0103;SAR=1;SRF=0;SRP=9.52472;SRR=3;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:5:3,2:3:150:2:67:-3.1308,0,-7.62253
+phiX174 2814 . T G 0.670623 . AB=0.4;ABP=3.44459;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=5;DPB=5;DPRA=0;EPP=3.0103;EPPR=3.73412;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=33;NS=1;NUMALT=1;ODDS=1.79322;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=57;QR=135;RO=3;RPL=1;RPP=3.0103;RPPR=3.73412;RPR=1;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=0;SRP=9.52472;SRR=3;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:5:3,2:3:135:2:57:-2.93048,0,-6.73807
+phiX174 2828 . T G 0.160629 . AB=0.25;ABP=7.35324;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=8;DPB=8;DPRA=0;EPP=7.35324;EPPR=3.0103;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=31;NS=1;NUMALT=1;ODDS=3.27866;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=100;QR=348;RO=6;RPL=0;RPP=7.35324;RPPR=4.45795;RPR=2;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=1;SRP=8.80089;SRR=5;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:8:6,2:6:348:2:100:-3.34336,0,-13.6819
+phiX174 2983 . T G 1.55467 . AB=0.4;ABP=3.44459;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=5;DPB=5;DPRA=0;EPP=3.0103;EPPR=3.73412;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=33;NS=1;NUMALT=1;ODDS=0.843116;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=57;QR=150;RO=3;RPL=1;RPP=3.0103;RPPR=3.73412;RPR=1;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=2;SRP=3.73412;SRR=1;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:5:3,2:3:150:2:57:-3.34311,0,-7.61298
+phiX174 3110 . T C 1.16857 . AB=0.333333;ABP=5.18177;AC=1;AF=0.5;AN=2;AO=3;CIGAR=1X;DP=9;DPB=9;DPRA=0;EPP=3.73412;EPPR=4.45795;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=35;NS=1;NUMALT=1;ODDS=1.17522;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=94;QR=361;RO=6;RPL=0;RPP=9.52472;RPPR=3.0103;RPR=3;RUN=1;SAF=1;SAP=3.73412;SAR=2;SRF=1;SRP=8.80089;SRR=5;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:9:6,3:6:361:3:94:-4.02963,0,-16.4039
+phiX174 3155 . T G 5.42491 . AB=0.285714;ABP=5.80219;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=7;DPB=7;DPRA=0;EPP=3.0103;EPPR=3.44459;GTI=0;LEN=1;MEANALT=1;MQM=37;MQMR=29.8;NS=1;NUMALT=1;ODDS=0.911202;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=75;QR=299;RO=5;RPL=1;RPP=3.0103;RPPR=3.44459;RPR=1;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=4;SRP=6.91895;SRR=1;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:7:5,2:5:299:2:75:-4.38484,0,-11.5688
+phiX174 3325 . A C 0.0705175 . AB=0.333333;ABP=4.45795;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=6;DPB=6;DPRA=0;EPP=7.35324;EPPR=5.18177;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=31;NS=1;NUMALT=1;ODDS=4.11232;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=58;QR=284;RO=4;RPL=2;RPP=7.35324;RPPR=3.0103;RPR=0;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=3;SRP=5.18177;SRR=1;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:6:4,2:4:284:2:58:-2.65031,0,-9.66348
+phiX174 3418 . A C 0.000937887 . AB=0.2;ABP=10.8276;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=10;DPB=10;DPRA=0;EPP=7.35324;EPPR=3.0103;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=34;NS=1;NUMALT=1;ODDS=8.44033;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=60;QR=526;RO=8;RPL=2;RPP=7.35324;RPPR=4.09604;RPR=0;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=5;SRP=4.09604;SRR=3;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:10:8,2:8:526:2:60:-1.49768,0,-21.0207
+phiX174 3729 . C T 0.295112 . AB=0.333333;ABP=4.45795;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=6;DPB=6;DPRA=0;EPP=7.35324;EPPR=5.18177;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=37;NS=1;NUMALT=1;ODDS=2.65478;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=63;QR=232;RO=4;RPL=2;RPP=7.35324;RPPR=5.18177;RPR=0;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=0;SRP=11.6962;SRR=4;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:6:4,2:4:232:2:63:-3.28331,0,-11.7007
+phiX174 4031 . T G 0.0784807 . AB=0.25;ABP=7.35324;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=8;DPB=8;DPRA=0;EPP=7.35324;EPPR=4.45795;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=27;NS=1;NUMALT=1;ODDS=4.00441;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=85;QR=322;RO=6;RPL=2;RPP=7.35324;RPPR=8.80089;RPR=0;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=3;SRP=3.0103;SRR=3;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:8:6,2:6:322:2:85:-3.02817,0,-11.6802
+phiX174 4502 . A C 0.232919 . AB=0.333333;ABP=4.45795;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=6;DPB=6;DPRA=0;EPP=3.0103;EPPR=5.18177;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=37;NS=1;NUMALT=1;ODDS=2.89868;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=55;QR=284;RO=4;RPL=1;RPP=3.0103;RPPR=11.6962;RPR=1;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=3;SRP=5.18177;SRR=1;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:6:4,2:4:284:2:55:-2.57533,0,-11.8832
+phiX174 4558 . C G 0.076247 . AB=0.333333;ABP=4.45795;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=6;DPB=6;DPRA=0;EPP=7.35324;EPPR=3.0103;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=28;NS=1;NUMALT=1;ODDS=4.03372;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=59;QR=204;RO=4;RPL=2;RPP=7.35324;RPPR=11.6962;RPR=0;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=2;SRP=3.0103;SRR=2;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:6:4,2:4:204:2:59:-2.68445,0,-8.4362
+phiX174 4655 . T G 0.399408 . AB=0.333333;ABP=4.45795;AC=1;AF=0.5;AN=2;AO=2;CIGAR=1X;DP=6;DPB=6;DPRA=0;EPP=3.0103;EPPR=3.0103;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=37;NS=1;NUMALT=1;ODDS=2.33999;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=70;QR=193;RO=4;RPL=1;RPP=3.0103;RPPR=3.0103;RPR=1;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=0;SRP=11.6962;SRR=4;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/1:6:4,2:4:193:2:70:-2.81796,0,-11.1489
diff -r 000000000000 -r d1116f934d3d test-data/freebayes-phix174-test2.vcf
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/freebayes-phix174-test2.vcf Thu Jan 25 06:52:44 2018 -0500
@@ -0,0 +1,78 @@
+##fileformat=VCFv4.2
+##reference=localref.fa
+##contig=
+##phasing=none
+##commandline="freebayes --region phiX174:0..5386 --bam b_0.bam --fasta-reference localref.fa --vcf ./vcf_output/part_phiX174:0..5386.vcf --haplotype-length 0 --min-alternate-count 1 --min-alternate-fraction 0 --pooled-continuous --report-monomorphic --standard-filters --min-coverage 14"
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT unknown
+phiX174 1134 . A . 0 . DP=14;DPB=14;EPPR=3.0103;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=916;RO=14;RPPR=8.59409 GT:DP:AD:RO:QR:AO:QA:GL 0/0:14:14:14:916:.:.:0
+phiX174 1135 . C . 0 . DP=14;DPB=14;EPPR=3.63072;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=957;RO=14;RPPR=5.49198 GT:DP:AD:RO:QR:AO:QA:GL 0/0:14:14:14:957:.:.:0
+phiX174 1136 . G . 0 . DP=14;DPB=14;EPPR=8.59409;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=919;RO=14;RPPR=3.0103 GT:DP:AD:RO:QR:AO:QA:GL 0/0:14:14:14:919:.:.:0
+phiX174 1137 . C . 0 . DP=14;DPB=14;EPPR=5.49198;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=927;RO=14;RPPR=3.63072 GT:DP:AD:RO:QR:AO:QA:GL 0/0:14:14:14:927:.:.:0
+phiX174 1138 . C . 0 . DP=14;DPB=14;EPPR=3.63072;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=876;RO=14;RPPR=5.49198 GT:DP:AD:RO:QR:AO:QA:GL 0/0:14:14:14:876:.:.:0
+phiX174 1139 . G . 0 . DP=14;DPB=14;EPPR=3.63072;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=946;RO=14;RPPR=5.49198 GT:DP:AD:RO:QR:AO:QA:GL 0/0:14:14:14:946:.:.:0
+phiX174 1140 . T . 0 . DP=15;DPB=15;EPPR=4.31318;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=927;RO=15;RPPR=4.31318 GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:15:15:927:.:.:0
+phiX174 1141 . T G 3.61953e-07 . AB=0;ABP=0;AC=0;AF=0;AN=2;AO=1;CIGAR=1X;DP=15;DPB=15;DPRA=0;EPP=5.18177;EPPR=3.63072;GTI=0;LEN=1;MEANALT=1;MQM=37;MQMR=37;NS=1;NUMALT=1;ODDS=16.3003;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=50;QR=880;RO=14;RPL=0;RPP=5.18177;RPPR=5.49198;RPR=1;RUN=1;SAF=0;SAP=5.18177;SAR=1;SRF=5;SRP=5.49198;SRR=9;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:14,1:14:880:1:50:0,-0.836684,-42.5367
+phiX174 1142 . G . 0 . DP=16;DPB=16;EPPR=5.18177;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=1082;RO=16;RPPR=3.55317 GT:DP:AD:RO:QR:AO:QA:GL 0/0:16:16:16:1082:.:.:0
+phiX174 1143 . G . 0 . DP=16;DPB=16;EPPR=5.18177;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=1029;RO=16;RPPR=3.55317 GT:DP:AD:RO:QR:AO:QA:GL 0/0:16:16:16:1029:.:.:0
+phiX174 1144 . C . 0 . DP=16;DPB=16;EPPR=5.18177;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=1114;RO=16;RPPR=3.55317 GT:DP:AD:RO:QR:AO:QA:GL 0/0:16:16:16:1114:.:.:0
+phiX174 1145 . G . 0 . DP=15;DPB=15;EPPR=6.62942;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=1042;RO=15;RPPR=3.15506 GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:15:15:1042:.:.:0
+phiX174 1146 . C . 0 . DP=15;DPB=15;EPPR=6.62942;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=1058;RO=15;RPPR=3.15506 GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:15:15:1058:.:.:0
+phiX174 1147 . T . 0 . DP=15;DPB=15;EPPR=4.31318;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=999;RO=15;RPPR=4.31318 GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:15:15:999:.:.:0
+phiX174 1148 . C . 0 . DP=15;DPB=15;EPPR=4.31318;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=1003;RO=15;RPPR=4.31318 GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:15:15:1003:.:.:0
+phiX174 1149 . T C 4.04114e-07 . AB=0;ABP=0;AC=0;AF=0;AN=2;AO=1;CIGAR=1X;DP=14;DPB=14;DPRA=0;EPP=5.18177;EPPR=3.17734;GTI=0;LEN=1;MEANALT=1;MQM=37;MQMR=37;NS=1;NUMALT=1;ODDS=16.1901;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=33;QR=830;RO=13;RPL=0;RPP=5.18177;RPPR=7.18621;RPR=1;RUN=1;SAF=0;SAP=5.18177;SAR=1;SRF=5;SRP=4.51363;SRR=8;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/0:14:13,1:13:830:1:33:0,-1.0599,-39.8344
+phiX174 1151 . C T 3.45599e-08 . AB=0;ABP=0;AC=0;AF=0;AN=2;AO=1;CIGAR=1X;DP=15;DPB=15;DPRA=0;EPP=5.18177;EPPR=8.59409;GTI=0;LEN=1;MEANALT=1;MQM=37;MQMR=37;NS=1;NUMALT=1;ODDS=18.6491;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=27;QR=946;RO=14;RPL=1;RPP=5.18177;RPPR=3.0103;RPR=0;RUN=1;SAF=1;SAP=5.18177;SAR=0;SRF=3;SRP=12.937;SRR=11;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:14,1:14:946:1:27:0,-1.85676,-43.833
+phiX174 1152 . G . 0 . DP=16;DPB=16;EPPR=5.18177;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=993;RO=16;RPPR=5.18177 GT:DP:AD:RO:QR:AO:QA:GL 0/0:16:16:16:993:.:.:0
+phiX174 1153 . T . 3.38668e-15 . DP=15;DPB=15;EPPR=4.31318;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=880;RO=15;RPPR=4.31318 GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:15:15:880:.:.:0
diff -r 000000000000 -r d1116f934d3d test-data/freebayes-phix174-test3.vcf
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/freebayes-phix174-test3.vcf Thu Jan 25 06:52:44 2018 -0500
@@ -0,0 +1,78 @@
+##fileformat=VCFv4.2
+##reference=localref.fa
+##contig=
+##phasing=none
+##commandline="freebayes --region phiX174:0..5386 --bam b_0.bam --fasta-reference localref.fa --vcf ./vcf_output/part_phiX174:0..5386.vcf --haplotype-length 0 --min-alternate-count 1 --min-alternate-fraction 0 --pooled-continuous --report-monomorphic --standard-filters --min-coverage 14"
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT unknown
+phiX174 1134 . A . 0 . DP=14;DPB=14;EPPR=3.0103;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=916;RO=14;RPPR=8.59409 GT:DP:AD:RO:QR:AO:QA:GL 0/0:14:14:14:916:.:.:0
+phiX174 1135 . C . 0 . DP=14;DPB=14;EPPR=3.63072;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=957;RO=14;RPPR=5.49198 GT:DP:AD:RO:QR:AO:QA:GL 0/0:14:14:14:957:.:.:0
+phiX174 1136 . G . 0 . DP=14;DPB=14;EPPR=8.59409;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=919;RO=14;RPPR=3.0103 GT:DP:AD:RO:QR:AO:QA:GL 0/0:14:14:14:919:.:.:0
+phiX174 1137 . C . 0 . DP=14;DPB=14;EPPR=5.49198;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=927;RO=14;RPPR=3.63072 GT:DP:AD:RO:QR:AO:QA:GL 0/0:14:14:14:927:.:.:0
+phiX174 1138 . C . 0 . DP=14;DPB=14;EPPR=3.63072;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=876;RO=14;RPPR=5.49198 GT:DP:AD:RO:QR:AO:QA:GL 0/0:14:14:14:876:.:.:0
+phiX174 1139 . G . 0 . DP=14;DPB=14;EPPR=3.63072;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=946;RO=14;RPPR=5.49198 GT:DP:AD:RO:QR:AO:QA:GL 0/0:14:14:14:946:.:.:0
+phiX174 1140 . T . 0 . DP=15;DPB=15;EPPR=4.31318;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=927;RO=15;RPPR=4.31318 GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:15:15:927:.:.:0
+phiX174 1141 . T G 3.61953e-07 . AB=0;ABP=0;AC=0;AF=0;AN=2;AO=1;CIGAR=1X;DP=15;DPB=15;DPRA=0;EPP=5.18177;EPPR=3.63072;GTI=0;LEN=1;MEANALT=1;MQM=37;MQMR=37;NS=1;NUMALT=1;ODDS=16.3003;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=50;QR=880;RO=14;RPL=0;RPP=5.18177;RPPR=5.49198;RPR=1;RUN=1;SAF=0;SAP=5.18177;SAR=1;SRF=5;SRP=5.49198;SRR=9;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:14,1:14:880:1:50:0,-0.836684,-42.5367
+phiX174 1142 . G . 0 . DP=16;DPB=16;EPPR=5.18177;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=1082;RO=16;RPPR=3.55317 GT:DP:AD:RO:QR:AO:QA:GL 0/0:16:16:16:1082:.:.:0
+phiX174 1143 . G . 0 . DP=16;DPB=16;EPPR=5.18177;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=1029;RO=16;RPPR=3.55317 GT:DP:AD:RO:QR:AO:QA:GL 0/0:16:16:16:1029:.:.:0
+phiX174 1144 . C . 0 . DP=16;DPB=16;EPPR=5.18177;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=1114;RO=16;RPPR=3.55317 GT:DP:AD:RO:QR:AO:QA:GL 0/0:16:16:16:1114:.:.:0
+phiX174 1145 . G . 0 . DP=15;DPB=15;EPPR=6.62942;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=1042;RO=15;RPPR=3.15506 GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:15:15:1042:.:.:0
+phiX174 1146 . C . 0 . DP=15;DPB=15;EPPR=6.62942;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=1058;RO=15;RPPR=3.15506 GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:15:15:1058:.:.:0
+phiX174 1147 . T . 0 . DP=15;DPB=15;EPPR=4.31318;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=999;RO=15;RPPR=4.31318 GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:15:15:999:.:.:0
+phiX174 1148 . C . 0 . DP=15;DPB=15;EPPR=4.31318;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=1003;RO=15;RPPR=4.31318 GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:15:15:1003:.:.:0
+phiX174 1149 . T C 4.04114e-07 . AB=0;ABP=0;AC=0;AF=0;AN=2;AO=1;CIGAR=1X;DP=14;DPB=14;DPRA=0;EPP=5.18177;EPPR=3.17734;GTI=0;LEN=1;MEANALT=1;MQM=37;MQMR=37;NS=1;NUMALT=1;ODDS=16.1901;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=33;QR=830;RO=13;RPL=0;RPP=5.18177;RPPR=7.18621;RPR=1;RUN=1;SAF=0;SAP=5.18177;SAR=1;SRF=5;SRP=4.51363;SRR=8;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/0:14:13,1:13:830:1:33:0,-1.0599,-39.8344
+phiX174 1151 . C T 3.45599e-08 . AB=0;ABP=0;AC=0;AF=0;AN=2;AO=1;CIGAR=1X;DP=15;DPB=15;DPRA=0;EPP=5.18177;EPPR=8.59409;GTI=0;LEN=1;MEANALT=1;MQM=37;MQMR=37;NS=1;NUMALT=1;ODDS=18.6491;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=27;QR=946;RO=14;RPL=1;RPP=5.18177;RPPR=3.0103;RPR=0;RUN=1;SAF=1;SAP=5.18177;SAR=0;SRF=3;SRP=12.937;SRR=11;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:14,1:14:946:1:27:0,-1.85676,-43.833
+phiX174 1152 . G . 0 . DP=16;DPB=16;EPPR=5.18177;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=993;RO=16;RPPR=5.18177 GT:DP:AD:RO:QR:AO:QA:GL 0/0:16:16:16:993:.:.:0
+phiX174 1153 . T . 3.38668e-15 . DP=15;DPB=15;EPPR=4.31318;GTI=0;MQMR=37;NS=1;NUMALT=0;ODDS=0;PAIREDR=0;PQR=0;PRO=0;QR=880;RO=15;RPPR=4.31318 GT:DP:AD:RO:QR:AO:QA:GL 0/0:15:15:15:880:.:.:0
diff -r 000000000000 -r d1116f934d3d test-data/freebayes-phix174-test4.vcf
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/freebayes-phix174-test4.vcf Thu Jan 25 06:52:44 2018 -0500
@@ -0,0 +1,85 @@
+##fileformat=VCFv4.2
+##reference=localref.fa
+##contig=
+##phasing=none
+##commandline="freebayes --region phiX174:0..5386 --bam b_0.bam --fasta-reference localref.fa --vcf ./vcf_output/part_phiX174:0..5386.vcf --theta 0.001 --ploidy 1"
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##INFO=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+##FORMAT=
+#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT unknown
+phiX174 311 . A G 5.72946e-15 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=9;DPB=9;DPRA=0;EPP=7.35324;EPPR=5.80219;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=33.5714;NS=1;NUMALT=1;ODDS=38.0687;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=70;QR=478;RO=7;RPL=0;RPP=7.35324;RPPR=3.32051;RPR=2;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=5;SRP=5.80219;SRR=2;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:9:7,2:7:478:2:70:0,-16.86
+phiX174 374 . T G 22.5297 . AB=0;ABP=0;AC=1;AF=1;AN=1;AO=2;CIGAR=1X;DP=4;DPB=4;DPRA=0;EPP=7.35324;EPPR=7.35324;GTI=0;LEN=1;MEANALT=1;MQM=37;MQMR=25;NS=1;NUMALT=1;ODDS=5.18205;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=75;QR=141;RO=2;RPL=1;RPP=3.0103;RPPR=7.35324;RPR=1;RUN=1;SAF=1;SAP=3.0103;SAR=1;SRF=0;SRP=7.35324;SRR=2;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 1:4:2,2:2:141:2:75:-1.64847,0
+phiX174 913 . A C 0.0178937 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=5;DPB=5;DPRA=0;EPP=3.0103;EPPR=3.73412;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=29;NS=1;NUMALT=1;ODDS=5.4898;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=67;QR=171;RO=3;RPL=1;RPP=3.0103;RPPR=3.73412;RPR=1;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=3;SRP=9.52472;SRR=0;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:5:3,2:3:171:2:67:0,-2.9351
+phiX174 1205 . A C 6.98431e-15 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=10;DPB=10;DPRA=0;EPP=7.35324;EPPR=7.35324;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=34;NS=1;NUMALT=1;ODDS=41.7843;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=67;QR=456;RO=8;RPL=2;RPP=7.35324;RPPR=7.35324;RPR=0;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=2;SRP=7.35324;SRR=6;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:10:8,2:8:456:2:67:0,-19.2238
+phiX174 1245 . G T 2.17569e-10 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=9;DPB=9;DPRA=0;EPP=7.35324;EPPR=3.32051;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=28.4286;NS=1;NUMALT=1;ODDS=23.7171;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=70;QR=389;RO=7;RPL=2;RPP=7.35324;RPPR=3.32051;RPR=0;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=0;SRP=18.2106;SRR=7;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:9:7,2:7:389:2:70:0,-11.7275
+phiX174 1249 . T G 9.08243e-14 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=9;DPB=9;DPRA=0;EPP=3.0103;EPPR=3.32051;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=30.1429;NS=1;NUMALT=1;ODDS=31.5341;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=74;QR=464;RO=7;RPL=1;RPP=3.0103;RPPR=5.80219;RPR=1;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=2;SRP=5.80219;SRR=5;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:9:7,2:7:464:2:74:0,-14.6241
+phiX174 1445 . C A 3.85961e-09 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=7;DPB=7;DPRA=0;EPP=7.35324;EPPR=6.91895;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=32.2;NS=1;NUMALT=1;ODDS=20.8413;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=76;QR=273;RO=5;RPL=1;RPP=3.0103;RPPR=6.91895;RPR=1;RUN=1;SAF=1;SAP=3.0103;SAR=1;SRF=2;SRP=3.44459;SRR=3;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:7:5,2:5:273:2:76:0,-9.96463
+phiX174 1577 . A C 7.60962e-15 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=9;DPB=9;DPRA=0;EPP=7.35324;EPPR=3.32051;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=35.2857;NS=1;NUMALT=1;ODDS=38.7543;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=60;QR=460;RO=7;RPL=1;RPP=3.0103;RPPR=10.7656;RPR=1;RUN=1;SAF=1;SAP=3.0103;SAR=1;SRF=4;SRP=3.32051;SRR=3;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:9:7,2:7:460:2:60:0,-18.0151
+phiX174 1631 . T G 8.14828e-15 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=10;DPB=10;DPRA=0;EPP=7.35324;EPPR=4.09604;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=28;NS=1;NUMALT=1;ODDS=36.3199;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=68;QR=500;RO=8;RPL=0;RPP=7.35324;RPPR=3.0103;RPR=2;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=3;SRP=4.09604;SRR=5;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:10:8,2:8:500:2:68:0,-15.8506
+phiX174 1772 . T G 2.89903e-14 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=10;DPB=10;DPRA=0;EPP=3.0103;EPPR=3.32051;GTI=0;LEN=1;MEANALT=2;MQM=31;MQMR=31.8571;NS=1;NUMALT=1;ODDS=32.7439;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=59;QR=425;RO=7;RPL=1;RPP=3.0103;RPPR=3.32051;RPR=1;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=1;SRP=10.7656;SRR=6;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:10:7,2:7:425:2:59:0,-15.4048
+phiX174 1945 . T G 2.84647e-08 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=6;DPB=6;DPRA=0;EPP=3.0103;EPPR=5.18177;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=37;NS=1;NUMALT=1;ODDS=18.8432;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=59;QR=263;RO=4;RPL=1;RPP=3.0103;RPPR=5.18177;RPR=1;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=2;SRP=3.0103;SRR=2;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:6:4,2:4:263:2:59:0,-8.60945
+phiX174 2230 . T G 1.46203e-15 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=10;DPB=10;DPRA=0;EPP=7.35324;EPPR=3.0103;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=35.5;NS=1;NUMALT=1;ODDS=46.0348;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=80;QR=491;RO=8;RPL=0;RPP=7.35324;RPPR=4.09604;RPR=2;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=3;SRP=4.09604;SRR=5;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:10:8,2:8:491:2:80:0,-20.0697
+phiX174 2699 . C A 0.0286554 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=4;DPB=4;DPRA=0;EPP=7.35324;EPPR=3.0103;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=37;NS=1;NUMALT=1;ODDS=5.01766;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=76;QR=109;RO=2;RPL=1;RPP=3.0103;RPPR=3.0103;RPR=1;RUN=1;SAF=1;SAP=3.0103;SAR=1;SRF=0;SRP=7.35324;SRR=2;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:4:2,2:2:109:2:76:0,-2.17914
+phiX174 2722 . T G 0.000497664 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=5;DPB=5;DPRA=0;EPP=7.35324;EPPR=3.73412;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=33;NS=1;NUMALT=1;ODDS=9.07408;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=67;QR=150;RO=3;RPL=1;RPP=3.0103;RPPR=3.73412;RPR=1;RUN=1;SAF=1;SAP=3.0103;SAR=1;SRF=0;SRP=9.52472;SRR=3;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:5:3,2:3:150:2:67:0,-4.49173
+phiX174 2814 . T G 0.00240428 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=5;DPB=5;DPRA=0;EPP=3.0103;EPPR=3.73412;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=33;NS=1;NUMALT=1;ODDS=7.49878;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=57;QR=135;RO=3;RPL=1;RPP=3.0103;RPPR=3.73412;RPR=1;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=0;SRP=9.52472;SRR=3;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:5:3,2:3:135:2:57:0,-3.80759
+phiX174 2828 . T G 4.53266e-10 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=8;DPB=8;DPRA=0;EPP=7.35324;EPPR=3.0103;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=31;NS=1;NUMALT=1;ODDS=22.9831;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=100;QR=348;RO=6;RPL=0;RPP=7.35324;RPPR=4.45795;RPR=2;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=1;SRP=8.80089;SRR=5;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:8:6,2:6:348:2:100:0,-10.3385
+phiX174 2983 . T G 0.000276493 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=5;DPB=5;DPRA=0;EPP=3.0103;EPPR=3.73412;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=33;NS=1;NUMALT=1;ODDS=9.66185;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=57;QR=150;RO=3;RPL=1;RPP=3.0103;RPPR=3.73412;RPR=1;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=2;SRP=3.73412;SRR=1;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:5:3,2:3:150:2:57:0,-4.26987
+phiX174 3110 . T C 4.70002e-12 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=3;CIGAR=1X;DP=9;DPB=9;DPRA=0;EPP=3.73412;EPPR=4.45795;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=35;NS=1;NUMALT=1;ODDS=27.5529;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=94;QR=361;RO=6;RPL=0;RPP=9.52472;RPPR=3.0103;RPR=3;RUN=1;SAF=1;SAP=3.73412;SAR=2;SRF=1;SRP=8.80089;SRR=5;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:9:6,3:6:361:3:94:0,-12.3743
+phiX174 3155 . T G 1.16474e-06 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=7;DPB=7;DPRA=0;EPP=3.0103;EPPR=3.44459;GTI=0;LEN=1;MEANALT=1;MQM=37;MQMR=29.8;NS=1;NUMALT=1;ODDS=15.1316;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=75;QR=299;RO=5;RPL=1;RPP=3.0103;RPPR=3.44459;RPR=1;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=4;SRP=6.91895;SRR=1;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:7:5,2:5:299:2:75:0,-7.18392
+phiX174 3325 . A C 2.80881e-07 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=6;DPB=6;DPRA=0;EPP=7.35324;EPPR=5.18177;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=31;NS=1;NUMALT=1;ODDS=16.5539;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=58;QR=284;RO=4;RPL=2;RPP=7.35324;RPPR=3.0103;RPR=0;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=3;SRP=5.18177;SRR=1;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:6:4,2:4:284:2:58:0,-7.01317
+phiX174 3418 . A C 0 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=10;DPB=10;DPRA=0;EPP=7.35324;EPPR=3.0103;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=34;NS=1;NUMALT=1;ODDS=44.7759;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=60;QR=526;RO=8;RPL=2;RPP=7.35324;RPPR=4.09604;RPR=0;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=5;SRP=4.09604;SRR=3;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:10:8,2:8:526:2:60:0,-19.523
+phiX174 3729 . C T 6.64464e-08 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=6;DPB=6;DPRA=0;EPP=7.35324;EPPR=5.18177;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=37;NS=1;NUMALT=1;ODDS=17.9954;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=63;QR=232;RO=4;RPL=2;RPP=7.35324;RPPR=5.18177;RPR=0;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=0;SRP=11.6962;SRR=4;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:6:4,2:4:232:2:63:0,-8.41737
+phiX174 4031 . T G 2.20192e-08 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=8;DPB=8;DPRA=0;EPP=7.35324;EPPR=4.45795;GTI=0;LEN=1;MEANALT=1;MQM=31;MQMR=27;NS=1;NUMALT=1;ODDS=19.0999;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=85;QR=322;RO=6;RPL=2;RPP=7.35324;RPPR=8.80089;RPR=0;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=3;SRP=3.0103;SRR=3;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:8:6,2:6:322:2:85:0,-8.65207
+phiX174 4502 . A C 3.4202e-08 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=6;DPB=6;DPRA=0;EPP=3.0103;EPPR=5.18177;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=37;NS=1;NUMALT=1;ODDS=18.6595;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=55;QR=284;RO=4;RPL=1;RPP=3.0103;RPPR=11.6962;RPR=1;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=3;SRP=5.18177;SRR=1;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:6:4,2:4:284:2:55:0,-9.30785
+phiX174 4558 . C G 1.36745e-05 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=6;DPB=6;DPRA=0;EPP=7.35324;EPPR=3.0103;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=28;NS=1;NUMALT=1;ODDS=12.6685;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=59;QR=204;RO=4;RPL=2;RPP=7.35324;RPPR=11.6962;RPR=0;RUN=1;SAF=2;SAP=7.35324;SAR=0;SRF=2;SRP=3.0103;SRR=2;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:6:4,2:4:204:2:59:0,-5.75175
+phiX174 4655 . T G 1.44134e-07 . AB=0;ABP=0;AC=0;AF=0;AN=1;AO=2;CIGAR=1X;DP=6;DPB=6;DPRA=0;EPP=3.0103;EPPR=3.0103;GTI=0;LEN=1;MEANALT=1;MQM=25;MQMR=37;NS=1;NUMALT=1;ODDS=17.2211;PAIRED=0;PAIREDR=0;PAO=0;PQA=0;PQR=0;PRO=0;QA=70;QR=193;RO=4;RPL=1;RPP=3.0103;RPPR=3.0103;RPR=1;RUN=1;SAF=0;SAP=7.35324;SAR=2;SRF=0;SRP=11.6962;SRR=4;TYPE=snp GT:DP:AD:RO:QR:AO:QA:GL 0:6:4,2:4:193:2:70:0,-8.33095
diff -r 000000000000 -r d1116f934d3d test-data/freebayes-phix174.bam
Binary file test-data/freebayes-phix174.bam has changed
diff -r 000000000000 -r d1116f934d3d test-data/freebayes-phix174.fasta
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/freebayes-phix174.fasta Thu Jan 25 06:52:44 2018 -0500
@@ -0,0 +1,2 @@
+>phiX174
+GAGTTTTATCGCTTCCATGACGCAGAAGTTAACACTTTCGGATATTTCTGATGAGTCGAAAAATTATCTTGATAAAGCAGGAATTACTACTGCTTGTTTACGAATTAAATCGAAGTGGACTGCTGGCGGAAAATGAGAAAATTCGACCTATCCTTGCGCAGCTCGAGAAGCTCTTACTTTGCGACCTTTCGCCATCAACTAACGATTCTGTCAAAAACTGACGCGTTGGATGAGGAGAAGTGGCTTAATATGCTTGGCACGTTCGTCAAGGACTGGTTTAGATATGAGTCACATTTTGTTCATGGTAGAGATTCTCTTGTTGACATTTTAAAAGAGCGTGGATTACTATCTGAGTCCGATGCTGTTCAACCACTAATAGGTAAGAAATCATGAGTCAAGTTACTGAACAATCCGTACGTTTCCAGACCGCTTTGGCCTCTATTAAGCTCATTCAGGCTTCTGCCGTTTTGGATTTAACCGAAGATGATTTCGATTTTCTGACGAGTAACAAAGTTTGGATTGCTACTGACCGCTCTCGTGCTCGTCGCTGCGTTGAGGCTTGCGTTTATGGTACGCTGGACTTTGTGGGATACCCTCGCTTTCCTGCTCCTGTTGAGTTTATTGCTGCCGTCATTGCTTATTATGTTCATCCCGTCAACATTCAAACGGCCTGTCTCATCATGGAAGGCGCTGAATTTACGGAAAACATTATTAATGGCGTCGAGCGTCCGGTTAAAGCCGCTGAATTGTTCGCGTTTACCTTGCGTGTACGCGCAGGAAACACTGACGTTCTTACTGACGCAGAAGAAAACGTGCGTCAAAAATTACGTGCaGAAGGAGTGATGTAATGTCTAAAGGTAAAAAACGTTCTGGCGCTCGCCCTGGTCGTCCGCAGCCGTTGCGAGGTACTAAAGGCAAGCGTAAAGGCGCTCGTCTTTGGTATGTAGGTGGTCAACAATTTTAATTGCAGGGGCTTCGGCCCCTTACTTGAGGATAAATTATGTCTAATATTCAAACTGGCGCCGAGCGTATGCCGCATGACCTTTCCCATCTTGGCTTCCTTGCTGGTCAGATTGGTCGTCTTATTACCATTTCAACTACTCCGGTTATCGCTGGCGACTCCTTCGAGATGGACGCCGTTGGCGCTCTCCGTCTTTCTCCATTGCGTCGTGGCCTTGCTATTGACTCTACTGTAGACATTTTTACTTTTTATGTCCCTCATCGTCACGTTTATGGTGAACAGTGGATTAAGTTCATGAAGGATGGTGTTAATGCCACTCCTCTCCCGACTGTTAACACTACTGGTTATATTGACCATGCCGCTTTTCTTGGCACGATTAACCCTGATACCAATAAAATCCCTAAGCATTTGTTTCAGGGTTATTTGAATATCTATAACAACTATTTTAAAGCGCCGTGGATGCCTGACCGTACCGAGGCTAACCCTAATGAGCTTAATCAAGATGATGCTCGTTATGGTTTCCGTTGCTGCCATCTCAAAAACATTTGGACTGCTCCGCTTCCTCCTGAGACTGAGCTTTCTCGCCAAATGACGACTTCTACCACATCTATTGACATTATGGGTCTGCAAGCTGCTTATGCTAATTTGCATACTGACCAAGAACGTGATTACTTCATGCAGCGTTACCgTGATGTTATTTCTTCATTTGGAGGTAAAACCTCTTATGACGCTGACAACCGTCCTTTACTTGTCATGCGCTCTAATCTCTGGGCATCTGGCTATGATGTTGATGGAACTGACCAAACGTCGTTAGGCCAGTTTTCTGGTCGTGTTCAACAGACCTATAAACATTCTGTGCCGCGTTTCTTTGTTCCTGAGCATGGCACTATGTTTACTCTTGCGCTTGTTCGTTTTCCGCCTACTGCGACTAAAGAGATTCAGTACCTTAACGCTAAAGGTGCTTTGACTTATACCGATATTGCTGGCGACCCTGTTTTGTATGGCAACTTGCCGCCGCGTGAAATTTCTATGAAGGATGTTTTCCGTTCTGGTGATTCGTCTAAGAAGTTTAAGATTGCTGAGGGTCAGTGGTATCGTTATGCGCCTTCGTATGTTTCTCCTGCTTATCACCTTCTTGAAGGCTTCCCATTCATTCAGGAACCGCCTTCTGGTGATTTGCAAGAACGCGTACTTATTCGCCACCATGATTATGACCAGTGTTTCCAGTCCGTTCAGTTGTTGCAGTGGAATAGTCAGGTTAAATTTAATGTGACCGTTTATCGCAATCTGCCGACCACTCGCGATTCAATCATGACTTCGTGATAAAAGATTGAGTGTGAGGTTATAACGCCGAAGCGGTAAAAATTTTAATTTTTGCCGCTGAGGGGTTGACCAAGCGAAGCGCGGTAGGTTTTCTGCTTAGGAGTTTAATCATGTTTCAGACTTTTATTTCTCGCCATAATTCAAACTTTTTTTCTGATAAGCTGGTTCTCACTTCTGTTACTCCAGCTTCTTCGGCACCTGTTTTACAGACACCTAAAGCTACATCGTCAACGTTATATTTTGATAGTTTGACGGTTAATGCTGGTAATGGTGGTTTTCTTCATTGCATTCAGATGGATACATCTGTCAACGCCGCTAATCAGGTTGTTTCTGTTGGTGCTGATATTGCTTTTGATGCCGACCCTAAATTTTTTGCCTGTTTGGTTCGCTTTGAGTCTTCTTCGGTTCCGACTACCCTCCCGACTGCCTATGATGTTTATCCTTTGAATGGTCGCCATGATGGTGGTTATTATACCGTCAAGGACTGTGTGACTATTGACGTCCTTCCCCGTACGCCGGGCAATAAtGTTTATGTTGGTTTCATGGTTTGGTCTAACTTTACCGCTACTAAATGCCGCGGATTGGTTTCGCTGAATCAGGTTATTAAAGAGATTATTTGTCTCCAGCCACTTAAGTGAGGTGATTTATGTTTGGTGCTATTGCTGGCGGTATTGCTTCTGCTCTTGCTGGTGGCGCCATGTCTAAATTGTTTGGAGGCGGTCAAAAAGCCGCCTCCGGTGGCATTCAAGGTGATGTGCTTGCTACCGATAACAATACTGTAGGCATGGGTGATGCTGGTATTAAATCTGCCATTCAAGGCTCTAATGTTCCTAACCCTGATGAGGCCGCCCCTAGTTTTGTTTCTGGTGCTATGGCTAAAGCTGGTAAAGGACTTCTTGAAGGTACGTTGCAGGCTGGCACTTCTGCCGTTTCTGATAAGTTGCTTGATTTGGTTGGACTTGGTGGCAAGTCTGCCGCTGATAAAGGAAAGGATACTCGTGATTATCTTGCTGCTGCATTTCCTGAGCTTAATGCTTGGGAGCGTGCTGGTGCTGATGCTTCCTCTGCTGGTATGGTTGACGCCGGATTTGAGAATCAAAAAGAGCTTACTAAAATGCAACTGGACAATCAGAAAGAGATTGCCGAGATGCAAAATGAGACTCAAAAAGAGATTGCTGGCATTCAGTCGGCGACTTCACGCCAGAATACGAAAGACCAGGTATATGCACAAAATGAGATGCTTGCTTATCAACAGAAGGAGTCTACTGCTCGCGTTGCGTCTATTATGGAAAACACCAATCTTTCCAAGCAACAGCAGGTTTCCGAGATTATGCGCCAAATGCTTACTCAAGCTCAAACGGCTGGTCAGTATTTTACCAATGACCAAATCAAAGAAATGACTCGCAAGGTTAGTGCTGAGGTTGACTTAGTTCATCAGCAAACGCAGAATCAGCGGTATGGCTCTTCTCATATTGGCGCTACTGCAAAGGATATTTCTAATGTCGTCACTGATGCTGCTTCTGGTGTGGTTGATATTTTTCATGGTATTGATAAAGCTGTTGCCGATACTTGGAACAATTTCTGGAAAGACGGTAAAGCTGATGGTATTGGCTCTAATTTGTCTAGGAAATAACCGTCAGGATTGACACCCTCCCAATTGTATGTTTTCATGCCTCCAAATCTTGGAGGCTTTTTTATGGTTCGTTCTTATTACCCTTCTGAATGTCACGCTGATTATTTTGACTTTGAGCGTATCGAGGCTCTTAAACCTGCTATTGAGGCTTGTGGCATTTCTACTCTTTCTCAATCCCCAATGCTTGGCTTCCATAAGCAGATGGATAACCGCATCAAGCTCTTGGAAGAGATTCTGTCTTTTCGTATGCAGGGCGTTGAGTTCGATAATGGTGATATGTATGTTGACGGCCATAAGGCTGCTTCTGACGTTCGTGATGAGTTTGTATCTGTTACTGAGAAGTTAATGGATGAATTGGCACAATGCTACAATGTGCTCCCCCAACTTGATATTAATAACACTATAGACCACCGCCCCGAAGGGGACGAAAAATGGTTTTTAGAGAACGAGAAGACGGTTACGCAGTTTTGCCGCAAGCTGGCTGCTGAACGCCCTCTTAAGGATATTCGCGATGAGTATAATTACCCCAAAAAGAAAGGTATTAAGGATGAGTGTTCAAGATTGCTGGAGGCCTCCACTATGAAATCGCGTAGAGGCTTTaCTATTCAGCGTTTGATGAATGCAATGCGACAGGCTCATGCTGATGGTTGGTTTATCGTTTTTGACACTCTCACGTTGGCTGACGACCGATTAGAGGCGTTTTATGATAATCCCAATGCTTTGCGTGACTATTTTCGTGATATTGGTCGTATGGTTCTTGCTGCCGAGGGTCGCAAGGCTAATGATTCACACGCCGACTGCTATCAGTATTTTTGTGTGCCTGAGTATGGTACAGCTAATGGCCGTCTTCATTTCCATGCGGTGCAtTTTATGCGGACACTTCCTACAGGTAGCGTTGACCCTAATTTTGGTCGTCGGGTACGCAATCGCCGCCAGTTAAATAGCTTGCAAAATACGTGGCCTTATGGTTACAGTATGCCCATCGCAGTTCGCTACACGCAGGACGCTTTTTCACGTTCTGGTTGGTTGTGGCCTGTTGATGCTAAAGGTGAGCCGCTTAAAGCTACCAGTTATATGGCTGTTGGTTTCTATGTGGCTAAATACGTTAACAAAAAGTCAGATATGGACCTTGCTGCTAAAGGTCTAGGAGCTAAAGAATGGAACAACTCACTAAAAACCAAGCTGTCGCTACTTCCCAAGAAGCTGTTCAGAATCAGAATGAGCCGCAACTTCGGGATGAAAATGCTCACAATGACAAATCTGTCCACGGAGTGCTTAATCCAACTTACCAAGCTGGGTTACGACGCGACGCCGTTCAACCAGATATTGAAGCAGAACGCAAAAAGAGAGATGAGATTGAGGCTGGGAAAAGTTACTGTAGCCGACGTTTTGGCGGCGCAACCTGTGACGACAAATCTGCTCAAATTTATGCGCGCTTCGATAAAAATGATTGGCGTATCCAACCTGCA
diff -r 000000000000 -r d1116f934d3d test-data/left-align-input.bam
Binary file test-data/left-align-input.bam has changed
diff -r 000000000000 -r d1116f934d3d test-data/left-align-output.bam
Binary file test-data/left-align-output.bam has changed
diff -r 000000000000 -r d1116f934d3d test-data/leftalign.fa
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/leftalign.fa Thu Jan 25 06:52:44 2018 -0500
@@ -0,0 +1,2 @@
+>phiX174
+GAGTTTTATCGCTTCCATGACGCAGAAGTTAACACTTTCGGATATTTCTGATGAGTCGAAAAATTATCTTGATAAAGCAGGAATTACTACTGCTTGTTTACGAATTAAATCGAAGTGGACTGCTGGCGGAAAATGAGAAAATTCGACCTATCCTTGCGCAGCTCGAGAAGCTCTTACTTTGCGACCTTTCGCCATCAACTAACGATTCTGTCAAAAACTGACGCGTTGGATGAGGAGAAGTGGCTTAATATGCTTGGCACGTTCGTCAAGGACTGGTTTAGATATGAGTCACATTTTGTTCATGGTAGAGATTCTCTTGTTGACATTTTAAAAGAGCGTGGATTACTATCTGAGTCCGATGCTGTTCAACCACTAATAGGTAAGAAATCATGAGTCAAGTTACTGAACAATCCGTACGTTTCCAGACCGCTTTGGCCTCTATTAAGCTCATTCAGGCTTCTGCCGTTTTGGATTTAACCGAAGATGATTTCGATTTTCTGACGAGTAACAAAGTTTGGATTGCTACTGACCGCTCTCGTGCTCGTCGCTGCGTTGAGGCTTGCGTTTATGGTACGCTGGACTTTGTGGGATACCCTCGCTTTCCTGCTCCTGTTGAGTTTATTGCTGCCGTCATTGCTTATTATGTTCATCCCGTCAACATTCAAACGGCCTGTCTCATCATGGAAGGCGCTGAATTTACGGAAAACATTATTAATGGCGTCGAGCGTCCGGTTAAAGCCGCTGAATTGTTCGCGTTTACCTTGCGTGTACGCGCAGGAAACACTGACGTTCTTACTGACGCAGAAGAAAACGTGCGTCAAAAATTACGTGCaGAAGGAGTGATGTAATGTCTAAAGGTAAAAAACGTTCTGGCGCTCGCCCTGGTCGTCCGCAGCCGTTGCGAGGTACTAAAGGCAAGCGTAAAGGCGCTCGTCTTTGGTATGTAGGTGGTCAACAATTTTAATTGCAGGGGCTTCGGCCCCTTACTTGAGGATAAATTATGTCTAATATTCAAACTGGCGCCGAGCGTATGCCGCATGACCTTTCCCATCTTGGCTTCCTTGCTGGTCAGATTGGTCGTCTTATTACCATTTCAACTACTCCGGTTATCGCTGGCGACTCCTTCGAGATGGACGCCGTTGGCGCTCTCCGTCTTTCTCCATTGCGTCGTGGCCTTGCTATTGACTCTACTGTAGACATTTTTACTTTTTATGTCCCTCATCGTCACGTTTATGGTGAACAGTGGATTAAGTTCATGAAGGATGGTGTTAATGCCACTCCTCTCCCGACTGTTAACACTACTGGTTATATTGACCATGCCGCTTTTCTTGGCACGATTAACCCTGATACCAATAAAATCCCTAAGCATTTGTTTCAGGGTTATTTGAATATCTATAACAACTATTTTAAAGCGCCGTGGATGCCTGACCGTACCGAGGCTAACCCTAATGAGCTTAATCAAGATGATGCTCGTTATGGTTTCCGTTGCTGCCATCTCAAAAACATTTGGACTGCTCCGCTTCCTCCTGAGACTGAGCTTTCTCGCCAAATGACGACTTCTACCACATCTATTGACATTATGGGTCTGCAAGCTGCTTATGCTAATTTGCATACTGACCAAGAACGTGATTACTTCATGCAGCGTTACCgTGATGTTATTTCTTCATTTGGAGGTAAAACCTCTTATGACGCTGACAACCGTCCTTTACTTGTCATGCGCTCTAATCTCTGGGCATCTGGCTATGATGTTGATGGAACTGACCAAACGTCGTTAGGCCAGTTTTCTGGTCGTGTTCAACAGACCTATAAACATTCTGTGCCGCGTTTCTTTGTTCCTGAGCATGGCACTATGTTTACTCTTGCGCTTGTTCGTTTTCCGCCTACTGCGACTAAAGAGATTCAGTACCTTAACGCTAAAGGTGCTTTGACTTATACCGATATTGCTGGCGACCCTGTTTTGTATGGCAACTTGCCGCCGCGTGAAATTTCTATGAAGGATGTTTTCCGTTCTGGTGATTCGTCTAAGAAGTTTAAGATTGCTGAGGGTCAGTGGTATCGTTATGCGCCTTCGTATGTTTCTCCTGCTTATCACCTTCTTGAAGGCTTCCCATTCATTCAGGAACCGCCTTCTGGTGATTTGCAAGAACGCGTACTTATTCGCCACCATGATTATGACCAGTGTTTCCAGTCCGTTCAGTTGTTGCAGTGGAATAGTCAGGTTAAATTTAATGTGACCGTTTATCGCAATCTGCCGACCACTCGCGATTCAATCATGACTTCGTGATAAAAGATTGAGTGTGAGGTTATAACGCCGAAGCGGTAAAAATTTTAATTTTTGCCGCTGAGGGGTTGACCAAGCGAAGCGCGGTAGGTTTTCTGCTTAGGAGTTTAATCATGTTTCAGACTTTTATTTCTCGCCATAATTCAAACTTTTTTTCTGATAAGCTGGTTCTCACTTCTGTTACTCCAGCTTCTTCGGCACCTGTTTTACAGACACCTAAAGCTACATCGTCAACGTTATATTTTGATAGTTTGACGGTTAATGCTGGTAATGGTGGTTTTCTTCATTGCATTCAGATGGATACATCTGTCAACGCCGCTAATCAGGTTGTTTCTGTTGGTGCTGATATTGCTTTTGATGCCGACCCTAAATTTTTTGCCTGTTTGGTTCGCTTTGAGTCTTCTTCGGTTCCGACTACCCTCCCGACTGCCTATGATGTTTATCCTTTGAATGGTCGCCATGATGGTGGTTATTATACCGTCAAGGACTGTGTGACTATTGACGTCCTTCCCCGTACGCCGGGCAATAAtGTTTATGTTGGTTTCATGGTTTGGTCTAACTTTACCGCTACTAAATGCCGCGGATTGGTTTCGCTGAATCAGGTTATTAAAGAGATTATTTGTCTCCAGCCACTTAAGTGAGGTGATTTATGTTTGGTGCTATTGCTGGCGGTATTGCTTCTGCTCTTGCTGGTGGCGCCATGTCTAAATTGTTTGGAGGCGGTCAAAAAGCCGCCTCCGGTGGCATTCAAGGTGATGTGCTTGCTACCGATAACAATACTGTAGGCATGGGTGATGCTGGTATTAAATCTGCCATTCAAGGCTCTAATGTTCCTAACCCTGATGAGGCCGCCCCTAGTTTTGTTTCTGGTGCTATGGCTAAAGCTGGTAAAGGACTTCTTGAAGGTACGTTGCAGGCTGGCACTTCTGCCGTTTCTGATAAGTTGCTTGATTTGGTTGGACTTGGTGGCAAGTCTGCCGCTGATAAAGGAAAGGATACTCGTGATTATCTTGCTGCTGCATTTCCTGAGCTTAATGCTTGGGAGCGTGCTGGTGCTGATGCTTCCTCTGCTGGTATGGTTGACGCCGGATTTGAGAATCAAAAAGAGCTTACTAAAATGCAACTGGACAATCAGAAAGAGATTGCCGAGATGCAAAATGAGACTCAAAAAGAGATTGCTGGCATTCAGTCGGCGACTTCACGCCAGAATACGAAAGACCAGGTATATGCACAAAATGAGATGCTTGCTTATCAACAGAAGGAGTCTACTGCTCGCGTTGCGTCTATTATGGAAAACACCAATCTTTCCAAGCAACAGCAGGTTTCCGAGATTATGCGCCAAATGCTTACTCAAGCTCAAACGGCTGGTCAGTATTTTACCAATGACCAAATCAAAGAAATGACTCGCAAGGTTAGTGCTGAGGTTGACTTAGTTCATCAGCAAACGCAGAATCAGCGGTATGGCTCTTCTCATATTGGCGCTACTGCAAAGGATATTTCTAATGTCGTCACTGATGCTGCTTCTGGTGTGGTTGATATTTTTCATGGTATTGATAAAGCTGTTGCCGATACTTGGAACAATTTCTGGAAAGACGGTAAAGCTGATGGTATTGGCTCTAATTTGTCTAGGAAATAACCGTCAGGATTGACACCCTCCCAATTGTATGTTTTCATGCCTCCAAATCTTGGAGGCTTTTTTATGGTTCGTTCTTATTACCCTTCTGAATGTCACGCTGATTATTTTGACTTTGAGCGTATCGAGGCTCTTAAACCTGCTATTGAGGCTTGTGGCATTTCTACTCTTTCTCAATCCCCAATGCTTGGCTTCCATAAGCAGATGGATAACCGCATCAAGCTCTTGGAAGAGATTCTGTCTTTTCGTATGCAGGGCGTTGAGTTCGATAATGGTGATATGTATGTTGACGGCCATAAGGCTGCTTCTGACGTTCGTGATGAGTTTGTATCTGTTACTGAGAAGTTAATGGATGAATTGGCACAATGCTACAATGTGCTCCCCCAACTTGATATTAATAACACTATAGACCACCGCCCCGAAGGGGACGAAAAATGGTTTTTAGAGAACGAGAAGACGGTTACGCAGTTTTGCCGCAAGCTGGCTGCTGAACGCCCTCTTAAGGATATTCGCGATGAGTATAATTACCCCAAAAAGAAAGGTATTAAGGATGAGTGTTCAAGATTGCTGGAGGCCTCCACTATGAAATCGCGTAGAGGCTTTaCTATTCAGCGTTTGATGAATGCAATGCGACAGGCTCATGCTGATGGTTGGTTTATCGTTTTTGACACTCTCACGTTGGCTGACGACCGATTAGAGGCGTTTTATGATAATCCCAATGCTTTGCGTGACTATTTTCGTGATATTGGTCGTATGGTTCTTGCTGCCGAGGGTCGCAAGGCTAATGATTCACACGCCGACTGCTATCAGTATTTTTGTGTGCCTGAGTATGGTACAGCTAATGGCCGTCTTCATTTCCATGCGGTGCAtTTTATGCGGACACTTCCTACAGGTAGCGTTGACCCTAATTTTGGTCGTCGGGTACGCAATCGCCGCCAGTTAAATAGCTTGCAAAATACGTGGCCTTATGGTTACAGTATGCCCATCGCAGTTCGCTACACGCAGGACGCTTTTTCACGTTCTGGTTGGTTGTGGCCTGTTGATGCTAAAGGTGAGCCGCTTAAAGCTACCAGTTATATGGCTGTTGGTTTCTATGTGGCTAAATACGTTAACAAAAAGTCAGATATGGACCTTGCTGCTAAAGGTCTAGGAGCTAAAGAATGGAACAACTCACTAAAAACCAAGCTGTCGCTACTTCCCAAGAAGCTGTTCAGAATCAGAATGAGCCGCAACTTCGGGATGAAAATGCTCACAATGACAAATCTGTCCACGGAGTGCTTAATCCAACTTACCAAGCTGGGTTACGACGCGACGCCGTTCAACCAGATATTGAAGCAGAACGCAAAAAGAGAGATGAGATTGAGGCTGGGAAAAGTTACTGTAGCCGACGTTTTGGCGGCGCAACCTGTGACGACAAATCTGCTCAAATTTATGCGCGCTTCGATAAAAATGATTGGCGTATCCAACCTGCA
diff -r 000000000000 -r d1116f934d3d tool-data/fasta_indexes.loc.sample
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/tool-data/fasta_indexes.loc.sample Thu Jan 25 06:52:44 2018 -0500
@@ -0,0 +1,29 @@
+#This is a sample file distributed with Galaxy that enables tools
+#to use a directory of Samtools indexed sequences data files. You will need
+#to create these data files and then create a fasta_indexes.loc file
+#similar to this one (store it in this directory) that points to
+#the directories in which those files are stored. The fasta_indexes.loc
+#file has this format (white space characters are TAB characters):
+#
+#
+#
+#So, for example, if you had hg19 Canonical indexed stored in
+#
+# /depot/data2/galaxy/hg19/sam/,
+#
+#then the fasta_indexes.loc entry would look like this:
+#
+#hg19canon hg19 Human (Homo sapiens): hg19 Canonical /depot/data2/galaxy/hg19/sam/hg19canon.fa
+#
+#and your /depot/data2/galaxy/hg19/sam/ directory
+#would contain hg19canon.fa and hg19canon.fa.fai files.
+#
+#Your fasta_indexes.loc file should include an entry per line for
+#each index set you have stored. The file in the path does actually
+#exist, but it should never be directly used. Instead, the name serves
+#as a prefix for the index file. For example:
+#
+#hg18canon hg18 Human (Homo sapiens): hg18 Canonical /depot/data2/galaxy/hg18/sam/hg18canon.fa
+#hg18full hg18 Human (Homo sapiens): hg18 Full /depot/data2/galaxy/hg18/sam/hg18full.fa
+#hg19canon hg19 Human (Homo sapiens): hg19 Canonical /depot/data2/galaxy/hg19/sam/hg19canon.fa
+#hg19full hg19 Human (Homo sapiens): hg19 Full /depot/data2/galaxy/hg19/sam/hg19full.fa
diff -r 000000000000 -r d1116f934d3d tool_data_table_conf.xml.sample
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_data_table_conf.xml.sample Thu Jan 25 06:52:44 2018 -0500
@@ -0,0 +1,8 @@
+
+
+
+
+ value, dbkey, name, path
+
+
+