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