while (my $align = $bam->read1) { my $seqid = $align->seq_id; my $start = $align->start; my $end = $align->end; my $strand = $align->strand; my $ref_dna= $align->dna; }