Skip to content

how to query? #39

Description

@shenweima

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

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions