telseq
telseq copied to clipboard
Error while using a list of paths
Hi,
Thank you for this helpful tool. I am running into a problem while I was trying to use a file containing list of paths as input. Every BAM file becomes unreadable except for the last row in the file.
My input command looks like this:
telseq -f "/my/own/path/telseq_test/bamlist.txt" > "/my/own/path/telseq_test/test_list.txt"
And the ouput looks like this:
> Start analysing BAM /my/own/path/first.bam
> Warning: can't find RG tag in the BAM header
> Warning: treat all reads in BAM as if they were from a same sample
> [scan] total reads in BAM scanned 1
> Completed scanning BAM
> Start analysing BAM /my/own/path/second.bam
> Warning: can't find RG tag in the BAM header
> Warning: treat all reads in BAM as if they were from a same sample
> [scan] total reads in BAM scanned 1
> Completed scanning BAM
> Start analysing BAM /my/own/path/third.bam
> Warning: can't find RG tag in the BAM header
> Warning: treat all reads in BAM as if they were from a same sample
> [scan] total reads in BAM scanned 1
> Completed scanning BAM
> Start analysing BAM /my/own/path/fourth.bam
> Warning: can't find RG tag in the BAM header
> Warning: treat all reads in BAM as if they were from a same sample
> [scan] total reads in BAM scanned 1
> Completed scanning BAM
> Start analysing BAM /my/own/path/fifth.bam
> Specified BAM has 4 read groups
> [scan] processed 10000000 reads
> [scan] processed 20000000 reads
> ...
I created the bamlist as instructed: one path per row. All the files works normally while they were used as a single input BAM file. How can I fix this? Thanks!