New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 453250712 (0.390)17 (0.610) 13/150 0.6 43.2% intergenic (+57/‑425) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 = 4532803 34 (1.230)intergenic (+353/‑129) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

ATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATG  >  NZ_CP009273/4532418‑4532595
                                                                                         |                                                                                        
gggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCa                                                                                          <  2:356916‑M1/9‑1 (MQ=255)
   cagtcgttctgtaccctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                       >  1:210590‑M1/79‑90 (MQ=255)
   cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                       >  2:90267‑M1/79‑90 (MQ=255)
   cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                       >  1:493483‑M1/79‑90 (MQ=255)
   cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                       >  2:513726‑M1/79‑90 (MQ=255)
   aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGATAATTGGGGCCAAAc                                                                                       <  2:595516/90‑1 (MQ=255)
   aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                                                       <  2:540236/90‑1 (MQ=255)
   aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                                                       <  1:711368/90‑1 (MQ=255)
       aaaaaTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgac                                                                                   >  2:686685‑M1/1‑83 (MQ=255)
            tgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCatat                                                                              >  2:842014‑M1/70‑90 (MQ=255)
             gtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATc                                                                             >  2:1107‑M1/69‑90 (MQ=255)
                  ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                       >  1:164147‑M1/64‑75 (MQ=255)
                  ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                       <  2:164147‑M1/12‑1 (MQ=255)
                  ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaa                                                                        <  1:320714‑M1/27‑1 (MQ=255)
                    ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaata                                                                      <  1:138237‑M1/29‑1 (MQ=255)
                       agagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGt                                                                   <  1:600741/90‑1 (MQ=255)
                                    aaCTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactt                                                      >  1:752421‑M1/1‑54 (MQ=255)
                                       gactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa                                                   >  2:451352‑M1/43‑90 (MQ=255)
                                        tATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttc                                                  <  1:649183‑M1/90‑41 (MQ=255)
                                                ataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatt                                          >  2:545855‑M1/1‑42 (MQ=255)
                                                   tgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAata                                       <  1:545855‑M1/60‑1 (MQ=255)
                                                    aGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgtt                                      <  2:392579‑M1/90‑53 (MQ=255)
                                                         gttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                       >  1:780318‑M1/25‑36 (MQ=255)
                                                         gttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                       <  2:780318‑M1/12‑1 (MQ=255)
                                                         aaaaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                                 <  1:39481‑M1/90‑58 (MQ=255)
                                                           aaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaa                               <  2:752421‑M1/90‑60 (MQ=255)
                                                               cAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGAc                           >  2:106561/1‑90 (MQ=255)
                                                                   taaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                       <  1:831639/90‑1 (MQ=255)
                                                                   taaataCAAGACAATTGGGGCCAAACTGCCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                       <  1:672983/90‑1 (MQ=255)
                                                                   aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                       >  2:448534‑M1/15‑90 (MQ=255)
                                                                   aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                       >  2:13002‑M1/15‑90 (MQ=255)
                                                                   aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                       >  1:666617‑M1/15‑90 (MQ=255)
                                                                         cAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttat                 <  1:190987‑M1/90‑74 (MQ=255)
                                                                         cAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATc                 >  2:716375/1‑90 (MQ=255)
                                                                           aGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                                 <  1:352208‑M1/72‑58 (MQ=255)
                                                                           aGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                                 >  2:352208‑M1/1‑15 (MQ=255)
                                                                           aGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa               >  2:408237/1‑90 (MQ=255)
                                                                            gACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctt              <  1:748846‑M1/90‑77 (MQ=255)
                                                                            gACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTaa              <  1:644536/90‑1 (MQ=255)
                                                                            aacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTaa              >  1:71349‑M1/6‑90 (MQ=255)
                                                                            aacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTaa              >  1:823499‑M1/6‑90 (MQ=255)
                                                                              cAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatcttta            <  2:763997‑M1/90‑79 (MQ=255)
                                                                              cAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                                 >  1:595400‑M1/1‑12 (MQ=255)
                                                                              cAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                                 <  2:595400‑M1/69‑58 (MQ=255)
                                                                                aTTGGGGCCAttttgactcatagaggaaagtatagcggacaaactttttgaggttatgtgttggcttaatatgatattgttatatttatt          >  1:137101‑M1/1‑10 (MQ=255)
                                                                                aTTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttc                                                  >  2:599994‑M1/1‑10 (MQ=255)
                                                                                aTTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttc                                                  <  1:599994‑M1/50‑41 (MQ=255)
                                                                                aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGt          >  2:59696/1‑90 (MQ=255)
                                                                                aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGt          >  1:875156/1‑90 (MQ=255)
                                                                                   ggggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttat       <  2:342949‑M1/90‑84 (MQ=255)
                                                                                   ggggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttat       <  2:717255‑M1/90‑84 (MQ=255)
                                                                                   ggggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttat       <  2:137101‑M1/90‑84 (MQ=255)
                                                                                     ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatag     >  1:491373‑M1/1‑5 (MQ=255)
                                                                                     ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatag     <  2:374191‑M1/90‑86 (MQ=255)
                                                                                     ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatag     <  2:183772‑M1/90‑86 (MQ=255)
                                                                                     ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttat                 >  1:895306‑M1/1‑5 (MQ=255)
                                                                                     ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttat                 <  2:895306‑M1/78‑74 (MQ=255)
                                                                                       ccAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagat   <  2:868348‑M1/90‑88 (MQ=255)
                                                                                        cAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatg  >  2:271268‑M1/1‑2 (MQ=255)
                                                                                        cAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaa                               <  2:570145‑M1/61‑60 (MQ=255)
                                                                                        cAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaa                               >  1:570145‑M1/1‑2 (MQ=255)
                                                                                         |                                                                                        
ATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATG  >  NZ_CP009273/4532418‑4532595

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ 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.