Create random 1000 INDEL in genome !
(base) ➜ dupStudy git:(master) ✗ perl ../simuG.pl -refseq simuSNP.simseq.genome.fa -indel_count 1000 -prefix sim...on: simuINDEL.refseq2simseq.map.txt Generating reference-based vcf file for genomic...1212 days ago
Create random 1000 CNVs in genome !
(base) ➜ dupStudy git:(master) ✗ perl ../simuG.pl -refseq simuINDEL.simseq.genome.fa -cnv_count 100 -prefix simuCN...tion: simuCNV.refseq2simseq.map.txt Generating reference-based vcf file for genomic...1212 days ago
Create random 5 inversions in genome !
(base) ➜ dupStudy git:(master) ✗ perl ../simuG.pl -refseq simuCNV.simseq.genome.fa -inversion_count 5 -prefix si...tion: simuINV.refseq2simseq.map.txt Generating reference-based vcf file for genomic...1212 days ago
Create random 2 translocations in genome !
(base) ➜ dupStudy git:(master) ✗ perl ../simuG.pl -refseq simuINV.simseq.genome.fa -translocation_count 2 -prefi...on: simuTRANS.refseq2simseq.map.txt Generating reference-based vcf file for genomic...1212 days ago
1204 days ago
Parse the NCBI taxonomy database with Perl !
use Bio::DB::Taxonomy; use warnings; my $sps="human"; # Get one from a NCBI taxonomy database my $dbh =...ancestor($nam); print "$ancestor_taxon->id ****"; #my $tree = $dbh->get_tree('human'...1173 days ago
Commandline for paired end reads simulation with BBMap !
(JitMetaENV) ➜ mixedSample git:(main) ✗ /home/urbe/Tool...mininsert=400, maxinsert=600, gaussian] Writing reference. Executing dna.FastaToChromA..., max=0, len=(0-0) nRate =0.0, max=0, len=(0-0) genome=1 PERFECT_READ_RATIO=0.0 AD...991 days ago
Tadpole is 250x faster than SPADes assembler !
lege@jit-Lenovo-ideapad-320-15ISK:~/Downloads/...mmended that mce=1 for low-depth metagenomes. mincountretain=0 (mincr...may yield more accurate genome size estimation. contigpasse...but they are very slow for large metagenomes. Overlap parameters (for...977 days ago
962 days ago
929 days ago