bam-readcount empty output
1
0
Entering edit mode
8 months ago
D_Jens • 0

Hi,

I am using bam-readcount. However, when I run the following input I get a blank output;

$ bam-readcount -f NCBI_ReferenceGenome.fna -l decomposed.vcf Bowtie2.bam [-b 20] > output.readcount

After running the script, I get the following;

Minimum mapping quality is set to 0
Expect library: Bowtie2 in BAM

The documentation says the output should be a tab-delimited file with column names, however output.readcount is empty. Assuming maybe my files are the problem, I should still get the column names in the output, right? I am a novice and am not sure where to look further in troubleshooting.

Thank you for your help

Protocol

bam-readcount Documentation

sam bam-readcount bam • 381 views
ADD COMMENT
1
Entering edit mode
8 months ago
ATpoint 64k

The file given to -l should be a tab-sep file (chr-start-end) not a VCF.

ADD COMMENT
0
Entering edit mode

I went back and converted my BAM files to BED files and it worked! Thank you!

ADD REPLY

Login before adding your answer.

Traffic: 678 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