#!/usr/bin/perl

$DNA_Fragment_1 = 'AAAGGCGCGCGCGCGCGCGTATAGCGCGCGCGCGAAAG';
$DNA_Fragment_1_Laenge = length($DNA_Fragment_1);
$DNA_Fragment_1_GC_Paare = 12;
$GC_Gehalt_Frag_1 = int((100/$DNA_Fragment_1_Laenge)*$DNA_Fragment_1_GC_Paare);

$DNA_Fragment_2 = 'attatttatatagcgcg';
$DNA_Fragment_2_Laenge = length($DNA_Fragment_2);
$DNA_Fragment_2_GC_Paare = 2;
$GC_Gehalt_Frag_2 = int((100/$DNA_Fragment_2_Laenge)*$DNA_Fragment_2_GC_Paare);

print "Sequenz 1:\t$DNA_Fragment_1\n";
$L_1 = length($DNA_Fragment_1);
print "Laenge der Sequenz 1:\t$L_1\n";
print "Anzahl der GC-Paare in Sequenz 1:\t$DNA_Fragment_1_GC_Paare\n";
print "GC-Gehalt in %:\t$GC_Gehalt_Frag_1\n\n";


print "Sequenz 2:\t$DNA_Fragment_2\n";
$L_2 = length($DNA_Fragment_2);
print "Laenge der Sequenz 2:\t$L_2\n";
print "Anzahl der GC-Paare in Sequenz 2:\t$DNA_Fragment_2_GC_Paare\n";
print "GC-Gehalt in %:\t$GC_Gehalt_Frag_2\n\n";


sub KONKAT {
$Konkat_Sequenz = $DNA_Fragment_1.$DNA_Fragment_2;
$Differenz_Frag_1_Frag_2 = $DNA_Fragment_1_Laenge - $DNA_Fragment_2_Laenge;
$Konkate_Frag_1_Frag_2 = $DNA_Fragment_1_Laenge + $DNA_Fragment_2_Laenge;
$Konkate_GC_Paare = $DNA_Fragment_1_GC_Paare + $DNA_Fragment_2_GC_Paare;
$GC_Gehalt_Konkat = int((100/$Konkate_Frag_1_Frag_2)*$Konkate_GC_Paare);

print "Konkatenierung Fragment1 - Fragment2:\n\n";
print "Sequenz: \t$Konkat_Sequenz\n";
print "Laenge: \t$Konkate_Frag_1_Frag_2\n";
print "Anzahl der GC Paare: \t$Konkate_GC_Paare\n";
print "GC-Gehalt (in % der Gesamtlaenge): \t$GC_Gehalt_Konkat\n";
print "Laengendifferenz Fragment1 - Fragment2: \t$Differenz_Frag_1_Frag_2\n\n";
}  ;