BOL

phyloXML (example) is an XML language designed to describe phylogenetic trees (or networks) and associated data. PhyloXML provides elements for commonly used features, such as taxonomic information, gene names and identifiers, branch...

https://github.com/crazyhottommy/getting-started-with-genomics-tools-and-resources/blob/master/Next%20Generation%20Sequencing%20spreadsheet%20-%20T.tsv #NGS #CheatSheet

#BlobToolKit can be used during #assembly to filter non-target DNA, helping researchers produce assemblies with high biological credibility. https://blobtoolkit.genomehubs.org/view #Genome #Blob #Con

Miniasm is a great long-read assembly tool: straight-forward, effective and very fast. However, it does not include a polishing step, so its assemblies have a high error rate – they are essentially made of stitched-together pieces of long...

Perl script to run in parallel https://perlmaven.com/speed-up-calculation-by-running-in-parallel #Perl #Parallel #Script #Core #Thread

RepeatModeler2 represents a valuable addition to the genome annotation toolkit that will enhance the identification and study of TEs in eukaryotic genome sequences. RepeatModeler2 is available as source code or a containerized package under an open...

Perl onliner to add string to each of the fasta header. $ perl -p -e 's/^(>.*)$/$1-origin/g' in.fasta > out.fasta #String #Header #Fasta

Simulated genomes with pre-defined and random genomic variants can be very useful for benchmarking genomic and bioinformatics analyses. Here we introduce simuG, a lightweight tool for simulating the full-spectrum of genomic variants (single...

The Environment for Tree Exploration (ETE) is a Python programming toolkit that assists in the recontruction, manipulation, analysis and visualization of phylogenetic trees (although clustering trees or any other tree-like data structure are also...

Multifasta 2 singlefasta cat toy.fasta| awk '{ if (substr($0, 1, 1)==">" {filename=(substr($0,2) ".fa")} print $0 > filename }' #Split #fasta #multifasta

One can evaluate the accuracy of a mapper or a SNP caller with wgsim_eval.pl that comes with the package. #Evaluate #Mapping #SNP

Deep_learning course https://www.slideshare.net/TessFerrandez/notes-from-coursera-deep-learning-courses-by-andrew-ng #AI #DeepLearn

Multiline to singleline covert. $ awk '/^>/ {printf("%s%s\t",(N>0?"\n":""),$0);N++;next;} {printf("%s",$0);} END {printf("\n");}' < input.fa #Multi #Single #Convert #awk

#Estimate the #insert #size with mapped reads $ head -10000 mappings.sam | awk '{ if ($9 > 0) { N+=1; S+=$9; S2+=$9*$9 }} END { M=S/N; print "n="N", mean="M", stdev="sqrt ((S2-M*M*N)/(N-1))}'

You need to add this file https://bioinformaticsonline.com/file/view/36952/getoptspl-file at requested location. Cheers

Cannot locate getopts.pl