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
Jump to navigation
Jump to search
Line 4: | Line 4: | ||
=Code fix/cleanup= | =Code fix/cleanup= | ||
#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 | ||
#check all getoptions stderr output. | #check all getoptions stderr output. |
Revision as of 17:49, 5 March 2014
Wiki doc
- write how on website: Thetas,Tajima,Neutrality test needs documentation and examples files
- user custom class missing
Code fix/cleanup
- if(doMajorMinor==4&&refToInt[pars->anc[s]]==4) in analysisMajorMinor
- don't quit program if a chr is not included in the -sites file
- check all getoptions stderr output.
- -samglfclean not documented
- soap usage is not documented
- 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
- make haploid 1dsfs
- simplefy filereading. from glf files
- fix -doFasta for single chromosomes.