New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AE015451 5827746 =29 (0.710)3 (0.080) 3/320 NT 10.1% coding (659/711 nt) mtgA Monofunctional biosynthetic peptidoglycan transglycosylase
?AE015451 5827764 = 26 (0.670)coding (677/711 nt) mtgA Monofunctional biosynthetic peptidoglycan transglycosylase

CCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCC  >  AE015451/5827554‑5827988
                                                                                                                                                                                                                  |                                                                                                                                                                                                                                
ccGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACtcctcc                                                                                                                                                                                                            <  2:464465/233‑1 (MQ=255)
     tGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACtcctcc                                                                                                                                                                                                            <  1:554414/228‑1 (MQ=255)
      ggggCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCtt                                                                                                                                                                                                   >  1:508570/1‑236 (MQ=255)
           aaGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGc                                                                                                                                                                                                     <  2:591127/229‑1 (MQ=255)
           aaGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGc                                                                                                                                                                                                     >  1:591128/1‑229 (MQ=255)
                       ggTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggcccggcc                                                                                                                                                                                  >  1:646419/1‑236 (MQ=255)
                                 cAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggcc                                                                                                                                                                                       >  2:594430/1‑221 (MQ=255)
                                    cggcgCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGc                                                                                                                                                                                                     <  2:682904/204‑1 (MQ=255)
                                      gcgcGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAAc                                                                                                                                                                   <  1:651112/236‑1 (MQ=255)
                                          gTTATCGCTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggcc                                                                                                                                                                                       <  1:233590/212‑1 (MQ=255)
                                                  tttGCCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGCGCGGGCCCAGCTAcc                                                                                                                                                                                                                                <  1:354322/163‑1 (MQ=255)
                                                   cctacaggggatgcggttgcctgtaggcgcgggtttgcccgcgaggaggccggggatgcggttgcccgtaggagcgggtttacccgcgaaagggccgggcctgcattcaaagctttctggaggagtccagctgcatcaggtagctgggcccgcccagctAcc                                                                                                                                                                                                                                <  1:713784‑M2/3‑1 (MQ=255)
                                                            aCGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAg                                                                                                                                                                                          <  2:197029/191‑1 (MQ=255)
                                                                            ccgccAGCGGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGa                                                                                                                                                                                                <  2:158263/168‑1 (MQ=255)
                                                                               ccAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCggg                                                                                                                           >  1:704510/1‑235 (MQ=255)
                                                                                              gcTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTAc                                                                                                            <  1:57856/235‑1 (MQ=255)
                                                                                                gaggccggggatgcggttgcccgtaggagcgggtttacccgcgaaagggccgggcctgcattcaaagctttctggaggagtccagctgcatcaggtAGCTGGGCGGGCCCCGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCCggcccg                                                                                                                                                                                     <  1:265844‑M1/64‑1 (MQ=255)
                                                                                                    ccggggatgcggtttcccgtaggagcgggtttacccgcgaaagggccgggcctgcattcaaagctttctggaggagtccagctgcatcaggtagctgggcccgcccagctACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAg                                                                                                                                                                                          <  1:162574‑M2/41‑1 (MQ=255)
                                                                                                       ggTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAAc                                                                                                     >  1:689411/1‑233 (MQ=255)
                                                                                                       ggTGCTGCCAAGCCAAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAAcc                                                                                                    <  1:392959/234‑1 (MQ=255)
                                                                                                                               gAGTGCCAGCCGACCGGGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCgg                                                                                 <  2:531614/229‑1 (MQ=255)
                                                                                                                                   gacccgcgaaagggccgggcctgcattcaaagctttctggaggagtccagctgcatcaggtAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTg                                                                                                                                                                                                 >  2:265843‑M1/62‑113 (MQ=255)
                                                                                                                                     cAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTcg                                                                                                                               >  2:755970/1‑177 (MQ=255)
                                                                                                                                                  cgggccggcattcaaagctttctggaggagtccagctgcatcaggtagcggggcccgcccagctACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGcctc                                                                                                                                   <  1:265844‑M2/96‑1 (MQ=255)
                                                                                                                                                   ccTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCaa                                                                                                                        <  1:98154/170‑1 (MQ=255)
                                                                                                                                                       ctgcattcaaagctttctggaggagtccagctgcatcaggtAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGAAACCGCATCCCCgg                                                                                                                                       <  1:162574‑M1/110‑1 (MQ=255)
                                                                                                                                                        gTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTa                                                                                        <  2:508569/197‑1 (MQ=255)
                                                                                                                                                          ggCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAg                                                                                                          <  1:518396/177‑1 (MQ=255)
                                                                                                                                                              caaagctttctggaggagtccagctgcatcaggtagctgggcccgcccagctACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGgtc                                                                                                                                                                      >  2:265843‑M2/53‑112 (MQ=255)
                                                                                                                                                                ccgggccgtCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGGCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCg                                                                                                                                                >  2:664305/1‑133 (MQ=255)
                                                                                                                                                                   ggccgaCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCCACCTGAGGCAGCTGGACTCCTCCCGAAATCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTAc                                                                                                                                                         <  2:361950/121‑1 (MQ=38)
                                                                                                                                                                                      cAGATGAGCCAGCCGGGCGGGCCCAGCTACCTGATGCAGCt                                                                                                                                                                                                                      <  2:492029/41‑1 (MQ=255)
                                                                                                                                                                                         aTGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTcgcg                                                                                                                             <  2:157943/127‑1 (MQ=255)
                                                                                                                                                                                            aGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGcc                                        <  2:588061/209‑1 (MQ=255)
                                                                                                                                                                                             ggtAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAgg                                                                                      <  1:713784‑M1/159‑1 (MQ=255)
                                                                                                                                                                                                 GCTGGGCGGGCCCCGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTGTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCTGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAAcc                                                                                                    >  2:247302‑M1/1‑144 (MQ=255)
                                                                                                                                                                                                     GGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGgtc                                                                                                                                                                      <  2:168402‑M1/74‑2 (MQ=255)
                                                                                                                                                                                                       CGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAAccc  >  1:256761‑M1/1‑236 (MQ=255)
                                                                                                                                                                                                          gcccagctACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCg                                       <  1:764047‑M2/188‑1 (MQ=255)
                                                                                                                                                                                                                  |                                                                                                                                                                                                                                
CCGAGTGGGGCAAGGGCGTGTTTGGTGCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCC  >  AE015451/5827554‑5827988

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ 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.