Skip to content

Commit

Permalink
Update RNA-Seq example: For a stranded paired-end library (#796)
Browse files Browse the repository at this point in the history
Merging temporary files with using a wrong name - rev1.bam instead of a.rev1.bam and similar
  • Loading branch information
janoppelt authored and dpryan79 committed Feb 18, 2019
1 parent d3bd4a5 commit 23b2c0b
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions docs/content/tools/bamCoverage.rst
Original file line number Diff line number Diff line change
Expand Up @@ -123,7 +123,7 @@ To get the file for transcripts that originated from the **forward strand**:
$ samtools view -b -f 80 a.bam > a.fwd2.bam
# combine the temporary files
$ samtools merge -f fwd.bam fwd1.bam fwd2.bam
$ samtools merge -f fwd.bam a.fwd1.bam a.fwd2.bam
# index the filtered BAM file
$ samtools index fwd.bam
Expand All @@ -147,7 +147,7 @@ To get the file for transcripts that originated from the **reverse strand**:
$ samtools view -b -f 64 -F 16 a.bam > a.rev2.bam
# merge the temporary files
$ samtools merge -f rev.bam rev1.bam rev2.bam
$ samtools merge -f rev.bam a.rev1.bam a.rev2.bam
# index the merged, filtered BAM file
$ samtools index rev.bam
Expand Down

0 comments on commit 23b2c0b

Please sign in to comment.