Blast script to index and extract sequence !!
# look at the file $ head EC4115.fa >NC_011353.1 Escherichia coli O157:H7 str. EC4115 chromosome,...s # retreive an entry by id $ blastdbcmd -db EC -entry 'NC_011353.1' | head >lcl|NC_011353.1 Escherichia...2857 days ago
Download the gff files from NCBI using bash script/command
...ate a directory with appending G mkdir "GFF$filename" cd "GFF$filename" # take action on each file. $f store current file name head -n 4 $f > $f.head wget --input $f.head gu...2528 days ago
Calculate Dinucleotide Frequency with Perl
#!/usr/bin/perl -w use strict; my ($genome, $head, $tail); my (%mono_nt, %di_nt); $/ = ">"; open my $fasta, '2343 days ago
2474 days ago
Perl script to read multi fasta sequence one by one
...die "couldn't open the file $outfile $!"; print OUT "$key\n$fastaSeq{$key}\n"; } sub readfasta { (my $file)=@_; my %sequence; my $header; my $temp_seq; #supp...2226 days ago
Perl subroutine for reading multifasta file !
sub readfasta { (my $file)=@_; my %sequence; my $header; my $temp_seq; #suppose fasta files contains multiple sequences; open (IN, "2149 days ago
To convert just one specific read group to fastq
...ep only valid SRR numbers. cat runinfo.txt | cut -f 1 -d , | grep SRR | head -$N > selected.txt # Store...mtools merge -f all.bam bam/*.bam # Investigate the readgroups in the header. echo "" echo "SAM file h...1547 days ago
1529 days ago
Pack a perl program with their dependencies on Ubuntu !
#Follow steps to create your own executable ./web jit@jit-HP-Pro-3335-MT:~/Downloads/autoConTAMPR/bin$ sudo apt install libpar-packer-perl Reading package lists......1516 days ago
Script to extract the cluster detail !
$ lsb_release -a No LSB modules are available. Distributor ID: Ubuntu Description: Ubuntu 18.04.1 LTS Release: 18.04 Codename: bionic $ cat /proc/cpuinfo | grep -i 'model name' | head -n 1 model name : Intel(R) Xeon(R) CPU E5-2620 v4 @ 2.10GHz1383 days ago