Perl script to insert sequence in contig !!
# sub signature: #insertSEQintoCONTIGatLOC( SEQ , CONTIG , LOC ) ; sub insertSEQintoCONTIGatLOC{ my ( $SEQ , $CON...ATAGATAG"; print "s = $s \n"; print 'insert = '.($s = insertSEQintoCO...GATAGATGATAGTAGATAGATAGATAGATAGATAG # insert = 'CCCC ' # (now, s...2723 days ago
Perl script to insert the DNA string in genome
...q; my $file = $ARGV[0]; # input fasta file (genome file) my $out = $ARGV[1]; # output fasta file my $chr="test"; #insertion chromosome my $pos=10; # position of the insertion my $seqI = "AAAA"; #seque...2674 days ago
Genetic Algorithms demonstration with word DNA in Perl
...ot defined $individual); push @parent_population, $individual; # insert the individual in the parent...ft @_; # for each individual foreach my $id (1 .. $pop_size) { # insert an anonymous hash reference i...2387 days ago
Insert the sequence at desire location in multi-fasta file with Perl
#!/usr/bin/perl use warnings; use strict; use Bio::SeqIO; use Bio::Seq; use File::Copy; #ARGV[0] should be in following format --- Keep the coordinate sorted...2329 days ago
2090 days ago
picard tools command to get some insert statistics
#picard tools to get some insert statistics to see whether our reads seem to be in the correct place #.../2.0.1 java -Xmx16g -XX:PermSize=8g -jar $PICARD_HOME/picard.jar CollectInsertSizeMetrics MINIMUM_PCT=0 HIST...1352 days ago
Commandline for paired end reads simulation with BBMap !
...paired interleaved reads=100k length=250 mininsert=400 maxinsert=600 gaussian java -ea -Xmx16...paired interleaved reads=100k length=250 mininsert=400 maxinsert=600 gaussian E...paired=true read length=250 reads=100000 insert size=400-600 Wrote reads_BBM...1005 days ago
Tadpole is 250x faster than SPADes assembler !
...d optionally perform error correction. correct: Error correct only. insert: Measure insert sizes. d...991 days ago