Kraken: ultrafast metagenomic sequence classification using exact alignments
...BLAST program. In its fastest mode, Kraken classifies 100 base pair reads at a rate of over 4.1 million reads per minute, 909 times faster than Megablast and 11 times faster than the...2832 days ago
2625 days ago
BASH script for SelfBLAST a genome
...blast" SEQ=$FASTAFILE else echo "Something went wrong $USER - Contact jitendra" fi echo "Doing alignments -- BLASting"; blastn -task megablast -query $SEQ -db $MYDB -evalue...2615 days ago
Read a tab delimited file and search with perl
use strict; use warnings; use Data::Dumper; use Text::CSV; use IO::Handle; my $file = "/home/urbe/Tools/Alienomics_v0.1/Alienomics/output/intermediate_files/rRNA/refGene.megablast"; open my $fh, "[0]\n"; warn Dumper $row; # To see the structure }2488 days ago
How to print only the first best hit from blastn tool?
I used the following command to print the blast hit. [Terminal] blastn -query aaa.fa -db blastdb/nt -evalue 10e-5 -num_threads 4 -max_target_seqs 1 -outfmt '6 qseqid staxid qstart qend sallseqi sstart send evalue length frames qcovs' -out testBlastSee But the outfile print many hits, I only...Tags: blastn, blastall, megablast, top, hits
2584 days ago
Converting BLAST output into CSV
Suppose we wanted to do something with all this BLAST output. Generally, that’s the case - you want to retrieve all matches, or do a reciprocal BLAST, or something. As with most programs that run on UNIX, the text output is in some specific format. If the program is popular enough, there w...Tags: Convert, BLAST, output, CSV, blastn, megablast
2300 days ago
2597 days ago