Questions tagged [mpileup]

The tag has no usage guidance.

Filter by
Sorted by
Tagged with
0
votes
0answers
32 views

bcftools mpileup double counting for overlapping reads

I have a dataset of paired-end reads, which can overlap. So when i run bcftools mpileup this overlap gives a double count in the depths so we see DP = 2 even though it should be DP = 1. The question ...
2
votes
1answer
67 views

How to generate a consensus sequence from a multi-reference BAM file?

I am trying to generate a consensus sequence from a BAM file that was generated by mapping reads to a reference FASTA containing multiple sequences. Usually, I generate consensus sequences from BAM ...
0
votes
0answers
9 views

How do I tell IsoMut2py what to input into mpileup?

Just a warning, this question is extremely specific to the indel-detection software IsoMut2py. The paper for isomut can be found here: https://bmcbioinformatics.biomedcentral.com/articles/10.1186/...
3
votes
1answer
120 views

Is there a tool that can perform a read-group-aware mpileup from a single file?

I would like to perform a samtools mpileup from a single file that contains thousands of read groups with different SM tags. I could split the bam by read group ...
0
votes
2answers
167 views

Generate VCF from different .bam files with different chromosome names

I have two resources of .bam files. One is generated by our lab (1 sample = 1 bam). One is downloaded online (again 1 sample = 1 bam). For the downloaded samples the chromosomes are labelled: chr1, ...
1
vote
0answers
19 views

Indel quality from mpileup

Could someone please confirm that I understand this correctly? I have the res and qual columns from mpileup and I would like to match them to get the qual per base. It seems that the indel initiation (...
2
votes
1answer
69 views

Why does samtools mpileup sometimes include ref bases (other than ',' or '.')?

This is my first post here. I can think of no way of giving an easily reproducible example, as per the stack ethos. SO apologies in advance and any feedback on question format appreciated. I have ...
3
votes
2answers
233 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 ...
5
votes
2answers
226 views

samtools mpileup skipping read

I run the command: samtools mpileup -O -s -q20-B -Q20 -f hg19.fa -r chr1:569929-569931 myFile.bam and get: ...
0
votes
1answer
164 views
7
votes
2answers
262 views

variant calling on ChIP-seq style data: samtools mpileup with minimal filters

I am running samtools mpileup (v1.4) on a bam file with very choppy coverage (ChIP-seq style data). I want to get a first-pass list of positions with SNVs and their frequency as reported by the read ...
7
votes
1answer
321 views

samtools mpileup empty when filtering out flags

I produced a bam file by aligning reads to a small set of synthetic sequences using bwa-mem. I am heavily filtering reads that are not paired and of a certain orientation. Applying the filtering, I ...