2057 days ago
Perl script to reverse complement a DNA sequence !
#!/usr/bin/perl -w $DNA = 'ACGGGAGGACGGGAAAATTACTACGG...print "Here is the starting DNA:\n\n"; print "$DNA\n\n"; $revcom = reverse $DNA; $revcom =~ s/A/T/g; $rev...ere is the reverse complement DNA: WRONG:\n\n"; print "$revc...ere is the reverse complement DNA:\n\n"; print "$revcom\n";...2040 days ago
Perl One-Liner to print only non-uppercase letters
#Go through file and only print words that do not have any uppercase letters. perl -ne 'print unless m/[A-Z]/' dna.fa > dnaOnlyLowercase.fa #To lowercase everything perl -pne 'tr/[A-Z]/[a-z]/' dnaUpperCase.fa >dnawithoutuppercase.fa;1380 days ago
1337 days ago
Commandline for paired end reads simulation with BBMap !
(JitMetaENV) ➜ mixedSample git:(main) ✗ /home/urbe/Tools/bbmap/randomreads.sh ref=m...h=250, mininsert=400, maxinsert=600, gaussian] Writing reference. Executing dna.FastaToChromArrays2 [mixed.fa...987 days ago
834 days ago
Perl script to find edit distance between two sequences !
#!/usr/bin/perl use strict; use warnings; sub edit_distance { my ($s1, $s2) = @_; my $len1 = length($s1); my $len2 = length($s2); my @dp; for (my $i = 0; $i421 days ago
Raku script to find palindrome in genomes !
sub is-palindrome(Str $str) returns Bool {....flip; } sub find-palindromes(Str $dna, Int $min-length, Int $max-le...-length -> $length { for 0..^$dna.chars - $length -> $pos {...} } } # Example usage my $dna = "GGATCCATGGCCTAGG"; # examp...421 days ago
Raku script to calculate GC content !
sub calculate-gc-content(Str $sequence) { my $gc-cou...return $gc-count / $total-bases * 100; } my $dna_sequence = "ATGCGCTAAAGCGCGCG...GCGCGCGCGC"; my $gc_content = calculate-gc-content($dna_sequence); say "DNA Sequen...108 days ago
91 days ago