samtools
Tools for handling high-throughput sequencing (genomics) data. tldr-pages
On several platforms
Several platforms
Tools for handling high-throughput sequencing (genomics) data. Used for reading/writing/editing/indexing/viewing of data in SAM/BAM/CRAM format.
Convert a SAM input file to BAM stream and save to file
samtools view -S --bam input.sam >output.bamIdentifyTake input from stdin (-) and print the SAM header and any reads overlapping a specific region to stdout
other_command | samtools view --with-header - chromosome:start-endIdentifySort file and save to BAM (the output format is automatically determined from the output file's extension)
samtools sort input --output output.bamIdentifyIndex a sorted BAM file (creates sorted_input.bam.bai)
samtools index sorted_input.bamIdentifyPrint alignment statistics about a file
samtools flagstat sorted_inputIdentifyCount alignments to each index (chromosome/contig)
samtools idxstats sorted_indexed_inputIdentifyMerge multiple files
samtools merge output input1 input2 ...IdentifySplit input file according to read groups
samtools split merged_inputIdentify