bam-readcount empty output
1
0
Entering edit mode
7 weeks 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 • 218 views
ADD COMMENT
1
Entering edit mode
7 weeks ago
ATpoint 57k

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: 2741 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