BOL

Results for "awk"

Wire posts

• • BioStar@biostar

  

  MultiLine fasta to oneline fasta converter awk '/^>/ {printf("\n%s\n",$0);next; } { printf("%s",$0);} END {printf("\n");}' < file.fa #Multiline #Oneline

  2044 days ago

• • Neelam Jha@neelam

  

  Multifasta to singleline fasta. $ awk '/^>/ {printf("\n%s\n",$0);next; } { printf("%s",$0);} END {printf("\n");}' < file.fa #Multi #Single #Fasta #Convert

  2035 days ago

• • Jit@jit.aber

  

  Tab file to FASTA file $ awk '{print ">"$1"\n"$2}' allSeq_tm.txt > allSeq_tm.fa #Fasta #Tab #Convert #Tab2Fasta

  2014 days ago

• • Abhimanyu Singh@abhimanyu

  

  #Intersect your bed data with coverage $ bedtools genomecov -ibam aln.bam -bga | awk '$4==0' | bedtools intersect -a cds.bed -b - -f 1.0

  1990 days ago

• • Rahul Nayak@rahul

  

  Add number in fasta header $ awk '/^>/{$0=$0"_"(++i)}1' infile.fa #fasta #header #add #number #count

  1955 days ago

• • Rahul Nayak@rahul

  

  Split a multi-FASTA file into individual FASTA files: awk '/^>/{s=++d".fa"} {print > s}' multi.fa #split #multifasta

  1930 days ago

• • Jit@jit.aber

  

  Get the secondary hists reads count from SAM awk '{ print $2 }' out.sam| grep "2048" | wc -l #SAM #BAM #Secondary #Hits

  1829 days ago

• • LEGE@lege

  

  #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))}'

  1615 days ago

• • LEGE@lege

  

  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

  1615 days ago

• • Neelam Jha@neelam

  

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

  1615 days ago

• « Previous
• 1
• 2
• 3
• 4
• 5
• Next »