That's a really good idea!
• 1 views
•
link
written
by
pjmaguire3
how to locate a gene sequence among fastq files containing short reads
Hello!
I've got a dozen different strains of bacteria for which we've sequenced the whole genomes of (we have paired end reads - forward and reverse - for each strain). I wish to find and locate a specific house keeping gene in each strain.
Could I convert the fastq files into fasta files, set up a blast database containing the fasta short read files and then blast the query gene sequence against those? Or would I need to assemble each genome first and then make a database out of the assemblies and then blast the query gene sequence against those?
Would appreciate your input, thanks :-)
• 4,208 views
•
link
updated
by
Ram
•
written
by
jerrybug109 •
•
written
by
jerrybug109 •
1 answer
Using blast that way is very inefficient. If you are really impatient to see quick results and if you already have a sequence of the house keeping gene from the same species, than you may take this sequence as reference sequence and map all your reads on it with 'bwa mem'. If you can afford to wait about 5 minutes longer, you should assemble your reads with SPades.
After assembly, there is also no need to blast. It is much easier to map the contigs to the sequence of the house keeping gene with 'bwa mem'. You can even fed the contigs from all of your isolates into 'bwa mem' in a single run and you will get a nice little BAM file showing a multi sequence alignment of all the isolates comprising the house keeping gene. However, if it is really a house keeping gene, than it will be present in all of the isolates.
Log in to answer this question.
More posts like this
-
Blast locally on multiple fasta files with multiple database
written by fec2 •
-
Better method to find specific enzymes in whole genome sequences?
written by Mel161803 •
-
Separate the Longest Strain in Bacteria
written by Shelle •
-
How to merge BLAST XML output from different database?
written by tcf.hcdg
-
How to check for contamination by one contaminant across multiple whole genome sequenced samples?
written by jerrybug109 •
-
how to BLAST search for known genes in local set of genomes?
written by jerrybug109 •
-
what should I use for reference guided assembly of bacterial genomes intended for comparative genom…
written by jerrybug109 •
-
How to discover new mutations in a derived strain relative to a reference genome?
written by vladimir.gritsenko
-
Extract fasta from BLAST hits
written by biotech
-
blasting genome contigs against local SILVA 16S RNA database
written by robjohn70000
Don't do any of that .. yet. Make a "genome" with the gene(s) (if known or choose examples from related strains) you need and then align with BBMap. Depending on how similar "different" strains in your pool are there is some risk that reads may multimap. It sounds like you are just looking to see if a specific gene is there so go ahead and use option
ambig=allwith BBMap to allow reads to multi-map at all possible locations.You could also try using BBSplit to bin the reads if you have the reference genomes for these strains.