Matt Shirley

Another example that replicates my issue: ``` curl -s ftp://platgene:[email protected]/IlluminaPlatinumGenomes_v7.0/merged_platinum/NA12878.vcf.gz | gzip -cdfq | grep -v "^#" | vawk '{ if (S$NA12878$GT==S$NA12879$GT) print }' | head -n2 chr1 116549 ....

And FYI, example1 works under linux, fails on my Mac. Example2 works on both linux and Mac.

Doesn't just plain "print" do that?

Good point. Matt Shirley http://mattshirley.com/ > On Oct 20, 2014, at 2:27 PM, Colby Chiang [email protected] wrote: > > at the moment, " { print } | cut -f 10-...

@naoyak We're planning a deployment in our data science group, and the author and date templating are features we discussed as "nice to have", so I'm surprised to find that...

It looks like they're starting work on scRNA support: https://github.com/COMBINE-lab/salmon/blob/a41c6b4e38fb23e51b59dc4a0a450071dc92c180/src/CollapsedCellOptimizer.cpp Seems like @k3yavi is doing most of the implementation so far. Our group would definitely be interested in this functionality,...

Can you confirm which pyfaidx version you're using?

@yang-dongxu could you provide some example code that I could use to reproduce the issue? If so I would be glad to take a look.

Thank you both for providing more details! I'll reproduce this and then figure out a solution.

I've just run your code and can confirm I can reproduce the behavior you're describing. I haven't had a chance to look further yet, but did reproduce using: Python 3.11.5...