Merge remote branch 'remotes/master/master' into fb-annots

backup commit

Merge branch 'master' of github.com:samtools/samtools

bcftools view: Exit with an error if the -l file cannot be read

Merge pull request #30 from leecbaker/master

Fixed several small memory leaks.

Fix memory leaks:
* When downloading a remote FASTA index
* When indexing a BAM, we can't open the FAI file
* When indexing a BCF, we can't open the index file
* Parsing a VCF header with no sample line
* Encountering multiple adjacent tabs in a SAM header

backup commit

bcftools view: always output the header, even when no other arguments are given

Adding -f flag to samtools sort to give output name directly.

Merge remote branch 'remotes/pierre/master' into develop

fix html

cont

works

cont

works with curses

cont

...

dump log-likelihood

bamcheck: new stats, number of filtered vs raw sequences

Merge branch 'vsbuffalo-master' into develop

fixed missing braces

big fix: using samtools view file.bam gibberish returns alignments; iter not checked against kh_end

mpileup: New --rf and --ff options for filtering reads

bamcheck: New stats - the number of QC failed reads and non-primary alignments

bam2depth: Enabled reading the list of BAMs from a file list

ksort.h: declared __ks_insertsort_##name as static to compile with C99. (On behalf of Hanspeter Niederstrasser.)

bam_plcmd.c: more descriptive -b usage help message

bam_plcmd.c: added sanity check in read_file_list, complain if file list is not a file list

bamcheck: New stats, number of pairs mapped to different chromosomes

misc/Makefile, on behalf of ghuls: changed the order of bamcheck rules with -lz last

bamcheck: max quality value increased to 255 to handle missing qual field; pairs with IS=0 can still count as paired; output IS=0 stats as well

bamcheck: Fixed a bug introduced by previous change igcd to uint: not a good idea to init with -1.

Merge remote branch 'remotes/pd3master/master' into develop

Drop empty bins from coverage distribution

Changed default reflen; more informative error message when GCD bin size too small or reflen too big

More user-friendly error output

Complain when BAM cannot be open. Severe bug fixed in -m haploid calling.

Merge remote branch 'remotes/pd3/master'

Conflicts:
misc/bamcheck.c

Fix in insertion cycle calculation of rev reads

GC-depth calculation now more robust to small GCD bin sizes

New --GC-depth option to bamcheck. Take into account not so frequent CIGAR codes. Fix in realloc_buffers. Prevent segfault with long command line options; Set xrange correctly in plot-bamcheck coverage.gp

Merge pull request #14 from peterjc/depad3

Update MPOS/PNEXT during 'samtools depad'

Fix 'samtools depad' seg-fault when MPOS/PNEXT given but MRNM/RNEXT is not

Tweak order free things (avoids an 'samtools depad' seg fault on error)

First attempt at updating MPOS in 'samtools depad'

r580: bugfix - bam_aux_get() error given 'f' type

Use int64_t to prevent count overflow in view -c; on behalf of Joshua Randall: Set y2range in quals2.png plot

Fix in output of missing haploid genotypes

Allow filtering by read group or sample name in bamcheck. Bug fix in rm_info.

VCF annotations with and without -m and more robust rm_info

Bug fixes

Alternative model for multiallelic and rare-variant calling. Proof-of-principle implementation.

New -p switch for increased sensititivy of indel-calling. Here -m and -F are used per-sample which enables to detect rare indels in low-coverage bams. Experimental.

Merge branch 'master' of github.com:samtools/samtools

New -p option to start tview at desired position

r579: after merging Peter's depad changes

fixed a couple of minor compiling warnings

Merge pull request #11 from peterjc/depad2

Depad2

Update @SQ lines during 'samtools depad'

Note does not (yet) preserve additional @SQ line tags (e.g. AS, SP, UR).

Calculate unpadded reference length, update BAM header (SAM header pending)

SAM/BAM input to 'samtools depad' not linked to FASTA reference

Wording of error messages

Catch missing reference error in faidx (avoids floating point exception)

Check embedded reference against FASTA if given in 'samtools depad'

Use FASTA file when no embedded reference in 'samtools depad'

Load FAI in 'samtools depad' (not used yet)

Add reference FASTA file to 'samtools depad' command line (not used yet)

Support SAM or BAM input/output in 'samtools depad' (default still BAM -> BAM)

Minimal 'samtools depad' command line parsing (-? only so far)

Adding a usage function to 'samtools depad' following 'samtools view' style

Add a tiny bit of help text to 'samtools depad' command.

Revert one of my earlier changes - Heng was right, CIGAR P not sensible in a padded reference

Fix buffer overflow in 'samtools depad' with CIGAR D operators

Fixed the next silly error on that line.

Fix silly typo

More elegant solution to avoid superfluous leading P ops

Refine logic for removing superfluous P ops (still testing ...)

Merge remote branch 'rq/getopt_mac'

Fix missing declaration for getopt on Mac OS X

Update bam1_core_t and bam1_t documentation

* Query strand is no longer represented by a field of its own.
* The order of segments in data appears to have changed.
* Added notes about encoding of cigar and seq data in the bam1_t.data byte array.

New -s option for multisample BAMs: include only reads from a given sample in the display

Fixed errors introduced by the 8c15f916dabce475febdf508a9cc0c708c5a7747 commit.

The -E switch to recalculate extended BAQ on the fly, ignoring existing BQs

Fix to prevent long stretches of Ns to be mistaken for indels; Always output the most likely ALT allele

Merge branch 'master' of github.com:samtools/samtools

Added -f, -F options. Count and plot indels per cycle for each mate pair separately.

change the default clevel to 1 for speed

added bamshuf from htslib

bugfix: bgzf-mt may write prematured output

added length threshold

Reverted to the previous behaviour, the unmapped reads are not removed. Added -r switch for backward compatibility.

init_regions: check for non-matching sequence names

Resolved merge conflict

Indexed access can be now combined with target regions. The number of mapped bases more accurate now with -t.

r572: added "bedcov" command

change getline() to mygetline()

updated to the latest kseq

Added check for out-of-bounds indels at the end of reads

On behalf of Bradford Powell: fixed typo in VCF header, output DV not DP

Removed glibc dependency (getline)