Skip to main content
Search type Search syntax
Tags [tag]
Exact "words here"
Author user:1234
user:me (yours)
Score score:3 (3+)
score:0 (none)
Answers answers:3 (3+)
answers:0 (none)
isaccepted:yes
hasaccepted:no
inquestion:1234
Views views:250
Code code:"if (foo != bar)"
Sections title:apples
body:"apples oranges"
URL url:"*.example.com"
Saves in:saves
Status closed:yes
duplicate:no
migrated:no
wiki:no
Types is:question
is:answer
Exclude -[tag]
-apples
For more details on advanced search visit our help page
Results tagged with
Search options questions only not deleted user 272

Questions specific to interacting with and post-processing sequence alignments using the SAMtools package. For questions specifically about formats SAM, BAM or CRAM use tags sam or bam,

4 votes
2 answers
655 views

What is a samtools mpileup reference skip?

The samtools documentation for mpileup states: At this column, a dot stands for a match to the reference base on the forward strand, a comma for a match on the reverse strand, a '>' or '<' for a reference …
mattm's user avatar
  • 784
15 votes
1 answer
2k views

Are soft-clipped bases used for variant calling in samtools + bcftools?

If there are soft clipped base pairs specified in the CIGAR string for a read in a SAM/BAM file, will these be used for variant calling in a samtools + bcftools workflow? … The samtools documentation does not mention clipped bases. …
mattm's user avatar
  • 784
4 votes
2 answers
342 views

Filter bam using SNP list in bed format with minimum mapping quality and base quality

Samtools seems to almost solve this problem, but not quite: samtools depth -b snp_list -q 20 -Q 30 example.bam allows me to count the number of alignments meeting my desired criteria, but not to produce … There is no -Q option for samtools view, and the -q option is for mapping quality Can this be achieved in samtools? Or is there another preexisting tool that can? …
mattm's user avatar
  • 784
10 votes
1 answer
4k views

What are all the reference files produced by bwa index, and are these dependent upon whether...

I have discovered that samtools does not take a gzipped reference, so I am planning to use an unzipped version of the reference for my workflow instead of dealing with two separate representations of the …
mattm's user avatar
  • 784