This is a test version of Biostars. For the public version, visit https://www.biostars.org.
bam-readcount fails to report on all alignments in bam file

Hi All

I’m using bam-readcount on a bam file with ~60000 alignments (stitched reads which span the 374 bp length of the amplicon reference).

bam-readcount   -w -1 -q 0 -b 0 -f my_ref.fa  sort.bam > bam_readcount.q0b0_w-1  2> bam_readcount.error

Here is a “head” of the output:

my_ref  1       G       7999    =:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00  A:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00  C:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00        G:7999:44.00:37.95:44.00:7999:0:0.00:0.00:0.00:0:0.00:0.00:0.00 T:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00  N:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00
my_ref  2       A       8000    =:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00  A:8000:44.00:37.83:44.00:8000:0:0.00:0.00:0.00:0:0.00:0.00:0.00 C:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00        G:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00  T:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00  N:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00
my_ref  3       T       8001    =:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00  A:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00  C:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00        G:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00  T:8001:44.00:37.96:44.00:8001:0:0.00:0.00:0.00:0:0.00:0.00:0.00 N:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00
my_ref  4       G       8002    =:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00  A:1:44.00:11.00:44.00:1:0:0.00:0.00:0.00:0:0.00:0.00:0.00       C:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00        G:8000:44.00:37.94:44.00:8000:0:0.00:0.00:0.00:0:0.00:0.00:0.00 T:1:44.00:38.00:44.00:1:0:0.00:0.00:0.00:0:0.00:0.00:0.00       N:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00

Note that only about 8000 reads are reported for each position and this is true for each position in the reference. Why are not all 60000 reads accounted for? Interestingly I get warnings for 8020 reads (about the number of alignments reported in the output) 

WARNING: In read M00645:234:000000000-AT4W2:1:1106:7492:11191: Couldn't find the generated tag.

This is repeated 374 times for each read with a warning.

Thanks for your help

Mark

bam-readcount error warning alignment bam
• 2,239 views • link
written by mark.rose • mark.rose

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • Varscan fpfilter error
    written by a01039012387 • a01039012387
    Hello. I'm new to Varscan. I want to call the SNPs from the bam file using varscan in the order of somatic - processSomatic-fpfilter but …
  • Bam-readcount output empty from certain regions (RNA seq bams)
    written by nils.rudqvist • nils.rudqvist
    Hi all, I am using bam-readcount to find reads in a RNA seq bam at the location of variants found using RNA seq. However, I …
  • bam-readcount : Reference_base is not the base with highest count
    written by AB AB
    Hi, I recently ran bam-readcount on my bamfile and I noticed something weird about my result chr1 201334382 G 58077 =:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00 A:57989:238.35:35.90:0.01:34667:23322:0.51:0.01:43.91:34670:0.52:122.12:0.49 C:8:255.00:33.00:0.00:6:2:0.37:0.02:56.38:6:0.76:120.62:0.62 G:61:242.61:34.13:0.00:41:20:0.37:0.01:20.82:41:0.77:116.07:0.60 T:17:255.00:35.24:0.00:15:2:0.25:0.02:50.18:15:0.76:117.06:0.69 …
  • Bam-readcounts min-base-quality not giving expected output
    written by sheenams sheenams
    I would like to limit the base-quality-score to phred 13 (similar to default of samtools mpileup (1.3.1)), using bam-readcount (0.8.0). Despite various command modifications, I'm …
  • How to convert bam-readcount output to BED-VAF format
    written by Haitao • Haitao
    Dear All, When I use bam-readcount to count each SNP position's depth for making %VAF, Especially for those pair samples, only happened mutation in one …
  • bam-readcount fails to report all positions
    written by mark.rose • mark.rose
    Hello I'm using bam-readcount (version 0.7.4-unstable-35-0f3b356-dirty). Though it's version name doesn't exactly instill confidence, it is apparently the latest version available on github. Anyway, my …
  • bam- readcount error
    written by mark.rose • mark.rose
    When running bam-readcount I get (minus the 8M previous warnings of the same type): WARNING: In read m151222_205817_42179_c100908792550000001823210904291644_s1_p0/4606/0_7297: Couldn't find the generated tag. WARNING: In …
  • unexpected bowtie2 unpaired alignment behavior
    written by mark.rose • mark.rose
    Hello Seemingly according to the manual, bowtie2 will attempt to align paired-end reads first as pairs but then as unpaired reads if necessary. If this …
  • bam-readcount reports all reference bases as N
    written by Chris F. • Chris F.
    Hi, I'm running bam-readcount (commit `6c3f3ae901`) on a few hundred bam files against a single reference fasta file (designated with `-f`). However, when I look …
  • How To Interpret Output From Bam-Readcount
    written by merbergd • merbergd
    <p>After reading the documentation at <a href='https://github.com/genome/bam-readcount/blob/master/README.textile'>https://github.com/genome/bam-readcount/blob/master/README.textile</a> there is one column in the output from <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a>-readcount that I still don't understand. </p> <p>Columns 1-4 …