Para contar o número de linhas que conte nucleotideos

perl -lane 'print $_ if ($.%4==2);' [nome do arquivo.fastq] | wc -l

Comentários

Postar um comentário

Postagens mais visitadas deste blog

Retirar SVG de página da web - exemplo Google Chart

Entrar pelo google ferramentas do desenvolvedor  When the SVG is integrated as   <svg ...></svg>   markup directly into the HTML page. Right click on the SVG to inspect it in developer tools Find the root of the  <svg>  element and right click to  "Copy element" Go to  https://jakearchibald.github.io/svgomg/  and  "Paste markup" Download your optimized SVG file and enjoy
Leia mais

Comandos comuns do PAUP - block

retirado de : http://www.peter.unmack.net/molecular/programs/paup.command.blocks.html Common PAUP analysis commands By Peter Unmack I usually paste the blocks below to the end of my nexus files, open the command line version (which runs faster than the GUI) and execute the file (I usually copy the datafile to my paup directory, or copy paup to where my datafile is).   Stuff highlighted in yellow is what usually gets changed by the user depending upon their dataset and needs. Parsimony Analysis begin paup; set autoclose=yes; set criterion=parsimony; set root=outgroup; set storebrlens=yes; set increase=auto; outgroup fish1 fish2 etc ; hsearch addseq=random nreps=1000 swap=tbr hold=1; savetrees file= datasetname .cb.pa.tree.nex format=altnex brlens=yes; pscores /tl ci ri rc; end; Parsimony Bootstrap Analysis begin paup; set autoclose=yes; set criterion=parsimony; set root=outgroup; set storebrlens=yes; set increase=auto; outgr...
Leia mais