#Find the reverse complement of a DNA string.
#Given: A DNA string Pattern.
#Return: Pattern, the reverse complement of Pattern.

use strict;
use warnings;

my $string="AAAACCCGGT";
my $finalString="";
my %hash = (
    "C" => "G",
    "A" => "T",
    "T" => "A",
    "G" => "C",
);

for (my $aa=0; $aa<=(length($string)-1); $aa++) {
    my $char=substr $string, $aa, 1;
    #print $hash{$char};
    $finalString="$hash{$char}"."$finalString";
}

print $finalString;
print "\n";

Script are moved to http://bioinformaticsonline.com/snippets/view/34633/clump-finding-problem-solved-with-perl

Address of the bookmark: http://engr.case.edu/li_jing/papers/00798gpattern.pdf

https://github.com/roblanf/minion_qc

Address of the bookmark: https://github.com/roblanf/minion_qc

• -e Makes the line of code be executed instead of a script
• -n Forces your line to be called in a loop. Allows you to take lines from the diamond operator (or stdin)
• -p Forces your line to be called in a loop. Prints $_ at the end

• This counts the number of quotation marks in each line and prints it

  perl -ne '$cnt = tr/"//;print "$cnt\n"' inputFileName.txt

• Adds string to each line, followed by tab

  perl -pe 's/(.*)/string\t$1/' inFile > outFile

• Append a new line to each line

  perl -pe 's//\n/' all.sent.classOnly > all.sent.classOnly.sep

• Replace all occurrences of pattern1 (e.g. [0-9]) with pattern2

  perl -p -i.bak -w -e 's/pattern1/pattern2/g' inputFile

• Go through file and only print words that do not have any uppercase letters.

  perl -ne 'print unless m/[A-Z]/' allWords.txt > allWordsOnlyLowercase.txt

• Go through file, split line at each space and print words one per line.

  perl -ne 'print join("\n", split(/ /,$_));print("\n")' someText.txt > wordsPerLine.txt

• or in other words, delete every character that is not a letter, white space or line end (replace with nothing)

  perl -pne 's/[^a-zA-Z\s]*//g' text_withSpecial.txt > text_lettersOnly.txt

• perl -pne 'tr/[A-Z]/[a-z]/' textWithUpperCase.txt > textwithoutuppercase.txt;

• Print only the second column of the data when using tabular as a separator

  perl -ne '@F = split("\t", $_); print "$F[1]";' columnFileWithTabs.txt > justSecondColumn.txt

• One-Liner: Sort lines by their length

  perl -e 'print sort {length $a <=> length $b} <>' textFile

• One-Liner: Print second column, unless it contains a number

  perl">perl -lane 'print $F[1] unless $F[1] =~ m/[0-9]/' wordCounts.txt

install and load required packages

  require(RColorBrewer)
  require(ComplexHeatmap)
  require(circlize)
  require(digest)
  require(cluster)

If all load successfully, proceed to Part 3. Otherwise, go through the following code chunks in order to ensure that each package is installed and loaded properly.

BiocManager (Morgan (2019))

Address of the bookmark: https://github.com/kevinblighe/E-MTAB-6141

~$ ./sviper -s short-reads.bam -l long-reads.bam -r ref.fa -c variants.vcf -o polished_variants

This will output a polished_variants.vcf file, that contains all the refined variants.

Sometimes it is helpful to look at the polished sequence, e.g. with the IGV browser. In that case you want SViper to output the polished and aligned sequences in a bam file via the option --output-polished-bam:

~$ ./sviper -s short-reads.bam -l long-reads.bam -r ref.fa -c variants.vcf -o polished_variants --output-polished-bam

Address of the bookmark: https://github.com/smehringer/SViper

https://genome.sph.umich.edu/wiki/Vt

https://github.com/atks/vt

Address of the bookmark: https://genome.sph.umich.edu/wiki/Vt

Address of the bookmark: https://bernatgel.github.io/karyoploter_tutorial/