mahulchak / svmu

A program to call variants from genome alignment
GNU General Public License v3.0
75 stars 18 forks source link

The coordinates of CNV in sv.prefix.txt #27

Open LiuShang-777 opened 1 year ago

LiuShang-777 commented 1 year ago

Hello, I used SVMU for SV detection based on mummer result and found that in the sv.prefix.txt file, the length of CNV in both reference and query genomes were the same like below: REF_CHROM REF_START REF_END SV_TYPE Q_CHROM Q_START Q_END ID LEN COV_REF COV_Q A01 356308 356314 CNV-R A01 245745 245751 0000000023 6 1 2

356314-356308=245751-245745=6, but the covR and covQ are different. Is this means that segments in query have 2 copies in reference genome? I think that (ref_end-ref_start)/(query_end-query_start) should be equal to 2 in this condition.

mahulchak commented 1 year ago

Based on the information provided, your initial interpretation seems to be correct (reference has 2 copies of that sequence, but query has one).

On Tue, May 30, 2023, 07:29 LiuShang-777 @.***> wrote:

Hello, I used SVMU for SV detection based on mummer result and found that in the sv.prefix.txt file, the length of CNV in both reference and query genomes were the same like below: REF_CHROM REF_START REF_END SV_TYPE Q_CHROM Q_START Q_END ID LEN COV_REF COV_Q A01 356308 356314 CNV-R A01 245745 245751 0000000023 6 1 2

356314-356308=245751-245745=6, but the covR and covQ are different. Is this means that segments in query have 2 copies in reference genome? I think that (ref_end-ref_start)/(query_end-query_start) should be equal to 2 in this condition.

— Reply to this email directly, view it on GitHub https://urldefense.com/v3/__https://github.com/mahulchak/svmu/issues/27__;!!CzAuKJ42GuquVTTmVmPViYEvSg!NFvaW4mwDSTna2JZJBTQlrScyxunqFsiNhuwMWWchILoDgkC_WqJfRHg1DVbtxkeLLBd-Rus-jRol2WZ1btKHJMj$, or unsubscribe https://urldefense.com/v3/__https://github.com/notifications/unsubscribe-auth/ABZQH2BJFQMX6436AZQ7IBDXIXR2VANCNFSM6AAAAAAYT6AFWI__;!!CzAuKJ42GuquVTTmVmPViYEvSg!NFvaW4mwDSTna2JZJBTQlrScyxunqFsiNhuwMWWchILoDgkC_WqJfRHg1DVbtxkeLLBd-Rus-jRol2WZ1Zdr6ZMI$ . You are receiving this because you are subscribed to this thread.Message ID: @.***>

LiuShang-777 commented 1 year ago

Based on the information provided, your initial interpretation seems to be correct (reference has 2 copies of that sequence, but query has one). On Tue, May 30, 2023, 07:29 LiuShang-777 @.> wrote: Hello, I used SVMU for SV detection based on mummer result and found that in the sv.prefix.txt file, the length of CNV in both reference and query genomes were the same like below: REF_CHROM REF_START REF_END SV_TYPE Q_CHROM Q_START Q_END ID LEN COV_REF COV_Q A01 356308 356314 CNV-R A01 245745 245751 0000000023 6 1 2 356314-356308=245751-245745=6, but the covR and covQ are different. Is this means that segments in query have 2 copies in reference genome? I think that (ref_end-ref_start)/(query_end-query_start) should be equal to 2 in this condition. — Reply to this email directly, view it on GitHub <https://urldefense.com/v3/https://github.com/mahulchak/svmu/issues/27;!!CzAuKJ42GuquVTTmVmPViYEvSg!NFvaW4mwDSTna2JZJBTQlrScyxunqFsiNhuwMWWchILoDgkC_WqJfRHg1DVbtxkeLLBd-Rus-jRol2WZ1btKHJMj$>, or unsubscribe <https://urldefense.com/v3/https://github.com/notifications/unsubscribe-auth/ABZQH2BJFQMX6436AZQ7IBDXIXR2VANCNFSM6AAAAAAYT6AFWI;!!CzAuKJ42GuquVTTmVmPViYEvSg!NFvaW4mwDSTna2JZJBTQlrScyxunqFsiNhuwMWWchILoDgkC_WqJfRHg1DVbtxkeLLBd-Rus-jRol2WZ1Zdr6ZMI$> . You are receiving this because you are subscribed to this thread.Message ID: @.>

Ok! Thank you for your kind reply!