From Genome Analysis Wiki
Jump to navigationJump to search
379 bytes added
, 13:46, 29 October 2010
Line 1: |
Line 1: |
− | == Trim BAM == | + | == Polish BAM == |
− | The <code>trimBam</code> program is released as part of the StatGen Library & Tools download. | + | The <code>polishBam</code> program is released as part of the StatGen Library & Tools download. |
| | | |
− | <code>trimBam</code> trims the end of reads in a SAM/BAM file, changing read ends to ‘N’ and quality to ‘!’. | + | <code>polishBam</code> trims the end of reads in a SAM/BAM file, changing read ends to ‘N’ and quality to ‘!’. |
| | | |
| | | |
| === Parameters === | | === Parameters === |
| <pre> | | <pre> |
− | Required Parameters:
| + | Required parameters: |
− | inFile : the SAM/BAM file to be read | + | -i/--in : input BAM file |
− | outFile : the SAM/BAM file to be written | + | -o/--out : output BAM file |
− | num-bases-to-trim-on-each-side : the number of bases/qualities to trim from each side | + | Optional parameters: |
| + | -v : turn on verbose mode |
| + | -l/--log : writes logfile. <outBamFile>.log will be used if value is unspecified |
| + | --HD : add @HD header line |
| + | --RG : add @RG header line |
| + | --PG : add @PG header line |
| + | -f/--fasta : fasta reference file to compute MD5sums and update SQ tags |
| + | --AS : AS tag for genome assembly identifier |
| + | --UR : UR tag for @SQ tag (if different from --fasta) |
| + | --SP : SP tag for @SQ tag |
| + | --checkSQ : check the consistency of SQ tags (SN and LN) with existing header lines. Must be used with --fasta option |
| </pre> | | </pre> |
| | | |
Line 22: |
Line 32: |
| === Example Output === | | === Example Output === |
| <pre> | | <pre> |
− | Arguments in effect:
| + | polishBAM (options) --in=<inBamFile> --out=<outBamFile> |
− | Input file : testFiles/testSam.sam
| |
− | Output file : results/trimSam.sam
| |
− | #TrimBases : 2
| |
− | | |
− | Number of records read = 10
| |
− | Number of records written = 10
| |
| </pre> | | </pre> |
| | | |