This is a test version of Biostars. For the public version, visit https://www.biostars.org.
is our reads needs indexing before alignment like refrence genome or not ? some one explain to me plz

im new in this matter. and read some papers and watch some tutorial videos. i am in the way of aligning my dara. but in igv my data dosent match with my reference genome even iv got 92.7 % overal alignment with hisat2. some one help me pz

rna-seq alignment
• 2,197 views • link
updated by h.mon h.mon • written by sabaghianamir70 • sabaghianamir70

but in igv my data dosent match with my reference genome

what is "my data" ?

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Saccharomyces cerevisiae. i think its some thing wrong with my samtools

• 0 views • link
written by sabaghianamir70 • sabaghianamir70

Probably you are loading wrong genome in IGV. Check reference version in IGV and it should match the reference you have used in alignment sabaghianamir70

• 0 views • link
written by cpad0112 cpad0112

no i dont think so,. in igv i chose Saccharomyces cerevisiae refrence genome. both versions are same.

• 0 views • link
written by sabaghianamir70 • sabaghianamir70

1 answer

Both reference genome in fasta format and mapped reads file in bam format need to be indexed, and the bam file needs to be sorted as well. For the fasta reference:

samtools faidx ref.fasta

For the mapped bam:

 samtools sort -o sorted.bam mapped.bam
 samtools index sorted.bam

In addition, the reference genome loaded into OGV need to be the same used for mapping.

• 0 views • link
written by h.mon h.mon

i think that is my problem , samtools wont work .im trying to install it for 2 days. its really annoying. i have the same problem as he had error installing samtools on ubuntu . but eventually is see this

libncurses5-dev : Depends: libtinfo5 (= 6.1-1ubuntu1) but 6.1-1ubuntu1.18.04 is to be installed Depends: libncurses5 (= 6.1-1ubuntu1) but 6.1-1ubuntu1.18.04 is to be installed Depends: libtinfo-dev (= 6.1-1ubuntu1) but it is not going to be installed

and the libtinfo-dev does not work.

• 0 views • link
written by sabaghianamir70 • sabaghianamir70

Log in to answer this question.

More posts like this