Posts by aiteteji

... Thanks for posting the link. I read it, but I didn't get how to remove the "@" in .bam and "##" in vcf file. I was hoping that someone who knows better can help. Here is my approach:
import pysam
pysamlist = []
with pysam.AlignmentFile("ex1.bam", "rb") as f:
tags = pysam. ...

... Hi Biostars community,
Is there an analogous way in Pysam to the code below to extract the header information in .bam and .vcf (or .vcf.gz) files? Help with a code snippet will great.
The code snippet below uses exifread library to get the meta-info (in key-value pair) in jpeg and tiff files.
imp ...

... Hello Biostars,
I have a list of chromosome and coordinates in my vcf file, see examples:
1:10327
1:10434
1:10440
1:10469
...
Is there anyway to get the rs number for each coordinate using plink or bcftools? I found some old answers here, but I'll like to know how this is curr ...

... I want to archive some .ped and associated .map files. But do not know if any, how to extract the embedded metadata in them. These metadata are expected to aid data discovery that helps in the querying of the database to find matching objects.
Anyone has idea if there are standards to use to curate ...