SAM/BAM Convert Sequence
From Genome Analysis Wiki
Jump to navigationJump to search
SAM/BAM support conversion between '=' and the base in a sequence
Request: To be able to call SamRecord::getSequence and have it return the sequence with '=' or with the correct bases.
Requester: Goncalo Abecasis & Hyun Min Kang
Date Requested: Nov 23, 2010?
Notes
It was requested to have this handled within the getSequence call so the user doesn't have to convert the sequence after calling getSequence.
Proposed Solution
....More details coming soon!
Added a class, SamQuerySeqWithRef
, with methods seqWithEquals
and seqWithoutEquals