New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb 2109295 =40 (0.710)4 (0.080) 3/90 NT 9.4% coding (69/1167 nt) wbbH O‑antigen polymerase
?W3110S.gb 2109315 = 42 (0.850)coding (49/1167 nt) wbbH O‑antigen polymerase

ATCTGACGTTATTTCATATCCCAATAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTA  >  W3110S.gb/2109226‑2109365
                                                                     |                                                                      
atCTGACGTTATTTCATATCCCAATAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAata                                                                       >  1:3142915/1‑71 (MQ=255)
atCTGACGTTATTTCATATCCCAATAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAat                                                                        <  1:3908938/70‑1 (MQ=255)
aatataatgagtacggattaatgatctatcttgtaattagtgtctttctcattacagcatttatctgttta                                                                       >  1:1051049‑M1/70‑71 (MQ=255)
  tataatgagtacggattaatgatctatcttgtaattagtgtctttctcattacagcatttatctgttta                                                                       <  1:1213527‑M1/2‑1 (MQ=255)
           tttCATATCCCAATAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAata                                                                       >  1:4233787/1‑60 (MQ=255)
           tttCATATCCCAATAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTa                                                            <  1:3579691/71‑1 (MQ=255)
             tCATATCCCAATAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCcttctt                                                             >  1:2723679/1‑68 (MQ=255)
                  tCCCAATAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGAtata                                                     <  1:4238674/71‑1 (MQ=255)
                     cAATAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAAc                                                  <  1:2418971/71‑1 (MQ=255)
                      aaTAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAatata                                                                     >  1:3985273/1‑51 (MQ=255)
                       aTAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAatat                                                                      <  1:420066/49‑1 (MQ=255)
                       aTAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATaaa                                                   >  1:3617274/1‑68 (MQ=255)
                       aTAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAg                                                <  1:2720310/71‑1 (MQ=255)
                         ctatcttgtaattagtgtctttctcattacagcatttatctgttTatat                                                                    <  1:671615‑M1/5‑1 (MQ=255)
                         aaGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGa                                               >  1:540732/1‑70 (MQ=255)
                          tatcttgtaattagtgtctttctcattacagcatttatctgtttat                                                                      >  1:4421315‑M1/44‑46 (MQ=255)
                          aGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATa                                             <  1:2199917/71‑1 (MQ=255)
                           gACCAGTGCGAAGATGATATTAACGCATAAGGCTGGATAAAATATATCCTTCTTAAGATATaa                                                    <  1:250592/63‑1 (MQ=255)
                            tcttgtaattagtgtctttctcattacagcatttatctgttt                                                                        >  1:4971064‑M1/42‑42 (MQ=255)
                            aCCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATaaa                                           <  1:4413216/71‑1 (MQ=255)
                                gTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCt                                       >  1:5149212/1‑71 (MQ=255)
                                 tGCGAAGATGATATTAACGCATACGGCTGGATAAAatata                                                                     <  1:1536421/40‑1 (MQ=255)
                                   attagtgtctttctcattacagcatttatctgttTATATCCTTCTTAAGATATaa                                                    <  1:2436248‑M1/21‑1 (MQ=255)
                                     aaGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAAt                                  <  1:2929565/71‑1 (MQ=255)
                                      aGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATaaa                                                   >  1:4878007/1‑53 (MQ=255)
                                      aGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTaa                                   >  1:4958494/1‑69 (MQ=255)
                                      aGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAAt                                  >  1:3841251/1‑70 (MQ=255)
                                      aGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAAt                                  >  1:1561209/1‑70 (MQ=255)
                                      aGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAAt                                  >  1:4394183/1‑70 (MQ=255)
                                         tgatATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGat                                                        >  1:1963393/1‑45 (MQ=255)
                                         tgatATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGat                                                        >  1:3009417/1‑45 (MQ=255)
                                         tgatATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATgaga                              <  1:2973464/71‑1 (MQ=255)
                                         tgatATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATgag                               >  1:2412919/1‑70 (MQ=255)
                                          gatATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAAc                                                  >  1:2012209/1‑50 (MQ=255)
                                            tatTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATgag                               >  1:3557520/1‑67 (MQ=255)
                                            tatTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAg                           <  1:2566342/71‑1 (MQ=255)
                                             atTAACGCATACGGCTGGATAAAATATATCcttctt                                                             <  1:3074934/36‑1 (MQ=255)
                                              ttAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATaaa                                                   >  1:4582828/1‑45 (MQ=255)
                                              ttAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGa                          <  1:3855587/70‑1 (MQ=255)
                                              ttAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGa                          >  1:4339163/1‑70 (MQ=255)
                                                 aCGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACt                      <  1:3004491/71‑1 (MQ=255)
                                                  attacagcatttatctgttTATATCTTAAGAAGGATATAAACAGATAAATGCTGTAATGAGAAAGACACt                      <  1:4172186‑M1/51‑1 (MQ=255)
                                                  attacagcatttatctgttTATATCCTTCTTAAGat                                                        <  1:1868489‑M1/17‑1 (MQ=255)
                                                  attacagcatttatctgttTATATACTTCTTAAGat                                                        <  1:11928‑M1/17‑1 (MQ=255)
                                                   gCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGa                                               >  1:5151712/1‑44 (MQ=255)
                                                     aTACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGa                          <  1:2490106/63‑1 (MQ=255)
                                                        cGGCTGGATAAAATATATCCTTCTTAAGATATaaa                                                   >  1:3647772/1‑35 (MQ=255)
                                                        cGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTAc                <  1:361447/70‑1 (MQ=255)
                                                        cGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTAc                <  1:4514452/70‑1 (MQ=255)
                                                         ggCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACaa              <  1:3568055/71‑1 (MQ=255)
                                                         ggCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACaa              <  1:2555154/71‑1 (MQ=255)
                                                         ggCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACaa              <  1:1892866/71‑1 (MQ=255)
                                                           cTGGATAAAATATATCCTTCTTAAGATATAAACAGAt                                              <  1:310234/37‑1 (MQ=255)
                                                            tGGATAAAATATATCCTTCTTAAGATATAAACAGATa                                             <  1:799697/37‑1 (MQ=255)
                                                             ggATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTAc                <  1:4912149/65‑1 (MQ=255)
                                                                tAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTAc                <  1:4393505/62‑1 (MQ=255)
                                                                tAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAaga            >  1:3830383/1‑66 (MQ=255)
                                                                tAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAaga            >  1:1550565/1‑66 (MQ=255)
                                                                   aaTATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTAc                <  1:4860218/59‑1 (MQ=255)
                                                                   aaTATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTAc                <  1:2911979/59‑1 (MQ=255)
                                                                   aaTATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCAt    >  1:4740725/1‑71 (MQ=255)
                                                                   aaTATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCAt    <  1:1030301/71‑1 (MQ=255)
                                                                    atatatCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGaaa                            >  1:1708156/1‑46 (MQ=255)
                                                                    atatatCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCAt    >  1:3954474/1‑70 (MQ=255)
                                                                     TATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAagatag         <  1:306850‑M1/64‑1 (MQ=255)
                                                                     TATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTa  >  1:4653668‑M1/1‑71 (MQ=255)
                                                                     |                                                                      
ATCTGACGTTATTTCATATCCCAATAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTA  >  W3110S.gb/2109226‑2109365

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.