Generate histogram Via BBmap tool
0
0
Entering edit mode
3.0 years ago
Fizzah ▴ 30

Hello there, I am trying to generate histogram for paired end reads via bbmap tool program. However I am getting this output:

root@cemb-ls1-pgs:/mnt/cemb/S017679/raw# /home/fizzah/bbmap/reformat.sh BU21R1.fastq qchist=qchist.txt
java -ea -Xms300m -cp /home/fizzah/bbmap/current/ jgi.ReformatReads BU21R1.fastq qchist=qchist.txt
Executing jgi.ReformatReads [BU21R1.fastq, qchist=qchist.txt]

Set qcount histogram output to qchist.txt
No output stream specified.  To write to stdout, please specify 'out=stdout.fq' or similar.
Input is being processed as unpaired
Input:                      37926174 reads              5726852274 bases
Output:                     37926174 reads (100.00%)    5726852274 bases (100.00%)

Time:                           50.960 seconds.
Reads Processed:      37926k    744.23k reads/sec
Bases Processed:       5726m    112.38m bases/sec

when I specify stream output, it says:

File qchist.txt exists and overwrite=false

Please guide about its solution

histogram bbmap • 956 views
ADD COMMENT
0
Entering edit mode

BBMap does not produce graphical output. What you will get is a text based histogram.

In general for all BBMap programs you have to explicitly add option ow=t to overwrite existing files. Be aware that this option instantly wipes the contents of existing files so if you are not careful you can lose data.

ADD REPLY
0
Entering edit mode

Ok. So can you name some tool that can give graphical outcome based histogram for paired end reads?

ADD REPLY

Login before adding your answer.

Traffic: 1765 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6