Public snipts » dgg32's snipts » filter the SSAHA2 result

• #to filter the SSAHA2 result like
  #===================================================
  #Matches For Query 20489 (913 bases): b201.1e11.f.b1
  #===================================================
  #Score     Q_Name             S_Name            Q_Start    Q_End  #S_Start    S_End Direction #Bases identity
  #ALIGNMENT::00 730   b201.1e11.f.b1 SC02       15      816    288648    #289465   F     802 97.51 913
  #ALIGNMENT::00 730   b201.1e11.f.b1 SC05       15      816     34214     #35031   F     802 97.51 913
  
  cat '/megx/home/shuang/result/ANME_vs_ANME1_SSAHA2/ANME_vs_ANME1_Endseq.txt' | egrep "ALIGNMENT" | sed 's/::/:/g' | sed 's/ \+/:/g' | cut -d: -f4,12
  

  bash

  1

  #to filter the SSAHA2 result like
  #===================================================
  #Matches For Query 20489 (913 bases): b201.1e11.f.b1
  #===================================================
  #Score     Q_Name             S_Name            Q_Start    Q_End  #S_Start    S_End Direction #Bases identity
  #ALIGNMENT::00 730   b201.1e11.f.b1 SC02       15      816    288648    #289465   F     802 97.51 913
  #ALIGNMENT::00 730   b201.1e11.f.b1 SC05       15      816     34214     #35031   F     802 97.51 913
  
  cat '/megx/home/shuang/result/ANME_vs_ANME1_SSAHA2/ANME_vs_ANME1_Endseq.txt' | egrep "ALIGNMENT" | sed 's/::/:/g' | sed 's/ \+/:/g' | cut -d: -f4,12
  

  copy | embed

0 Comments