how to count the number of reads align per viruses in my sample
2
0
Entering edit mode
3 months ago
viviane20 • 0

Hello,

I would like to count the number of reads aligning per viruses in my sample based on the output of minimap2 in a bam file.

Does anyone how to proceed?

By using the command: samtools view -f 2 ech22ReadsViralSorted.bam, I am able to see the alignments, but I am unsure about how to process the file.

Attached is the link to download the sample ech22ReadsViralSorted.bam

Thank you!

RNAseq samtools illumina • 593 views
ADD COMMENT
0
Entering edit mode

All virus sequences are in a single file, which means that in the minimap2 output, I need to know the aligned viruses and the number of reads aligned for each virus

ADD REPLY
0
Entering edit mode

All virus sequences are in a single file

If the reference is a multi-fasta file in following format

> virus_1
Sequence
> virus_2
sequence
....
>virus_N
sequence

then run the command I posted below. It will show you all reference sequence names and reads aligned to each.

ADD REPLY
4
Entering edit mode
3 months ago
GenoMax 139k

count the number of reads aligning per viruses

Assuming the virii were in the multi-fasta reference as independent entries you should use

samtools idxstats ech22ReadsViralSorted.bam
ADD COMMENT
0
Entering edit mode

yes the command line samtools idxstats ech22ReadsViralSorted.bam worked. Thanks you

ADD REPLY
0
Entering edit mode

Please accept the answer (green checkmark) to provide closure to this thread.

ADD REPLY
0
Entering edit mode
3 months ago
$ samtools view

(...)
  -c, --count                Print only the count of matching records
(....)
ADD COMMENT

Login before adding your answer.

Traffic: 1088 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