Advertisement
jhangyu

get_rid_of_VEP's_problem.sh

Jan 31st, 2018
125
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Bash 0.89 KB | None | 0 0
  1. ### Clear the Parent info ###
  2. sed -i 's/;Parent=.*//g' mushroom_remove_parent.gff3
  3. ### Clear the scaffold size info incase the error in the output html report file ###
  4. sed -i 's/|size[0-9]*//g' ${ID}.vcf
  5. sed -i 's/|size[0-9]*//g' mushroom_remove_parent.gff3
  6. sed -i 's/|size[0-9]*//g' mushroom_clear_scaffold.fasta
  7. ### Add additional "0" to avoid the sort problem in html report file ###
  8. sed -r 's/scaffold[0-9]{1}\b/scaffold0&/g' ${ID}.vcf | sed -r 's/0scaffold/0/g' > ${ID}_add_zero
  9. sed -r 's/scaffold/scaffold\t/g' mushroom_remove_parent.gff3 | sed -r 's/\t[0-9]{1}\t/0&/g' > mushroom_add_zero.gff3
  10. sed -r 's/scaffold[0-9]{1}\b/scaffold0&/g' mushroom_clear_scaffold.fasta | sed -r 's/0scaffold/0/g' > mushroom_add_zero.fasta
  11.  
  12. ### Finally we can run VEP to output report ###
  13. vep -i ${ID}_add_zero -gff ${genomeDir}/${gffName}.gff.gz -fasta ${genomeDir}/${genomeFile}.gz -o ${ID}_variant_effect.txt
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement