New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb = 20185830 (0.000)51 (1.340) 5/118 0.2 100% coding (579/687 nt) fliH flagellar biosynthesis protein
?W3110S.gb = 2521059 NA (NA)intergenic (+170/+30) insL/yfeA predicted transposase/predicted diguanylate cyclase

TAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCA  >  W3110S.gb/2521001‑2521116
                                                          |                                                         
tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATa                                                     >  1:1255156/1‑65 (MQ=31)
tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATa                                                     >  1:2981627/1‑65 (MQ=31)
           cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtt                                          <  1:805859/65‑1 (MQ=255)
           cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtt                                          <  1:2703171/65‑1 (MQ=255)
           cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtt                                          <  1:2333441/65‑1 (MQ=255)
           cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtt                                          <  1:2281180/65‑1 (MQ=255)
           cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtt                                          <  1:2065607/65‑1 (MQ=255)
           cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtt                                          <  1:319017/65‑1 (MQ=255)
           cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtt                                          <  1:1470084/65‑1 (MQ=255)
           cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtt                                          <  1:389222/65‑1 (MQ=255)
           cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtt                                          <  1:1207211/65‑1 (MQ=255)
                             aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTcc                                                >  1:2914286/1‑41 (MQ=255)
                             aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGtttt                                        >  1:2338893/1‑49 (MQ=255)
                             aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAAt                         >  1:1342001/1‑64 (MQ=255)
                             aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                        >  1:3055527/1‑65 (MQ=255)
                             aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                        >  1:692253/1‑65 (MQ=255)
                             aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                        >  1:28048/1‑65 (MQ=255)
                             aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                        >  1:2340414/1‑65 (MQ=255)
                             aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                        >  1:2649583/1‑65 (MQ=255)
                             aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                        >  1:1138122/1‑65 (MQ=255)
                             aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                        >  1:1659808/1‑65 (MQ=255)
                             aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                        >  1:1319347/1‑65 (MQ=255)
                             aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATc                        >  1:1204627/1‑65 (MQ=255)
                                                    cgcTTATGGGATAATTCCCCGGTTTTTTTAGGCCTGTTAATCAGCTAATGGCTGCGGGCGACCa  <  1:54712/64‑1 (MQ=255)
                                                    cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCa  <  1:2764289/64‑1 (MQ=255)
                                                    cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCa  <  1:2768286/64‑1 (MQ=255)
                                                    cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCa  <  1:2084477/64‑1 (MQ=255)
                                                    cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCa  <  1:199604/64‑1 (MQ=255)
                                                    cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCa  <  1:1581141/64‑1 (MQ=255)
                                                    cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCa  <  1:449637/64‑1 (MQ=255)
                                                    cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCa  <  1:505001/64‑1 (MQ=255)
                                                    cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCa  <  1:1312084/64‑1 (MQ=255)
                                                    cgcTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCa  <  1:785154/64‑1 (MQ=255)
                                                          |                                                         
TAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATAATTCCCCGGTTTTTTTACGCCTGTTAATCAGCTAATGGCTGCGGGCGACCA  >  W3110S.gb/2521001‑2521116

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.