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...2643 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 }2516 days ago
Script to rapid genome clustering based on pairwise ANI
First, create a blast+ database: makeblastdb -in -dbtype nucl -out Next, use megablast from blast+ package to perform all-vs-all blastn of sequences: blastn -query -db -outfmt '...627 days ago