Skip to content

GitLab

  • Menu
Projects Groups Snippets
    • Loading...
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
  • Sign in
  • B Bio-DB-HTS
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
    • Locked Files
  • Issues 6
    • Issues 6
    • List
    • Boards
    • Service Desk
    • Milestones
    • Iterations
  • Jira
    • Jira
  • Merge requests 0
    • Merge requests 0
  • Requirements
    • Requirements
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Schedules
    • Test Cases
  • Deployments
    • Deployments
    • Environments
    • Releases
  • Monitor
    • Monitor
    • Incidents
  • Packages & Registries
    • Packages & Registries
    • Package Registry
    • Container Registry
    • Infrastructure Registry
  • Analytics
    • Analytics
    • CI/CD
    • Code review
    • Insights
    • Issue
    • Repository
    • Value stream
  • Wiki
    • Wiki
  • Snippets
    • Snippets
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
Collapse sidebar
  • ensembl-gh-mirror
  • Bio-DB-HTS
  • Issues
  • #54

Closed
Open
Created Jul 11, 2017 by Marek Szuba@mks

'missing SAM header' when opening SAM files

Created by: AlexanderDilthey

Hi,

There seems to be an issue with iterating over alignments in a SAM file. I get the errors

[E::sam_parse1] missing SAM header
[W::sam_read1] parse error at line 4

although there are headers.

Code for reproducing the problem:

use Bio::DB::HTS;
my $sam =  Bio::DB::HTS->new(-bam => 'test.sam');
my $alignment_iterator = $sam->features(-iterator => 1);	
$alignment_iterator->next_seq;

and test.sam:

@HD	VN:1.5	SO:coordinate
@SQ	SN:ref	LN:100
r001	99	ref	1	30	1M3D2M	=	37	3	TAT	*
r002	99	ref	1	30	1M2D2M	=	37	3	TAT	*

samtools can view the SAM without problems.

Assignee
Assign to
Time tracking