Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

BioAlignments issue #7

Open
code12412f opened this issue Aug 2, 2020 · 3 comments
Open

BioAlignments issue #7

code12412f opened this issue Aug 2, 2020 · 3 comments

Comments

@code12412f
Copy link

Hello,

I'm getting this issue for some of my Fastq files.

[ Info: 2020-08-02 07:06:27] Program start
[ Info: 2020-08-02 07:06:28] Handing Paired-end sample: /projects/b1042/MeeksLab/abacus/6/62C51828A4DD56C3C570218C79F140DB_1_val_
1.fq.gz
[ Info: 2020-08-02 07:06:28] Aligning
[ Info: 2020-08-02 07:35:48] Reading aligned results
ERROR: LoadError: BioAlignments.SAM.Reader file format error on line 1 ~>"D00760:3"
Stacktrace:
[1] error(::Type, ::String, ::Int64, ::String, ::String) at ./error.jl:42
[2] _read!(::BioAlignments.SAM.Reader, ::BioCore.Ragel.State{BufferedStreams.BufferedInputStream{IOStream}}, ::BioAlignments.SAM.Record) at /home/kam9898/.julia/packages/BioCore/YBJvb/src/ReaderHelper.jl:164
[3] read!(::BioAlignments.SAM.Reader, ::BioAlignments.SAM.Record) at /home/kam9898/.julia/packages/BioCore/YBJvb/src/ReaderHelper.jl:134
[4] tryread!(::BioAlignments.SAM.Reader, ::BioAlignments.SAM.Record) at /home/kam9898/.julia/packages/BioCore/YBJvb/src/IO.jl:73
[5] iterate at /home/kam9898/.julia/packages/BioCore/YBJvb/src/IO.jl:84 [inlined]
[6] read_alignrs(::Dict{String,Any}, ::Array{Any,1}, ::Array{Any,1}, ::String) at /projects/b1042/MeeksLab/catt/CATT/catt.jl:230
[7] mainflow(::Dict{String,Any}, ::Array{Any,1}, ::Array{Any,1}, ::String, ::String) at /projects/b1042/MeeksLab/catt/CATT/catt.jl:700
[8] proc(::Dict{String,Any}) at /projects/b1042/MeeksLab/catt/CATT/catt.jl:764
[9] top-level scope at /projects/b1042/MeeksLab/catt/CATT/catt.jl:825
[10] include at ./boot.jl:328 [inlined]
[11] include_relative(::Module, ::String) at ./loading.jl:1105
[12] include(::Module, ::String) at ./Base.jl:31
[13] exec_options(::Base.JLOptions) at ./client.jl:287
[14] _start() at ./client.jl:460
in expression starting at /projects/b1042/MeeksLab/catt/CATT/catt.jl:824

This is the command I used:
for prefix in $(ls /projects/b1042/MeeksLab/abacus/6/62C*gz | rev | cut -c 14- | rev | uniq)
do
/projects/b1042/MeeksLab/catt/CATT/./catt --f1 ${prefix}1_val_1.fq.gz --f2 ${prefix}2_val_2.fq.gz -o $prefix -t 4
done

@aqzas
Copy link
Collaborator

aqzas commented Aug 3, 2020

Hello,@code12412f
May I ask what version of BioAligment you are using ?
And could you mind give me a test sample of your input file, which could help find the problem

@code12412f
Copy link
Author

I'm using version 1.0.1. It's a paired sample
1.fq.txt
2.fq.txt

@aqzas
Copy link
Collaborator

aqzas commented Aug 5, 2020

@code12412f
Sorry, I can't replicate this bugs. Please provide the the detail information of your Julia environment, by

using Pkg
Pkg.status()

By the way, it is a better choice to use Docker (The recommended way). The Julia environment is in a rapidly changing environment and there is a poor compatibility between package in different versions.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants