-
Notifications
You must be signed in to change notification settings - Fork 37
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Add subsample mt #508
Add subsample mt #508
Conversation
|
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Small comments but looks good otherwise
docs/output.md
Outdated
@@ -115,6 +116,18 @@ The pipeline is built using [Nextflow](https://www.nextflow.io/) and processes d | |||
- `*.metrics`: Text file containing the dedup metrics. | |||
</details> | |||
|
|||
#### Subsample mitochondrial alignments | |||
|
|||
[Samtools view](https://www.htslib.org/doc/samtools-view.html) is used by the pipeline to subsample mitochondrial alignments to a user specified coverage. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
[Samtools view](https://www.htslib.org/doc/samtools-view.html) is used by the pipeline to subsample mitochondrial alignments to a user specified coverage. | |
[Samtools view](https://www.htslib.org/doc/samtools-view.html) is used by the pipeline to subsample mitochondrial alignments to a user specified coverage. The file is mainly intended to be used for visualization of MT alignments in IGV. The non-subsampled bam file is used for variant calling and other downstream analysis steps. |
subworkflows/local/subsample_mt.nf
Outdated
) | ||
.csv | ||
.join(ch_mt_bam_bai, failOnMismatch:true) | ||
.map{meta, seedfrac, bam, bai -> | ||
return [meta + [seedfrac: file(seedfrac).text.readLines()[0]], bam, bai] | ||
} | ||
.set { ch_subsample_in } |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Could you fix the indentation of this part. Think it needs to be left shifted by one level
PR checklist
Addresses #497
nf-core lint
).nextflow run . -profile test,docker --outdir <OUTDIR>
).nextflow run . -profile test_one_sample,docker --outdir <OUTDIR>
).docs/usage.md
is updated.docs/output.md
is updated.CHANGELOG.md
is updated.README.md
is updated (including new tool citations and authors/contributors).