Perl script to remove the duplicate sequences from multifasta file
use strict; use Bio::SeqIO; my %unique; my $file = "myseqs.fa"; my $seqio = Bio::SeqIO...display_id; my $seq = $seqs->seq; unless(exists($unique{$seq})) { $outseq->write_seq($seqs); $unique{$seq...2680 days ago
2065 days ago
Installing ggplot2 and its dependencies on Ubuntu !
jit@jit-HP-Pro-3335-MT:~/Downloads/MitoHunter/minidot/bin$ sudo R [sudo] password for jit: R version 3.4.4...ng -Wformat -Werror=format-security -Wdate-time -D_FORTIFY_SOURCE=2 -g -c unique.c -o unique.o g...1540 days ago
1392 days ago
Commandline for paired end reads simulation with BBMap !
(JitMetaENV) ➜ mixedSample git:(main) ✗ /home/urbe/Tools/bbmap/randomreads.sh ref=mixed.fa ou...rue Writing chunk 1 Waiting for writing to finish. Finished. snpRate=0.0, max=0, unique=true...980 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 su..._id.fa filtered_without_omi.fasta > omi_single_id_plus_all.fa #Extract unique kmer ./uniquekmer...818 days ago
784 days ago
Raku script to find microsatellites in DNA fragments !
sub find-microsatellites($sequence, $min-repeat-length = 2, $max-repeat-length = 6, $min-repea...es.push($substring); } } } return @microsatellites.unique; }...85 days ago