Skip Menu |

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

Report information
The Basics
Id: 66138
Status: new
Priority: 0/
Queue: Bio-SamTools

Owner: Nobody in particular
Requestors: Yi-Shiou.Chen [...]

Bug Information
Severity: (no value)
Broken in: (no value)
Fixed in: (no value)

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 [...]" <bug-Bio-SamTools [...]>
From: Yi-Shiou Chen <Yi-Shiou.Chen [...]>
Download (untitled) / with headers
text/plain 2.2k
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] 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. [cid:image001.gif@01CBD422.D952A600] [cid:image002.gif@01CBD422.D952A600] Yi-Shiou Chen Principal Engineer Accelrys, Inc. 10188 Telesis Court, Suite 100 San Diego, CA 92121, USA +858 799 5610 [cid:image003.gif@01CBD422.D952A600]<> [cid:image004.gif@01CBD422.D952A600] <> [cid:image005.gif@01CBD422.D952A600] <> Stay Connected with Accelrys and the Scientific & Engineering Community Blog<> | Forum<> | Website<> | Support<>
Download (untitled) / with headers
text/html 14.3k

Message body is not shown because it is too large.

Download image001.gif
image/gif 1.1k
Download image002.gif
image/gif 934b
Download image003.gif
image/gif 454b
Download image004.gif
image/gif 493b
Download image005.gif
image/gif 547b

This service is sponsored and maintained by Best Practical Solutions and runs on infrastructure.

Please report any issues with to