https://github.com/poisonalien/basecounts
A tool to extract ATGC/Indel read counts at user specific loci from BAM files
https://github.com/poisonalien/basecounts
bam basecounts bioinformatics nucleotide-frequency readcounts sequencing wgs-data
Last synced: 8 days ago
JSON representation
A tool to extract ATGC/Indel read counts at user specific loci from BAM files
- Host: GitHub
- URL: https://github.com/poisonalien/basecounts
- Owner: PoisonAlien
- Created: 2020-05-24T14:02:15.000Z (almost 5 years ago)
- Default Branch: master
- Last Pushed: 2022-02-22T09:54:02.000Z (about 3 years ago)
- Last Synced: 2025-04-01T00:41:22.657Z (about 2 months ago)
- Topics: bam, basecounts, bioinformatics, nucleotide-frequency, readcounts, sequencing, wgs-data
- Language: C
- Homepage:
- Size: 38.1 KB
- Stars: 5
- Watchers: 3
- Forks: 1
- Open Issues: 1
-
Metadata Files:
- Readme: README.md
Awesome Lists containing this project
README
## basecounts
`basecounts` is a simple tool to extract nucleotide counts from BAM files for the user specific loci. This tool can also be used to generate VAF table for known somatic varaints without the proper variant calling.
### Installation
Clone the repository and run `make`. Requires [htslib](https://github.com/samtools/htslib).
```
git clone https://github.com/PoisonAlien/basecounts
cd basecounts
make
```
```
$ basecounts
-------------------------------------------------------------------------
basecounts: A tool to extract nucleotide counts
of genomic loci from the BAM file. Version: 0.1.0
USAGE:
basecounts [OPTIONS]
e.g; basecounts target_loci.tsv Tumor.bam
OPTIONS:
-f Indexed fasta file. If provided, extracts and adds reference base to the ouput tsv
-q Min. mapping quality threshold. Default 10 [Read filter]
-F Exclude reads with FLAGS >= INT. Default 3852 [Read filter]
-o Output file basename. Default parses from basename of BAM file
ARGS:
loci file (chr\tpos)
Indexed BAM file
OUPUT:
.tsv TSV file with nucletide counts for all variants
-------------------------------------------------------------------------
```
### Usage
`basecounts` takes a 2 column tsv file with loci and the sorted/indexed BAM file as positional arguments.
First two columns of loci file should be chromsome and position.
```
$ head loci.tsv
1 6304615
1 10407830
1 46269731
1 47133825
1 50884996
1 84662414
1 152328587
1 216019313
1 226252155
```
```
#basic usage
$ basecounts loci.tsv tumor.bam
#providing fasta file will add the reference base at the loci
$ basecounts -f GRCh37.fa loci.tsv tumor.bam
#Output is a TSV file with readcounts suporting all 4 nucleotide bases and INDELS
$ head tumor.tsv
loci fa_ref A T G C Ins Del
1:6304615 A 70 0 0 19 0 0
1:10407830 C 0 0 39 50 0 0
1:46269731 A 50 3 1 21 0 0
1:47133825 T 3 37 14 4 0 0
1:50884996 A 51 0 0 15 0 0
1:84662414 G 0 2 81 0 0 0
1:152328587 A 85 3 0 0 0 0
1:216019313 T 0 53 17 0 0 0
1:226252155 G 0 41 52 0 0 0
```