Question: bedtools multicov - tab-delimited bam files?
0
gravatar for susibing
13 days ago by
susibing20
susibing20 wrote:

Dear all,

I am trying to use bedtools multicov using six bam files as input and a gff file to compare with and unfortunately receive the error message: "It looks as though you have less than 3 columns at line: 1. Are you sure your files are tab-delimited?"

Does anyone know how to find out whether a bam file is delimited or not? And how to change that?

This is the first 5 lines of one of the bam files:

SRR2960489.11401.1      0       chr1    629572  255     31M     *       0       0       AGTAAGGTCAGCTAATTAAGCTATCGGGCCC BBBFFFFFFFFFFIIIIIIIIIIIIIIIIII XA:i:0  MD:Z:15A15      NM:i:1
SRR2960489.33430.1      0       chr1    629572  255     31M     *       0       0       AGTAAGGTCAGCTAATTAAGCTATCGGGCCC BBBFFFFFFFFFFIIIIIIIIIIIIIIIIII XA:i:0  MD:Z:15A15      NM:i:1
SRR2960489.48476.1      0       chr1    629572  255     31M     *       0       0       AGTAAGGTCAGCTAATTAAGCTATCGGGCCC BBBFFFFFFFFFFFIIIIFFFFFFFFFFFFI XA:i:0  MD:Z:15A15      NM:i:1
SRR2960489.50759.1      0       chr1    629572  255     31M     *       0       0       AGTAAGGTCAGCTAATTAAGCTATCGGGCCA BBBFFFFFFFFFFIIIIIIIIIIIIIIIIII XA:i:0  MD:Z:15A14C0    NM:i:2
SRR2960489.71848.1      0       chr1    629572  255     31M     *       0       0       AGTAAGGTCAGCTAATTAAGCTATCGGGCCC BBBFFFFFFFFFFIIIIIIIIIIIIIIIIII XA:i:0  MD:Z:15A15      NM:i:1

I am new to both samtools and bedtools and happy to hear about any suggestions from you.

Thanks in advance

ADD COMMENTlink modified 12 days ago • written 13 days ago by susibing20
2

What is the command of your bedtools multicov? Did you sort your bam and indexed?

ADD REPLYlink written 13 days ago by Ankit60
1

Looks like the error message is not referring to the BAM files but the input BED file specifying the interval/loci for coverage calculation. Can you post the first few lines of the input BED file?

ADD REPLYlink written 13 days ago by Vitis1.7k
1
gravatar for susibing
12 days ago by
susibing20
susibing20 wrote:

Thanks! It was indeed the gff input file that was not tab-delimited.. Thank you very much. (Feeling kind of stupid now..)

ADD COMMENTlink written 12 days ago by susibing20
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1648 users visited in the last hour