[Bioperl-l] adding tags

[Jonas Zierer]

Hi everybody,

I have a bam file of aligned transcriptome reads and use bioperl
(Bio::DB::Sam) to work with it.
Now i wanted to ask if there is any possibility to add tags to
alignments (=reads) in this File?


I use 
	my $iterator     = $bam->features(-iterator => 1, -flags    => {M_UNMAPPED=>0});
	while (my $align = $iterator->next_seq) { 
        	...
	}
to loop through the aligned reads and I am searching for anthing like 
	$align->addTag(key=>value)

and i have one more question: is it possible to get the coverage of a
segment using only the read start positions? (that means that each read
is treated as it would have Length 1)

Thx
bye