This queue is for tickets about the Bio-SamTools CPAN distribution.

Report information
The Basics
Id:
66138
Status:
new
Priority:
Low/Low
Queue:

People
Owner:
Nobody in particular
Requestors:
Yi-Shiou.Chen [...] accelrys.com
Cc:
AdminCc:



Subject: _coverage function should not count reads that skip the input region
Date: Thu, 24 Feb 2011 13:27:47 -0800
To: "bug-Bio-SamTools@rt.cpan.org" <bug-Bio-SamTools@rt.cpan.org>
From: Yi-Shiou Chen <Yi-Shiou.Chen@accelrys.com>

Hi Bio-Samtools Team,

 

I have noticed a problem of getting coverage of RNA-Seq reads. These reads usually have long skips in middle (i.e. long “N” in CIGAR string) to cross introns. However, if I use the _coverage function on an intron, I still get very high coverage because the function counts the reads that cross the introns to map to exons at both ends . The solution, as you already mentioned in the function, is to use pileup instead in the _coverage function. The latest samtools pileup function returns bam_pileup1_t data structure with is_refskip and is_del flags to be used to ignore reads that skip the region.

 

I am using Bio-Samtools 1.27 by Perl 5.8.9 on Linux64.

 

Thanks,

 

Yi-Shiou Chen

Accelrys

 

 

 

 

If you do not wish to register with Bitcard, please report your bug by sending email to: bug-Bio-SamTools [at] rt.cpan.org.

Try to include any information you think might help the developer isolate, reproduce and fix your problem.

Be sure to include at least the following information:

  • Distribution name and version (For example, "DBIx-SearchBuilder-0.46")
  • Perl version (Find this by running the command perl -v)
  • Operating System vendor and version (Find this by running the command uname -a)

Other things that are often helpful:

  • Details about your operating environment that might be related to the issue being described
  • Exact cut and pasted error or warning messages
  • The shortest, clearest code you can manage to write which reproduces the bug described.
  • A patch against the latest released version of this distribution which fixes this bug.

 

 

 

sig-accl-logo

sig-accl-end

 

Yi-Shiou Chen
Principal Engineer

Accelrys, Inc.
10188 Telesis Court, Suite 100
San Diego, CA 92121, USA
+858 799 5610

twitter   Linked in   Blog

Stay Connected with Accelrys and the
Scientific & Engineering Community

Blog   |   Forum   |   Website   |   Support

 

 

Image displayed inline above

Image displayed inline above

Image displayed inline above

Image displayed inline above

Image displayed inline above



This service runs on Request Tracker, is sponsored by The Perl Foundation, and maintained by Best Practical Solutions.

Please report any issues with rt.cpan.org to rt-cpan-admin@bestpractical.com.