Perl script to check modules installed !
#use strict; use warnings; use ExtUtils::Installed; # # First, check if all the required modules have been installed in the system # BEGIN { my @import_modules...1207 days ago
1106 days ago
Install grabseqs using conda !
vik@vik-Lenovo-ideapad-320-15ISK:~/Downloads/setu/setu$ conda install grabseqs -c louiejtaylor -c bioconda -c conda-forge Collecting package metadata (current_repodata.json): done Solving environment: done ==> WARNING: A newer version of conda exists.1106 days ago
Perl onliner to print fasta headers !
#Save all your fasta in seq.fa and run the following ... perl -ne 'print if /^>/' seq.fa #Print header with line number perl -ne 'print "$. $_" if /^>/ ' seq.fa1048 days ago
Oneliner to convert lower-case to sequence masked with Ns
perl -pe '/^[^>]/ and $_=~ s/[a-z]/N/g' genomic.fna > genomic.N-masked.fna awk '{if(/^[^>]/)gsub(/[a-z]/,"N");print $0}' genomic.fna > genomic.N-masked.fna1002 days ago
Perl script for Smith-Waterman Algorithm
# Smith-Waterman Algorithm # usage statement die "usage: $0 \n" unless @ARGV == 2; # get sequences from command line my ($seq1, $seq2) = @ARGV; # scoring scheme...1002 days ago
Perl onliner to check the ids in two files !
perl -lane 'BEGIN{open(A,"ids2.txt"); while(){chomp; $k{$_}++}} if (defined($k{$F[0]})) {print "$_\t$F[0]\t1"} else {print "$_\tNA\t0"}; ' ids1.txt > aaa.xls979 days ago
Installing Covid19 Environment !
(base) vikas@vikas-Lenovo-ideapad-320-15ISK:~/vinodLab/Genepi$ conda env create -f covid19-environment.yml Collecting package metadata (repodata.json): done Solving enviro...975 days ago
Perl script to rename the fasta file !
#Script #1 #!/usr/bin/perl -w use strict; #USAGE #perl extractPattern.pl kmerfasta > uniref100_result_broad my %kHash; local $/ = '>'; my $infile2 = "$ARGV[0]"; # Kmer fasta open( FH2, '910 days ago
Extract all fasta sequences except ids !
awk 'BEGIN{while((getline0)l[">"$1]=1}/^>/{f=!l[$1]}f' genomic.fna > filtered_without_omi.fasta #extract subseq seqtk subseq omi_ids.fa omi_single_id.txt > omi_single_id...880 days ago