Advertisement
Guest User

Untitled

a guest
Jun 17th, 2019
77
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.65 KB | None | 0 0
  1. rs36810213_HET rs2438689 rs70927523570_HET rs54666437 ...
  2. 1 0 2 0
  3. 0 1 0 1
  4. 2 0 1 1
  5. ... ... ... ...
  6.  
  7. rs2438689 rs54666437 ...
  8. 0 0
  9. 1 1
  10. 0 1
  11. ... ...
  12.  
  13. #!/bin/bash
  14.  
  15. path="/data/folder"
  16.  
  17. awk -v OFS='t' '
  18.  
  19. NR==1{
  20. for (i=1;i<=NF;i++)
  21. if ($i=="_HET") {
  22. n=i-1
  23. m=NF-(i==NF)
  24. }
  25. }
  26.  
  27. {
  28. for(i=1;i<=NF;i+=1+(i==n))
  29. printf "%s%s",$i,i==m?ORS:OFS
  30. }
  31.  
  32. ' $path/input.txt >> $path/output.txt
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement