New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb = 149941331 (1.360)3 (0.150) 3/92 NT 9.5% coding (814/1626 nt) mdoD glucan biosynthesis protein, periplasmic
?W3110S.gb = 1499426 29 (1.410)coding (827/1626 nt) mdoD glucan biosynthesis protein, periplasmic

GAAAATCACCTGTATGCGCGCAAAGACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAATTCATGACTCTGATCGTCTGTC  >  W3110S.gb/1499347‑1499483
                                                                  |                                                                      
gAAAATCACCTGTATGCGCGCAAAGACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTg                                                                    <  1:834623/71‑1 (MQ=255)
  aaaTCACCTGTATGCGCGCAAAGACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTGCg                                                                  <  1:2059220/71‑1 (MQ=255)
   aaTCACCTGTATGCGCGCGAAGACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTGCg                                                                  >  1:1543439/1‑70 (MQ=255)
       aCCTGTATGCGCGCAAAGACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTGCGGTAc                                                              >  1:1887786/1‑70 (MQ=255)
                  cgcAAAGACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCg                                                  <  1:1321081/71‑1 (MQ=255)
                  cgcAAAGACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCg                                                  <  1:1323582/71‑1 (MQ=255)
                     aaaGACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGAt                                               >  1:1385832/1‑71 (MQ=255)
                      aaGACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGAtg                                              >  1:442833/1‑71 (MQ=255)
                       aGACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCAttagtaccgcagctgaac                                                      >  1:1507302‑M1/1‑44 (MQ=255)
                       aGACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGAtgt                                             >  1:567411/1‑71 (MQ=255)
                        gACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTGCgg                                                                 <  1:1066483/50‑1 (MQ=255)
                         aCATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGc                                           >  1:429703/1‑71 (MQ=255)
                              aaaCAGCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATa                                       >  1:1582231/1‑70 (MQ=255)
                                 cAGCTGGGCATTGCGCCGATGACCAGTATGTTCAttagtaccg                                                               >  1:160412‑M1/1‑34 (MQ=255)
                                  aGCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGa                                         <  1:434328/64‑1 (MQ=255)
                                   gCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGa                                                <  1:1749366/56‑1 (MQ=255)
                                      gggCATTGCGCCGATGACCAGTATGTTCAttagtaccgc                                                              <  1:432036‑M1/39‑11 (MQ=255)
                                       ggCATTGCGCCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCAt                               <  1:1742905/69‑1 (MQ=255)
                                       ggCATTGCGCCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCAt                               <  1:2100104/69‑1 (MQ=255)
                                       ggCATTGCGCCGATGACCAGTATGTTCACCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCAt                               <  1:1056773/69‑1 (MQ=255)
                                          aTTGCGCCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCa                          <  1:1170493/71‑1 (MQ=255)
                                             gcgcCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAAt                       >  1:1370976/1‑71 (MQ=255)
                                             gcgcCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAAt                       >  1:1363303/1‑71 (MQ=255)
                                               gcCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAAtt                      <  1:193531/70‑1 (MQ=255)
                                                ccGATGACCAGTATGTTCAGCTGCGGTACTAATGAACg                                                     <  1:1732001/38‑1 (MQ=255)
                                                ccGATGACCAGTATGTTCAGCTGCGGTACTAATGAACg                                                     <  1:1393191/38‑1 (MQ=255)
                                                ccGATGACCAGTATGTTCAGCTGCGGTACTAATGAACg                                                     <  1:1372818/38‑1 (MQ=255)
                                                 cGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAAttaa                    >  1:1336557/1‑68 (MQ=255)
                                                  gATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAATTc                     >  1:243111/1‑68 (MQ=255)
                                                   aTGACCAGTATGTTCAGCTGCTGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAATTCATg                  >  1:1799806/1‑70 (MQ=255)
                                                   aTGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAATTCATGa                 >  1:2035560/1‑71 (MQ=255)
                                                    tGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAATTCAt                   >  1:1048378/1‑68 (MQ=255)
                                                          gTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAATTCATGACTCTGAt          >  1:1830848/1‑71 (MQ=255)
                                                           tATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCAt                               <  1:779505/49‑1 (MQ=255)
                                                           tATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATc                              >  1:1428138/1‑50 (MQ=255)
                                                              gTTCAGCTGCGGTACTAAtgaacatactggtcatcggcgcaatgcccagctgtttaatgtct               >  1:1507302‑M2/1‑18 (MQ=255)
                                                               ttCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAATTCATGACTCTGATCGTCt     <  1:1386711/71‑1 (MQ=255)
                                                               ttCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAATTCATGACTCTGATCGTCt     >  1:720404/1‑71 (MQ=255)
                                                               ttCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAATTCATGACTCTGATCGTCt     <  1:1384384/71‑1 (MQ=255)
                                                                  aGCTGCGGTACTAATGAACGTCGGATGTGCGATACAAt                                   >  1:1066018/1‑38 (MQ=255)
                                                                  aGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAATTCATGACTCTGATCGTCTGt   <  1:988790/70‑1 (MQ=255)
                                                                  aGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAATTCATGACTCTGATCGTCTGTc  <  1:50513/71‑1 (MQ=255)
                                                                  |                                                                      
GAAAATCACCTGTATGCGCGCAAAGACATTAAACAGCTGGGCATTGCGCCGATGACCAGTATGTTCAGCTGCGGTACTAATGAACGTCGGATGTGCGATACAATTCATCCGCAAATTCATGACTCTGATCGTCTGTC  >  W3110S.gb/1499347‑1499483

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ 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.