good first issuegood next issueplugin_todo
倉庫指標
- 星標
- (555 顆星)
- PR 合併指標
- (平均合併 4小時 7分鐘) (30 天內合併 8 個 PR)
描述
yesterday, I spent some time digging into a sourmash use case with @shannonekj, and a few different reasonably generic utility script needs emerged.
the code for this is in a private repository so I'll try to describe things here - we were looking for differential presence of hashes in a genome between two samples (specifically, looking for hashes that correlated with male vs female genomes).
to do this, we needed the following new functionality -
- code to export hashes, together with their abundances, from a signature computed with track-abund; viz #1098
- code to intersect one track-abund signature with another signature, without flattening the abundances in the first signature. (note that
sourmash sig intersectflattens all signatures) - this could be maybe be done by updatingsourmash sig intersect - code to select sequences in a FASTA/FASTQ file that have some number of overlapping hashes with a signature (this has been a repeatedly useful utility that I've implemented a dozen times in various contexts 😆 )
- code to estimate the abundance of sequences based on median hash abundance from a signature (i.e. estimate sequence abundance in a FASTA/FASTQ file using abundances from a track-abund sourmash signature) - this may be too niche to implement in sourmash directly, but I feel like it has come in handy.
I implemented all of this in a Jupyter notebook fairly easily, but it'd nice to have this in the sourmash CLI.
since code exists for all of this and I can make it available upon request, I'll label this as a good first issue...