More posts like this
-
How to download fasta sequence of my interest from NCBI?
written by KumarHi all, I am trying to download the protein fasta sequence of "HSP listeria monocytogenes" from NCBI. Please suggest how I can download these all …
-
Number of Entries in Non Redundant DB
written by 6schulte •Hello, Has anyone an idea of wether there is some link to a `website` telling me `how many sequences` (as in entries) are currently to …
-
Discarding duplicate reads from fastq file on the basis of random barcode
written by naeem40thju •Say, I have used an oligo containing random barcode (Ns) for excluding the duplicated reads (PCR biased reads) Oligo: CCCTACACGACGCTCTTCCGATCTNNNNNNNNNNNNNNNNNNTGTAAAACGACGGCCAGT So, how can I exclude/remove …
-
Different results between remote and online blast
written by anasofiamoreira94I perform remote blast with this command line: blastn -task megablast -num_threads 6 -outfmt "6 std ssciname" -max_target_seqs 10 -db nt_v5 -query input_file -taxidlist Plants_DB …
-
Simple Tutorial failing
written by biodreck •Hi there, I found this tutorial: https://www.ncbi.nlm.nih.gov/books/NBK52637/ See this figure: https://www.ncbi.nlm.nih.gov/books/NBK52637/figure/blast_setup_pc.F5/?report=objectonly The database refseq_rna.00 is installed correctly When I do this command: blastdbcmd -db refseq_rna.00 …
-
how can I remove duplicated variants from vcf file?
written by kk.mahsaHow can I remove duplicated variants from vcf file? I googled and searched in biostars history but I did not find any way to do …
-
comparing two lists using command line tools
written by Farbod**Dear Friends**, Hi (***Sorry** if this question is simple or duplicated*) I have two lists of IDs, list A.txt and list B.txt and B>A (A=blast …
-
Search a proteome file to blast database (command line usage)
written by MehmetDear All: I have a proteome file generated from Next Gen. Sequencing. I want to search my proteome file against to blast database using command …
-
Command line BLAST
written by KumarHow to align Accession ids with NCBI database using command line BLAST. Thanks
-
How To Show More Blast Results Using Biopython?
written by Niek De Klein<p>Hi,</p> <p>I'm using biopython to BLAST over the internet. However, it only saves 30 results (there are more than 30 results that are under the …
define "duplicatedentries". Hsp ? Hit ?
thanks,
Could you describe this in command format.
Sure.
Yes, i mean duplicate hits. as an example, like
>gi|10001|,>gi|1000|such as duplicate entries..Are these duplicate hits seen for the same query? If they are seen for different queries, they are not duplicates - just asking so we are closer to defining the problem.
Yes..these are the duplicate hits that generated by same hits..is any way at command line BLAST to remove these type of all duplicated hits from result file.
I'm sorry, duplicate hits by same hits? Are you referring to multiple HSPs perhaps? Remember, each query can have multiple hits. Each hit can have multiple HSPs. Can you give us a sample maybe?
I don't think BLAST would give you a multiple hits where the subject as well as the query are the same sequence. BLAST would just specify them as a single hit with multiple HSPs, unless we are dealing with repeat sequences, in which case repeatmasker could be a prudent pre-processing step.
Hey RamRS.
I have an HSPs table. There are some cases in which queries of different length hit regions of the subject that overlap. The table has a column with the coordinates per HPS. I am interested in removing the redundant HPSs, is there a script available I can use. Thanks in advance for your reply.