SAMtools-Teaching: Difference between revisions

From Research Computing Center Wiki
Jump to navigation Jump to search
No edit summary
No edit summary
Line 9: Line 9:


=== Version ===
=== Version ===
0.1.20, 1.9, 1.10
0.1.20, 1.6, 1.14, 1.16.1
   
   
=== Author / Distributor ===
=== Author / Distributor ===
Line 23: Line 23:
'''Version 0.1.20'''
'''Version 0.1.20'''


Installed in /apps/eb/SAMtools/0.1.20-GCC-8.3.0
Installed in /apps/eb/SAMtools/0.1.20-GCC-11.2.0


To use this version, please load the module with
To use this version, please load the module with
<pre class="gscript">
<pre class="gscript">
ml SAMtools/0.1.20-GCC-8.3.0
ml SAMtools/0.1.20-GCC-11.2.0
</pre>  
</pre>  


'''Version 1.9'''
'''Version 1.6'''


Installed in /apps/eb/SAMtools/1.9-GCC-8.3.0
Installed in /apps/eb/SAMtools/1.6-GCC-11.2.0


To use this version, please load the module with
To use this version, please load the module with
<pre class="gscript">
<pre class="gscript">
ml SAMtools/1.9-GCC-8.3.0
ml SAMtools/1.6-GCC-11.2.0
</pre>  
</pre>  


'''Version 1.10'''
'''Version 1.14'''


Installed in /apps/eb/SAMtools/1.10-GCC-8.3.0
Installed in /apps/eb/SAMtools/1.14-GCC-11.2.0


To use this version, please load the module with
To use this version, please load the module with
<pre class="gscript">
<pre class="gscript">
ml SAMtools/1.10-GCC-8.3.0
ml SAMtools/1.14-GCC-11.2.0
</pre>
 
'''Version 1.16.1'''
 
Installed in /apps/eb/SAMtools/1.16.1-GCC-11.3.0
 
To use this version, please load the module with
<pre class="gscript">
ml SAMtools/1.16.1-GCC-11.3.0
</pre>  
</pre>  


Line 63: Line 72:
   
   
cd $SLURM_SUBMIT_DIR<br>
cd $SLURM_SUBMIT_DIR<br>
ml SAMtools/1.10-GCC-8.3.0<br>     
ml SAMtools/1.16.1-GCC-11.3.0<br>     
samtools <u>[options]</u><br>   
samtools <u>[options]</u><br>   
</div>
</div>
Line 79: Line 88:
   
   
<pre  class="gcommand">
<pre  class="gcommand">
ml SAMtools/1.10-GCC-8.3.0
ml SAMtools/1.16.1-GCC-1.3.0


samtools --help
samtools --help


Program: samtools (Tools for alignments in the SAM format)
Program: samtools (Tools for alignments in the SAM format)
Version: 1.10 (using htslib 1.10)
Version: 1.16.1 (using htslib 1.16)


Usage:  samtools <command> [options]
Usage:  samtools <command> [options]
Line 102: Line 111:
     addreplacerg  adds or replaces RG tags
     addreplacerg  adds or replaces RG tags
     markdup        mark duplicates
     markdup        mark duplicates
    ampliconclip  clip oligos from the end of reads


   -- File operations
   -- File operations
     collate        shuffle and group alignments by name
     collate        shuffle and group alignments by name
     cat            concatenate BAMs
     cat            concatenate BAMs
    consensus      produce a consensus Pileup/FASTA/FASTQ
     merge          merge sorted alignments
     merge          merge sorted alignments
     mpileup        multi-way pileup
     mpileup        multi-way pileup
Line 113: Line 124:
     fastq          converts a BAM to a FASTQ
     fastq          converts a BAM to a FASTQ
     fasta          converts a BAM to a FASTA
     fasta          converts a BAM to a FASTA
    import        Converts FASTA or FASTQ files to SAM/BAM/CRAM
    reference      Generates a reference from aligned data


   -- Statistics
   -- Statistics
Line 122: Line 135:
     phase          phase heterozygotes
     phase          phase heterozygotes
     stats          generate stats (former bamcheck)
     stats          generate stats (former bamcheck)
    ampliconstats  generate amplicon specific stats


   -- Viewing
   -- Viewing
     flags          explain BAM flags
     flags          explain BAM flags
    head          header viewer
     tview          text alignment viewer
     tview          text alignment viewer
     view          SAM<->BAM<->CRAM conversion
     view          SAM<->BAM<->CRAM conversion
     depad          convert padded BAM to unpadded BAM
     depad          convert padded BAM to unpadded BAM
    samples        list the samples in a set of SAM/BAM/CRAM files


  -- Misc
    help [cmd]    display this help message or help for [cmd]
    version        detailed version information


</pre>
</pre>

Revision as of 14:39, 5 September 2023

Category

Bioinformatics

Program On

Teaching

Version

0.1.20, 1.6, 1.14, 1.16.1

Author / Distributor

SAMtools

Description

"SAM Tools provide various utilities for manipulating alignments in the SAM format, including sorting, merging, indexing and generating alignments in a per-position format." More details are at SAMtools

Running Program

Version 0.1.20

Installed in /apps/eb/SAMtools/0.1.20-GCC-11.2.0

To use this version, please load the module with

ml SAMtools/0.1.20-GCC-11.2.0

Version 1.6

Installed in /apps/eb/SAMtools/1.6-GCC-11.2.0

To use this version, please load the module with

ml SAMtools/1.6-GCC-11.2.0

Version 1.14

Installed in /apps/eb/SAMtools/1.14-GCC-11.2.0

To use this version, please load the module with

ml SAMtools/1.14-GCC-11.2.0

Version 1.16.1

Installed in /apps/eb/SAMtools/1.16.1-GCC-11.3.0

To use this version, please load the module with

ml SAMtools/1.16.1-GCC-11.3.0

Here is an example of a shell script, sub.sh, to run on the batch queue:

#!/bin/bash
#SBATCH --job-name=j_SAMtools
#SBATCH --partition=batch
#SBATCH --mail-type=ALL
#SBATCH --mail-user=username@uga.edu
#SBATCH --ntasks=1
#SBATCH --mem=10gb
#SBATCH --time=08:00:00
#SBATCH --output=SAMtools.%j.out
#SBATCH --error=SAMtools.%j.err

cd $SLURM_SUBMIT_DIR
ml SAMtools/1.16.1-GCC-11.3.0
samtools [options]

In the real submission script, at least all the above underlined values need to be reviewed or to be replaced by the proper values.

Please refer to Running_Jobs_on_the_teaching_cluster, Run X window Jobs and Run interactive Jobs for more details of running jobs at Teaching cluster.


Here is an example of job submission command:

sbatch ./sub.sh 

Documentation

ml SAMtools/1.16.1-GCC-1.3.0

samtools --help

Program: samtools (Tools for alignments in the SAM format)
Version: 1.16.1 (using htslib 1.16)

Usage:   samtools <command> [options]

Commands:
  -- Indexing
     dict           create a sequence dictionary file
     faidx          index/extract FASTA
     fqidx          index/extract FASTQ
     index          index alignment

  -- Editing
     calmd          recalculate MD/NM tags and '=' bases
     fixmate        fix mate information
     reheader       replace BAM header
     targetcut      cut fosmid regions (for fosmid pool only)
     addreplacerg   adds or replaces RG tags
     markdup        mark duplicates
     ampliconclip   clip oligos from the end of reads

  -- File operations
     collate        shuffle and group alignments by name
     cat            concatenate BAMs
     consensus      produce a consensus Pileup/FASTA/FASTQ
     merge          merge sorted alignments
     mpileup        multi-way pileup
     sort           sort alignment file
     split          splits a file by read group
     quickcheck     quickly check if SAM/BAM/CRAM file appears intact
     fastq          converts a BAM to a FASTQ
     fasta          converts a BAM to a FASTA
     import         Converts FASTA or FASTQ files to SAM/BAM/CRAM
     reference      Generates a reference from aligned data

  -- Statistics
     bedcov         read depth per BED region
     coverage       alignment depth and percent coverage
     depth          compute the depth
     flagstat       simple stats
     idxstats       BAM index stats
     phase          phase heterozygotes
     stats          generate stats (former bamcheck)
     ampliconstats  generate amplicon specific stats

  -- Viewing
     flags          explain BAM flags
     head           header viewer
     tview          text alignment viewer
     view           SAM<->BAM<->CRAM conversion
     depad          convert padded BAM to unpadded BAM
     samples        list the samples in a set of SAM/BAM/CRAM files

  -- Misc
     help [cmd]     display this help message or help for [cmd]
     version        detailed version information

Back to Top

Installation

Source code is obtained from SAMtools

System

64-bit Linux