def _read(self, arr, start, end, strand): # get bitarray b = BitArrayTrack._read(self, arr, start, end) # decode to sequence codec_dict = _decoders[self.bpb] seq = ''.join(b.decode(codec_dict)) if strand == NEG_STRAND: seq = dna_reverse_complement(seq) return seq