Bioinformatics/base count: Difference between revisions

=={{header|VBScript}}==

<lang vb>

b=_

"CGTAAAAAATTACAACGTCCTTTGGCTATCTCTTAAACTCCTGCTAAATG" &_

"CTCGTGCTTTCCAATTATGTAAGCGTTCCGAGACGGGGTGGTCGATTCTG" &_

"AGGACAAAGGTCAAGATGGAGCGCATCGAACGCAATAAGGATCATTTGAT" &_

"GGGACGTTTCGTCGACAAAGTCTTGTTTCGAGAGTAACGGCTACCGTCTT" &_

"CGATTCTGCTTATAACACTATGTTCTTATGAAATGGATGTTCTGAGTTGG" &_

"TCAGTCCCAATGTGCGGGGTTTCTTTTAGTACGTCGGGAGTGGTATTATA" &_

"TTTAATTTTTCTATATAGCGATCTGTATTTAAGCAATTCATTTAGGTTAT" &_

"CGCCGCGATGCTCGGTTCGGACCGCCAAGCATCTGGCTCCACTGCTAGTG" &_

"TCCTAAATTTGAATGGCAAACACAAATAAGATTTAGCAATTCGTGTAGAC" &_

"GACCGGGGACTTGCATGATGGGAGCAGCTTTGTTAAACTACGAACGTAAT"

s="SEQUENCE:"

acnt=0:ccnt=0:gcnt=0:tcnt=0

for i=0 to len(b)-1

if (i mod 30)=0 then s = s & vbcrlf & right(" "& i+1,3)&": "

if (i mod 5)=0 then s=s& " "

m=mid(b,i+1,1)

s=s & m

select case m

case "A":acnt=acnt+1

case "C":ccnt=ccnt+1

case "G":gcnt=gcnt+1

case "T":tcnt=tcnt+1

case else

wscript.echo "error at ",i+1, m

end select

next

wscript.echo s & vbcrlf

wscript.echo "Count: A="&acnt & " C=" & ccnt & " G=" & gcnt & " T=" & tcnt

</lang>

{{out}}

<pre>

SEQUENCE:

1: CGTAA AAAAT TACAA CGTCC TTTGG CTATC

31: TCTTA AACTC CTGCT AAATG CTCGT GCTTT

61: CCAAT TATGT AAGCG TTCCG AGACG GGGTG

91: GTCGA TTCTG AGGAC AAAGG TCAAG ATGGA

121: GCGCA TCGAA CGCAA TAAGG ATCAT TTGAT

151: GGGAC GTTTC GTCGA CAAAG TCTTG TTTCG

181: AGAGT AACGG CTACC GTCTT CGATT CTGCT

211: TATAA CACTA TGTTC TTATG AAATG GATGT

241: TCTGA GTTGG TCAGT CCCAA TGTGC GGGGT

271: TTCTT TTAGT ACGTC GGGAG TGGTA TTATA

301: TTTAA TTTTT CTATA TAGCG ATCTG TATTT

331: AAGCA ATTCA TTTAG GTTAT CGCCG CGATG

361: CTCGG TTCGG ACCGC CAAGC ATCTG GCTCC

391: ACTGC TAGTG TCCTA AATTT GAATG GCAAA

421: CACAA ATAAG ATTTA GCAAT TCGTG TAGAC

451: GACCG GGGAC TTGCA TGATG GGAGC AGCTT

481: TGTTA AACTA CGAAC GTAAT

Count: A=129 C=97 G=119 T=155

</pre>