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.
Allele Counts: Difference between revisions
(Created page with "==dumping counts== -dumpCounts 1, print overall depth in the .pos file -dumpCounts 2, 1) plus .counts with overall sample depth, each line nInd long -dumpCounts 3,...") |
|||
| Line 1: | Line 1: | ||
==dumping counts== | ==dumping counts== | ||
; -dumpCounts [int] | |||
1: print overall depth in the .pos file. This depth is the sum of reads covering a sites for all individuals. The first colum is the chromosome, the second it the position the third is the total depth | |||
<pre> | |||
1 13999959 3 | |||
1 13999960 3 | |||
1 13999961 3 | |||
1 13999962 3 | |||
1 13999963 4 | |||
1 13999964 5 | |||
1 13999965 6 | |||
1 13999966 6 | |||
1 13999967 6 | |||
1 13999968 6 | |||
</pre> | |||
2: prints the depth of each individual. Example of the depth of 10 individuals. Each line corresponce to the same line in the postion file. | |||
<pre> | |||
0 0 0 0 0 0 0 1 0 2 | |||
0 0 0 0 0 0 0 1 0 2 | |||
0 0 0 0 0 0 0 1 0 2 | |||
0 0 0 0 0 0 0 1 0 2 | |||
0 0 0 0 0 0 0 1 0 3 | |||
0 0 0 0 0 1 0 1 0 3 | |||
1 0 0 0 0 1 0 1 0 3 | |||
1 0 0 0 0 1 0 1 0 3 | |||
1 0 0 0 0 1 0 1 0 3 | |||
1 0 0 0 0 1 0 1 0 3 | |||
</pre> | |||
3: Prints the depth for each of the four bases for each indivial for each site. Example with the first four column belonging to the first individuals the counts of the number of A C G and Ts. Only two indivduals are shown. Each line corresponce to the same line in the postion file. | |||
<pre> | |||
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 1 0 0 0 0 0 0 ... | |||
0 0 1 0 0 0 0 0 ... | |||
1 0 0 0 0 0 0 0 ... | |||
0 0 1 0 0 0 0 0 ... | |||
</pre> | |||
===requred=== | |||
In order to print the counts the options '-doCounts' have to be used and the input data needs to be sequence data. | |||
==Example== | |||
Print the individuals depth from bam files | |||
<pre> | |||
./angsd -out out -doCounts 1 -dumpCounts 2 -bam bam.filelist | |||
</pre> | |||
Revision as of 18:53, 18 June 2012
dumping counts
- -dumpCounts [int]
1: print overall depth in the .pos file. This depth is the sum of reads covering a sites for all individuals. The first colum is the chromosome, the second it the position the third is the total depth
1 13999959 3 1 13999960 3 1 13999961 3 1 13999962 3 1 13999963 4 1 13999964 5 1 13999965 6 1 13999966 6 1 13999967 6 1 13999968 6
2: prints the depth of each individual. Example of the depth of 10 individuals. Each line corresponce to the same line in the postion file.
0 0 0 0 0 0 0 1 0 2 0 0 0 0 0 0 0 1 0 2 0 0 0 0 0 0 0 1 0 2 0 0 0 0 0 0 0 1 0 2 0 0 0 0 0 0 0 1 0 3 0 0 0 0 0 1 0 1 0 3 1 0 0 0 0 1 0 1 0 3 1 0 0 0 0 1 0 1 0 3 1 0 0 0 0 1 0 1 0 3 1 0 0 0 0 1 0 1 0 3
3: Prints the depth for each of the four bases for each indivial for each site. Example with the first four column belonging to the first individuals the counts of the number of A C G and Ts. Only two indivduals are shown. Each line corresponce to the same line in the postion file.
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 1 0 0 0 0 0 0 ... 0 0 1 0 0 0 0 0 ... 1 0 0 0 0 0 0 0 ... 0 0 1 0 0 0 0 0 ...
requred
In order to print the counts the options '-doCounts' have to be used and the input data needs to be sequence data.
Example
Print the individuals depth from bam files
./angsd -out out -doCounts 1 -dumpCounts 2 -bam bam.filelist