Commit 79a51d31 authored by Matteo Barcella's avatar Matteo Barcella
Browse files

Update README.md

parent 2815127d
...@@ -20,14 +20,11 @@ Fastq files were aligned with [BWA aligner](https://github.com/lh3/bwa) (v0.7.17 ...@@ -20,14 +20,11 @@ Fastq files were aligned with [BWA aligner](https://github.com/lh3/bwa) (v0.7.17
- gatk BaseRecalibrator --input SampleID.dedup_reads_mouse_recal.bam --reference $genome --known-sites $vreference --output SampleID_post_recal_data_mouse.table - gatk BaseRecalibrator --input SampleID.dedup_reads_mouse_recal.bam --reference $genome --known-sites $vreference --output SampleID_post_recal_data_mouse.table
``` ```
## Disambiguation ## ## Disambiguation ##
We perform alignment using both human and mouse reference genomes in order to perform reads disambiguation using [disambiguate](https://pubmed.ncbi.nlm.nih.gov/27990269/) and discard reads from human mapping that belong to mouse cells. We perform alignment using both human and mouse reference genomes in order to perform reads disambiguation using [disambiguate](https://pubmed.ncbi.nlm.nih.gov/27990269/) and discard reads from human mapping that belong to mouse cells.
<details><summary>Code</summary> `python disambiguate.py -a bwa -s SampleID 02-Alignment_human/SampleID.dedup_reads_rehead_recal.bam 02-Alignment_mouse/SampleID.dedup_reads_mouse_recal.bam`
`python disambiguate.py -a bwa -s Sample_10_395_3 02-Alignment_human/SampleID.dedup_reads_rehead_recal.bam 02-Alignment_mouse/SampleID.dedup_reads_mouse_recal.bam`
</details>
# Variant calling # # Variant calling #
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment