Return to list

View alignment of MAF (12-files) vs LASTZ top hits (20-flies)

View all possible LASTZ hits (20-flies)

dme-mir-287

chr2L:17574610-17574702 +

Known_mirbase

[View on miRBase] [View on UCSC Genome Browser {Cornell Mirror}]

intergenic


Legend:

   mature strand at position 17574674-17574694
   mutations
   read mismatches

Notes: droBia1 folds better without the flanking regions


20-way Fast Statistical Aligner (FSA) re-alignment

Visually pleasing best hits from LASTZ and MAF


Species Coordinate Source Confidence Alignment
dm3 chr2L:17574595-17574717 + 1 AA-GCC--------CAC---------TCAC-----------------AGAGGACGCCGGG--GATGTATGGGTGTGTAGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACG-ACTGTGAGCAA--AACAC---G-TCGGA
droSim1 chr2L:17277816-17277938 + LASTZ 5 AA-GCT--------CAC---------TCAC-----------------AGAGGACGCCGTG--GATGTATGGGTGTGTAGGGTCGGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACG-ACTGTGAGCAA--AACAC---G-TCGGA
droSec1 super_7:1215433-1215555 + LASTZ 5 AA-GCT--------CAC---------TCAC-----------------AGAGGACGCCGTG--GATGTATGGGTGTGTAGGGTCGGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACG-ACTGTGAGCAA--AACAC---G-TCGGA
droYak2 chr2R:3525115-3525236 - LASTZ 5 AC-A--------T----ATATACAC---AC-----------------AGAGCAAGCCGTG--GATGTATGGGTGT---GGGTCTGAAATTTTGCACACATTTACAGTAATTGTAAATGTGTTGAAAATCGTTTGCATG-ACTGTGAGCAA--AACAC---T-TCGGA
droEre2 scaffold_4845:6676077-6676186 - LASTZ 5 GT-GTG-------------------------------------------GA--CGC--TG--GATGTATGGGTGTGTGGGCTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACG-ACTGCGAGCAA--AACAC---G-TCGGA
droEug1 scf7180000409016:198057-198168 + LASTZ 5 AATA--------------------------------A----------AGTG--CTC--TA--GATGTATGGGTGTATAGGGTCTGAAATTTTGCACACATTTACAATTATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AGCAC---G-TCGGT
droBia1 scf7180000302422:6132274-6132381 - LASTZ 5 CG----------------------------------------------GGG--CGC--TG--GTTGTGTGGGTGTGTAGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGA
droTak1 scf7180000415617:9738-9849 - LASTZ 5 CA-C----------------------T--------------------GGGG--CGC--TG--GTTGTATGGGTGTGTGGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGA
droEle1 scf7180000491046:278921-279029 + LASTZ 5 GC-A--------C----------------------------------GGGG--CAC--TG--GGTGTATGGGT--GTAAGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGA
droRho1 scf7180000777709:31207-31296 + LASTZ 5 G---------------------------------------------------------------------GTT--GTAAGGTCTGAAATTTTGCACACATTTACAATTATTGTAAATGTGTTGAAAATCGCTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGA
droFic1 scf7180000453855:358907-359016 + LASTZ 5 AC-A--------C----------------------------------TGGG--CGC--TG--GGTGTATGGGT--GTAAGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCCGA
droKik1 scf7180000302246:110146-110268 - LASTZ 5 AG-GAC--------GAC---------GACA-ACG--A----------AGGC--CAG--TGCGTTTGTGTGGGT--GTAGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--GACAT---G-TCGGA
droAna3 scaffold_12916:3302442-3302546 - LASTZ 5 AT-A--------------------------------------------GG----GC--TG--TGTGTATTGGT--AGTGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AGCAC---G-TCGGA
droBip1 scf7180000396737:1470812-1470916 + LASTZ 5 AT-A--------------------------------------------GG----GC--TG--TGTGTATTGGT--AGTGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AGCAC---G-TCGGA
dp4 chr4_group3:10858860-10858990 - LASTZ 5 AC-A--------C----ACACACATA----------AATATATTCACGGAGACAGTC--T--TATGTGTGAGT--GTGGGGCCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGA
droPer1 super_8:2069946-2070076 - LASTZ 5 AC-A--------C----ACACACATA----------AATATATTCACGGAG--ACC--TGCTTATGTGTGAGT--GTGGGGCCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGA
droWil1 scaffold_180764:928810-928950 - LASTZ 5 AC-T--------CCCAC---------ATACATTGGAAAGAGCTTCTCCGTAATATG--TA--TATGTATTGCT--GTAAAGTCTGAAATTTTGCACACATTTACAATTATTGTAAATGTGTTGAAAATCGTTTGCACAAACTGTGAGCAA--AATCCTGTAGCTGGG
droVir3 scaffold_12963:16502531-16502646 + LASTZ 5 AG-ATT--------GAC---------TCCC-----------------ACGCTTTCGCACA--CACACAAACGTGCATAATGTTTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA-------------CCGA
droMoj3 scaffold_6500:21428983-21429089 + LASTZ 5 AC-T--------CCCAC---------GCAC---------------------------------ACACACACGCGCATAATGTTTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--CCAAC---A---GCA
droGri2 scaffold_15252:13917089-13917209 + LASTZ 5 CA-GTTTTCACTCCCAC---------ACA-----------------------------CT--CACACACACATGCATAATGCTTGAAATTTTGCACACATTTACAATATTTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAAGCAACAC---A-CATGT

Legend

0- The ortholog is absent. A better synthenic sequence candidate could not be found.
1- Reads w/ good pileup on mature and star strand & folds into hairpin structure.
2- Reads w/ good pileup on mature and star strand BUT does not fold into hairpin structure.
3- Sporadic reads & folds into a hairpin structure.
4- No reads but folds into a hairpin structure.
5- Ortholog is confident but miRNA appears absent (no reads and does not fold into a hairpin)

Re-alignment after removal of non-confident orthologs

Species Alignment
dm3 AA-GCC--------CAC---------TCAC-----------------AGAGGACGCCGGG--GATGTATGGGTGTGTAGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACG-ACTGTGAGCAA--AACAC---G-TCGGA
droSim1 AA-GCT--------CAC---------TCAC-----------------AGAGGACGCCGTG--GATGTATGGGTGTGTAGGGTCGGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACG-ACTGTGAGCAA--AACAC---G-TCGGA
droSec1 AA-GCT--------CAC---------TCAC-----------------AGAGGACGCCGTG--GATGTATGGGTGTGTAGGGTCGGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACG-ACTGTGAGCAA--AACAC---G-TCGGA
droYak2 AC-A--------T----ATATACAC---AC-----------------AGAGCAAGCCGTG--GATGTATGGGTGT---GGGTCTGAAATTTTGCACACATTTACAGTAATTGTAAATGTGTTGAAAATCGTTTGCATG-ACTGTGAGCAA--AACAC---T-TCGGA
droEre2 GT-GTG-------------------------------------------GA--CGC--TG--GATGTATGGGTGTGTGGGCTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACG-ACTGCGAGCAA--AACAC---G-TCGGA
droEug1 AATA--------------------------------A----------AGTG--CTC--TA--GATGTATGGGTGTATAGGGTCTGAAATTTTGCACACATTTACAATTATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AGCAC---G-TCGGT
droBia1 CG----------------------------------------------GGG--CGC--TG--GTTGTGTGGGTGTGTAGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGA
droTak1 CA-C----------------------T--------------------GGGG--CGC--TG--GTTGTATGGGTGTGTGGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGA
droEle1 GC-A--------C----------------------------------GGGG--CAC--TG--GGTGTATGGGT--GTAAGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGA
droRho1 G---------------------------------------------------------------------GTT--GTAAGGTCTGAAATTTTGCACACATTTACAATTATTGTAAATGTGTTGAAAATCGCTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGA
droFic1 AC-A--------C----------------------------------TGGG--CGC--TG--GGTGTATGGGT--GTAAGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCCGA
droKik1 AG-GAC--------GAC---------GACA-ACG--A----------AGGC--CAG--TGCGTTTGTGTGGGT--GTAGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--GACAT---G-TCGGA
droAna3 AT-A--------------------------------------------GG----GC--TG--TGTGTATTGGT--AGTGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AGCAC---G-TCGGA
droBip1 AT-A--------------------------------------------GG----GC--TG--TGTGTATTGGT--AGTGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AGCAC---G-TCGGA
dp4 AC-A--------C----ACACACATA----------AATATATTCACGGAGACAGTC--T--TATGTGTGAGT--GTGGGGCCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGA
droPer1 AC-A--------C----ACACACATA----------AATATATTCACGGAG--ACC--TGCTTATGTGTGAGT--GTGGGGCCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGA
droWil1 AC-T--------CCCAC---------ATACATTGGAAAGAGCTTCTCCGTAATATG--TA--TATGTATTGCT--GTAAAGTCTGAAATTTTGCACACATTTACAATTATTGTAAATGTGTTGAAAATCGTTTGCACAAACTGTGAGCAA--AATCCTGTAGCTGGG
droVir3 AG-ATT--------GAC---------TCCC-----------------ACGCTTTCGCACA--CACACAAACGTGCATAATGTTTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA-------------CCGA
droMoj3 AC-T--------CCCAC---------GCAC---------------------------------ACACACACGCGCATAATGTTTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--CCAAC---A---GCA
droGri2 CA-GTTTTCACTCCCAC---------ACA-----------------------------CT--CACACACACATGCATAATGCTTGAAATTTTGCACACATTTACAATATTTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAAGCAACAC---A-CATGT

20-way Secondary Structure

dm3 droSim1 droSec1 droYak2 droEre2 droEug1 droBia1 droTak1 droEle1 droRho1 droFic1 droKik1 droAna3 droBip1 dp4 droPer1 droWil1 droVir3 droMoj3 droGri2
RNAfold Image
UNAFold Image

Small RNA-seq read pileup

Species Alignment
droSim1
AA-GCT--------CAC---------TCAC-----------------AGAGGACGCCGTG--GATGTATGGGTGTGTAGGGTCGGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACG-ACTGTGAGCAA--AACAC---G-TCGGASizeHit CountTotal NormV044
Embryo
M023
Head
M024
Male body
M025
Embryo
.......................................................................................................................TTTGAAAATGGTTTGCAC..............................1811.000.000.001.000.00
................................................................TGTATGTGTGTGTAGGTTC....................................................................................1911.000.000.001.000.00
.........................................................................................................................TGAAAATGGTTTGCACG-C...........................1811.001.000.000.000.00
.......................................................................................................CAAAAATTGTACATGTGT..............................................1820.500.500.000.000.00
droSec1
AA-GCT--------CAC---------TCAC-----------------AGAGGACGCCGTG--GATGTATGGGTGTGTAGGGTCGGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACG-ACTGTGAGCAA--AACAC---G-TCGGASizeHit CountTotal NormV113
Male body
V114
Embryo
V115
Head
.......................................................................................................CAAAAATTGTACATGTGT..............................................1821.000.001.000.00
.........................................................GTG--TATGTGTGGGTGTGT..........................................................................................1811.000.001.000.00
.........................................................................................................................TGAAAATGGTTTGCACG-C...........................1811.000.001.000.00
droYak2
AC-A--------T----ATATACAC---AC-----------------AGAGCAAGCCGTG--GATGTATGGGTGT---GGGTCTGAAATTTTGCACACATTTACAGTAATTGTAAATGTGTTGAAAATCGTTTGCATG-ACTGTGAGCAA--AACAC---T-TCGGASizeHit CountTotal NormV046
Embryo
V052
Head
V120
Male body
V058
Head
M026
Head
..........................................................TG--GGTGGATGGGTGT---GGG......................................................................................1811.000.000.001.000.000.00
.........................................................GTG--GGTGTGTGGGTGT---GGGT.....................................................................................2030.330.000.000.330.000.00
droEre2
GT-GTG-------------------------------------------GA--CGC--TG--GATGTATGGGTGTGTGGGCTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACG-ACTGCGAGCAA--AACAC---G-TCGGASizeHit CountTotal NormV040
Embryo
V060
Head
V047
Embryo
V108
Head
.................................................................................................TATTTACAACAATTGTAA....................................................1811.000.000.000.001.00
.........................................................................................................AAAATTGTACATGTGTTGA...........................................1911.001.000.000.000.00
.................................................GA--CGC--xG--GATGGATGGGTG.............................................................................................1811.001.000.000.000.00
.......................................................................................................CAAAAATTGTACATGTGT..............................................1811.001.000.000.000.00
droEug1
AATA--------------------------------A----------AGTG--CTC--TA--GATGTATGGGTGTATAGGGTCTGAAATTTTGCACACATTTACAATTATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AGCAC---G-TCGGTSizeHit CountTotal Norm
droBia1
CG----------------------------------------------GGG--CGC--TG--GTTGTGTGGGTGTGTAGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGASizeHit CountTotal Norm
droTak1
CA-C----------------------T--------------------GGGG--CGC--TG--GTTGTATGGGTGTGTGGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGASizeHit CountTotal Norm
droEle1
GC-A--------C----------------------------------GGGG--CAC--TG--GGTGTATGGGT--GTAAGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGASizeHit CountTotal Norm
droRho1
G---------------------------------------------------------------------GTT--GTAAGGTCTGAAATTTTGCACACATTTACAATTATTGTAAATGTGTTGAAAATCGCTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGASizeHit CountTotal Norm
droFic1
AC-A--------C----------------------------------TGGG--CGC--TG--GGTGTATGGGT--GTAAGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCCGASizeHit CountTotal Norm
droKik1
AG-GAC--------GAC---------GACA-ACG--A----------AGGC--CAG--TGCGTTTGTGTGGGT--GTAGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--GACAT---G-TCGGASizeHit CountTotal Norm
droAna3
AT-A--------------------------------------------GG----GC--TG--TGTGTATTGGT--AGTGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AGCAC---G-TCGGASizeHit CountTotal NormV105
Male body
V106
Head
V055
Unknown
V039
Embryo
..............................................................................................................TGAAAATGCGTTGAAAAT.......................................1820.500.000.000.000.50
.......................................................................................................................GTTGAAAATAGTATGCACA-............................1920.500.500.000.000.00
...........................................................................................................AAGTGTAAATxTGTTGAAAA........................................1920.500.000.500.000.00
droBip1
AT-A--------------------------------------------GG----GC--TG--TGTGTATTGGT--AGTGGGTCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AGCAC---G-TCGGASizeHit CountTotal Norm
dp4
AC-A--------C----ACACACATA----------AATATATTCACGGAGACAGTC--T--TATGTGTGAGT--GTGGGGCCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGASizeHit CountTotal NormV112
Male body
V043
Embryo
V051
Head
M022
Male body
droPer1
AC-A--------C----ACACACATA----------AATATATTCACGGAG--ACC--TGCTTATGTGTGAGT--GTGGGGCCTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--AACAC---G-TCGGASizeHit CountTotal NormV111
Male body
V042
Embryo
V050
Head
V057
Head
M021
Embryo
.......................................................C--TGCTTGTGTGTGTGT--GTG.........................................................................................1980.250.000.000.000.000.25
.......................................................C--TGTGTATGTGTGAGT--GT..........................................................................................1870.140.140.000.000.000.00
.............................................................................................................GTGTAAATxTGTTGAAAAT.......................................1830.330.000.000.000.000.33
droWil1
AC-T--------CCCAC---------ATACATTGGAAAGAGCTTCTCCGTAATATG--TA--TATGTATTGCT--GTAAAGTCTGAAATTTTGCACACATTTACAATTATTGTAAATGTGTTGAAAATCGTTTGCACAAACTGTGAGCAA--AATCCTGTAGCTGGGSizeHit CountTotal NormV117
Male body
V118
Embryo
V119
Head
M020
Head
....................................................................................................TTACAATTxTTGGAAATGTG...............................................1911.000.001.000.000.00
...........................................................................................................GATTGTAAATGTGTGGAA..........................................1820.500.000.500.000.00
...........................................................................................................................ATAATCGTTTGCACAAAA..........................1811.000.001.000.000.00
.................................................................................TCTGATATTTTGCACACC....................................................................1811.000.000.000.001.00
..........................................................................................................................AATAATCGTTTGCACAAA...........................1811.000.000.000.001.00
................................................................................................................AAAATGTGTTGATAATCGT....................................1911.001.000.000.000.00
............................................................................................................ATTGTAAATGTATTGAAT.........................................1811.000.001.000.000.00
droVir3
AG-ATT--------GAC---------TCCC-----------------ACGCTTTCGCACA--CACACAAACGTGCATAATGTTTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA-------------CCGASizeHit CountTotal NormV116
Male body
V053
Unknown
V047
Embryo
M027
Male body
M028
Head
........................................................CACA--CACACAACCGGGCA...........................................................................................1820.500.000.000.000.000.50
droMoj3
AC-T--------CCCAC---------GCAC---------------------------------ACACACACGCGCATAATGTTTGAAATTTTGCACACATTTACAATAATTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAA--CCAAC---A---GCASizeHit CountTotal NormV110
Male body
V056
Unknown
V049
Head
V041
Embryo
............................................................................................................................................................C---A---ACA2030.330.000.000.330.00
...............AT---------ACAC---------------------------------ACACACACGCGC............................................................................................18290.030.030.000.000.00
.....................................................................................................................GTGTTGAAAATxGGTTGCAC..............................1920.500.000.500.000.00
droGri2
CA-GTTTTCACTCCCAC---------ACA-----------------------------CT--CACACACACATGCATAATGCTTGAAATTTTGCACACATTTACAATATTTGTAAATGTGTTGAAAATCGTTTGCACA-ACTGTGAGCAAGCAACAC---A-CATGTSizeHit CountTotal NormV109
Male body

Generated: 04/20/2012 at 05:20:55 EDT