kmtricks icon indicating copy to clipboard operation
kmtricks copied to clipboard

how to query?

Open shenweima opened this issue 6 months ago • 0 comments

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

shenweima avatar Aug 02 '24 02:08 shenweima