New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb = 210832620 (0.360)3 (0.060) 3/88 NT 12.7% coding (1038/1167 nt) wbbH O‑antigen polymerase
?W3110S.gb = 2108333 24 (0.490)coding (1031/1167 nt) wbbH O‑antigen polymerase

TGCTAATATTAGTCATGAAGCTTTCATGATAAAAAATGAAAGAAAAGGTATAAATAAAATATGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAACACCTGAAATACAGCCATGAATAAC  >  W3110S.gb/2108265‑2108400
                                                                    |                                                                   
tGCTAATATTAGTCATGAAGCTTTCATGATAAAAAATGAAAGAAAAGGTATAAATAAAATATGAATAAAa                                                                    <  1:1744475‑M2/70‑2 (MQ=255)
 gCTAATATTAGTCATGAAGCTTTCATGATAAAAAATGAAAGAAAAGGTATAAATAAAATATGAATAAAata                                                                  <  1:1150321‑M2/71‑4 (MQ=255)
  cTAATATTAGTCATGAAGCTTTCATGATAAAAAATGAAAGAAAAGGTATAAATAAAATATGAATAAAatat                                                                 <  1:4933179‑M2/71‑5 (MQ=255)
   tAATATTAGTCATGAAGCTTTCATGATAAAAAATGAAAGAAAAGGTATAAATAAAATATGAATAAAatatt                                                                <  1:1060981‑M2/71‑6 (MQ=255)
   tAATATTAGTCATGAAGCTTTCATGATAAAAAATGAAAGAAAAGGTATAAATAAAATATGAATAAAat                                                                   >  1:1160147‑M2/1‑66 (MQ=255)
     atatTAGTCATGAAGCTTTCATGATAAAAAATGAAAGAAAAGGTATAAATAAAATATGAATAAAatatttt                                                              <  1:2821199‑M2/71‑8 (MQ=255)
     atatTAGTCATGAAGCTTTCATGATAAAAAATGAAAGAAAAGGTATAAATAAAATATGAATAAAatatttt                                                              >  1:2735378‑M2/1‑64 (MQ=255)
       atTAGTCATGAAGCTTTCATGATAAAAAATGAAAGAAAAGGTATAAATAAAATATGAATAAAa                                                                    <  1:4398963‑M2/63‑2 (MQ=255)
              aTGAAGCTTTCATGATAAAAAATGAAAGAAAAGGTATAAATAAAATATGAATAAAATATTTTc                                                             >  1:2494394/1‑63 (MQ=255)
                 aaGCTTTCATGATAAAAAATGAAAGAAAAGGTATAAATAAAATATGAATAAAATATTTTCACAGATATGt                                                   >  1:3141582/1‑70 (MQ=255)
                 aaGCTTTCATGATAAAAAATGAAAGAAAAGGTATAAATAAAATATGAATAAAATATTTTCACAGATATGt                                                   >  1:5174395/1‑70 (MQ=255)
                  aGCTTTCATGATAAAAAATGAAAGAAAAGGTATAAATAAAATATGAATAAAATATTTTCACAGATATGTaa                                                 >  1:1915537/1‑71 (MQ=255)
                            aTAAAAAATGAAAGAAAAGGTATAAATAAAATATGAATAAAATATTTTCACAGATATGTAATTTCGAGACa                                       <  1:340487/71‑1 (MQ=255)
                              aaaaaaTGAAAGAAAAGGTATAAATAAAATATGAATAAAATATTTTCACAGATATGTAATTTCGAGACAAt                                     >  1:4886037/1‑71 (MQ=255)
                                     gaaagaaaAGGTATAAATAAAATATGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCAt                               >  1:2827940/1‑70 (MQ=255)
                                        agaaaAGGTATAAATAAAATATGAATAAAATATTTTCACAGATATGTAATTTCGAGACAAtct                                   <  1:1029156/63‑1 (MQ=255)
                                         gaaaAGGTATAAATAAAATATGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATaaaa                           >  1:3556271/1‑70 (MQ=255)
                                           aaaGGTATAAATAAAATATGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAacac                        >  1:2424792/1‑71 (MQ=255)
                                           aaaGGTATAAATAAAATATGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAacac                        >  1:474697/1‑71 (MQ=255)
                                           aaaGGTATAAATAAAATATGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAacac                        >  1:2646733/1‑71 (MQ=255)
                                           aaaGGTATAAATAAAATATGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAaca                         <  1:4124657/70‑1 (MQ=255)
                                                  taaataaaATATGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCAt                               <  1:2912124/57‑1 (MQ=255)
                                                  taaataaaATATGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAACACCTGaaa                  <  1:5145204/70‑1 (MQ=255)
                                                    aataaaataTGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAACACCTGAAATa                >  1:2038199/1‑70 (MQ=255)
                                                    aataaaataTGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAACACCTGAAATa                >  1:4935190/1‑70 (MQ=255)
                                                     ataaaataTGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAACACCTGAAATa                >  1:4278243/1‑69 (MQ=255)
                                                        aaataTGAATAAAatattttatttataccttttctttcattttttatcatgaaagcttcat                     >  1:3962146‑M2/1‑13 (MQ=255)
                                                        aaataTGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAACACCt                      >  1:2446310/1‑60 (MQ=255)
                                                           tatGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCAt                               <  1:534138/48‑1 (MQ=255)
                                                           tatGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAACACCTGa                    <  1:858481/59‑1 (MQ=255)
                                                           tatGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAACACCTGAAATACAGCCAt         >  1:2259488/1‑70 (MQ=255)
                                                            atGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAACAcc                       <  1:2306490/55‑1 (MQ=255)
                                                            atGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAACACCTGAAATACAGCCATg        >  1:3219564/1‑70 (MQ=255)
                                                             tGAATAAAatattttatttataccttttctttcattttttatcatgaaagcttcatgactaatattagcag      >  1:3734508‑M2/1‑8 (MQ=255)
                                                             tGAATAAAatattttatttataccttttctttcatttt                                       >  1:1951895‑M2/1‑8 (MQ=255)
                                                              gAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAACACCTGAAATACAGCCATGaa      >  1:3071451/1‑70 (MQ=255)
                                                                 tAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAac                          <  1:4512783/47‑1 (MQ=255)
                                                                  aaaaTATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAACACCTGAAATACAGCCATGAATAAc  >  1:3465813/1‑70 (MQ=255)
                                                                    aaTATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAacac                        >  1:1347067/1‑46 (MQ=255)
                                                                    aaTATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAACACCTGAAATACa              >  1:3304989/1‑56 (MQ=255)
                                                                    |                                                                   
TGCTAATATTAGTCATGAAGCTTTCATGATAAAAAATGAAAGAAAAGGTATAAATAAAATATGAATAAAATATTTTCACAGATATGTAATTTCGAGACAATCTCCATAAAACACCTGAAATACAGCCATGAATAAC  >  W3110S.gb/2108265‑2108400

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.