Skip to content

Commit 6599dc6

Browse files
committed
fixed formatting and links for deduplication step
1 parent 189c835 commit 6599dc6

File tree

1 file changed

+8
-2
lines changed

1 file changed

+8
-2
lines changed

Methyl-Seq/Pipeline_GL-DPPD-7113_Versions/GL-DPPD-7113.md

Lines changed: 8 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -43,6 +43,8 @@ Jonathan Galazka (GeneLab Project Scientist)
4343
- [4c. Sort Alignment Files](#4c-sort-alignment-files)
4444
- [5. Alignment QC](#5-alignment-qc)
4545
- [6. Deduplicate (skip if data are RRBS)](#6-deduplicate-skip-if-data-are-rrbs)
46+
- [6a. Deduplicate](#6a-deduplicate)
47+
- [6b. Sort Deduplicated Alignment Files](#6b-sort-deduplicated-alignment-files)
4648
- [7. Extract methylation calls](#7-extract-methylation-calls)
4749
- [8. Generate individual sample report](#8-generate-individual-sample-report)
4850
- [9. Generate combined summary report](#9-generate-combined-summary-report)
@@ -559,10 +561,14 @@ qualimap bamqc -bam sample-1_bismark_bt2_sorted.bam \
559561

560562
---
561563

562-
## 6a. Deduplicate (skip if data are RRBS)
564+
## 6. Deduplicate (skip if data are RRBS)
563565
> **NOTE**
564566
> This step should **not** be done if the data are RRBS (reduced representation bisulfite sequencing; see e.g., [bismark documentation](https://felixkrueger.github.io/Bismark/bismark/deduplication/)).
565567
568+
<br>
569+
570+
### 6a. Deduplicate
571+
566572
```bash
567573
deduplicate_bismark sample-1_bismark_bt2.bam
568574
```
@@ -664,7 +670,7 @@ bismark_methylation_extractor --parallel 4 \
664670

665671
**Input data:**
666672

667-
* sample-1_bismark_bt2_sorted*.bam (bismark bowtie2 alignment bam file sorted by chromosomal coordinates, output from [Step 4c](#4c-sort-alignment-files) above if data are RRBS, or deduplicated bam file from [step 6](#6b-sort-deduplicated-alignment-files) if data are not RRBS and the bam file was deduplicated (e.g., sample-1_bismark_bt2_sorted.deduplicated.bam from above))
673+
* sample-1_bismark_bt2_sorted*.bam (bismark bowtie2 alignment bam file sorted by chromosomal coordinates, output from [Step 4c](#4c-sort-alignment-files) above if data are RRBS, or deduplicated bam file from [step 6b](#6b-sort-deduplicated-alignment-files) if data are not RRBS and the bam file was deduplicated (e.g., sample-1_bismark_bt2_sorted.deduplicated.bam from above))
668674
* a directory holding the reference genome in fasta format (this pipeline version uses the Ensembl fasta file indicated in the `fasta` column of the [GL-DPPD-7110_annotations.csv](../../GeneLab_Reference_Annotations/Pipeline_GL-DPPD-7110_Versions/GL-DPPD-7110/GL-DPPD-7110_annotations.csv) GeneLab Annotations file))
669675

670676

0 commit comments

Comments
 (0)