Page 1 of 7 123 ... LastLast
Results 1 to 10 of 67

Thread: BAM Analysis Kit (Windows and Linux). New release.

  1. #1
    Registered Users
    Posts
    165
    Sex
    Nationality
    Finnish
    Y-DNA (P)
    R1b-Z142
    mtDNA (M)
    H10g

    BAM Analysis Kit (Windows and Linux). New release.

    I have released a new version of BAM Analysis kit. Please report any bugs you find!

    Highlights:

    -Update SNP filters to include more SNPs (23andMe V3 to V5, FTDNA old and new, Ancestry)

    -Replace samtools and bcftools with native versions. Htslib has fixed a bug that prevented using native samtools/bcftools. Using native version of samtools speeds up handling large BAM files.

    -Remove old mtDNA haplogroup prediction

    -Add Yleaf v2 and use it to predict Y-haplogroup

    -Update haplogrep to 2.1.20

    Download the zip file or use a Git client.

    https://github.com/teepean/BAM-Analysis-Kit/releases

    I have also released an experimental Linux version of BAM Analysis Kit. Required software: bedtools, samtools 1.7 or higher, bcftools 1.7 or higher, gawk, Python 3.6 or higher with pandas installed

    https://github.com/teepean/BAM-Analy...ases/tag/v2.09
    Last edited by teepean47; 06-03-2019 at 03:58 AM.

  2. The Following 10 Users Say Thank You to teepean47 For This Useful Post:

     AbdoNumen (06-03-2019),  ArmandoR1b (06-02-2019),  DMXX (06-02-2019),  evon (01-03-2020),  pmokeefe (08-01-2019),  poi (06-02-2019),  Robert1 (06-04-2019),  shadowhite (11-02-2019),  Tonev (06-02-2019),  tsunami (08-10-2019)

  3. #2
    Registered Users
    Posts
    2,298
    Sex

    Windows BAK v2.091 fails for me. I tested it with a random file. Here is the first error with v2.091.

    Code:
    Input BAM : "C:\BAM\Crusades\SI-41.bam"
    
    Pre-execution Cleanup ...
    
    Sorting ...
    sort: invalid option -- '@'
    open: No such file or directory
    [bam_sort_core] fail to open file 4
    BAM Analysis Kit 1.8 works fine for the same file. The BAM is in a different directory than the programs.
    Last edited by ArmandoR1b; 06-02-2019 at 04:31 PM.

  4. The Following 2 Users Say Thank You to ArmandoR1b For This Useful Post:

     Robert1 (06-04-2019),  Tonev (06-02-2019)

  5. #3
    Registered Users
    Posts
    38
    Sex
    Location
    Bulgaria
    Ethnicity
    Bulgarian
    Nationality
    Bulgarian
    Y-DNA (P)
    G-FT9612 G2a2b2a1a..
    mtDNA (M)
    H-a* (C16291T)

    Bulgaria
    Quote Originally Posted by ArmandoR1b View Post
    Windows BAK v2.091 fails for me. I tested it with a random file. Here is the first error with v2.091.

    Code:
    Input BAM : "C:\BAM\Crusades\SI-41.bam"
    
    Pre-execution Cleanup ...
    
    Sorting ...
    sort: invalid option -- '@'
    open: No such file or directory
    [bam_sort_core] fail to open file 4
    BAM Analysis Kit 1.8 works fine for the same file. The BAM is in a different directory than the programs.
    Same errors and behaviour for me too...

  6. The Following User Says Thank You to Tonev For This Useful Post:

     ArmandoR1b (06-02-2019)

  7. #4
    Registered Users
    Posts
    165
    Sex
    Nationality
    Finnish
    Y-DNA (P)
    R1b-Z142
    mtDNA (M)
    H10g

    I have released a fix, version 2.0.92. For some reason Git did not update the old executables.

    https://github.com/teepean/BAM-Analysis-Kit/releases

  8. #5
    Registered Users
    Posts
    115
    Sex

    Were the issues Generalissimo brought up ever addressed?

  9. #6
    Registered Users
    Posts
    165
    Sex
    Nationality
    Finnish
    Y-DNA (P)
    R1b-Z142
    mtDNA (M)
    H10g

    Quote Originally Posted by xenus View Post
    Were the issues Generalissimo brought up ever addressed?
    If you mean base alignment quality then yes. BAQ is not enabled in samtools or GATK.

    I have been working on a new tool based on pileupCaller. The problem with that tool at the moment is that I haven't figured out a decent way to annotate the results.

    https://github.com/teepean/FADAK

  10. #7
    Registered Users
    Posts
    115
    Sex

    Quote Originally Posted by teepean47 View Post
    If you mean base alignment quality then yes. BAQ is not enabled in samtools or GATK.

    I have been working on a new tool based on pileupCaller. The problem with that tool at the moment is that I haven't figured out a decent way to annotate the results.

    https://github.com/teepean/FADAK
    What annotations do you need?
    Are you just trying to work pileupCaller into a pipeline?

  11. #8
    Registered Users
    Posts
    165
    Sex
    Nationality
    Finnish
    Y-DNA (P)
    R1b-Z142
    mtDNA (M)
    H10g

    Quote Originally Posted by xenus View Post
    What annotations do you need?
    Are you just trying to work pileupCaller into a pipeline?
    The aim is to replace GATK with pileUpcaller or any other program (snpAD, ANGSD, FreeBayes) suitable for ancient DNA. One aim is to produce Gedmatch compatible files and therefore annotating missing rsids is a necessity. The pipeline has to work on Windows but so far pileupCaller is the only one that is working and is stable for everyday usage. The problem with pileupCaller is that is not multi-threading so could be a lot faster but I have no experience with Haskell.

  12. #9
    Registered Users
    Posts
    115
    Sex

    Quote Originally Posted by teepean47 View Post
    The aim is to replace GATK with pileUpcaller or any other program (snpAD, ANGSD, FreeBayes) suitable for ancient DNA. One aim is to produce Gedmatch compatible files and therefore annotating missing rsids is a necessity. The pipeline has to work on Windows but so far pileupCaller is the only one that is working and is stable for everyday usage. The problem with pileupCaller is that is not multi-threading so could be a lot faster but I have no experience with Haskell.
    If the problem is just parallelizing without shared resources it should be easy. My Haskell ability is rudimentary at best and Haskell's Cabal package manager has always wanted to lock me in dependency hell which I wouldn't imagine to be any nicer on windows.

    Seeing that your repo includes Cygwin i'm wondering if that's just because it works better for your use case than WSL or does WSL just not work at all?

  13. #10
    Registered Users
    Posts
    165
    Sex
    Nationality
    Finnish
    Y-DNA (P)
    R1b-Z142
    mtDNA (M)
    H10g

    Quote Originally Posted by xenus View Post
    If the problem is just parallelizing without shared resources it should be easy. My Haskell ability is rudimentary at best and Haskell's Cabal package manager has always wanted to lock me in dependency hell which I wouldn't imagine to be any nicer on windows.
    Cabal works on Windows. Most of the time. I am currently experimenting with parallel options (RTS+) as pileupCaller is compiled with the -threaded option.

    Seeing that your repo includes Cygwin i'm wondering if that's just because it works better for your use case than WSL or does WSL just not work at all?
    I am using Windows/Cygwin combination so that the program is more accessible to average users. Linux test version is working fine using WSL.

Page 1 of 7 123 ... LastLast

Similar Threads

  1. Replies: 166
    Last Post: 01-25-2018, 04:55 PM
  2. little script ( windows + R) converting to 23&me
    By anglesqueville in forum DNA.Land
    Replies: 2
    Last Post: 08-16-2017, 11:56 AM
  3. Replies: 3
    Last Post: 04-28-2017, 03:18 PM

Posting Permissions

  • You may not post new threads
  • You may not post replies
  • You may not post attachments
  • You may not edit your posts
  •