Pysam is a Python package for reading, manipulating, and writing genomics data such as SAM/BAM/CRAM and VCF/BCF files. It's a lightweight wrapper of the HTSlib API, the same one that powers samtools, bcftools, and tabix.
I was looking at using pysam for handling some FASTA file data for the Tree of Life project, and creating an output assembly. Is there a reverse-complement function available via pysam?
I thought pysam would be a nice thing to use since it would also handle bgzipped FASTA.
Also wondering about very long sequences. We now have some chromosomes > 2^31 bp long.
I was looking at using pysam for handling some FASTA file data for the Tree of Life project, and creating an output assembly. Is there a reverse-complement function available via pysam? I thought pysam would be a nice thing to use since it would also handle bgzipped FASTA. Also wondering about very long sequences. We now have some chromosomes > 2^31 bp long.