New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 45325071 (0.050)18 (0.920) 12/152 0.4 85.9% intergenic (+57/‑425) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 = 4532803 5 (0.260)intergenic (+353/‑129) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

AAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGT  >  NZ_CP009273/4532420‑4532596
                                                                                       |                                                                                         
gcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCaaa                                                                                         <  2:271705‑M1/11‑1 (MQ=255)
 cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                        >  2:423324‑M1/79‑90 (MQ=255)
 cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                        >  2:352499‑M1/79‑90 (MQ=255)
 cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                        >  1:219851‑M1/79‑90 (MQ=255)
  agtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACt                                                                                       >  1:366847‑M1/78‑90 (MQ=255)
   gtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTg                                                                                      <  2:69051‑M1/14‑1 (MQ=255)
      gttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTcc                                                                                   >  2:347922‑M1/74‑90 (MQ=255)
        tctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCat                                                                                 <  2:122415‑M1/19‑1 (MQ=255)
         ataaaaaagagAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcat                                                                                <  2:360592‑M1/87‑12 (MQ=255)
         aTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcat                                                                                <  1:122743‑M1/90‑12 (MQ=255)
         aTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGACAttttgactcat                                                                                <  1:130108‑M1/90‑12 (MQ=255)
          tgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCatat                                                                               >  1:122415‑M1/70‑90 (MQ=255)
           tAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatag                                                                              <  2:236663‑M1/90‑14 (MQ=255)
           gtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATc                                                                              >  2:349380‑M1/69‑90 (MQ=255)
               actttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataa                                                                          >  1:602363‑M1/65‑90 (MQ=255)
                ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATc                                                                              <  1:315159‑M1/22‑1 (MQ=255)
                ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATc                                                                              >  2:315159‑M1/64‑85 (MQ=255)
                ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaa                                                                         <  1:393694‑M1/27‑1 (MQ=255)
                  ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCa                                                                                  >  2:226095‑M1/62‑79 (MQ=255)
                  ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCa                                                                                  <  1:226095‑M1/18‑1 (MQ=255)
                          tttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGt                                                               <  2:319147/90‑1 (MQ=255)
                           ttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcgga                                                              >  2:188529‑M1/1‑61 (MQ=255)
                           ttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcgga                                                              >  2:436850‑M1/1‑61 (MQ=255)
                               tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaa                                                          <  2:49269‑M1/90‑34 (MQ=255)
                                      tATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttc                                                   >  2:164168‑M1/1‑50 (MQ=255)
                                                 tgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAata                                        >  2:432581‑M1/31‑90 (MQ=255)
                                                  aGTTAAAAAAACAAATTAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgtt                                       <  1:148951‑M1/90‑53 (MQ=255)
                                                  aGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgtt                                       <  1:461915‑M1/90‑53 (MQ=255)
                                                  aGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgtt                                       <  2:549928‑M1/90‑53 (MQ=255)
                                                       aaaaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                                  <  2:203988‑M1/90‑58 (MQ=255)
                                                         aaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaa                                <  2:405674‑M1/90‑60 (MQ=255)
                                                         aaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaa                                <  1:583509‑M1/90‑60 (MQ=255)
                                                            cccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGa                             <  2:184991‑M1/71‑1 (MQ=255)
                                                            cccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGa                             <  2:494522‑M1/71‑1 (MQ=255)
                                                            aCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatat                             <  2:428814‑M1/90‑63 (MQ=255)
                                                              aaataaatACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattc                           <  2:323887‑M1/90‑65 (MQ=255)
                                                                 aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                        >  1:483030‑M1/15‑90 (MQ=255)
                                                                   aagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAg                      >  1:158276‑M1/13‑90 (MQ=255)
                                                                         aaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa                >  1:293984‑M1/7‑90 (MQ=255)
                                                                         aGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatct                <  1:316880‑M1/90‑76 (MQ=255)
                                                                         aGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatct                >  1:43512‑M1/1‑15 (MQ=255)
                                                                                   ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatag      <  2:392334‑M1/90‑86 (MQ=255)
                                                                                   ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatag      <  1:213768‑M1/90‑86 (MQ=255)
                                                                                   ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatag      <  1:116932‑M1/90‑86 (MQ=255)
                                                                                   ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgtt                                       >  1:603682‑M1/1‑5 (MQ=255)
                                                                                   ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgtt                                       <  2:603682‑M1/57‑53 (MQ=255)
                                                                                      cAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatg   >  1:214365‑M1/1‑2 (MQ=255)
                                                                                       AAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGt  <  1:97167‑M1/90‑1 (MQ=255)
                                                                                       |                                                                                         
AAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGT  >  NZ_CP009273/4532420‑4532596

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