New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? USA300TCH1516_ALE 34173 =73 (1.250)9 (0.160) 5/248 NT 12.1% coding (462/480 nt) rlmH Ribosomal RNA large subunit methyltransferase H
?USA300TCH1516_ALE = 57858 61 (1.090)intergenic (+304/+613) USA300HOU_RS00240/USA300TCH1516_00049 hypothetical protein/hypothetical protein

TGTCAAAAATCATGAACCTCATTACTTATGATAAGCTTCTTAAAAACATAACAGCAATTCACATAAACCTCATATGTTCTGATACATTCAAAATCCCTTTATGAAGCGGCTGAAAAAACCGCATCATTTATGATATGCTTCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  USA300TCH1516_ALE/34314‑34173
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcttctGCTTATCAGTTAATGATGCGGTTTTTGAAAAATAACTACAGAAATTTATAGAAGAACAACAAAACTCAAACGCTATATCATATAATAAGTTCTTAGAAGCTTACAAAGTAGACACGTCTTC  <  USA300TCH1516_ALE/57858‑57738
                                                                                                                                                                                                                                                                       
TGTCAAAAATCATGAACCTCATTACTTATGATAAGCTTCTTAAAAACATAACAGCAATTCACATAAACCTCATATGTTCTGATACATTCAAAATCCCTTTATGAAGCGGCTGAAAAAACCGCATCATTTATGATATGCTTC                                                                                                                            <  2:553776/141‑1
  TCAAAAATCATGAACCTCATTACTTATGATAAGCTTCTTAAAAACATAACAGCAATTCACATAAACCTCATATGTTCTGATACATTCAAAATCCCTTTATGAAGCGGCTGAAAAAACCGCATCATTTATGATATGCTTCTG                                                                                                                          <  2:148793/141‑1
                         TTATGATAAGCTTCTTAAAAACATAACAGCAATTCACATAAACCTCATATGTTCTGATACATTCAAAATCCCTTTATGAAGCGGCTGAAAAAACCGCATCATTTATGATATGCTTCTGCTTATCAGTTAATGATGCGGTTT                                                                                                   <  2:534022/141‑1
                               TAAGCTTCTTAAAAACATAACAGCCATTCACATAAACCTCATATGTTCTGATACATTCAAAATCCCTTTATGAAGCGGCTGAAAAAACCTCATCATTTATGATATGCTTCTGCGTATCGGTTAATGATGCGGTTTTTaatt                                                                                             <  2:43339/141‑5
                                                           CACATAAACCTCATATGTTCTGATACATTCAAAATCCCTTTATGAAGCGGCTGAAAAAACCGCATCATTTATGATATGCTTCTGCTTATCAGTTAATGATGCGGTTTTTGAAAAATAACTACAGAAATTTATAGAAGAACA                                                                 >  1:494811/1‑141
                                                           CACATAAACCTCATATGTTCTGATACATTCAAAATCCCTTTATGAAGCGGCTGAAAAAACCGCATCATTTATGATATGCTTCTGCTTATCAGTTAATGATGCGGTTTTTGAAAAATAACTACAGAAATTTATAGAAGAACA                                                                 >  1:94472/1‑141
                                                           CACATAAACCTCATATGTTCTGATACATTCAAAATCCCTTTATGAAGCGGCTGAAAAAACCGCATCATTTATGATATGCTTCTGCTTATCAGTTAATGATGCGGTTTTTGAAAAATAACTACAGAAATTTATAGAAGAACA                                                                 >  2:269236/1‑141
                                                                              CTGATACATTCAAAATCCCTTTATGAAGCGGCTGAAAAAACCGCATCATTTATGATATGCTTCTGCTTATCAGTTAATGATGCGGTTTTTGAAAAATAACTACAGAAATTTATAGAAGAACAACAAAACTCAAACGCTATt                                              >  1:79998/1‑140
                                                                              CTGATACATTCAAAATCCCTTTATGAAGCGGCTGAAAAAACCGCATCATTTATGATATGCTTCTGCTTATCAGTTAATGATGCGGTTTTTGAAAAATAACTACAGAAATTTATAGAAGAACAACAAAACTCAAACGCTATA                                              >  1:356618/1‑141
                                                                                                             CTGAAAAAACCGCATCATTTATGATATGCTTCTGCTTATCAGTTAATGATGCGGTTTTTGAAAAATAACTACAGAAATTTATAGAAGAACAACAAAACTCAAACGCTATATCATATACTAAGTTCTTAGAAGCTTACAAAG               >  2:162799/1‑141
                                                                                                                          ATCATTTATGATATGCTTCTGCTTATCAGTTAATGATGCGGTTTTTGAAAAATAACTACAGAAATTTATAGAAGAACAACAAAACTCAAACGCTATATCATATAATAAGTTCTTAGAAGCTTACAAAGTAGACACGTCTTC  <  1:281766/141‑1
                                                                                                                                                                                                                                                                       
TGTCAAAAATCATGAACCTCATTACTTATGATAAGCTTCTTAAAAACATAACAGCAATTCACATAAACCTCATATGTTCTGATACATTCAAAATCCCTTTATGAAGCGGCTGAAAAAACCGCATCATTTATGATATGCTTCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  USA300TCH1516_ALE/34314‑34173
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcttctGCTTATCAGTTAATGATGCGGTTTTTGAAAAATAACTACAGAAATTTATAGAAGAACAACAAAACTCAAACGCTATATCATATAATAAGTTCTTAGAAGCTTACAAAGTAGACACGTCTTC  <  USA300TCH1516_ALE/57858‑57738

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.