ANGSD: Analysis of next generation Sequencing Data

Latest tar.gz version is (0.938/0.939 on github), see Change_log for changes, and download it here.

Depth

From angsd
Revision as of 00:15, 12 October 2012 by Thorfinn (talk | contribs) (→‎Output)
Jump to navigation Jump to search

Depth

This method will find depth distribution for every sample and for all samples jointly.

-doDepth 1

which requires -doCount 1.

Options

- maxDepth [int]

Defaults to 100. Sites with more than maxDepth reads are counted as having 100 reads.


Example

Find the depth for chromsome 7 for the first 5 samples in smallBam.filelist

angsd -bam smallBam.filelist -doDepth 1 -out CAN_DEL_ME -doCounts 1 -r 7: -nInd 5


Output

1147 	4240 	8357 	13017 	16077 	15481 	12871 	10738 	7288 	4537 	2998 	1638 	889 	451 	211 	92 	50 	
1831 	6725 	12572 	17090 	17451 	15359 	11588 	8274 	4538 	2551 	1199 	545 	209 	83 	20 	33 	13 	9 	
12749 	23974 	25261 	18900 	10294 	5388 	2326 	850 	241 	83 	28 	2 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 
15366 	25983 	24257 	16854 	9152 	4832 	2585 	791 	214 	42 	20 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	
13131 	24070 	25347 	18949 	10514 	5090 	1991 	688 	213 	82 	21 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 


The outputfile .depth has one line per individual, with maxDepth entries.