This is a test version of Biostars. For the public version, visit https://www.biostars.org.
quantify how many reads detect ORF9 and ORF10

Hi there,

Is there any way to quantify how many reads detect ORF9? I have a word doc with the different parts annotated and sequences of ORF9. I don't have any bam file or other reference file.

Any kind of help regarding how I can get reads by using these information would be awesome.

Thank you so much in advance! enter image description here

reads orf9 annotated

you "don't have any bam" but you want "to count he number of reads" ?(!)

1 answer

Yes. Remove ORF9 from word and place it in a text file in fasta format (google it). Convert this file to an alignment index for your favorite aligner. Use your favorite aligner and the index you just made and align your reads to generate a BAM file. Use samtools to count the reads in your BAM file that align to your ORF9 reference.

Thank you so much seidel for your response. I want to use igv aligner (need reference and Bam fileenter image description here) and get the image (attached) regarding the alignment. I already made the reference file based on my ORF9 sequences. Could you please describe briefly how i can get the BAM file using these annotated parts?

Thanks again for your help! Khatun

enter image description here

Log in to answer this question.