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: Difference between revisions

From angsd
Jump to navigation Jump to search
Line 20: Line 20:
==Output==
==Output==
<pre>
<pre>
1147 4240 8357 13017 16077 15481 12871 10738 7288 4537 2998 1638 889 451 211 92 50
#.depthSample
1147 4240 8357 13017 16077 15481 12871 10738 7288 4537 2998 1638 889 451 211 92 50     14
1831 6725 12572 17090 17451 15359 11588 8274 4538 2551 1199 545 209 83 20 33 13 9
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  
12749 23974 25261 18900 10294 5388 2326 850 241 83 28 2 0 0 0 0 0 0 0 0 0 0 0  
Line 26: Line 27:
13131 24070 25347 18949 10514 5090 1991 688 213 82 21 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  
</pre>
</pre>
 
<pre>
#.depthGlobal
150500 21708 39928 65820 87290 97337 105127 115961 123207 126721 129377 127351 125596 118588 107096 93171 78612 66905 53952 43224 34777 25672 19393 14431 9771 6958 4837 3149 2027 1362 919 471 308 198 79 48 41 8 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
</pre>


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

Revision as of 00:16, 12 October 2012

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

#.depthSample
1147 	4240 	8357 	13017 	16077 	15481 	12871 	10738 	7288 	4537 	2998 	1638 	889 	451 	211 	92 	50     14 	
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 
#.depthGlobal
150500 	21708 	39928 	65820 	87290 	97337 	105127 	115961 	123207 	126721 	129377 	127351 	125596 	118588 	107096 	93171 	78612 	66905 	53952 	43224 	34777 25672 	19393 	14431 	9771 	6958 	4837 	3149 	2027 	1362 	919 	471 	308 	198 	79 	48 	41 	8 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0 	0

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