How to call CNV(copy number variant) in Pacbio ?
0
1
Entering edit mode
6.1 years ago
qincai158 ▴ 10

Hi, I have some data on rice.BWA software is used for sequence alignment, Be based for alignment result,I get many structural variants called by sniffles ,Now I want to get copy number variant. I don't know what to do.Can you recommend some software that can detect CNV(copy number variant ) ?

 I'll appreciate your help!
alignment sequence gene • 2.3k views
ADD COMMENT
0
Entering edit mode

CPV is type of SV specially deletion or duplication which both identified by Sniffles . Their is exact thread for what you asked (specially that you used sniffles ):

How to group CNV

ADD REPLY
0
Entering edit mode

Thank for your reply ! How to judge copy number variants from deletion or duplication or vcf file?

ADD REPLY
0
Entering edit mode

You can count both as CNV

ADD REPLY
0
Entering edit mode

Thank for your reply! But Sniffles software get duplication which contains Tandem duplication and Dispersed duplication, CNV refers to Tandem duplication. How to abstract tandem duplication from duplication? Do you agree my idea ?

ADD REPLY
0
Entering edit mode

CNV refers to Tandem duplication

I don't think I agree here. An element which is duplicated will have a higher copy number, regardless whether the duplication is in tandem or dispersed.

ADD REPLY
0
Entering edit mode

Sniffles detects structural variants, which include copy number variants. So you already have copy number variants.

ADD REPLY
0
Entering edit mode

Thank for your reply ! I read an artical named "Detection and assessment of copy number variation using PacBio long-read and Illumina sequencing in New Zealand dairy cattle". The author grouped DUP and DEL called by Sniffles as CNV . Here is the original words in the artical:

Sniffles software identified a total of 38,709 putative SV, of which 19,797 were CNV (deletions, n = 18,577; duplications, n = 1,220) in a PacBio sequence generated from Esteem’s DNA. Of the 3,532 CNV (deletions, n = 3,055; duplications, n = 477) that remained after filtering, 899 (deletions, n = 740; duplications, n = 159) and 126 (deletions, n = 76; duplications, n = 50) overlapped by at least 1 bp with Ensembl (build 84; http:// mar2016.archive.ensembl.org/index.html) annotated genes and exons, respectively (Table 1).

I am not sure whether he is right.

ADD REPLY
0
Entering edit mode

A copy number variant is an element of which the amount varies in the genome. You have for example more or less copies of gene x. A deletion or a duplication is a copy number variant: you have more or less of a feature.

ADD REPLY

Login before adding your answer.

Traffic: 1289 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6