Skip to content

Commit

Permalink
Update fastq_to_bam.rst
Browse files Browse the repository at this point in the history
  • Loading branch information
michalsakin committed Jan 24, 2021
1 parent a2d56a4 commit 264d8f9
Showing 1 changed file with 1 addition and 4 deletions.
5 changes: 1 addition & 4 deletions docs/source/fastq_to_bam.rst
Original file line number Diff line number Diff line change
Expand Up @@ -28,9 +28,6 @@ From fastq to final valid pairs bam file
bwa mem -5SP -T0 -t16 hg38.fasta OmniC_2M_R1.fastq OmniC_2M_R2.fastq| pairtools parse --min-mapq 40 --walks-policy 5unique --max-inter-align-gap 30 --nproc-in 8 --nproc-out 8 --chroms-path hg38.genome | pairtools sort --tmpdir=/home/ubuntu/ebs/temp/ --nproc 16|pairtools dedup --nproc-in 8 --nproc-out 8 --mark-dups --output-stats stats.txt|pairtools split --nproc-in 8 --nproc-out 8 --output-pairs mapped.pairs --output-sam -|samtools view -bS -@16 | samtools sort -@16 -o mapped.PT.bam;samtools index mapped.PT.bam
|clock| The full command above, with 2M read pairs on an Ubuntu 18.04 machine with 16 CPUs and 64GiB was completed in less than 5 minutes.
On the same machine type.

|clock| The full command above, with 2M read pairs on an Ubuntu 18.04 machine with 16 CPUs and 64GiB was completed in less than 5 minutes.
On the same machine type.

Expand Down Expand Up @@ -253,4 +250,4 @@ Index the bam file:
The `mapped.PT.bam` is the final bam file that will be used downstream steps.


The above steps resulted in multiple intermediate files, to simplify the process and avoid intermediate files, you can pipe the steps as in the example above (:ref: `fastq to final valid pairs bam file - for the impatient <impatient>`)
The above steps resulted in multiple intermediate files, to simplify the process and avoid intermediate files, you can pipe the steps as in the example above (:ref: `fastq to final valid pairs bam file - for the impatient <impatient>`)

0 comments on commit 264d8f9

Please sign in to comment.