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)

Merge remote branch 'main/master'

Conflicts:
misc/Makefile
misc/bamcheck.c

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

compiling errors on some machines

Temporary fix to prevent segfault in kpa_glocal

shrink bam1_t::m_data when sorting

Otherwise most of allocated memory may not be used

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

Conflicts:
misc/Makefile

Also, killed warnings on Mac

Merge pull request #2 from peterjc/depad

Improvments to 'samtools depad' command

Compile cleanly from top-level Makefile. Check for out of bounds tids.

Compile cleanly from top-level Makefile. Check for out of bounds tids.

Added rule for bamcheck creation

Added rule for bamcheck creation

Initial release of bamcheck and plot-bamcheck

Initial release of bamcheck and plot-bamcheck

Check reference length in BAM header against embedded reference in 'samtools depad'

Cope with trailing CIGAR H (or S and H) in 'samtools depad'

Cope with leading CIGAR H (or H and S) in 'samtools depad'

Tolerate CIGAR H operators in 'samtools depad'

Include leading CIGAR P operators in 'samtools depad'

Record actual sequence in unpad_seq (fixes apparently harmless bug)

filter on query sequence length (not aln len)

move seqtk to a separate project

Correct my error message; assed some comments

Add TODO note about reference lengths

Cope with CIGAR D or P to mark gaps in an embedded padded reference

Check current reference sequence matches RNAME in samtools depad

Support CIGAR X/= in depad

Helpful error message on bad CIGAR in unpad_seq

change SO

a little more accurate mem estimate

control compression level at command line

optionally surpress /[12] in bam2fq

Merge branch 'master' into mt

bugfix: downsampling not always working

name sorting checkes READ1 and READ2

sorting also considers strand

Previously, merging considers strand but sorting not. This is a bug and breaks
stability. Now it seems that sorting is stable and it gives the same output no
matter what parameters in use.

change the default block size in sorting