my @all_reads = $segment->features( -type => 'match', -seq_id => $seq_ids[0], -start => $window_start, -end => $window_end, -filter => sub { shift->strand == -1 }, ); #### for my $read (@all_reads) { print $read->qname,',',$read->aux_get("MD"), ',',$read->cigar_array,',', $read->strand, "\n"; };