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

how to query? #39

Open
shenweima opened this issue Aug 2, 2024 · 0 comments
Open

how to query? #39

shenweima opened this issue Aug 2, 2024 · 0 comments

Comments

@shenweima
Copy link

I have a kmer which has known present or absent status, like below:

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAGCCTCCAT 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 1 0 0 1 0 0 1 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 1 0 0 1 0 1 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0

I bulid index using the same input fastq files.

kmtricks pipeline --file rice.blast.files.list.fof --kmer-size 51 -t 48 --run-dir ./blast.index --mode hash:bft:bin --hard-min 2

Then I want to search the kmer AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAGCCTCCAT in this index. Below is a part of the result. All samples were same.

# FORMAT:
# * [query name]
# For each target, 3 or 4 fields:
#   [taget name]
#   (unless --no-detail option) string in {+-} showing positions covered (+) by at least a shared kmer, else (-)
#   Ratio of kmers of the query shared with the target
#   Ratio of positions of the query covered by at least a kmer shared with the target
* [test.txt1] 
[ro1-1] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[rb22] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[Z8] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[Z5] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[Z1] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[YN187] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[Y9] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[Y8] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[Y7] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[Y6] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[Y5] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[Y4] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[Y3] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[Y2] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[Y12] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[Y11] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[Y10] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[Y1] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[T8] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[T6] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[T4] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[T3] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[T2] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[T1] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[SRR6669224] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[SRR6669223] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[SRR6669222] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[SRR6669221] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[SRR6669220] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[SRR6669219] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[SRR6669218] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[SRR6669217] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[SRR6669216] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[SRR6669215] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[SRR6669214] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
[SRR6669213] +++++++++++++++++++++++++++++++++++++++++++++++++++ 1.00 1.00
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

1 participant