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 4: Line 4:


=Code fix/cleanup=
=Code fix/cleanup=
#fix offset in perchunk printout to stderr
#if(doMajorMinor==4&&refToInt[pars->anc[s]]==4) in analysisMajorMinor
#if(doMajorMinor==4&&refToInt[pars->anc[s]]==4) in analysisMajorMinor
# don't quit program if a chr is not included in the -sites file
# don't quit program if a chr is not included in the -sites file

Revision as of 19:05, 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. fix offset in perchunk printout to stderr
  2. if(doMajorMinor==4&&refToInt[pars->anc[s]]==4) in analysisMajorMinor
  3. don't quit program if a chr is not included in the -sites file
  4. check all getoptions stderr output.
  5. -samglfclean not documented
  6. soap usage is not documented
  7. 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.