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

CAATAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTACTCATTATATTT  >  AP009048/2109247‑2109384
                                                                    |                                                                     
cAATAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAAc                                                                     <  1:2418971/71‑1 (MQ=255)
  aTAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATaaa                                                                      >  1:3617274/1‑68 (MQ=255)
  aTAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAg                                                                   <  1:2720310/71‑1 (MQ=255)
    aaGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGa                                                                  >  1:540732/1‑70 (MQ=255)
     aGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATa                                                                <  1:2199917/71‑1 (MQ=255)
      gACCAGTGCGAAGATGATATTAACGCATAAGGCTGGATAAAATATATCCTTCTTAAGATATaa                                                                       <  1:250592/63‑1 (MQ=255)
       aCCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATaaa                                                              <  1:4413216/71‑1 (MQ=255)
           gTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCt                                                          >  1:5149212/1‑71 (MQ=255)
              attagtgtctttctcattacagcatttatctgttTATATCCTTCTTAAGATATaa                                                                       <  1:2436248‑M1/21‑1 (MQ=255)
                aaGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAAt                                                     <  1:2929565/71‑1 (MQ=255)
                 agtgtctttctcattacagcatttatctgtttatatccttcttaagatataAACAGATAAATGCTGTAAt                                                     <  1:4172186‑M2/19‑1 (MQ=255)
                 aGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATaaa                                                                      >  1:4878007/1‑53 (MQ=255)
                 aGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTaa                                                      >  1:4958494/1‑69 (MQ=255)
                 aGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAAt                                                     >  1:4394183/1‑70 (MQ=255)
                 aGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAAt                                                     >  1:3841251/1‑70 (MQ=255)
                 aGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAAt                                                     >  1:1561209/1‑70 (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)
                         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)
                              gCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGa                                                                  >  1:5151712/1‑44 (MQ=255)
                                aTACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGa                                             <  1:2490106/63‑1 (MQ=255)
                                   cGGCTGGATAAAATATATCCTTCTTAAGATATaaa                                                                      >  1:3647772/1‑35 (MQ=255)
                                   cGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTAc                                   <  1:4514452/70‑1 (MQ=255)
                                   cGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTAc                                   <  1:361447/70‑1 (MQ=255)
                                    ggCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACaa                                 <  1:1892866/71‑1 (MQ=255)
                                    ggCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACaa                                 <  1:3568055/71‑1 (MQ=255)
                                    ggCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACaa                                 <  1:2555154/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:2911979/59‑1 (MQ=255)
                                              aaTATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTAc                                   <  1:4860218/59‑1 (MQ=255)
                                              aaTATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCAt                       <  1:1030301/71‑1 (MQ=255)
                                              aaTATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCAt                       >  1:4740725/1‑71 (MQ=255)
                                               ttatatcttaagaaggatataAACAGATAAATGCTGTAATGAGAAAGACACTAAt                                      <  1:2436248‑M2/34‑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)
                                                 ATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTaa                    <  1:2506331‑M1/71‑1 (MQ=255)
                                                  TATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAagatag                            >  1:4267201‑M1/1‑62 (MQ=255)
                                                  TATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAAt                   >  1:3664604‑M1/1‑71 (MQ=255)
                                                   atcttaagaagtatataAACAGATAAATGCTGTAAt                                                     <  1:11928‑M2/19‑1 (MQ=255)
                                                   atcttaagaaggatataAACAGATAAATGCTGTAAt                                                     <  1:1868489‑M2/19‑1 (MQ=255)
                                                      cttcttAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAagatagat                          <  1:4500324/60‑1 (MQ=255)
                                                      cttcttAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCg                <  1:2272087/70‑1 (MQ=255)
                                                      cttcttAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCg                <  1:291429/70‑1 (MQ=255)
                                                       ttcttAAGATATAAACAGATAAATGCTGTAATGAGAAAGAAACTAATTACAAGATAGATCAtt                      >  1:1484256/1‑63 (MQ=255)
                                                        tcttAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAagat                              >  1:2439572/1‑54 (MQ=255)
                                                        tcttAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCa                        >  1:3919575/1‑60 (MQ=255)
                                                        tcttAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTAc             >  1:281515/1‑71 (MQ=255)
                                                         cttAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACt                                         <  1:2409068/42‑1 (MQ=255)
                                                         cttAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTa              <  1:1757058/69‑1 (MQ=255)
                                                         cttAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTACt            >  1:1733036/1‑71 (MQ=255)
                                                         cttAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTACt            >  1:2216033/1‑71 (MQ=255)
                                                            aaGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAAtt                                     <  1:2901256/43‑1 (MQ=255)
                                                            aaGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTACtga          >  1:3198769/1‑68 (MQ=255)
                                                               atataAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAagatag                            <  1:671615‑M2/44‑1 (MQ=255)
                                                                 ataAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAagata                             >  1:4421315‑M2/4‑46 (MQ=255)
                                                                  taAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTACTCATtata     <  1:1213527‑M2/67‑1 (MQ=255)
                                                                  taAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTACTCATTATAtt   >  1:1051049‑M2/3‑71 (MQ=255)
                                                                   aAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAaga                               >  1:4971064‑M2/2‑42 (MQ=255)
                                                                    aaCAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTACTCATTATAttt  <  1:4615808/70‑1 (MQ=255)
                                                                    aaCAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTACTCATTATAttt  <  1:4035925/70‑1 (MQ=255)
                                                                    aaCAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTACTCATTATAttt  <  1:1262759/70‑1 (MQ=255)
                                                                    aaCAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTACTCATTATAttt  <  1:729645/70‑1 (MQ=255)
                                                                    |                                                                     
CAATAAGACCAGTGCGAAGATGATATTAACGCATACGGCTGGATAAAATATATCCTTCTTAAGATATAAACAGATAAATGCTGTAATGAGAAAGACACTAATTACAAGATAGATCATTAATCCGTACTCATTATATTT  >  AP009048/2109247‑2109384

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.