Hi Perlmonks
This is for a bioinformatics assignment that I have to complete. I'm aware there are forums for this and I have posted there but still have no reply.
I need help with getting the output to look the way I want it. I am comparing two strings (sequences) with one across the top and the other down the side with dots printed everywhere there's a match.
The problem is that the sequences are too long to fit on one line. Is there a way, whether it be in the code itself or in the UNIX terminal, to get it to display within the terminal space. The letters don't necessarily need to be legible. It's the pattern of dots that is important. I would appreciate it if you could take a look at my code and let me know what changes I could make to achieve this.
#!/usr/bin/perl
use strict;use warnings;
my $i;my $j;my $a;
open (FH1,'<', "./CAM3UCA.txt");
open (FH2,'<', "./TCH3A.txt");
my $seq1 = do { local $/; <FH1> };
my $seq2 = do { local $/; <FH2> };
my @s1=split('',$seq1);
my @s2=split('',$seq2);
my @matrix=();
for($i=0;$i<scalar(@s2);$i++)
{
for($j=0;$j<scalar(@s1);$j++)
{
if($s2[$i] eq $s1[$j])
{
$matrix[$i][$j]=".";
}
else
{
$matrix[$i][$j]=" ";
}
}
}
#Printing matrix of dot plot
open (FHOUT,">dotplot.txt") or die "cannot open outfile\n";
for($a=0;$a<scalar(@s1);$a++)
{
print FHOUT "$s1[$a]";
}
print FHOUT "\n";
#print"\n\n";
for($i=0;$i<scalar(@s2);$i++)
{
print FHOUT "$s2[$i]";
for($j=0;$j<scalar(@s1);$j++)
{
print FHOUT "$matrix[$i][$j]";
}
print FHOUT "\n";
#print"\n\n";
}
close FH1;
close FH2;
close FHOUT;
Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!
Titles consisting of a single word are discouraged, and in most cases are disallowed outright.
Read Where should I post X? if you're not absolutely sure you're posting in the right place.
Please read these before you post! —
Posts may use any of the Perl Monks Approved HTML tags:
- a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr
You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)
| |
For: |
|
Use: |
| & | | & |
| < | | < |
| > | | > |
| [ | | [ |
| ] | | ] |
Link using PerlMonks shortcuts! What shortcuts can I use for linking?
See Writeup Formatting Tips and other pages linked from there for more info.