Then, as you can see, we have now as many nuc
files in directory S4/NMalign/ than the number of studied chromosomes. These .nuc
files are text files organized as a table which columns are seperated by a single whitespace. From left to right, the meaning of each column is:
- the reference sequence,
- the reference chunk index,
- the ignore status of the reference nucleosome (y=yes/n=no),
- the outlier status of the reference nucleosome (0=ok,1=outlier),
- the starting position of the nucleosome on the reference,
- the ending position of the nucleosome on the reference,
- the status of the nucleosome on the reference (1=fuzzy, 2=well-localized),
- the shift in bp between the midpoints of the reference and the corresponding query nucleosomes,
- the fraction of the insertion/deletion spanning the nucleosome w.r.t the total length of the reference nucleosome,
- the average probability that the probes covering the reference nucleosome fall into a nucleosome (either fuzzy or well-localized),
- the query chunk index,
- the ignore status of the query nucleosome (y=yes/n=no),
- the outlier status of the query nucleosome (0=ok,1=outlier),
- the starting position of the query nucleosome on the reference,
- the ending position of the query nucleosome on the reference,
- the status of the nucleosome on the query (1=fuzzy, 2=well-localized),
- the shift in bp between the midpoints of the query and the corresponding reference nucleosomes,
- the fraction of the insertion/deletion spanning the nucleosome w.r.t the total length of the query nucleosome,
- the average probability that the probes covering the query nucleosome fall into a nucleosome (either fuzzy or well-localized).
- if aligned, the p-value of the nucleosome (see next section),
- if aligned, the effect size corresponding to the p-value (see next section),
- if aligned, the status of the aligned nucleosome:
- -1 = insertion/deletion,
- 0 = fuzzy/fuzzy,
- 1 = fuzzy/well-localized or well-localized/fuzzy,
- 2 = well-localized/well-localized,
- the probability that the region corresponds in fact to a linker / linker configuration,
- the probability that the region corresponds in fact to a nucleosome (reference) / linker (query) configuration (insertion),
- the probability that the region corresponds in fact to a linker (query) / nucleosome (reference) configuration (deletion),
- the probability that the region corresponds in fact to a nucleosome / nucleosome configuration.
Now, we can use the function NMnxcplot
of the nucleominer R package to visualize the content
of these .nuc
files. An exemple is depicted in Figure 8.
Jean-Baptiste Veyrieras
2010-05-28