UPGRADE NOTICE:
The University Wiki Service was upgraded to Confluence 7.13.7. Please refer to the University Wiki Service Help Pages for a list of changes.
If you encounter any issues with the new version, submit a ticket by emailing us at help@wikis.utexas.edu.
Page tree
Skip to end of metadata
Go to start of metadata

You are viewing an old version of this page. View the current version.

Compare with Current View Page History

« Previous Version 3 Next »

Samtools

Setup output directory.

mkdir -p 03_Output/variant_calling/samtools
 

If you do not have a alignment file in the SAM format you may want to start with Introduction to mapping.

cp 03_Output/mapping/bowtie/REL606.5.sam 03_Output/variant_calling/samtools/
 
cp 03_Output/mapping/bowtie/REL606.5.fasta 03_Output/variant_calling/samtools/
 

Prepare reference file.

samtools faidx 03_Output/variant_calling/samtools/REL606.5.fasta
 

Prepare alignment file.

Convert SAM to BAM format.

samtools view -bS -o 03_Output/variant_calling/samtools/REL606.5.bam 03_Output/mapping/bowtie/REL606.5.sam
[samopen] SAM header is present: 1 sequences.

Variant call output.

samtools mpileup -uf 03_Output/variant_calling/samtools/REL606.5.fasta 03_Output/variant_calling/samtools/REL606.5.bam
[mpileup] 1 samples in 1 input files
<mpileup> Set max per-file depth to 8000
[bam_pileup_core] the input is not sorted (reads out of order)
[afs] 0:0.000 1:0.000 2:0.000
[bam_plp_destroy] memory leak: 1. Continue anyway.
  • No labels