Hello everyone,
Basically, I have clustered and annotated the barcodes in R, then I wanted to look at reads in several particular clusters in IGV. I generated a barcode list following the 10X tutorial as the picture below shown, briefly subset the clusters in R then tagged the barcodes with CB:Z: and saved as filter.txt. While I could generate a 2 Kb header file, the later step always gives me a 0 Kb filtered_SAM_body file.
I tried:
samtools view -D CB:FILTER.txt BAMFILE.bam -h -b -o FILTEREDBAM.bam
===
(following the a third thread...[lost the page])
samtools view BAMFILE.bam | LC_ALL=C grep -F -f FILTER.txt > filtered_SAM_body
===
(modified 10X)
samtools view -H BAMFILE.bam > SAM_header
samtools view BAMFILE.bam -X INDEX.bai | LC_ALL=C grep -F -f FILTER.txt > filtered_SAM_body
===
(10X)
samtools view -H BAMFILE.bam > SAM_header
samtools view BAMFILE.bam | LC_ALL=C grep -F -f FILTER.txt > filtered_SAM_body
===
I also tried adding CB: prior to filter path or removing CB:Z from the filter.txt file... none worked all gave me the 0 Kb file.
Any Suggestions will be very very very much appreciated.
Thank you.
0 answers
No answers yet.
Log in to answer this question.
-
Samtools merge bam issue with header tags
written by sam •
-
How to filter the BAM file using barcode list produced by 10x pipelines?
written by Info.shi •
-
paired-end reads, idxstats, and samtools and grep
written by Bogdan
-
Output of samtools view, what does the third column actually represent?
written by DNAngel
-
Disappearing CB, the bam tag after samtools sort -t CB
written by akh22
-
Set bam read group from existing tag
written by jbalberge
-
samtools fixmate: Couldn't write to output file: Broken pipe
written by Ric
-
How to format input data for Picard CollectRnaSeqMetrics
written by Lina F
-
Split SAM file into smaller SAM files by line count
written by cacampbell •
-
How to remove reads with hard/soft clipping along with its mate?
written by komal.rathi
Just double checking, but for the first approach, did you ensure your samtools is <=1.17?
Yes. took me sometime.... I deleted 1.19, suppose I unzipped 1.17 (wherein I couldn't find sourceme.bash). Eventually I figured out I was using the 1.10 version originally by typing just "samtools", I then changed to use the CellRanger built in samtools, which is 1.12.
The results were the same, blank filtered bam files.
This code is weird:
It should just be:
Also we don't need to see the file names. You should use placeholders like I have above. It makes it easier to debug the code.
Are you sure the
LC_ALL=Cis not required or is that change defeating the point of adding that modifier in the first place?Unless the tools are outputting some weird character formats, the answer is no don't set
LC_ALL. See here for more details: https://unix.stackexchange.com/questions/87745/what-does-lc-all-c-do I think OP copy pasted some code from somewhere without really understanding whyI think that might be relevant here (users on the forum don't go out of their way to add that setting unless something weird happened without that setting). But you could be right - maybe the OP on the original thread with this setting had some special requirement. I'll see if I can find that thread.
EDIT: I could not find that thread. I'd recommend trying it both with and without the
LC_ALL=C(maybe a limitedgrep -m100) to see if anything changes.Changed that, thank you very much and sorry for the inconvenience.