Secondary structure and small RNA-seq evidence for dme-mir-1

Legend:maturestarmismatch

12 Drosophila species corrected alignment

Species Coordinate Alignment
dm3 chr2L:20487426-20487546 + AAAGACC-----------------ATCCCAT--ATTCAGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGT-TAT-ATTCAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAGACATCG----------------------------------G----AGTTGAAA
droSim1 chr2L:20052694-20052814 + AAAGACC-----------------ATCCCAT--CTTCAGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGT-TAT-ACTCAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAGACATCG----------------------------------G----AGTTGAAA
droSec1 super_18:365651-365771 + AAGGACC-----------------ATCCCAT--CTTCAGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGT-TAT-ACTCAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAGACATCG----------------------------------G----AGTTGAAA
droYak2 chr2R:6842709-6842830 + AAAGAAC-----------------ATCCCAC--CATCAGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGTATAT-ACTCAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAGTCATCG----------------------------------G----AGTTGAAA
droEre2 scaffold_4845:3954268-3954372 - AAGGACC---------------------------TCCAGCCT-GTGAGAGTTCCATGCTTCCTTGCATTCAATAGT-TAT-ACTCAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAGACATCG----------------------------------G------------
droAna3 scaffold_12916:6985813-6985928 + AAATATC-------------------------AAGAGAGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGTATAA-TATCAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCAAGACTTCA----------------------------------A----GTTTAAAT
dp4 chr4_group3:7845311-7845422 - GGAAA--------------------------TCAGAAGGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGTATAA-CATAAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCAAGACAGCG-----------------------------------------TCTACA
droPer1 super_1:9292009-9292126 - AAAGACGG-A-------------A-------TCAGAAGGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGTATAA-CATAAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCAAGACATCG----------------------------------C------CTACAG
droWil1 scaffold_180708:1845911-1846049 + ATAGATGGTAGATATATATGTAG--------CTACTTGGCCTCTTGAGAGTTCCATGCTTCCTTGCATTCAATAGTATATTAAATAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGTCCGACATCA----------------------------------TCTACAATATAAT
droVir3 scaffold_12963:13613413-13613527 - TGTGC--------------------------TTACTTGGCCT-TTGAGAGTTCCATGCTTCCTTACATTCAATAGTATTC-AAATAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAAACATCA----------------------------------C----AAAAACAA
droMoj3 scaffold_6500:24555167-24555325 + TGGGTCTT-AAA--------TTGT-------TCATTGCGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGTATTT-ATATAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAAGCATCATACCAAACTGTAAAAAAAAACACATTATGAAATTG----TGTTTAAA
droGri2 scaffold_15126:6250403-6250508 + ATTGT--------------------------TTACTTGGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGTATAT-AAATAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAAACATC------------------------------------------------

Reduced alignment with Hairpin orthologs

Species Coordinate Alignment
dm3 chr2L:20487426-20487546 + AAAGACC-----------------ATCCCAT--ATTCAGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGT-TAT-ATTCAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAGACATCG----------------------------------G----AGTTGAAA
droSim1 chr2L:20052694-20052814 + AAAGACC-----------------ATCCCAT--CTTCAGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGT-TAT-ACTCAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAGACATCG----------------------------------G----AGTTGAAA
droSec1 super_18:365651-365771 + AAGGACC-----------------ATCCCAT--CTTCAGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGT-TAT-ACTCAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAGACATCG----------------------------------G----AGTTGAAA
droYak2 chr2R:6842709-6842830 + AAAGAAC-----------------ATCCCAC--CATCAGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGTATAT-ACTCAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAGTCATCG----------------------------------G----AGTTGAAA
droEre2 scaffold_4845:3954268-3954372 - AAGGACC---------------------------TCCAGCCT-GTGAGAGTTCCATGCTTCCTTGCATTCAATAGT-TAT-ACTCAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAGACATCG----------------------------------G------------
droAna3 scaffold_12916:6985813-6985928 + AAATATC-------------------------AAGAGAGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGTATAA-TATCAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCAAGACTTCA----------------------------------A----GTTTAAAT
dp4 chr4_group3:7845311-7845422 - GGAAA--------------------------TCAGAAGGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGTATAA-CATAAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCAAGACAGCG-----------------------------------------TCTACA
droPer1 super_1:9292009-9292126 - AAAGACGG-A-------------A-------TCAGAAGGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGTATAA-CATAAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCAAGACATCG----------------------------------C------CTACAG
droWil1 scaffold_180708:1845911-1846049 + ATAGATGGTAGATATATATGTAG--------CTACTTGGCCTCTTGAGAGTTCCATGCTTCCTTGCATTCAATAGTATATTAAATAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGTCCGACATCA----------------------------------TCTACAATATAAT
droVir3 scaffold_12963:13613413-13613527 - TGTGC--------------------------TTACTTGGCCT-TTGAGAGTTCCATGCTTCCTTACATTCAATAGTATTC-AAATAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAAACATCA----------------------------------C----AAAAACAA
droMoj3 scaffold_6500:24555167-24555325 + TGGGTCTT-AAA--------TTGT-------TCATTGCGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGTATTT-ATATAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAAGCATCATACCAAACTGTAAAAAAAAACACATTATGAAATTG----TGTTTAAA
droGri2 scaffold_15126:6250403-6250508 + ATTGT--------------------------TTACTTGGCCT-TTGAGAGTTCCATGCTTCCTTGCATTCAATAGTATAT-AAATAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAAACATC------------------------------------------------

Small RNA-sequencing read pileup

Species Read pileup
droSim1
AAAGACCATCCCATCTTCA----GCCTTTGAGAGTTCCATGCTTCCTTGCATTCAATAGTTATA--CTCAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCGAGACATCGGAGTTGAAASizeHit CountTotal NormTotalGSM343915
Embryo
............................................................................TGGAATGTAAAGAAGTATGGAG.............................221383333.00383333383333
............................................................................TGGAATGTAAAGAAGTATGGA..............................2112370.0023702370
............................................................................TGGAATGTAAAGAAGTATGGAGC............................231484.00484484
.............................................................................GGAATGTAAAGAAGTATGGAG.............................211263.00263263
............................................................................TGGAATGTAAAGAAGTATGG...............................201166.00166166
..............................................................................GAATGTAAAGAAGTATGGAG.............................201119.00119119
................................................................................ATGTAAAGAAGTATGGAG.............................18113.001313
dp4
GGAAA-----------TCAGAAGGCCTTTGAGAGTTCCATGCTTCCTTGCATTCAATAGT-ATAACATAAAGCATATGGAATGTAAAGAAGTATGGAGCGAAATCTGGCAAGACA---GCGTCTACASizeHit CountTotal NormTotalGSM343916
Embryo
GSM444067
Head
............................................................................TGGAATGTAAAGAAGTATGGAG.............................221635117.00635117246218388899
............................................................................TGGAATGTAAAGAAGTATGGA..............................21119848.0019848217517673
............................................................................TGGAATGTAAAGAAGTATGG...............................2012776.0027761962580
............................................................................TGGAATGTAAAGAAGTATGGAGC............................231611.00611349262
.............................................................................GGAATGTAAAGAAGTATGGAG.............................211411.00411145266
............................................................................TGGAATGTAAAGAAGTATG................................191334.0033427307
..............................................................................GAATGTAAAGAAGTATGGAG.............................201235.0023594141
...............................................................................AATGTAAAGAAGTATGGAG.............................19135.00352015
...........................................................................ATGGAATGTAAAGAAGTATGGAG.............................23126.00261214
............................................................................TGGAATGTAAAGAAGTAT.................................18124.0024123

Secondary Structure

dm3

dG=-29.0,  p-value=0.009901 dG=-28.3,  p-value=0.009901

droSim1

dG=-29.0,  p-value=0.009901 dG=-28.3,  p-value=0.009901

droSec1

dG=-29.0,  p-value=0.009901 dG=-28.3,  p-value=0.009901

droYak2

dG=-27.9,  p-value=0.009901 dG=-27.4,  p-value=0.009901 dG=-27.2,  p-value=0.009901

droEre2

dG=-29.0,  p-value=0.009901 dG=-28.3,  p-value=0.009901

droAna3

dG=-27.9,  p-value=0.009901 dG=-27.4,  p-value=0.009901 dG=-27.2,  p-value=0.009901

dp4

dG=-27.9,  p-value=0.009901 dG=-27.4,  p-value=0.009901 dG=-27.2,  p-value=0.009901

droPer1

dG=-27.9,  p-value=0.009901 dG=-27.4,  p-value=0.009901 dG=-27.2,  p-value=0.009901

droWil1

dG=-28.5,  p-value=0.009901 dG=-27.8,  p-value=0.009901 dG=-27.8,  p-value=0.009901 dG=-27.7,  p-value=0.009901

droVir3

dG=-27.9,  p-value=0.009901 dG=-27.4,  p-value=0.009901 dG=-27.2,  p-value=0.009901 dG=-26.9,  p-value=0.009901

droMoj3

dG=-27.9,  p-value=0.009901 dG=-27.8,  p-value=0.009901 dG=-27.4,  p-value=0.009901 dG=-27.2,  p-value=0.009901 dG=-27.1,  p-value=0.009901

droGri2

dG=-27.9,  p-value=0.009901 dG=-27.8,  p-value=0.009901 dG=-27.4,  p-value=0.009901 dG=-27.2,  p-value=0.009901 dG=-27.1,  p-value=0.009901

Generated: 03/07/2013 at 04:44 PM