I received a text file which has following entries and I assume this is a Bowtie aligned map files and I would like to visualize the same on a local genome browser.
HWI-ST863:202:C2014ACXX:4:1101:2197:1894 1:N:0:NTTTTTCG + NM_011363 3354 TGTTGTTGTTGTTTTAAACAAAATGGAAAAGCATAA IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII 1
HWI-ST863:202:C2014ACXX:4:1101:3133:1962 1:N:0:NAATAAGA + NM_175190 729 GTGTGGTTGTCCCTTTTGTTAATAAACATATGAGCA IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII 0
HWI-ST863:202:C2014ACXX:4:1101:5688:1889 1:N:0:NTATCTGG + NM_020606 4378 CTGTCTTGATCCATTTCTTCTGCATGATTCCAGAAA IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII 0
HWI-ST863:202:C2014ACXX:4:1101:9703:1907 1:N:0:NTTAAAGG + NM_026030 397 GCTTCTGATGACTTAGATGATTTGAACTTCTTTAAT IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII 0
I would like to validate if the reads are enriched towards the 3' end of the mRNA (genes). I would also like a local genome browser like tool for this purpose. How can I see this in a genome browser like IGV?
•
written
by
k.nirmalraman
2 answers
Convert your sequences to FASTA and run them through, for example, a command-line instance of UCSC BLAT to get a PSL file. Run the PSL file through BEDOPS psl2bed to get a BED file. You can turn this into a BED-like BedGraph file pretty easily for uploading to and visualization with a, for example, local UCSC Genome Browser instance.
Log in to answer this question.
-
How to split Illumina Hiseq 2000 paired-end fastq files?
written by Yuka Takemon •
-
Java Error "Concurrent.Executionexception" While Running Trimmomatic V0.32
written by arnstrm
-
Bwa Run With Fastq File With Different Read Length
written by ttom
-
Fastx-Toolkit: Fastx_Barcode_Splitter
written by k.nirmalraman
-
Split Fastq Files With A Pattern
written by Assa Yeroslaviz
-
How To Extract Reads With Wrong Insert Size With Samtools
written by Assa Yeroslaviz
-
Manipulate Sequences In Fastq Files
written by lsvijfhuizen •
-
Working With Fastq Files In Gmap
written by Assa Yeroslaviz
-
Parsing Fastq Files
written by deepthitheresa •
-
Unmapped Reads And Sequence Name
written by Stefano Berri
I guess I should some how convert this file to
bedformat...I don't quite understand the question. But if you just want to visualize the bam file, you can load the BAM file in IGV browser. And you can load your reference genome as well and just focus on the region you are interested in.
I tried to load this file directly into IGV from file-> load file with mouse genome as reference, but I cannot see anything beyond that
Does your BAM file have an index file? I think IGV usually needs it.
You also need to have the same chromosome name between the bam and your reference sequence
I don't have a index file, but here it says, I should have a sorted file. I tried IGV tools for sorting but no succes.
File type or sorting not supportederror... which leaves me wondering, if the file is in bam format....you can use samtools for sorting the bam file.
then, index the bam file (Note: to index the bam file, the file needs to be sorted beforehand).