Is there a command to output the kmers of each sequence in a multifasta file?
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Trouble parsing header
Dear BBMap team:
I tried to use filterbytile.sh to remove the reads with low quality, but I encountered an error message saying that there was a trouble parsing the header. I've read the description of the script and Brian Bushnell said that was possible when the reads were renamed (such as in SRA) and to contact him if such error happened.
I downloaded the sequencing data (SRA) from ncbi and used fastq-dump to get the fastq files. I wonder if there is a solution to this?
Thank you very much!
Rose
Comment
-
BBsketch alltoall is incomplete
Can I ask a question about bbsketch?
I want to compare the ANI between many genomes (1000+) to each other.
I did
Code:bbsketch.sh perfile genome_folder/*.fasta out=sketch.gz k=31,24 threads=16 comparesketch.sh alltoall sketch.gz k=31,24 prealloc=0.75 format=3 threads=16 out=table.tsv
Code:Set threads to 16 Loading sketches. Loaded 1157 sketches in 59.541 seconds. Total Time: 59.784 seconds.
Code:Set threads to 16 Loading sketches. Executing kmer.KmerTableSet [ways=31, tabletype=10, prealloc=0.75] Initial size set to 45218398 Initial: Ways=31, initialSize=45218398, prefilter=f, prealloc=0.75 Memory: max=91268m, total=91268m, free=90848m, used=420m 3.713 seconds. Indexed 2880884 unique and 10513099 total hashcodes. Loaded 1157 sketches in 8.457 seconds. Ran 1225005 comparisons in 9.344 seconds. Total Time: 17.801 seconds.
- Genomes are highly similar.
#Query Ref ANI QSize RefSize QBases RBases QTaxID RTaxID KID WKID SSU
genome1.fasta genome2.fasta 94.223 1984118 1796930 1987598 1797650 -1 -1 24.952 27.523 .
- It is not simply due to the naming: I neither find "genome1 vs genome2" nor "genome 2 vs genome1"
Any idea?
Comment
-
I'm trying to use BBmap to find all perfect hits or hits with an indel length 1.
Code:bbmapskinner.sh in=kmer.fasta out=result.sam ambiguous=all strictmaxindel=1
Is there something that I am doing wrong?
Comment
Latest Articles
Collapse
-
by seqadmin
Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.
Somatic Genomics
“We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...-
Channel: Articles
05-24-2024, 01:16 PM -
-
by seqadmin
The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...-
Channel: Articles
05-06-2024, 07:48 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 05-24-2024, 07:15 AM
|
0 responses
16 views
0 likes
|
Last Post
by seqadmin
05-24-2024, 07:15 AM
|
||
Started by seqadmin, 05-23-2024, 10:28 AM
|
0 responses
18 views
0 likes
|
Last Post
by seqadmin
05-23-2024, 10:28 AM
|
||
Started by seqadmin, 05-23-2024, 07:35 AM
|
0 responses
22 views
0 likes
|
Last Post
by seqadmin
05-23-2024, 07:35 AM
|
||
Started by seqadmin, 05-22-2024, 02:06 PM
|
0 responses
11 views
0 likes
|
Last Post
by seqadmin
05-22-2024, 02:06 PM
|
Comment