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: README.md
+5-4Lines changed: 5 additions & 4 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -12,7 +12,7 @@ A Nextflow DSL2 pipeline to perform allele-specific expression analysis from RNA
12
12
## Features
13
13
14
14
- Genotype-aware RNA-seq alignment with **STAR + WASP correction**.
15
-
-Allele-specific read counting (REF/ALT) from RNA-seq BAMs.
15
+
-Optional allele-specific read counting (REF/ALT) from RNA-seq BAMs, performed only if a regions VCF is provided and exists.
16
16
- Reports total depth and variant allele fraction.
17
17
18
18
## Requirements
@@ -68,9 +68,9 @@ reference/
68
68
└── gencode.v38.annotation.gtf
69
69
```
70
70
71
-
### 4. Define SNPs of interest (VCF)
71
+
### 4. Define SNPs of interest (VCF) (optional)
72
72
73
-
The pipeline expects a VCF so REF and ALT alleles are explicit.
73
+
If you want to perform allele-specific counting, provide a VCF file with the SNPs of interest. If this file is not provided or does not exist, the allele counting step will be skipped.
74
74
75
75
`query_sites.vcf`
76
76
@@ -123,7 +123,8 @@ The pipeline consists of the following main steps:
123
123
- Uses per-sample phased variants via `--varVCFfile`.
124
124
- Output: `sampleID.Aligned.sortedByCoord.out.bam` (coordinate sorted BAM file with WASP tags)
0 commit comments