New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb = 20185834 (0.070)41 (0.790) 31/96 NT 91.4% coding (579/687 nt) fliH flagellar biosynthesis protein
?W3110S.gb = 2521059 NA (NA)intergenic (+170/+30) insL/yfeA predicted transposase/predicted diguanylate cyclase

TCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCT  >  W3110S.gb/2520998‑2521128
                                                             |                                                                     
tCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAAtt                                                              <  1:1527855/71‑1 (MQ=35)
tCTTAACGAGGCACCGAGGAGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAAtt                                                              <  1:4779361/71‑1 (MQ=25)
 cTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTc                                                             >  1:3014271/1‑71 (MQ=38)
   tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTccc                                                           >  1:4035694/1‑71 (MQ=39)
         ggCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGttt                                                     <  1:2237079/71‑1 (MQ=255)
            aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                                  <  1:351912/71‑1 (MQ=255)
                     tCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGttttt                                                   <  1:1029757/61‑1 (MQ=255)
                       gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGtt                                        <  1:3975608/70‑1 (MQ=255)
                        cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTaa                                      <  1:395162/71‑1 (MQ=255)
                         aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTAc                                               <  1:523547/61‑1 (MQ=255)
                          ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAAt                                     >  1:1088887/1‑70 (MQ=255)
                            cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATaa                                                                >  1:4994859/1‑41 (MQ=35)
                            cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATaa                                                                >  1:3297051/1‑41 (MQ=35)
                             ttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATaa                                                                >  1:2125739/1‑40 (MQ=37)
                             ttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                                    <  1:2407539/68‑1 (MQ=255)
                                aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCAGGTTTTTTTACGCCTGTTAATCAGCTa                               <  1:3755364/70‑1 (MQ=255)
                                    ggTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCg                                                         >  1:1797097/1‑40 (MQ=255)
                                      ttCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTg                        <  1:4509602/71‑1 (MQ=255)
                                       tCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGc                       >  1:3437323/1‑71 (MQ=255)
                                           cccTTAAGTTAGCGCTTATGGGATAATTCCCCGGttttttt                                                 >  1:4035797/1‑41 (MQ=255)
                                           cccTTAAGTTAGCGCTTATGGGATAATTCCCCGGttttttt                                                 >  1:1336738/1‑41 (MQ=255)
                                           cccTTAAGTTAGCGCTTATGGGATAATTCCCCGGttt                                                     <  1:3946485/37‑1 (MQ=255)
                                           cccTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGc                   >  1:3456144/1‑71 (MQ=255)
                                              ttAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGc                                             >  1:3966683/1‑42 (MQ=255)
                                              ttAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGAc                >  1:886023/1‑71 (MQ=255)
                                               tAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGAc                >  1:1928984/1‑70 (MQ=255)
                                                    tAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCgg                     >  1:3875927/1‑60 (MQ=255)
                                                    tAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCgg                     >  1:3602266/1‑60 (MQ=255)
                                                     aGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCa          <  1:2330865/70‑1 (MQ=255)
                                                       cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTa                               >  1:1562239/1‑47 (MQ=255)
                                                       cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCaaa        >  1:1671023/1‑70 (MQ=255)
                                                          ttATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGc                   <  1:3608024/56‑1 (MQ=255)
                                                           tATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCaaa        >  1:1080787/1‑66 (MQ=255)
                                                            aTGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCt  >  1:5179162/1‑71 (MQ=255)
                                                            aTGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCt  >  1:1193689/1‑71 (MQ=255)
                                                             |                                                                     
TCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAATACCCT  >  W3110S.gb/2520998‑2521128

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.