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.

User:Thorfinn: Difference between revisions

From angsd
Jump to navigation Jump to search
No edit summary
Line 21: Line 21:
h1 h2 h3 = 9 10 6  has less than 3 blocks. skipping
h1 h2 h3 = 9 10 6  has less than 3 blocks. skipping
h1…
h1…
</pre>
fix alignement of output
<pre>
-> [analysisKeepList.cpp] -sites is still beta, use at own risk...
Filter file contains major/minor information to use these in analysis supper '-doMajorMinor 3'
-> Reading fasta: /space/genomes/refgenomes/hg19/gatk/hg19.fa
Will call genotypes using sample allele frequencies
-> Parsing 33 number of samples
-> Printing at chr: 1 pos:566303 chunknumber 1800Segmentation fault (core dumped)
</pre>
</pre>



Revision as of 17:39, 5 March 2014

Wiki doc

  1. write how on website: Thetas,Tajima,Neutrality test needs documentation and examples files
  2. user custom class missing

Code fix/cleanup

  1. check all getoptions stderr output.
  2. -samglfclean not documented
  3. soap usage is not documented
  1. from morten shitoutput
h1 h2 h3 = 7 10 5  has less than 3 blocks. skipping
h1 h2 h3 = 8 10 5  has less than 3 blocks. skipping
h1 h2 h3 = 10 11 5  has less than 3 blocks. skipping
h1 h2 h3 = 7 12 5  has less than 3 blocks. skipping
h1 h2 h3 = 10 12 5  has less than 3 blocks. skipping
h1 h2 h3 = 10 13 5  has less than 3 blocks. skipping
h1 h2 h3 = 10 15 5  has less than 3 blocks. skipping
h1 h2 h3 = 7 10 6  has less than 3 blocks. skipping
h1 h2 h3 = 8 10 6  has less than 3 blocks. skipping
h1 h2 h3 = 9 10 6  has less than 3 blocks. skipping
h1…

fix alignement of output

	-> [analysisKeepList.cpp] -sites is still beta, use at own risk...
Filter file contains major/minor information to use these in analysis supper '-doMajorMinor 3'
	-> Reading fasta: /space/genomes/refgenomes/hg19/gatk/hg19.fa
Will call genotypes using sample allele frequencies
	-> Parsing 33 number of samples 
	-> Printing at chr: 1 pos:566303 chunknumber 1800Segmentation fault (core dumped)

Addtional methods and functionality

  1. make haploid 1dsfs
  2. simplefy filereading. from glf files
  3. fix -doFasta for single chromosomes.