Getting nucleotide frequencies per position
0
0
Entering edit mode
9 months ago

Hi,

I am very new to bioinformatics and was given the task to analyze some total RNA sequencing results. So far I have converted my .fastaq files to indexed .bam using samtools index. I found the following in literature:

# Generate nucleotide counts per position
samtools index -b ./sorted_bbmerge_norm_$sample.bam

pysamstats -f /Users/kasenriemersma/Bioinformatics/Reference/$reference.fasta -t variation -D 1000000 --format csv --output $outputdir/ntcounts_$sample.csv ./sorted_bbmerge_norm_$sample.bam

I tried running this but I received an error that I was missing the BAM file operand. Is this there any easy way to get past this? Or is there another method that would be better?

Thanks for the help!

pysamstats • 235 views
ADD COMMENT

Login before adding your answer.

Traffic: 1878 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6