WeiWei060512 / NUMTs-detection

Detecting NUMTs from WGS
11 stars 2 forks source link

Breakpoints output and circos plot #7

Open notthathime opened 1 year ago

notthathime commented 1 year ago

Thank you very much for your great work. I have a question about the breakpoint output and how to use it for circos plot. For example, I got seperate tsv results from [searchBreakpoint_fromblatoutputs.py: here is the example of one tsv file named 10T_chr4.86138150.86139914.Breakpoints.tsv

nu_Tend_Bleft nuleft chr4 86138985 + 9 -1 10T nu_Tstart_Bright nuright chr4 -1 + 1 86138929 10T nu_Tstart_Bright nuright chr4 -1 + 1 86138933 10T nu_Tstart_Bright nuright chr4 -1 + 1 86138935 10T nu_Tstart_Bright nuright chr4 -1 + 1 86138936 10T nu_Tstart_Bright nuright chr4 -1 + 1 86138938 10T nu_Tstart_Bright nuright chr4 -1 + 1 86138942 10T nu_Tstart_Bright nuright chr4 -1 + 1 86138943 10T nu_Tstart_Bright nuright chr4 -1 + 1 86138944 10T nu_Tstart_Bright nuright chr4 -1 + 1 86138945 10T nu_Tstart_Bright nuright chr4 -1 + 1 86138950 10T nu_Tstart_Bright nuright chr4 -1 + 1 86138955 10T nu_Tstart_Bright nuright chr4 -1 + 2 86138999 10T mt_Tend mtLeft chrM 6217 - 8 -1 10T mt_Tstart mtRight chrM -1 - 2 6147 10T

My question is how do you sort the breakpoints and make the matched pair as the output is just in the order of NU first and then MT And for Circosplot, I assumed that you need to prapare input as the paired breakpoints. May I ask how did you generate that?

Thank you very much!