Skip to content

error with coverage table #340

Description

@KonstantinaRi

after running concoct_coverage_table.py concoct_input.bed S*.sorted.bam > concoct_coverage_table.tsv i keep getting this error :
Errors in BED line 'k141_12453 0 3079 k141_12453.concoct_part_0'
Errors in BED line 'k141_50396 0 4349 k141_50396.concoct_part_0'
Errors in BED line 'k141_3487 0 4389 k141_3487.concoct_part_0'
Errors in BED line 'k141_41388 0 3616 k141_41388.concoct_part_0'
Errors in BED line 'k141_26985 0 6944 k141_26985.concoct_part_0'
Errors in BED line 'k141_41407 0 3213 k141_41407.concoct_part_0'
Errors in BED line 'k141_30590 0 10364 k141_30590.concoct_part_0'
Errors in BED line 'k141_3558 0 5350 k141_3558.concoct_part_0'
Errors in BED line 'k141_55705 0 4501 k141_55705.concoct_part_0'
Errors in BED line 'k141_56953 0 3195 k141_56953.concoct_part_0'
Errors in BED line 'k141_56959 0 4069 k141_56959.concoct_part_0'
Errors in BED line 'k141_56976 0 4505 k141_56976.concoct_part_0'
Errors in BED line 'k141_57027 0 3888 k141_57027.concoct_part_0'
Errors in BED line 'k141_57267 0 3275 k141_57267.concoct_part_0'
Traceback (most recent call last):
File "/home/s/sgeni/anaconda3/envs/concoct_env/bin/concoct_coverage_table.py", line 91, in
generate_input_table(args.bedfile, args.bamfiles, samplenames=samplenames)
File "/home/s/sgeni/anaconda3/envs/concoct_env/bin/concoct_coverage_table.py", line 61, in generate_input_table
df = pd.read_table(fh, header=None)
File "/home/s/sgeni/anaconda3/envs/concoct_env/lib/python3.7/site-packages/pandas/util/_decorators.py", line 311, in wrapper
return func(*args, **kwargs)
File "/home/s/sgeni/anaconda3/envs/concoct_env/lib/python3.7/site-packages/pandas/io/parsers/readers.py", line 683, in read_table
return _read(filepath_or_buffer, kwds)
File "/home/s/sgeni/anaconda3/envs/concoct_env/lib/python3.7/site-packages/pandas/io/parsers/readers.py", line 482, in _read
parser = TextFileReader(filepath_or_buffer, **kwds)
File "/home/s/sgeni/anaconda3/envs/concoct_env/lib/python3.7/site-packages/pandas/io/parsers/readers.py", line 811, in init
self._engine = self._make_engine(self.engine)
File "/home/s/sgeni/anaconda3/envs/concoct_env/lib/python3.7/site-packages/pandas/io/parsers/readers.py", line 1040, in _make_engine
return mapping[engine](self.f, **self.options) # type: ignore[call-arg]
File "/home/s/sgeni/anaconda3/envs/concoct_env/lib/python3.7/site-packages/pandas/io/parsers/c_parser_wrapper.py", line 69, in init
self._reader = parsers.TextReader(self.handles.handle, **kwds)
File "pandas/_libs/parsers.pyx", line 549, in pandas._libs.parsers.TextReader.cinit
i have created a new environment, i have used cut_up and created the bam files with bowtie2.
My bam files look like this:
VN:1.0 SO:coordinate
@sq SN:k141_1.concoct_part_0 LN:3071
@sq SN:k141_38980.concoct_part_0 LN:3213
@sq SN:k141_80949.concoct_part_0 LN:4443
@sq SN:k141_29985.concoct_part_0 LN:4156

the concoct_input.fa like this:

k141_1.concoct_part_0

the bed file like this:
k141_1 0 3071 k141_1.concoct_part_0
k141_38980 0 3213 k141_38980.concoct_part_0
k141_80949 0 4443 k141_80949.concoct_part_0
k141_29985 0 4156 k141_29985.concoct_part_0
k141_92941 0 4210 k141_92941.concoct_part_0
k141_32984 0 5353 k141_32984.concoct_part_0
k141_86952 0 4232 k141_86952.concoct_part_0

i run this to fix the columns just in case:
cut_up_fasta.py concoct_input.fa -c 10000 -o 0 --merge_last -b contigs_10K.bed > /dev/null and i get:
Errors in BED line 'k141_55705 0 4501 k141_55705.concoct_part_0'
Errors in BED line 'k141_56953 0 3195 k141_56953.concoct_part_0'
Errors in BED line 'k141_56959 0 4069 k141_56959.concoct_part_0'
Errors in BED line 'k141_56976 0 4505 k141_56976.concoct_part_0'
Errors in BED line 'k141_57027 0 3888 k141_57027.concoct_part_0'
Errors in BED line 'k141_57267 0 3275 k141_57267.concoct_part_0'
Traceback (most recent call last):
File "/home/anaconda3/envs/concoct_env/bin/concoct_coverage_table.py", line 91, in
generate_input_table(args.bedfile, args.bamfiles, samplenames=samplenames)
File "/home/anaconda3/envs/concoct_env/bin/concoct_coverage_table.py", line 61, in generate_input_table
df = pd.read_table(fh, header=None)
File "/home/anaconda3/envs/concoct_env/lib/python3.7/site-packages/pandas/util/_decorators.py", line 311, in wrapper
return func(*args, **kwargs)
File "/home/anaconda3/envs/concoct_env/lib/python3.7/site-packages/pandas/io/parsers/readers.py", line 683, in read_table
return _read(filepath_or_buffer, kwds)
File "/home/anaconda3/envs/concoct_env/lib/python3.7/site-packages/pandas/io/parsers/readers.py", line 482, in _read
parser = TextFileReader(filepath_or_buffer, **kwds)
File "/home/anaconda3/envs/concoct_env/lib/python3.7/site-packages/pandas/io/parsers/readers.py", line 811, in init
self._engine = self._make_engine(self.engine)
File "/home/anaconda3/envs/concoct_env/lib/python3.7/site-packages/pandas/io/parsers/readers.py", line 1040, in _make_engine
return mapping[engine](self.f, **self.options) # type: ignore[call-arg]
File "/home/anaconda3/envs/concoct_env/lib/python3.7/site-packages/pandas/io/parsers/c_parser_wrapper.py", line 69, in init
self._reader = parsers.TextReader(self.handles.handle, **kwds)
File "pandas/_libs/parsers.pyx", line 549, in pandas._libs.parsers.TextReader.cinit
pandas.errors.EmptyDataError: No columns to parse from file

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions