You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardExpand all lines: src/quick-start.md
+8-5Lines changed: 8 additions & 5 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -6,19 +6,18 @@
6
6
7
7
HiFi kinetics are required for predicting [m6A](glossary.md#m6a) with `fibertools`. **Check with your sequencing provider prior to sequencing** to ensure that the output file will have kinetics. Additionally, many of `fibertools` commands are compatible with CpG methylation which can be completed on instrument (if requested) or later with [Jasmine](https://github.com/PacificBiosciences/jasmine), e.g. `jasmine --keep-kinetics input.ccs.bam output.ccs.bam`. This command should be run prior to `fibertools` if CpG methylation information is desired as Jasmine will overwrite the m6A predictions in the MM and ML tags.
8
8
9
-
10
9
### Predict m6A and infer nucleosomes
11
10
12
11
#### Your PacBio data uses the SPQR chemistry or latter
13
12
14
13
As of the SPQR chemistry with PacBio their base modification caller `jasmine` can make m6A predictions on instrument in addition to 5mC.
15
-
This removes the need for `ft predict` and PacBio BAMs that have the kinetics information. However, after you must still run `ft add-nucleosomes` which is required for downstream analysis.
14
+
This removes the need for `ft predict` and PacBio BAMs that have the kinetics information. However, after you must still run `ft add-nucleosomes` which is required for downstream analysis.
16
15
17
16
```bash
18
17
ft add-nucleosomes -t 16 input.pacbio.bam output.fiberseq.bam
19
18
```
20
19
21
-
#### Your PacBio data predates the SPQR chemistry
20
+
#### Your PacBio data predates the SPQR chemistry
22
21
23
22
To create useable Fiber-seq data you must first call m6A base-mods on the PacBio CCS bam using `fibertools`. First [install fibertools](fibertools/install.md) and then process your bam file using the prediction command.
24
23
@@ -40,7 +39,8 @@ After this point, you will have a Fiber-seq BAM file that is compatible with all
40
39
41
40
### Validate your Fiber-seq BAM file
42
41
43
-
We have a quick validation tool which can test your BAM file for the desigered Fiber-seq features. At this point you should have m6A and nucleosome calls.
42
+
We have a quick validation tool which can test your BAM file for the desired Fiber-seq features.
43
+
At this point you should have m6A and nucleosome calls.
44
44
45
45
```bash
46
46
ft validate output.fiberseq.bam
@@ -55,7 +55,10 @@ Fraction with kinetics: 0.00%
55
55
56
56
### Fiber-seq peaks and UCSC browser tracks (FIRE)
57
57
58
-
Once you have a phased bam file, you can identify [Fiber-seq inferred regulatory elements (FIREs)](glossary.md#fires) to call Fiber-seq peaks and make a UCSC trackHub. Please see the FIRE page for more details [FIRE](fire/fire.md).
58
+
Once you have a phased bam file, you can identify [Fiber-seq inferred regulatory elements (FIREs)](glossary.md#fires) to call Fiber-seq peaks and make a UCSC trackHub.
59
+
60
+
You can find more details in on installing and running FIRE here:
61
+
[Running the FIRE pipeline](/fire/run.md).
59
62
60
63
# Fiber-seq starting with Oxford Nanopore Technologies (ONT)
0 commit comments