PennCNV Software ; Not reporting CNVs in some samples
0
0
Entering edit mode
3.4 years ago
zkafi1370 • 0

Hello

I am working with PennCNV software to analyze 300 samples (chickens) and I used below command for this analysis:

./detect_cnv.pl -test -hmm lib/hh550.hmm -pfb gallus.pfb sample1.txt sample2.txt sample3.txt -log sampleall.adjusted.log -out sampleall.adjusted.rawcnv -gcmodel gallus.gcmodel

but the software does not report any CNVs for a large number of samples (125 samples) without any errors or warnings.

For example, one of these samples is sample1.txt , which after running, the software gives the following notices but the output file is empty :

NOTICE: All program notification/warning messages that appear in STDERR will be also written to log file sampleall.adjusted.log
NOTICE: Reading marker coordinates and population frequency of B allele (PFB) from out.pfb ... Done with 51348 records (3980 records in chr W,0,Z,LGE64,LGE22C19W28_E50C23 were discarded)
NOTICE: Reading LRR and BAF values for from sample1.txt ... Done with 51348 records in 29 chromosomes (3981 records are discarded due to lack of PFB information for the markers)
NOTICE: Adjusting LRR by GC model: WF changes from -0.0105 to 0.0111, GCWF changes from -0.0024 to 0.0001</code>
NOTICE: Data from chromosome 32 will not be used in analysis
NOTICE: Median-adjusting LRR values for all autosome markers from sample1.txt by 0.0039
NOTICE: Median-adjusting BAF values for all autosome markers from sample1.txt by -0.0303
NOTICE: quality summary for sample1.txt: LRR_mean=-0.0063 LRR_median=0.0000 LRR_SD=0.1207 BAF_mean=0.4980 BAF_median=0.5000 BAF_SD=0.0399 BAF_DRIFT=0.000000 WF=0.0111 GCWF=0.0001

I hope you help me

Thank you

SNP snp CNV • 666 views
ADD COMMENT
0
Entering edit mode

ah something's gone afowl. hh550.hg18.pfb looks like a human reference to me. shouldn't you be using something that is more like galGal6?

ADD REPLY
0
Entering edit mode

I had made these changes before but there was no change in the result ...

ADD REPLY

Login before adding your answer.

Traffic: 1434 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