New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb = 20185830 (0.000)58 (1.560) 19/98 NT 100% 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
                                                              |                                                             
gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAAt                                                       >  1:1055374/1‑71 (MQ=34)
gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAAt                                                       >  1:12647/1‑71 (MQ=34)
gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAAt                                                       >  1:2458115/1‑71 (MQ=34)
                gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTa                                        >  1:1690506/1‑70 (MQ=255)
                 aGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACg                                      <  1:2945321/71‑1 (MQ=255)
                     gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                                  <  1:1182798/71‑1 (MQ=255)
                     gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                                  <  1:3183261/71‑1 (MQ=255)
                       cGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTc                                                     >  1:2603822/1‑50 (MQ=255)
                       cGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTc                                                     >  1:3029686/1‑50 (MQ=255)
                       cGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTc                                                     >  1:3345468/1‑50 (MQ=255)
                        gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTa                               <  1:3617662/71‑1 (MQ=255)
                        ccATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGtt                                >  1:370728/2‑70 (MQ=255)
                             cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTg                                  >  1:1656411/1‑63 (MQ=255)
                                    tGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                          >  1:739354/1‑48 (MQ=255)
                                    tGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATgg                   >  1:3467096/1‑71 (MQ=255)
                                      gTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                        <  1:2147328/64‑1 (MQ=255)
                                      gTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                        <  1:451417/64‑1 (MQ=255)
                                      gTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                        <  1:3124118/64‑1 (MQ=255)
                                      gTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                        <  1:1884070/64‑1 (MQ=255)
                                        tCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGc               >  1:408044/1‑71 (MQ=255)
                                                 aaGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTaa                              >  1:3487810/1‑47 (MQ=255)
                                                 aaGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTaa                              >  1:3088714/1‑47 (MQ=255)
                                                      aGCGCTTATTGGATAATTCCCCGGTTTTTTTACGCCTGTTa                               <  1:2563593/41‑1 (MQ=255)
                                                      aGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCa  <  1:147069/70‑1 (MQ=255)
                                                      aGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCa  >  1:1752515/1‑70 (MQ=255)
                                                              |                                                             
GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCA  >  W3110S.gb/2520997‑2521120

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