This is a test version of Biostars. For the public version, visit https://www.biostars.org.
why JGIL no output

I am using JGIL on sequence data of fruit fly. Command as following:

./JGIL/jgil-1.6/jgil_m -M SRR1.bam SRR1.bam SRR3.bam -n 15 -f ref/D_simulans_ref.fasta -m 1 -o new_file

merge three bam files; inbreeding for 15 generation; using reference genome; output format as type 1; output prefix as "new_file" There is no error on the running log, but no output files!!

log:

BAMs to be processed are merged.

Done parsing header text.

Sample names:

JGIL v1.6

Doing SNP calling.

Number of generations for JGIL: 15

SNP quality threshold = 0 (this corresponds to a p-value of 1)

Mapping quality threshold = 0

Base quality threshold = 0

SNP output formats: tab

Prefix for output files = new_file

Reference sequence FASTA file(s): ref/D_simulans_ref.fasta

Running in unthreaded mode.

Processing 3 BAM files.

Processing reference: Scf_2L

Processing reference: Scf_2R

Processing reference: Scf_3L

Processing reference: Scf_3R

Processing reference: Scf_X

Reads considered = 19358407

Reads rejected - flagged unmapped 3

Reads rejected - low quality 0

Reads rejected - flagged multihit 0

Reads rejected - anomalous pair 0

Does anyone know there is the problem? Thanks very much!

sequencing genome alignment

0 answers

No answers yet.

Log in to answer this question.