code:
  #!/usr/bin/bash
len=$(wc -l s7.fa  | awk '{print($1)}' )
(( len++ ))
for ((i=1; i<$len; i+=1))
do
        echo '>SEQUENCE_'${i}'_length_35' >> s7_fasted.fa
        head -n $i forSOAP_s7_final.fa | tail -n 1 >> s7_fasted.fa

done

Quote:

---

for ((i=1; i<$len; i+=1))
do
     echo '>SEQUENCE_'${i}'_length_35' >> s7_fasted.fa
     head -n $i forSOAP_s7_final.fa | tail -n 1 >> s7_fasted.fa

done

---

code:

#!/usr/bin/tclsh

proc Numerate {fileName} {
  set h [open $fileName r]
  set lineNumber 0
  while {[gets $h line] >= 0} {
    incr lineNumber
    puts ">SEQUENCE_${lineNumber}_length_35"
    puts $line
  }
  close $h
}

if {$argc >= 1} {
  Numerate [lindex $argv 0]
} else {
  puts "Usage: tclsh script.tcl filename"
}

code:
perl -pe 'print ">SEQUENCE_", ++$i, "_length_35\n"' < s7.fa

code:
 % i=1; cat forSOAP_s7_final.fa | while read line; do echo '>SEQUENCE_'${i}'_length_35'; echo $line; i=$((i+1)); done > s7_fasted.fa

code:

sed -e "=" | sed -e "s/^/>SEQUENCE_/;s/$/_length_35/;N"

code:
awk '{print "SEQUENCE_" NR "_length_35"} 1'  

code:

import System (getArgs)

parseFile :: String -> IO()
parseFile filename = do
  content <- readFile filename
  mapM_ putStrLn $ map (\(x, y) -> concat [">SEQUENCE_", show x, "_length_35\n", y]) $ zip [1..] (lines content)

main = do
  args <- getArgs
  case args of
    [] -> error "Usage: program filename"
    (filename:xs) -> parseFile filename 

Quote:

---

пейтон это

---

Quote:

---

Не, не получится

---

Quote:

---

зы: возможно, это уже нечестно сравнивать со скомпилированным скриптом =\

---

Quote:

---

Это да, чтобы по-честному, надо через runghc запускать.

---

Quote:

---

зы: возможно, это уже нечестно сравнивать со скомпилированным скриптом =\

---