36 lines
1.6 KiB
Plaintext
36 lines
1.6 KiB
Plaintext
SeedRandom[13122345];
|
|
seq = BioSequence["DNA", "ATAAACGTACGTTTTTAGGCT"];
|
|
randompos = RandomInteger[seq["SequenceLength"]];
|
|
seq = StringReplacePart[seq, RandomChoice[{"A", "T", "C", "G"}], {randompos, randompos}];
|
|
randompos = RandomInteger[seq["SequenceLength"]];
|
|
seq = StringReplacePart[seq, "", {randompos, randompos}];
|
|
randompos = RandomInteger[seq["SequenceLength"]];
|
|
seq = StringInsert[seq, RandomChoice[{"A", "T", "C", "G"}], randompos];
|
|
seq = BioSequence["DNA", StringJoin@RandomChoice[{"A", "T", "C", "G"}, 250]];
|
|
size = 50;
|
|
parts = StringPartition[seq["SequenceString"], UpTo[size]];
|
|
begins = Most[Accumulate[Prepend[StringLength /@ parts, 1]]];
|
|
ends = Rest[Accumulate[Prepend[StringLength /@ parts, 0]]];
|
|
StringRiffle[MapThread[ToString[#1] <> "-" <> ToString[#2] <> ": " <> #3 &, {begins, ends, parts}], "\n"]
|
|
Tally[Characters[seq["SequenceString"]]]
|
|
Do[
|
|
type = RandomChoice[{1, 2, 3}];
|
|
Switch[type, 1,
|
|
randompos = RandomInteger[seq["SequenceLength"]];
|
|
seq = StringReplacePart[seq, RandomChoice[{"A", "T", "C", "G"}], {randompos, randompos}];
|
|
, 2,
|
|
randompos = RandomInteger[seq["SequenceLength"]];
|
|
seq = StringReplacePart[seq, "", {randompos, randompos}];
|
|
, 3,
|
|
randompos = RandomInteger[seq["SequenceLength"]];
|
|
seq = StringInsert[seq, RandomChoice[{"A", "T", "C", "G"}], randompos];
|
|
]
|
|
,
|
|
{10}
|
|
]
|
|
parts = StringPartition[seq["SequenceString"], UpTo[size]];
|
|
begins = Most[Accumulate[Prepend[StringLength /@ parts, 1]]];
|
|
ends = Rest[Accumulate[Prepend[StringLength /@ parts, 0]]];
|
|
StringRiffle[MapThread[ToString[#1] <> "-" <> ToString[#2] <> ": " <> #3 &, {begins, ends, parts}], "\n"]
|
|
Tally[Characters[seq["SequenceString"]]]
|