Selecting sequences from a multi-fasta with certain kmers
Hi,
I have a fasta file with about 30 9-mers. I also have some metagenomes in fasta format which are about 15 GB (~100 million reads, 125 bp) I would like to use my kmer file to filter out sequences from my metagenomes that only have a match to one of the kmers.
There's a lot of kmer counting software out there like jellyfish, tallymer, meryl, but as far as I can tell, none of these have the utility to select or filter sequences based on the presence of kmers.
Does anyone know of software that would do this efficiently?
Thanks
• 2,728 views
•
link
written
by
michberr8 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
fasta sequences into non-overlapping k-mers using python
written by smrutimayipanda •
-
Interpreting Tallymer search output
written by jmartin •
-
Genome Size Estimation with Jellyfish and Genome Scope is Unexpectedly Small
written by dk0319 •
-
Have a question about run jellyfish in many directories with different number of zipped fasta files
written by schlogl
-
Kmer enrichment analysis
written by chipolino •
-
Can we use JELLYFISH for amino acid seqeunces?
written by marksingh1982 •
-
minia assembler - counting kmers 'gatb::core::system::Exception'
written by thackl
-
Normalization of certain gene abundance in a metagenome
written by bioinfo
-
Generate kmer profiles from a bunch of peptide sequences
written by Owen S.
-
Tool: Jellyfish - Fast, Parallel K-Mer Counting For Dna
written by Istvan Albert
BBduk from BBMap should be able to do this.
do you mean
ONE hit match => discard SAME kmer found twice => keep TWO different kmer => keep
?
Sorry, i meant one or more matches to my set of 30 kmers