X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=samtools.1;fp=samtools.1;h=4b3f75a1926f0c4248f228f6378d2ff06e3e0c77;hb=125ac895854cf760a04192a257a4279f2c541164;hp=718ac392016705cb758d4cded7eeded2755faff2;hpb=e08e3faf4cc66727674ecacb6914e69b9a5df8c9;p=samtools.git diff --git a/samtools.1 b/samtools.1 index 718ac39..4b3f75a 100644 --- a/samtools.1 +++ b/samtools.1 @@ -30,7 +30,7 @@ bcftools index in.bcf .PP bcftools view in.bcf chr2:100-200 > out.vcf .PP -bcftools view -vc in.bcf > out.vcf 2> out.afs +bcftools view -Nvm0.99 in.bcf > out.vcf 2> out.afs .SH DESCRIPTION .PP @@ -140,17 +140,38 @@ to another samtools command. .TP .B tview -samtools tview [ref.fasta] +samtools tview +.RB [ \-p +.IR chr:pos ] +.RB [ \-s +.IR STR ] +.RB [ \-d +.IR display ] +.RI +.RI [ref.fasta] Text alignment viewer (based on the ncurses library). In the viewer, press `?' for help and press `g' to check the alignment start from a region in the format like `chr10:10,000,000' or `=10,000,000' when viewing the same reference sequence. +.B Options: +.RS +.TP 14 +.BI -d \ display +Output as (H)tml or (C)urses or (T)ext +.TP +.BI -p \ chr:pos +Go directly to this position +.TP +.BI -s \ STR +Display only reads from this sample or read group +.RE + .TP .B mpileup -.B samtools mpileup -.RB [ \-EBug ] +samtools mpileup +.RB [ \-EBugp ] .RB [ \-C .IR capQcoef ] .RB [ \-r @@ -297,6 +318,10 @@ Phred-scaled gap open sequencing error probability. Reducing .I INT leads to more indel calls. [40] .TP +.BI -p +Apply -m and -F thresholds per sample to increase sensitivity of calling. +By default both options are applied to reads pooled from all samples. +.TP .BI -P \ STR Comma dilimited list of platforms (determined by .BR @RG-PL ) @@ -570,6 +595,8 @@ Minimum base quality to be used in het calling. [13] .IR mutRate ] .RB [ \-p .IR varThres ] +.RB [ \-m +.IR varThres ] .RB [ \-P .IR prior ] .RB [ \-1 @@ -652,6 +679,12 @@ Call per-sample genotypes at variant sites (force -c) .BI -i \ FLOAT Ratio of INDEL-to-SNP mutation rate [0.15] .TP +.BI -m \ FLOAT +New model for improved multiallelic and rare-variant calling. Another +ALT allele is accepted if P(chi^2) of LRT exceeds the FLOAT threshold. The +parameter seems robust and the actual value usually does not affect the results +much; a good value to use is 0.99. This is the recommended calling method. [0] +.TP .BI -p \ FLOAT A site is considered to be a variant if P(ref|D)