Hi, I am using the Fragment store for the read mapping problem. I append the matches using appendAlignedRead() and then
perform global alignment using convertMatchesToGlobalAlignment(),
but when I print the scores using fragStore.alignQualityStore[i].score
it shows 0 ("fragStore" is the FragmentStore object and "i" is
the alignedRead ID). I have generated the SAM file CIGAR
string and the edit distance (NM:i:) are perfectly OK.
-- Nauman Ahmed PhD candidate, Deparment of Software and Computer Technology Computer Engineering Laboratory, EWI, TU Delft |