FEATURE REQUEST: Improve seqinfo(<BamFile>) to return the genome assembly associated with each sequence

Bioconductor / Rsamtools

Binary alignment (BAM), FASTA, variant call (BCF), and tabix file import

https://bioconductor.org/packages/Rsamtools

Other

27 stars 27 forks source link

FEATURE REQUEST: Improve seqinfo(<BamFile>) to return the genome assembly associated with each sequence #61

Open hpages opened 8 months ago

hpages commented 8 months ago

Right now seqinfo(<BamFile>) returns a Seqinfo object with no genome information, even though that information is available in the file (the genome assembly is specified via the AS tag found in the @SQ header lines).

See https://github.com/Bioconductor/GenomeInfoDb/issues/105#issuecomment-1942470165 for some context.

I'll try to come up with a PR that fixes that.