]>
git.donarmstrong.com Git - samtools.git/log
Heng Li [Tue, 21 Dec 2010 17:16:33 +0000 (17:16 +0000)]
* samtools-0.1.12-9 (r891)
* allow to call SNPs from a subset of samples
Heng Li [Wed, 15 Dec 2010 16:28:16 +0000 (16:28 +0000)]
* samtools-0.1.12-12 (r889)
* set mapQ as 20 if it equals 255
Heng Li [Wed, 15 Dec 2010 03:41:09 +0000 (03:41 +0000)]
When -B is applied to mpileup, still use paired reads only unless -A is flagged.
Heng Li [Wed, 15 Dec 2010 03:37:05 +0000 (03:37 +0000)]
* samtools-0.1.12-6 (r887)
* added a hidden option -E to mpileup/calmd. -E triggers an alternative way to apply BAQ.
Heng Li [Tue, 14 Dec 2010 17:51:03 +0000 (17:51 +0000)]
(Arguably) improved the indel caller a tiny bit for lowCov data.
Petr Danecek [Tue, 14 Dec 2010 09:55:46 +0000 (09:55 +0000)]
Fixed the VCF header to pass validation
Heng Li [Mon, 13 Dec 2010 04:02:19 +0000 (04:02 +0000)]
* samtools-0.1.12-4 (r884)
* fixed a long-existing flaw in the INDEL calling model
Heng Li [Sun, 12 Dec 2010 01:05:42 +0000 (01:05 +0000)]
compute max SP and max GQ from sample genotypes
Heng Li [Fri, 10 Dec 2010 15:50:54 +0000 (15:50 +0000)]
drop bcf-fix.pl as it is redundant by the latest changes
Heng Li [Fri, 10 Dec 2010 15:50:29 +0000 (15:50 +0000)]
* fixed a minor issue in printing VCFs
* write bcftools specific INFO and FORMAT in the header
Heng Li [Fri, 10 Dec 2010 15:09:14 +0000 (15:09 +0000)]
Make sure that the GT genotype field is the first
Heng Li [Wed, 8 Dec 2010 22:27:05 +0000 (22:27 +0000)]
* samtools-0.1.12-2 (r877)
* allow to fine control the selection of indel candidates. The current
setting is okay for lowCov and highCov with ~100 samples, but it
skips too many indels for highCov with >250 samples.
Heng Li [Wed, 8 Dec 2010 03:40:35 +0000 (03:40 +0000)]
a spelling error..
Heng Li [Wed, 8 Dec 2010 03:39:57 +0000 (03:39 +0000)]
* samtools-0.1.12-1 (r873)
* added a switch to allow anomalous read pairs in calling
Heng Li [Tue, 7 Dec 2010 19:43:54 +0000 (19:43 +0000)]
fixed a bug in vcf2fq
Heng Li [Sun, 5 Dec 2010 06:18:06 +0000 (06:18 +0000)]
added a warning for the Windows version
Heng Li [Sun, 5 Dec 2010 06:05:51 +0000 (06:05 +0000)]
In ksprintf(), change "%lf" and "%lg" to "%f" and "%g", respectively.
According to the manual page, this change is valid. However, MinGW seems
to interpret "%lf" as "%Lf".
Heng Li [Sun, 5 Dec 2010 05:35:43 +0000 (05:35 +0000)]
bring back the windows support
Heng Li [Sun, 5 Dec 2010 04:33:51 +0000 (04:33 +0000)]
Fixed a compiling error when knetfile is not used.
Heng Li [Sat, 4 Dec 2010 05:13:22 +0000 (05:13 +0000)]
vcf->fastq
Heng Li [Fri, 3 Dec 2010 22:12:30 +0000 (22:12 +0000)]
* remove "-f". Instead always compute consensus quality
* increase the upper limit of quality
Heng Li [Fri, 3 Dec 2010 20:28:15 +0000 (20:28 +0000)]
more informative error message
Heng Li [Thu, 2 Dec 2010 21:16:08 +0000 (21:16 +0000)]
Release samtools-0.1.12a
Heng Li [Thu, 2 Dec 2010 20:55:06 +0000 (20:55 +0000)]
a possible fix to DP4=0,0,0,0; have not tested, but should have no side-effect
Heng Li [Thu, 2 Dec 2010 16:39:57 +0000 (16:39 +0000)]
Release samtools-0.1.12
Heng Li [Thu, 2 Dec 2010 16:24:41 +0000 (16:24 +0000)]
* samtools-0.1.11-1 (r858)
* fixed a bug in mpileup which causes segfaults
* bcftools: do not segfault when BCF contains errors
Heng Li [Wed, 1 Dec 2010 04:52:50 +0000 (04:52 +0000)]
fixed a memory leak in bam_fetch()
Heng Li [Fri, 26 Nov 2010 05:07:31 +0000 (05:07 +0000)]
* fixed a memory violation
* added splitchr to vcfutils.pl
Heng Li [Tue, 23 Nov 2010 14:05:08 +0000 (14:05 +0000)]
fixed a typo/bug in r^2 computation
Heng Li [Mon, 22 Nov 2010 03:20:20 +0000 (03:20 +0000)]
forget to change the version information
Heng Li [Mon, 22 Nov 2010 03:16:52 +0000 (03:16 +0000)]
Release samtools-0.1.11
Heng Li [Sat, 20 Nov 2010 04:16:08 +0000 (04:16 +0000)]
* samtools-0.1.10-9 (r844)
* added the "folded" or reference-free mode for variant calling
Heng Li [Sat, 20 Nov 2010 03:26:36 +0000 (03:26 +0000)]
In merging, if -R is specified, do not abort if the sequence dictionary is different.
On behalf of John Marshall [Sat, 20 Nov 2010 02:24:28 +0000 (02:24 +0000)]
When merging BAM headers, compare the list of target reference sequences
strictly (and fail/abort if there is a mismatch), but allow one list to be a
prefix of the other. (i.e., check that the lists are identical up until the
shorter runs out, and add the excess targets from the longer to the output.)
Heng Li [Fri, 19 Nov 2010 19:49:27 +0000 (19:49 +0000)]
* samtools-0.1.10 (r841)
* fixed a bug in pileup when the first CIGAR operation is D
* fixed a bug in view with range query
Heng Li [Fri, 19 Nov 2010 18:45:51 +0000 (18:45 +0000)]
* samtools-0.1.10-4 (r840)
* drop the MNP caller. It is slow while does not diliver too much
benefit. Possibly I will work on it in future given more time.
* there is a segfault in pileup
* someone has reported segfault from view/index/sort
Heng Li [Thu, 18 Nov 2010 22:30:11 +0000 (22:30 +0000)]
* samtools-0.1.10-6 (r839)
* call MNPs without realignment because it seems to me that it is not
worthwhile to significantly slow down SNP calling.
* the result looks quite different from the previous version. I have
work to do...
Heng Li [Thu, 18 Nov 2010 16:26:09 +0000 (16:26 +0000)]
Apply a patch by Rob Davis, which improves fault detection.
Heng Li [Thu, 18 Nov 2010 16:09:23 +0000 (16:09 +0000)]
* samtools-r836
* initiate MNP realignment when the MNP has at least 0.2% frequency (otherwise too slow)
Heng Li [Thu, 18 Nov 2010 05:25:13 +0000 (05:25 +0000)]
* modify the filtering rule: also filter SNPs around filtered indels
* added MNP filter
Heng Li [Thu, 18 Nov 2010 04:13:52 +0000 (04:13 +0000)]
* samtools-0.1.10-4 (r834)
* fixed a silly bug in printing MNP
* restrict to at most 1 alternative allele
Heng Li [Thu, 18 Nov 2010 02:58:58 +0000 (02:58 +0000)]
fixed a bug in printing MNPs
Heng Li [Thu, 18 Nov 2010 02:47:20 +0000 (02:47 +0000)]
minor change to how seqQ is applied
Heng Li [Thu, 18 Nov 2010 02:41:12 +0000 (02:41 +0000)]
* samtools-0.1.10 (r831)
* initial MNP caller
Heng Li [Wed, 17 Nov 2010 04:14:15 +0000 (04:14 +0000)]
Release samtools-0.1.10 (r829)
Heng Li [Wed, 17 Nov 2010 01:48:49 +0000 (01:48 +0000)]
update version information: samtools-0.1.9-20 (r828)
Heng Li [Tue, 16 Nov 2010 20:32:50 +0000 (20:32 +0000)]
bcftools: allow to skip indels
Heng Li [Tue, 16 Nov 2010 19:11:58 +0000 (19:11 +0000)]
remove ZQ if both BQ and ZQ are present
Heng Li [Tue, 16 Nov 2010 18:51:33 +0000 (18:51 +0000)]
* samtools-0.1.9-18 (r825)
* change the behaviour of calmd such that by default it does not change the base quality
Heng Li [Tue, 16 Nov 2010 04:31:53 +0000 (04:31 +0000)]
* samtools-0.1.9-17 (r824)
* added command line options to change the default parameters in indel calling
* update the manual
Heng Li [Mon, 15 Nov 2010 17:20:13 +0000 (17:20 +0000)]
* samtools-0.1.9-r823
* the BQ tag is now 64 shifted, not 33 shifted
Heng Li [Mon, 15 Nov 2010 05:30:18 +0000 (05:30 +0000)]
* samtools-0.1.9-16 (r822)
* keep the raw depth because in indel calling, DP4 may be way off the true depth
Heng Li [Sat, 13 Nov 2010 06:18:31 +0000 (06:18 +0000)]
* samtools-0.1.9-15 (r821)
* calmd: write BQ
* skip realignment if BQ is present
Heng Li [Sat, 13 Nov 2010 06:08:26 +0000 (06:08 +0000)]
* samtools-0.1.9-14 (r820)
* penalize reads with excessive differences in indel calling
Heng Li [Sat, 13 Nov 2010 02:36:27 +0000 (02:36 +0000)]
* samtools-0.1.9-13 (r819)
* fixed a bug in pileup given refskip
Heng Li [Fri, 12 Nov 2010 18:04:53 +0000 (18:04 +0000)]
* samtools-r818
* for indel calling, do two rounds of probabilistic realignments
Heng Li [Fri, 12 Nov 2010 01:04:07 +0000 (01:04 +0000)]
* samtools-0.1.19-11 (r817)
* only initiate indel calling when 0.2% of reads contain a gap
Heng Li [Thu, 11 Nov 2010 06:22:59 +0000 (06:22 +0000)]
* samtools-0.1.9-10 (r816)
* I know why the forward method fails. it is because of zero base
qualities. when that is fixed, the forward method seems to give
better results than Viterbi, as it should be. I am tired...
Heng Li [Thu, 11 Nov 2010 05:57:15 +0000 (05:57 +0000)]
effectively revert to the viterbi version. The forward realignment gives too many false positives.
Heng Li [Thu, 11 Nov 2010 05:18:02 +0000 (05:18 +0000)]
* samtools-0.1.9-9 (r810)
* use forward, instead of viterbi, for realignment
* realignment is now quality aware
Heng Li [Thu, 11 Nov 2010 03:45:24 +0000 (03:45 +0000)]
* prepare to replace kaln with kprobaln in realignment
Heng Li [Wed, 10 Nov 2010 22:28:50 +0000 (22:28 +0000)]
fixed a typo
Heng Li [Wed, 10 Nov 2010 21:54:46 +0000 (21:54 +0000)]
use zlib for direct reading when BCF_LITE is in use
Heng Li [Wed, 10 Nov 2010 21:32:13 +0000 (21:32 +0000)]
* do not use reads containing too many mismatches for indel calling
* fixed a trivial bug in case of multi-allelic indels
Heng Li [Wed, 10 Nov 2010 18:23:02 +0000 (18:23 +0000)]
* samtools-0.1.9-8 (r809)
* fixed a bug in the indel caller
Heng Li [Wed, 10 Nov 2010 17:24:10 +0000 (17:24 +0000)]
minor change to makefile
Heng Li [Wed, 10 Nov 2010 17:10:21 +0000 (17:10 +0000)]
* samtools-0.1.9-8 (r807)
* collect indel candidates only from specified platforms (@RG-PL)
* merge varFilter and filter4vcf in vcfutils.pl
Heng Li [Wed, 10 Nov 2010 03:05:46 +0000 (03:05 +0000)]
bcftools: compute equal-tail (Bayesian) credible interval
Heng Li [Tue, 9 Nov 2010 21:28:39 +0000 (21:28 +0000)]
added a double-hit filter to avoid overestimated indel likelihood
Heng Li [Tue, 9 Nov 2010 19:12:06 +0000 (19:12 +0000)]
* samtools-0.1.9-7 (r804)
* fixed a bug in the gap caller
Heng Li [Tue, 9 Nov 2010 15:45:33 +0000 (15:45 +0000)]
* samtools-0.1.9-6 (r803)
* mpileup: apply homopolymer correction when calculating GL, instead of before
* bcftools: apply a different prior to indels
Heng Li [Tue, 9 Nov 2010 04:53:15 +0000 (04:53 +0000)]
* samtools-0.1.9-5 (r802)
* relax tandem penalty. this will be made a command-line option in future.
Heng Li [Tue, 9 Nov 2010 04:35:52 +0000 (04:35 +0000)]
* samtools-0.1.9-4 (r801)
* fixed a minor issue in printing indel VCF
Heng Li [Mon, 8 Nov 2010 20:28:14 +0000 (20:28 +0000)]
fixed another silly bug in mpileup's indel caller
Heng Li [Mon, 8 Nov 2010 19:28:27 +0000 (19:28 +0000)]
fixed a silly bug in the indel caller
Heng Li [Mon, 8 Nov 2010 19:07:33 +0000 (19:07 +0000)]
Incorporate patches by Marcel Martin for read counting.
Heng Li [Mon, 8 Nov 2010 18:39:52 +0000 (18:39 +0000)]
* samtools-0.1.9-2 (r797)
* mpileup: indel calling seems to be working
Heng Li [Mon, 8 Nov 2010 15:54:46 +0000 (15:54 +0000)]
indel calling is apparently working, but more information needs to be collected
Heng Li [Mon, 8 Nov 2010 05:39:18 +0000 (05:39 +0000)]
fixed a few bugs in the indel caller. Probably there are more.
Heng Li [Mon, 8 Nov 2010 03:23:16 +0000 (03:23 +0000)]
prepare for the indel caller. It is not ready yet.
Heng Li [Fri, 5 Nov 2010 15:28:23 +0000 (15:28 +0000)]
Revert to r790. The recent changes are not good...
Heng Li [Fri, 5 Nov 2010 04:19:14 +0000 (04:19 +0000)]
* this revision is UNSTABLE
* indel caller seems working, but it is very insensitive and has
several things I do not quite understand.
Heng Li [Fri, 5 Nov 2010 02:58:43 +0000 (02:58 +0000)]
for backup. no effective changes
Heng Li [Wed, 3 Nov 2010 19:51:24 +0000 (19:51 +0000)]
fixed a minor problem in the example coming with kprobaln.c
Heng Li [Tue, 2 Nov 2010 19:41:27 +0000 (19:41 +0000)]
Separate kaln and kprobaln as I am preparing further changes. At
present, the results should be identical to the previous.
Petr Danecek [Tue, 2 Nov 2010 16:19:04 +0000 (16:19 +0000)]
Added -b option: read file names from a file
Heng Li [Sat, 30 Oct 2010 03:17:22 +0000 (03:17 +0000)]
* samtools-0.1.9-2 (r787)
* Allow to set a maximum per-sample depth to reduce memory. However,
BAQ computation is still applied to every read. The speed is not
improved.
Heng Li [Fri, 29 Oct 2010 16:10:40 +0000 (16:10 +0000)]
* samtools-0.1.9-1 (r786)
* samtools: optionally perform exact test for each sample
Heng Li [Fri, 29 Oct 2010 13:42:25 +0000 (13:42 +0000)]
Optionally output "DP", the individual read depth
Heng Li [Thu, 28 Oct 2010 03:10:27 +0000 (03:10 +0000)]
acknowledge Petr and John who have greatly contributed to the project.
Heng Li [Thu, 28 Oct 2010 02:47:47 +0000 (02:47 +0000)]
Release samtools-0.1.9 (r783)
Heng Li [Wed, 27 Oct 2010 23:58:54 +0000 (23:58 +0000)]
fixed a silly bug in pileup
Heng Li [Wed, 27 Oct 2010 18:39:48 +0000 (18:39 +0000)]
* samtools-0.1.8-22 (r781)
* made BAQ the default behavior of mpileup
* updated manual
* in merge, force to exit given inconsistent header when "-R" is not in use.
Heng Li [Wed, 27 Oct 2010 15:01:11 +0000 (15:01 +0000)]
* samtools-0.1.8-21 (r780)
* minor speedup to pileup
Heng Li [Wed, 27 Oct 2010 13:58:56 +0000 (13:58 +0000)]
improve pileup a little bit
Heng Li [Wed, 27 Oct 2010 04:14:43 +0000 (04:14 +0000)]
* samtools-0.1.8-20 (r778)
* speed up pileup, although I do not know how much is the improvement
Heng Li [Tue, 26 Oct 2010 21:26:04 +0000 (21:26 +0000)]
* samtools-0.1.8-19 (r777)
* integrate mpileup features to pileup: min_baseQ, capQ, prob_realn, paired-only and biased prior
Heng Li [Tue, 26 Oct 2010 19:27:46 +0000 (19:27 +0000)]
remove local realignment (probabilistic realignment is still there)
On behalf of John Marshall [Thu, 21 Oct 2010 10:52:38 +0000 (10:52 +0000)]
Add the relevant filename or region to error messages, and cause a failure
exit status where appropriate. Based on a patch provided by Marcel Martin.
Heng Li [Tue, 19 Oct 2010 23:44:31 +0000 (23:44 +0000)]
* Minor code changes. No real effect.
* change quality to 30 in toy.sam