Question: Can IGV be used to visualize Mitogenome readpools
1
gravatar for F.I.2011
6 months ago by
F.I.201110
F.I.201110 wrote:

Hello everyone,

I am struggling to use IGV to visualize my readspools against my reference Mitogenome. It is avian mitogenome. I upload new genome by downloading FASTA and GFF3 file.

For the workflow:

I convert paired end reads (R1 and R2) into bam file by this command:

java -jar ----picard.jar' FastqToSam F1=FQ-R1.fastq F2=FQ-R2.fastq O=.bam SM=GK-2 RG=null

sorted by :

SortSam \I=GK-2.bam \ O=GK-sorted.sam  \ SORT_ORDER=coordinate

and build Index by BuildBamIndex.

But when I uploaded that bam file against my own-made genome file, i had no reads and zero coverage. Here is the image of IGV after uploading bam files

Then I sorted bam with queryname SORT ORDER by this command

SortSam \I=GK-2.bam \ O=GK-sorted.bam  \ SORT_ORDER=queryname.

but then BuildBamIndex gave error and need to sort bam file with coordinate SORT ORDER not queryname and IGV also denied to upload file without index file.

Please need help me to fix it. Can I use IGV for small 17000 genome that dont have inbuilt reference genome.

Thanks

igv • 180 views
ADD COMMENTlink modified 6 months ago by b.nota5.5k • written 6 months ago by F.I.201110
0
gravatar for b.nota
6 months ago by
b.nota5.5k
Netherlands
b.nota5.5k wrote:

The FastqToSam tool makes an unaligned BAM file. For IGV you need aligned reads, so align (with e.g., bwa) your paired-end reads to your genome first, and use that bam file in IGV.

ADD COMMENTlink modified 6 months ago • written 6 months ago by b.nota5.5k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1090 users visited in the last hour