I have an illumina sequenced genome that is being mapped paired-end. The other 4 parts that I merge with it to create the fully mapped genome all look normal. When at the sampe step of BWA for this last part of the genome, I noticed it was taking an EXTREMELY long time. Here is the output concerning isize:
[infer_isize] (25, 50, 75) percentile: (30225, 52418, 70153)
[infer_isize] low and high boundaries: 45 and 150009 for estimating avg and std
[infer_isize] inferred external isize from 129 pairs: 49974.163 +/- 26814.967
[infer_isize] skewness: -0.056; kurtosis: -0.892; ap_prior: 1.00e-05
[infer_isize] inferred maximum insert size: 199065 (5.56 sigma)
The inferred size is definitely not what it should be (~250). Has anyone had this kind of problem before? Is the sequencing data for this sample usable?
[infer_isize] (25, 50, 75) percentile: (30225, 52418, 70153)
[infer_isize] low and high boundaries: 45 and 150009 for estimating avg and std
[infer_isize] inferred external isize from 129 pairs: 49974.163 +/- 26814.967
[infer_isize] skewness: -0.056; kurtosis: -0.892; ap_prior: 1.00e-05
[infer_isize] inferred maximum insert size: 199065 (5.56 sigma)
The inferred size is definitely not what it should be (~250). Has anyone had this kind of problem before? Is the sequencing data for this sample usable?
Comment