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

Bioinformatics Asked on September 5, 2021

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 called samtools mpileup (Version: 1.9 (using htslib 1.9)) on some oxford nanopore data (mapped to GRCh38).

I am seeing some instances like below (grouped with pythons Counter), here reference is T at this position:

 Counter({',': 34, '.': 33, 'c': 9, 'g': 6, 'a': 4, 't': 2, '-:ACTC': 1, '-:ACTCGG': 1, '-:5': 1, '$': 1, '*': 1})

Why is ‘t’: 2′ in there?

flags used: -A -B -Q 1 -R -x

One Answer

Why is 't': 2' in there?

How else do you think the pileup should show you that two reads (in the reverse direction) have a T in that position?

Answered by swbarnes2 on September 5, 2021

Add your own answers!

Related Questions

Downloading SRA Files from AWS

1  Asked on December 6, 2020


NCT search in PubMed via Entrez (python)

0  Asked on December 5, 2020 by nutarelli-federico


using snakemake shadow rules to store temp files on local nodes

1  Asked on December 2, 2020 by kamil-s-jaron


Software for taxonomic assignment?

2  Asked on December 1, 2020


Detecting Allelic Imbalance

0  Asked on November 26, 2020 by krizbomb


Mapping statistics from bam file using bbtools and sambamba

2  Asked on November 25, 2020 by bioinfonext


Preprocessing Affymetrix SNP Array Matrix

0  Asked on November 17, 2020 by thanh-nguyen


Ask a Question

Get help from others!

© 2023 All rights reserved. Sites we Love: PCI Database, MenuIva, UKBizDB, Menu Kuliner, Sharing RPP, SolveDir