Return to list

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

View all possible LASTZ hits (20-flies)

dme-mir-4983

chr2R:15524918-15525026 +

Known_mirbase

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

intergenic


Legend:

   mature strand at position 15524986-15525007
   star strand at position 15524936-15524957
   mutations
   read mismatches

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

Visually pleasing best hits from LASTZ and MAF


Species Coordinate Source Confidence Alignment
dm3 chr2R:15524903-15525041 + 1 A--ACA-CATCTGCC---CGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-CG--TTTGCATTT-------CTG------AAAA---ATGTATACGCAA-AAA-CCAAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TAC
droSim1 chr2R:14228631-14228769 + LASTZ 4 A--ACA-CATCTGCC---CGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-CG--TTTGCATTT-------CTG------AAAA---ATGTATACGCAA-AAA-CCAAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TAC
droSec1 super_1:13039564-13039702 + LASTZ 1 A--ACA-CATCTGCC---CGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTG------AAAA---ATGTATACGCAA-AAA-CCAAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TAC
droYak2 chr2R:13649050-13649195 - LASTZ 5 A--ACA-CATCTGCC---CGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-CG--TTTGCATTTGTATTTTCTG------AAAA---ATGTATACGCAA-AAA-CCAAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TAC
droEre2 scaffold_4845:9723500-9723642 + LASTZ 5 A--ACA-CATCTGCC---CGCTGG-CTGGCTG----------GCT----GT-------CTTGACCCACTTGCT-CG--CTTGCATTG-------CTG------AAAA---ATGTATACGCAA-AAA-CCAAAATCA------GAAATTCCTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TAC
droEug1 scf7180000409474:2425232-2425371 - LASTZ 5 A--ACA-CATCTGCC--CCGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTG------AAAA---ATGTATACGCAA-AAA-CCAAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TAC
droBia1 scf7180000302143:3024933-3025070 + LASTZ 5 A--ACA-CATCTGCC---CGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTG------AAAA---ATGTATACGCAAA--A-CCGAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TAC
droTak1 scf7180000415384:129223-129366 - LASTZ 5 A--ACA-CATCTGCC---CGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTGAATCTGAAAA---ATGTATACGCAAA--A-CCAAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TAC
droEle1 scf7180000491201:1980229-1980378 + LASTZ 5 A--ACA-CATCTGCC---CGCTGG-CTGGCTGGCTGGTTG--GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTG------AAAT---ATGTATACGCAAA--A-CCAAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TAC
droRho1 scf7180000780066:106697-106838 - LASTZ 5 A--ACA-CATCTGCC---CGCTGG-CTGGCTG----------GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTG------AAAT---ATGTATACGCAAA--A-CCAAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TAC
droFic1 scf7180000454039:817245-817390 + LASTZ 5 A--ACA-CATCTGCC---CGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-CGCATTTGCATTT-------CTCAAAAAAAAAA---ATGTATACGCAAA--A-CCAAAATCC------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TAC
droKik1 scf7180000302679:61967-62124 + LASTZ 5 A--ACAGCATCTGTGTGCTGCTGG-CTGG-----TGGCTGGCGCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTTCAAAAA--TG------AAAA---ATGTATACGCAA-AAACCCGAAAACA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TAC
droAna3 scaffold_13266:9754397-9754555 - LASTZ 5 AAAACA-CAAACGCC---CGCTG-------------------GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTCG-----CAAAAATATTCATATGCAAG--A-CCAAAATCAGAATCAGAAATTCTTTTGG-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCTCTATATGTA----CATCTAGTAC
droBip1 scf7180000396427:901727-901882 - LASTZ 5 A--CCA-CAAACGCC---CGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTG------AAAAAATATTTATATGCAAG--A-CCAAAATCA------GAAATTCTTTTGG-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATATGTACATACATCTAGTAC
dp4 chr3:12652733-12652881 + LASTZ 0 A--ACA-CGTCTCTC---CGTTTGTCAGTCTTTCAG------TCT----GTCAGCATGCTTGACCCACTTGATGTGCTTTTGAATTT-------TAG------AAAT---ATGTA---------AA-TCAAATTCT------GAAACATTTCGGGGACAAGTTTATGTGTCATTAAATTCCCCA----ATTTTATACT--AAA------------------
droPer1 super_2:8929657-8929759 - LASTZ 0 A--ACA-CTTCTCTC---CGTTTGTCAGTCTTTCAG------TCT----GTCAGCATGCTTGACCCACTTGATGTGCTTTTGAATTT-------TAG------AAAT---ATGTA---------AA-TCAAATTCT------GAAACATTT----------------------------------------------------------------------
droWil1 Unknown 0 -----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
droVir3 scaffold_12875:4153222-4153295 - LASTZ 0 A--ACA-CTTTTGAG---TGCCAG-CTT--------------TCTGAGCAT-------ATCAAA-CGCTCGA--CA--TTTGTCTTT-------TTG------AGGG---GGGTTCTCGCAC-A-------------------------------------------------------------------------------------------------
droMoj3 scaffold_6496:20896701-20896727 + LASTZ 0 A--ACA-CTTTTGAC---TGCTGC-CAG--------------CCT----GC-------C------------------------------------------------------------------------------------------------------------------------------------------------------------------
droGri2 Unknown 0 -----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

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 A--ACA-CATCTGCC---CGCTGGCT----G-----------GCTGTCTTGACCCACTTGCTCG--TTTGCATTT-------CTG------AAA--AATGTATACGCAAAAA-CCAAAA------TCAGAAATTCTTTTGAGCAAGTGTGTCATGAAATTGCCATTCGTTTCTGCTCT--ATA---------------TAC
droSim1 A--ACA-CATCTGCC---CGCTGGCT----G-----------GCTGTCTTGACCCACTTGCTCG--TTTGCATTT-------CTG------AAA--AATGTATACGCAAAAA-CCAAAA------TCAGAAATTCTTTTGAGCAAGTGTGTCATGAAATTGCCATTCGTTTCTGCTCT--ATA---------------TAC
droSec1 A--ACA-CATCTGCC---CGCTGGCT----G-----------GCTGTCTTGACCCACTTGCTTG--TTTGCATTT-------CTG------AAA--AATGTATACGCAAAAA-CCAAAA------TCAGAAATTCTTTTGAGCAAGTGTGTCATGAAATTGCCATTCGTTTCTGCTCT--ATA---------------TAC
droYak2 A--ACA-CATCTGCC---CGCTGGCT----G-----------GCTGTCTTGACCCACTTGCTCG--TTTGCATTTGTATTTTCTG------AAA--AATGTATACGCAAAAA-CCAAAA------TCAGAAATTCTTTTGAGCAAGTGTGTCATGAAATTGCCATTCGTTTCTGCTCT--ATA---------------TAC
droEre2 A--ACA-CATCTGCC---CGCTGGCTGGCTG-----------GCTGTCTTGACCCACTTGCTCG--CTTGCATTG-------CTG------AAA--AATGTATACGCAAAAA-CCAAAA------TCAGAAATTCCTTTGAGCAAGTGTGTCATGAAATTGCCATTCGTTTCTGCTCT--ATA---------------TAC
droEug1 A--ACA-CATCTGCC--CCGCTGGCT----G-----------GCTGTCTTGACCCACTTGCTTG--TTTGCATTT-------CTG------AAA--AATGTATACGCAAAAA-CCAAAA------TCAGAAATTCTTTTGAGCAAGTGTGTCATGAAATTGCCATTCGTTTCTGCTCT--ATA---------------TAC
droBia1 A--ACA-CATCTGCC---CGCTGGCT----G-----------GCTGTCTTGACCCACTTGCTTG--TTTGCATTT-------CTG------AAA--AATGTATACGCAA-AA-CCGAAA------TCAGAAATTCTTTTGAGCAAGTGTGTCATGAAATTGCCATTCGTTTCTGCTCT--ATA---------------TAC
droTak1 A--ACA-CATCTGCC---CGCTGGCT----G-----------GCTGTCTTGACCCACTTGCTTG--TTTGCATTT-------CTGAATCTGAAA--AATGTATACGCAA-AA-CCAAAA------TCAGAAATTCTTTTGAGCAAGTGTGTCATGAAATTGCCATTCGTTTCTGCTCT--ATA---------------TAC
droEle1 A--ACA-CATCTGCC---CGCTGGCTGGCTGG---CTGGTTGGCTGTCTTGACCCACTTGCTTG--TTTGCATTT-------CTG------AAA--TATGTATACGCAA-AA-CCAAAA------TCAGAAATTCTTTTGAGCAAGTGTGTCATGAAATTGCCATTCGTTTCTGCTCT--ATA---------------TAC
droRho1 A--ACA-CATCTGCC---CGCTGGCTGGCTG-----------GCTGTCTTGACCCACTTGCTTG--TTTGCATTT-------CTG------AAA--TATGTATACGCAA-AA-CCAAAA------TCAGAAATTCTTTTGAGCAAGTGTGTCATGAAATTGCCATTCGTTTCTGCTCT--ATA---------------TAC
droFic1 A--ACA-CATCTGCC---CGCTGGCT----G-----------GCTGTCTTGACCCACTTGCTCGCATTTGCATTT-------CTCAAAAAAAAA--AATGTATACGCAA-AA-CCAAAA------TCCGAAATTCTTTTGAGCAAGTGTGTCATGAAATTGCCATTCGTTTCTGCTCT--ATA---------------TAC
droKik1 A--ACAGCATCTGTGTGCTGCTGGCT----GGTGGCTGGC--GCTGTCTTGACCCACTTGCTTG--TTTGCATTTCAAAAA--TG------AAA--AATGTATACGCAAAAACCCGAAA------ACAGAAATTCTTTTGAGCAAGTGTGTCATGAAATTGCCATTCGTTTCTGCTCT--ATA---------------TAC
droAna3 AAAACA-CAAACGCC---CGCTG-------------------GCTGTCTTGACCCACTTGCTTG--TTTGCATTT-------CTCG----CAAAAATATTCATATGCAA-GA-CCAAAATCAGAATCAGAAATTCTTTTGGGCAAGTGTGTCATGAAATTGCCATTCGTTTCTGCTCTCTATATGTA----CATCTAGTAC
droBip1 A--CCA-CAAACGCC---CGCTGGCT----G-----------GCTGTCTTGACCCACTTGCTTG--TTTGCATTT-------CTGAA-----AAAATATTTATATGCAA-GA-CCAAAA------TCAGAAATTCTTTTGGGCAAGTGTGTCATGAAATTGCCATTCGTTTCTGCTCT--ATATGTACATACATCTAGTAC

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 None None
UNAFold Image None None

Small RNA-seq read pileup

Species Alignment
droSim1
A--ACA-CATCTGCC---CGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-CG--TTTGCATTT-------CTG------AAAA---ATGTATACGCAA-AAA-CCAAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TACSizeHit CountTotal NormV044
Embryo
M023
Head
M024
Male body
M025
Embryo
.......................................................................................................................CAA-AAA-CGAAAACCA------GAA............................................................................1811.001.000.000.000.00
droSec1
A--ACA-CATCTGCC---CGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTG------AAAA---ATGTATACGCAA-AAA-CCAAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TACSizeHit CountTotal NormV113
Male body
V114
Embryo
V115
Head
..............................................................................................................................................GAAATTCTTTTGA-GCAAGT----GTG....................................................2212.001.000.001.00
...............................................................................................................................................AAATTCTTTTGA-GCAAGT----GTG....................................................2111.000.000.001.00
..............................................................................................................................................GAAATTCTTTTGA-GCAAGT----GTGT...................................................2311.001.000.000.00
..............................................................................................................................................GAAATTCTTTTGA-GCAAGT----GT.....................................................2111.000.000.001.00
...........................................................................................................................AAA-CCAAAACCA------GAAGTT.........................................................................1811.000.001.000.00
...............................................................................................................................................AAATTCTTTTGA-GCAAGT----GT.....................................................2011.001.000.000.00
................................................................................................................................................AATTCTTTTGA-GCAAGT----GTGTCA.................................................2311.000.000.001.00
..........................................................................................................................................................................GATGAAATTGCC-ATTGGT................................1811.000.001.000.00
..................................................................................................................................................TTCTTTTGA-GCAAGT----GTGTCA.................................................2111.001.000.000.00
..............................................................................................................................................GAAATCCTTTTGA-GCAAGT----GT.....................................................2111.000.000.001.00
...............................................................................................................................................AAATTCTTTTGA-GCAACT----GTG....................................................2111.001.000.000.00
droYak2
A--ACA-CATCTGCC---CGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-CG--TTTGCATTTGTATTTTCTG------AAAA---ATGTATACGCAA-AAA-CCAAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TACSizeHit CountTotal NormV046
Embryo
V052
Head
V120
Male body
V058
Head
M026
Head
......................................................................................................................................AA------GAAATTCCTTTGA-GCA..............................................................1820.500.000.000.000.500.00
droEre2
A--ACA-CATCTGCC---CGCTGG-CTGGCTG----------GCT----GT-------CTTGACCCACTTGCT-CG--CTTGCATTG-------CTG------AAAA---ATGTATACGCAA-AAA-CCAAAATCA------GAAATTCCTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TACSizeHit CountTotal NormV040
Embryo
V060
Head
V047
Embryo
V108
Head
..................................................................ACTTGCT-CG--CTCGGATTG-------...............................................................................................................................1811.000.000.000.001.00
droEug1
A--ACA-CATCTGCC--CCGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTG------AAAA---ATGTATACGCAA-AAA-CCAAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TACSizeHit CountTotal Norm
droBia1
A--ACA-CATCTGCC---CGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTG------AAAA---ATGTATACGCAAA--A-CCGAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TACSizeHit CountTotal Norm
droTak1
A--ACA-CATCTGCC---CGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTGAATCTGAAAA---ATGTATACGCAAA--A-CCAAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TACSizeHit CountTotal Norm
droEle1
A--ACA-CATCTGCC---CGCTGG-CTGGCTGGCTGGTTG--GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTG------AAAT---ATGTATACGCAAA--A-CCAAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TACSizeHit CountTotal Norm
droRho1
A--ACA-CATCTGCC---CGCTGG-CTGGCTG----------GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTG------AAAT---ATGTATACGCAAA--A-CCAAAATCA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TACSizeHit CountTotal Norm
droFic1
A--ACA-CATCTGCC---CGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-CGCATTTGCATTT-------CTCAAAAAAAAAA---ATGTATACGCAAA--A-CCAAAATCC------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TACSizeHit CountTotal Norm
droKik1
A--ACAGCATCTGTGTGCTGCTGG-CTGG-----TGGCTGGCGCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTTCAAAAA--TG------AAAA---ATGTATACGCAA-AAACCCGAAAACA------GAAATTCTTTTGA-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATA---------------TACSizeHit CountTotal Norm
droAna3
AAAACA-CAAACGCC---CGCTG-------------------GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTCG-----CAAAAATATTCATATGCAAG--A-CCAAAATCAGAATCAGAAATTCTTTTGG-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCTCTATATGTA----CATCTAGTACSizeHit CountTotal NormV105
Male body
V106
Head
V055
Unknown
V039
Embryo
droBip1
A--CCA-CAAACGCC---CGCTGG-CTG--------------GCT----GT-------CTTGACCCACTTGCT-TG--TTTGCATTT-------CTG------AAAAAATATTTATATGCAAG--A-CCAAAATCA------GAAATTCTTTTGG-GCAAGT----GTGTCATGAAATTGCC-ATTCGTTTCTGCTCT--ATATGTACATACATCTAGTACSizeHit CountTotal Norm
dp4
A--ACA-CGTCTCTC---CGTTTGTCAGTCTTTCAG------TCT----GTCAGCATGCTTGACCCACTTGATGTGCTTTTGAATTT-------TAG------AAAT---ATGTA---------AA-TCAAATTCT------GAAACATTTCGGGGACAAGTTTATGTGTCATTAAATTCCCCA----ATTTTATACT--AAA------------------SizeHit CountTotal NormV112
Male body
V043
Embryo
V051
Head
M022
Male body
droPer1
A--ACA-CTTCTCTC---CGTTTGTCAGTCTTTCAG------TCT----GTCAGCATGCTTGACCCACTTGATGTGCTTTTGAATTT-------TAG------AAAT---ATGTA---------AA-TCAAATTCT------GAAACATTT----------------------------------------------------------------------SizeHit CountTotal NormV111
Male body
V042
Embryo
V050
Head
V057
Head
M021
Embryo
droWil1
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------SizeHit CountTotal NormV117
Male body
V118
Embryo
V119
Head
M020
Head
droVir3
A--ACA-CTTTTGAG---TGCCAG-CTT--------------TCTGAGCAT-------ATCAAA-CGCTCGA--CA--TTTGTCTTT-------TTG------AGGG---GGGTTCTCGCAC-A-------------------------------------------------------------------------------------------------SizeHit CountTotal NormV116
Male body
V053
Unknown
V047
Embryo
M027
Male body
M028
Head
droMoj3
A--ACA-CTTTTGAC---TGCTGC-CAG--------------CCT----GC-------C------------------------------------------------------------------------------------------------------------------------------------------------------------------SizeHit CountTotal NormV110
Male body
V056
Unknown
V049
Head
V041
Embryo
droGri2
-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------SizeHit CountTotal NormV109
Male body

Generated: 04/20/2012 at 05:16:46 EDT