Advertisement
jhangyu

Circos karyotype dataformat

Jun 21st, 2018
102
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 2.50 KB | None | 0 0
  1. Go to http://genome-test.soe.ucsc.edu/cgi-bin/hgTables to get Tilapia genome table
  2. Clade>Vertebrate    Genome>Nile Tilapia     Assembly>Nov. 2016(ASM185804v2/oreNil3)
  3. group>Mapping and Sequencing    track>Assemple
  4. And get the 'output' file like this:
  5.  
  6. #bin    chrom   chromStart  chromEnd    ix  type    frag    fragStart   fragEnd strand
  7. 585 chrM    0   16627   1   O   NC_013663.1 0   16627   +
  8. 1   chrLG1  0   38372991    1   O   MKQE01000001.1  0   38372991    +
  9. 1   chrLG2  0   35256741    1   O   MKQE01000012.1  0   35256741    +
  10. 1   chrLG4  0   38038224    1   O   MKQE01000018.1  0   38038224    +
  11. 1   chrLG5  0   34628617    1   O   MKQE01000019.1  0   34628617    +
  12. 1   chrLG6  0   44571662    1   O   MKQE01000020.1  0   44571662    +
  13. 1   chrLG7  0   62059223    1   O   MKQE01000021.1  0   62059223    +
  14. 1   chrLG8  0   30802437    1   O   MKQE01000022.1  0   30802437    +
  15. 1   chrLG9  0   27519051    1   O   MKQE01000023.1  0   27519051    +
  16. 1   chrLG10 0   32426571    1   O   MKQE01000002.1  0   32426571    +
  17. 1   chrLG11 0   36466354    1   O   MKQE01000003.1  0   36466354    +
  18. 1   chrLG12 0   41232431    1   O   MKQE01000004.1  0   41232431    +
  19. 1   chrLG13 0   32337344    1   O   MKQE01000005.1  0   32337344    +
  20. 1   chrLG14 0   39264731    1   O   MKQE01000006.1  0   39264731    +
  21. 1   chrLG15 0   36154882    1   O   MKQE01000007.1  0   36154882    +
  22. 1   chrLG16 0   43860769    1   O   MKQE01000008.1  0   43860769    +
  23. 1   chrLG17 0   40919683    1   O   MKQE01000009.1  0   40919683    +
  24. 1   chrLG18 0   37007722    1   O   MKQE01000010.1  0   37007722    +
  25. 1   chrLG19 0   31245232    1   O   MKQE01000011.1  0   31245232    +
  26. 1   chrLG20 0   36767035    1   O   MKQE01000013.1  0   36767035    +
  27. 1   chrLG22 0   37011614    1   O   MKQE01000014.1  0   37011614    +
  28. 1   chrLG23 0   44097196    1   O   MKQE01000015.1  0   44097196    +
  29. 1   chrLG3a 0   14041792    1   O   MKQE01000016.1  0   14041792    +
  30. 1   chrLG3b 0   54508961    1   O   MKQE01000017.1  0   54508961    +
  31.  
  32. Then you have to turn this genome table file into this format:(or just use my format)
  33.  
  34. chr - M M   0   16627       chrM
  35. chr - LG1   1   0   38372991        chr1
  36. chr - LG2   2   0   35256741        chr2
  37. chr - LG3a  3a  0   14041792        chr3
  38. chr - LG3b  3b  0   54508961        lum70chr3
  39. chr - LG4   4   0   38038224        chr4
  40. chr - LG5   5   0   34628617        chr5
  41. chr - LG6   6   0   44571662        chr6
  42. chr - LG7   7   0   62059223        chr7
  43. chr - LG8   8   0   30802437        chr8
  44. chr - LG9   9   0   27519051        chr9
  45. chr - LG10  10  0   32426571        chr10
  46. chr - LG11  11  0   36466354        chr11
  47. chr - LG12  12  0   41232431        chr12
  48. chr - LG13  13  0   32337344        chr13
  49. chr - LG14  14  0   39264731        chr14
  50. chr - LG15  15  0   36154882        chr15
  51. chr - LG16  16  0   43860769        chr16
  52. chr - LG17  17  0   40919683        chr17
  53. chr - LG18  18  0   37007722        chr18
  54. chr - LG19  19  0   31245232        chr19
  55. chr - LG20  20  0   36767035        chr20
  56. chr - LG22  22  0   37011614        chr22
  57. chr - LG23  23  0   44097196        chr23
  58.  
  59. Save these text as karyotype.nile.tilapia.oreNil3.txt
  60. put it into /$anaconda_install_path/data/karyotype/
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement