Finding all positions for all kmers belonging to a (not so short) list of kmers in a large genome
Hi,
Is there a tool for not only enumerating the counts of kmers (like jellyfish) but also will list their positions? I know its much more exhausting memory-wise but I'm looking for the best way to do this, even if a tool doesn't exist currently.
Thanks!
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Get VCF-format positions from coding positions
written by captainlabman •
-
Java -xmx option does not limit memory usage in the ImputePipelinePlugin?
written by twrl8 •
-
How can I get the positions where a sequence was aligned as output of a pair-wise alignment?
written by João •
-
Is there a way to do GO terms enrichment analysis and keep track of gene symbols?
written by DN99 •
-
Best way to obtain reference and non-reference k-mer
written by michele.schimd •
-
samtools sort by something other than position/name?
written by anderspitman •
-
Tools for finding gene clusters in RNA-seq differential expression data?
written by am •
-
ConsensusPathDB enrichment analysis: input data normalization?
written by abe •
-
alternatives to meme suite
written by kws15 •
-
Does RefSeq DB have a CDS file for its transcripts?
written by pwg46
Take a look at Finding 16 mer not present in GRCh38. In this a suggestion was to use
bowtieto align the kmers against the genome. I would do the alignment and then filter for matches with 100% sequence identity. It might help to set gap opening and mismatch penalties to like 10000 to only retain perfect matches.Is that really faster than for example implementing a search trie?
How large your Kmers?, all combinations?, all occurrences? I used to code some scripts in Perl for kmer counting (8-12 kmers) with their position for cis-regulatory elements in some plant genomes, so it is not hard to do, even on the 2 GB RAM machine I had.