Results for "Length"

Wire posts

  • perl -nle 'if (/^>/) { $sl = 0; print; next } next if ($sl >= 200); $_ = substr($_, 0, 200-$sl) if ($sl + length($_) > 200); $sl += length($_); print;' file.fa >new.fa #Perl #Extract #Fasta

    3633 days ago

  • Remove the blank space from a fasta file. perl -nlwe 'tr/ //d; print if length' fileName #Delete #Space #Fasta #Perl #Oneliner

    3633 days ago

  • Count the contigs length in index fasta file with samtools. $ awk '{s+=$2}END{print s}' scaffolds.fasta.fai #Fasta #Index #Length #Count #Size #Samtools

    2156 days ago

  • Sort FASTA sequence by length: cat f1.fasta| awk '{if($0 ~ /^>/){printf "%st", $0}else{print length($0)"t"$0}}'| sort -k2n,2| awk '{print $1"n"$3}' > sequence.fasta #Sort #Fasta #Trick #Onliner

    3200 days ago

  • #Extract #tabbed info from #BAM file with samtools: samtools view test.bam|awk '{print $1 "\t" $2 "\t" $3 "\t" $4 "\t" $4+length($10)-1}' > file.xls

    2885 days ago

  • #Extract #Chromosome #Number And Start Position #Reads: samtools view bamfile.bam|awk '{print $3 "\t" $4 "\t" $4+length($10)-1}' > newfile.tab

    2879 days ago

  • Fasta sequence length: awk '/^>/ {if (seqlen){print seqlen}; print ;seqlen=0;next; } { seqlen += length($0)}END{print seqlen}' file.fa #Length #Linux #OneLiner

    2671 days ago

  • Run length coding https://www.rosettacode.org/wiki/Run-length_encoding#Perl #Perl #String

    2527 days ago

  • #N50 = contig length so that half of the contigs are longer and 1/2 of contigs are shorter

    2191 days ago

  • perl -e '$/="\n>"; while (<>) { s/>//g; my ($id, $seq) = split /\n/; print ">$_" if length $seq; }' < all_p_ctg.fasta > all_p_ctg_corrected.fasta #Convert #Fasta #Zero

    2151 days ago