Hi

If the masked width is the total number of nucleotide positions that are masked, then why the maskedwidth in 'all masks together:" is not equal to the sum of maskedwidth of all 3 masks?
chr2L=Dmelanogaster$chr2L
chr2L
23011544-letter "MaskedDNAString" instance (# for masking) seq: CGACAATGCACGACAGAGGAAGCAGAACAGATATTTAGATTGCCTCTCATTTTCTCTCCCATATTATAGGGAGAAATAT...CAATCAAACTGTGTTCGAAAAAGAGAAAACTAACATTTTTTTGGCATATTTGCAAATTTTGATGAACCCCCCTTTCAAA
masks:
maskedwidth maskedratio active names 1 200 8.691290e-06 FALSE assembly gaps 2 1966561 8.545976e-02 FALSE RepeatMasker 3 61603 2.677048e-03 FALSE Tandem Repeats Finder [period<=12]
all masks together:
  maskedwidth maskedratio
      1988181  0.08639929
all active masks together:
  maskedwidth maskedratio
            0           0


sum(200, 1966561, 61603)
[1] 2028364




Regards,
Joseph

Joseph M. Dhahbi, PhD
Childrens Hospital Oakland Research Institute
5700 Martin Luther King Jr. Way
Oakland, CA 94609
USA
Ph.(510)428-3885 EXT.5743
Cell.(702)335-0795
Fax (510)450-7910
[EMAIL PROTECTED]
The email message (and any attachments) is for the sole...{{dropped:3}}

_______________________________________________
Bioc-sig-sequencing mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing

Reply via email to