Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

match to a query is missing from branchwater Web site #25

Open
ctb opened this issue May 3, 2024 · 3 comments
Open

match to a query is missing from branchwater Web site #25

ctb opened this issue May 3, 2024 · 3 comments

Comments

@ctb
Copy link

ctb commented May 3, 2024

The accession BK010471 is for a crAssphage that is ubiquitous in human gut metagenomes (link), and in particular is found in the 454 data set SRR073439.

When I do a containment search, I see:

% sourmash search --containment BK010471.fa.sig SRR073439.sig -k 31

selecting specified query k=31
loaded query: BK010471.fa... (k=31, DNA)
--
loaded 3 total signatures from 1 locations.
after selecting signatures compatible with search, 1 remain.

1 matches above threshold 0.080:
similarity   match
----------   -----
 59.0%       SRR073439

and the Venn diagram is pleasing:

venn2

However, the FASTA sequence does not have any matches when searched at https://branchwater.jgi.doe.gov/. Any ideas?

thanks!

SRR073439.k31.sig.zip
BK010471.k31.sig.zip
BK010471.fa.zip

@luizirber
Copy link
Member

luizirber commented May 3, 2024

This dataset is not in the index, likely due to it having "amplicon" in the metadata (abstract):

... datasets of sequenced bacterial 16S rRNA gene amplicons and total fecal ...

https://trace.ncbi.nlm.nih.gov/Traces/index.html?view=run_browser&acc=SRR073439&display=metadata

which connects with the discussion in #24 (comment) =]

@ctb
Copy link
Author

ctb commented May 3, 2024

It's in wort,

/group/ctbrowngrp/irber/data/wort-data/wort-sra/sigs/SRR073439.sig

@luizirber
Copy link
Member

It's in wort,

/group/ctbrowngrp/irber/data/wort-data/wort-sra/sigs/SRR073439.sig

Yes, see footnotes 1 and 2 here:
#24 (comment)
#24 (comment)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants