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

GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAA  >  W3110S.gb/2520997‑2521122
                                                              |                                                               
gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATaa                                                          <  1:2587763/70‑1 (MQ=33)
gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATaa                                                          <  1:2713948/70‑1 (MQ=33)
gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATaa                                                          <  1:3381367/70‑1 (MQ=33)
         aGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtt                                                <  1:1933320/71‑1 (MQ=255)
             aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGttttt                                             >  1:3317763/1‑70 (MQ=255)
                gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTc                                                       >  1:2601775/1‑57 (MQ=40)
                gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTc                                                       >  1:2952185/1‑57 (MQ=40)
                gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTAc                                         >  1:3728628/1‑71 (MQ=255)
                gAGGCGTCGCATTCTTCAGATGGTTAAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCg                                                   >  1:970013/1‑61 (MQ=255)
                   gCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGcc                                      >  1:2924564/1‑71 (MQ=255)
                        gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTa                                 <  1:1743050/71‑1 (MQ=255)
                           ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTcccc                                                    <  1:102071/49‑1 (MQ=255)
                             cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttttt                                            >  1:1710353/1‑55 (MQ=255)
                                cAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                          <  1:1099055/70‑1 (MQ=255)
                                  gATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGcc                                      <  1:2308636/56‑1 (MQ=255)
                                   aTGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCt                          <  1:3212518/67‑1 (MQ=255)
                                      gTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGc                    <  1:3529115/70‑1 (MQ=255)
                                        tCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCt                                     >  1:631262/1‑51 (MQ=255)
                                           aCCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCggg              >  1:2762680/1‑71 (MQ=255)
                                            cccTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGc             <  1:3807998/71‑1 (MQ=255)
                                               ttAAGTTAGCGCTTATGGGATAATTCCCCGGttttt                                             <  1:2928104/36‑1 (MQ=255)
                                                      aGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCa    <  1:3072146/70‑1 (MQ=255)
                                                       gcgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCaaa  <  1:2954474/71‑1 (MQ=255)
                                                              |                                                               
GTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCAATCAAA  >  W3110S.gb/2520997‑2521122

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.