Bioinformatics/base count: Difference between revisions
Content added Content deleted
No edit summary |
|||
Line 2,004: | Line 2,004: | ||
Base counts: A:129, C:97, G:119, T:155, total:500 |
Base counts: A:129, C:97, G:119, T:155, total:500 |
||
</pre> |
</pre> |
||
=={{header|Picat}}== |
|||
<lang Picat>main => |
|||
dna(DNA, ChunkSize), |
|||
Count = 0, |
|||
println("Sequence:"), |
|||
Map = new_map(['A'=0,'C'=0,'G'=0,'T'=0]), |
|||
foreach(Chunk in DNA.chunks_of(ChunkSize)) |
|||
printf("%4d: %s\n", Count, Chunk), |
|||
Count := Count + Chunk.len, |
|||
foreach(C in Chunk) |
|||
Map.put(C,Map.get(C)+1) |
|||
end |
|||
end, |
|||
println("\nBase count:"), |
|||
foreach(C in "ACGT") |
|||
printf("%5c: %3d\n", C, Map.get(C)) |
|||
end, |
|||
printf("Total: %d\n", Count), |
|||
nl. |
|||
dna(DNA,ChunkSize) => |
|||
DNA = "CGTAAAAAATTACAACGTCCTTTGGCTATCTCTTAAACTCCTGCTAAATG |
|||
CTCGTGCTTTCCAATTATGTAAGCGTTCCGAGACGGGGTGGTCGATTCTG |
|||
AGGACAAAGGTCAAGATGGAGCGCATCGAACGCAATAAGGATCATTTGAT |
|||
GGGACGTTTCGTCGACAAAGTCTTGTTTCGAGAGTAACGGCTACCGTCTT |
|||
CGATTCTGCTTATAACACTATGTTCTTATGAAATGGATGTTCTGAGTTGG |
|||
TCAGTCCCAATGTGCGGGGTTTCTTTTAGTACGTCGGGAGTGGTATTATA |
|||
TTTAATTTTTCTATATAGCGATCTGTATTTAAGCAATTCATTTAGGTTAT |
|||
CGCCGCGATGCTCGGTTCGGACCGCCAAGCATCTGGCTCCACTGCTAGTG |
|||
TCCTAAATTTGAATGGCAAACACAAATAAGATTTAGCAATTCGTGTAGAC |
|||
GACCGGGGACTTGCATGATGGGAGCAGCTTTGTTAAACTACGAACGTAAT".delete_all('\n'), |
|||
ChunkSize = 50.</lang> |
|||
{{out}} |
|||
<pre>Sequence: |
|||
0: CGTAAAAAATTACAACGTCCTTTGGCTATCTCTTAAACTCCTGCTAAATG |
|||
50: CTCGTGCTTTCCAATTATGTAAGCGTTCCGAGACGGGGTGGTCGATTCTG |
|||
100: AGGACAAAGGTCAAGATGGAGCGCATCGAACGCAATAAGGATCATTTGAT |
|||
150: GGGACGTTTCGTCGACAAAGTCTTGTTTCGAGAGTAACGGCTACCGTCTT |
|||
200: CGATTCTGCTTATAACACTATGTTCTTATGAAATGGATGTTCTGAGTTGG |
|||
250: TCAGTCCCAATGTGCGGGGTTTCTTTTAGTACGTCGGGAGTGGTATTATA |
|||
300: TTTAATTTTTCTATATAGCGATCTGTATTTAAGCAATTCATTTAGGTTAT |
|||
350: CGCCGCGATGCTCGGTTCGGACCGCCAAGCATCTGGCTCCACTGCTAGTG |
|||
400: TCCTAAATTTGAATGGCAAACACAAATAAGATTTAGCAATTCGTGTAGAC |
|||
450: GACCGGGGACTTGCATGATGGGAGCAGCTTTGTTAAACTACGAACGTAAT |
|||
Base count: |
|||
A: 129 |
|||
C: 97 |
|||
G: 119 |
|||
T: 155 |
|||
Total: 500</pre> |
|||
=={{header|PicoLisp}}== |
=={{header|PicoLisp}}== |