Bio::Tradis::Parser::Bam(3pm)

Bio::Tradis::Parser::Bam3pm

Very basic bam parser. limited functionality.

Usage

Parses BAM files and gives access to basic info in them.

   use Bio::Tradis::Parser::Bam;
   
   my $pipeline = Bio::Tradis::Parser::Bam->new(file => 'abc');
   $pipeline->read_info;
   $pipeline->next_read;
   $pipeline->seq_info;
   $pipeline->is_mapped;
   $pipeline->is_reverse;

=seq_info Reads BAM header and returns a hash (keys are sequence ids, values are hash refs with keys as tags (like LN and M5))

=next_read Moves _currentread to the next entry in the BAM. Returns 0 if EOF.

=read_info Returns info from _currentread = hash reference with field name as key. Standard fields are named as per the SAM format specification: 1 : QNAME 2 : FLAG 3 : RNAME 4 : POS 5 : MAPQ 6 : CIGAR 7 : RNEXT 8 : PNEXT 9 : TLEN 10 : SEQ 11 : QUAL Additional fields will use their tag names. Complete line is returned with key READ

=is_mapped Parses the flag for the current read and determines if mapped. Returns 0 or 1.

=is_reverse Parses the flag for the current read and determines if reverse complemented. Returns 0 or 1.

Version

version 1.4.5

Releases

3 releases since Debian 11 (Bullseye); the text changed 1 time.

  1. 11–12
  2. 13

Author

Carla Cummins <[email protected]>

From bio-tradis 1.4.5+dfsg2-2 in Debian 13 (Trixie). Parsed by carta from the page's man source.