New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb = 20185833 (0.040)75 (0.960) 45/94 NT 96.3% coding (579/687 nt) fliH flagellar biosynthesis protein
?W3110S.gb = 2521059 NA (NA)intergenic (+170/+30) insL/yfeA predicted transposase/predicted diguanylate cyclase

GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCA  >  W3110S.gb/2520997‑2521120
                                                              |                                                             
gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATaa                                                        <  1:716984/70‑1 (MQ=33)
  cTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTc                                                     >  1:7416954/1‑71 (MQ=38)
       cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCgg                                                <  1:561075/71‑1 (MQ=40)
          ggCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGttt                                             <  1:2367008/71‑1 (MQ=255)
             aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                          <  1:5665597/71‑1 (MQ=255)
             aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                          <  1:2458563/71‑1 (MQ=255)
             aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                          <  1:2215213/71‑1 (MQ=255)
               cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTa                                        >  1:2641518/1‑71 (MQ=255)
                   gCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGcc                                    <  1:313649/71‑1 (MQ=255)
                       cGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTa                                        >  1:411954/1‑63 (MQ=255)
                         cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTaa                              <  1:2198635/71‑1 (MQ=255)
                          aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGGTTTTTTACGCCTGTTAAt                             <  1:2426727/71‑1 (MQ=255)
                           ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                          >  1:6859132/1‑57 (MQ=255)
                           ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                            <  1:1208055/71‑1 (MQ=255)
                           ttcttcAGATGGTTCAACCCTTAAGGTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                                  <  1:1865089/65‑1 (MQ=255)
                            tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                            <  1:3852801/70‑1 (MQ=255)
                             cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCa                           <  1:2200102/70‑1 (MQ=255)
                             cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCa                           <  1:5835551/70‑1 (MQ=255)
                              ttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGc                         >  1:4054464/1‑71 (MQ=255)
                                 aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTa                       <  1:2677716/70‑1 (MQ=255)
                                  gATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAAt                     <  1:937614/71‑1 (MQ=255)
                                   aTGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                        >  1:1867471/1‑67 (MQ=255)
                                   aTGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                        >  1:5032340/1‑67 (MQ=255)
                                   aTGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATg                    >  1:972015/1‑71 (MQ=255)
                                       ttCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTg                <  1:5216563/71‑1 (MQ=255)
                                       ttCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTg                <  1:6058764/71‑1 (MQ=255)
                                         cAACCCTTAAGTTAGCGCTTATGGGATAATTCGCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGc               <  1:6432868/69‑1 (MQ=37)
                                          aaCCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCgg             >  1:7302611/1‑71 (MQ=255)
                                          aaCCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCg              <  1:4125203/70‑1 (MQ=255)
                                           aCCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGttttt                                           >  1:5069693/1‑40 (MQ=255)
                                           aCCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGttttt                                           >  1:1419246/1‑40 (MQ=255)
                                               ttAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGc                                     >  1:4213588/1‑42 (MQ=255)
                                                tAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGa         >  1:1116736/1‑69 (MQ=255)
                                                   gTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCa                           <  1:7397480/48‑1 (MQ=255)
                                                    ttAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATc   <  1:6423634/71‑1 (MQ=255)
                                                      aGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCa  >  1:1400021/1‑70 (MQ=255)
                                                              |                                                             
GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCA  >  W3110S.gb/2520997‑2521120

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