[Bioc-sig-seq] Getting optional fields from BAM file using Rsamtools

Steve Lianoglou mailinglist.honeypot at gmail.com
Mon Mar 8 17:04:33 CET 2010


Hi Martin,

On Mon, Mar 8, 2010 at 4:12 AM, Martin Morgan <mtmorgan at fhcrc.org> wrote:
> In v. 0.1.40, ScanBamParam(tag=c("XX", "YY", "ZZ")) tells scanBam to
> parse those fields, if they are encountered. Additional small test cases
> welcome -- I do not have access to a very rich set of real tags.

Just tested it out, and it seems like it will work nicely. I also have
a small set of flags I need, so have limited testing scenarios, but
everything seems to be fine.

If I run into any problems down the road, I'll provide specific examples.

Thanks,
-steve

-- 
Steve Lianoglou
Graduate Student: Computational Systems Biology
 | Memorial Sloan-Kettering Cancer Center
 | Weill Medical College of Cornell University
Contact Info: http://cbio.mskcc.org/~lianos/contact



More information about the Bioc-sig-sequencing mailing list