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
No edit summary |
|||
Line 4: | Line 4: | ||
=Code fix/cleanup= | =Code fix/cleanup= | ||
#-> If you really want angsd to exit uncleanly ctrl+c: 0 more times | |||
#fix ./angsd0.591.smartcount/angsd -i angsd0.591.smartcount/DNK02.hg19.flt.sort.rmdup.realign.md.bam -out tmp -minQ 20 -minMapQ 30 -P 2 does nothing but crashes | #fix ./angsd0.591.smartcount/angsd -i angsd0.591.smartcount/DNK02.hg19.flt.sort.rmdup.realign.md.bam -out tmp -minQ 20 -minMapQ 30 -P 2 does nothing but crashes | ||
# check that first line in .arg file is the command used | # check that first line in .arg file is the command used |
Revision as of 08:28, 2 April 2014
Wiki doc
- write how on website: Thetas,Tajima,Neutrality test needs documentation and examples files
- user custom class missing
Code fix/cleanup
- -> If you really want angsd to exit uncleanly ctrl+c: 0 more times
- fix ./angsd0.591.smartcount/angsd -i angsd0.591.smartcount/DNK02.hg19.flt.sort.rmdup.realign.md.bam -out tmp -minQ 20 -minMapQ 30 -P 2 does nothing but crashes
- check that first line in .arg file is the command used
- check domajorminor for -domaf complains about it lagging
- fix -r problem if chromosome contains :
- 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…
Addtional methods and functionality
- make haploid 1dsfs
- simplefy filereading. from glf files
- fix -doFasta for single chromosomes.