New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4532499 =1 (0.050)22 (1.130) 17/152 0.2 88.2% intergenic (+49/‑433) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 4532813 = 5 (0.260)intergenic (+363/‑119) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

ATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTA  >  NZ_CP009273/4532413‑4532582
                                                                                      |                                                                                   
gacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgg                                                                                  <  2:219851‑M1/4‑1 (MQ=255)
gacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgg                                                                                  >  1:71600‑M1/87‑90 (MQ=255)
  catgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgggg                                                                                >  2:357792‑M1/85‑90 (MQ=255)
       gcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCaaa                                                                           <  2:271705‑M1/11‑1 (MQ=255)
        cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                          >  2:423324‑M1/79‑90 (MQ=255)
        cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                          >  1:219851‑M1/79‑90 (MQ=255)
        cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                          >  2:352499‑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                                                                >  2:315159‑M1/64‑85 (MQ=255)
                       ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATc                                                                <  1:315159‑M1/22‑1 (MQ=255)
                       ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaa                                                           <  1:393694‑M1/27‑1 (MQ=255)
                         ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGc                                                                               >  2:222686‑M1/62‑68 (MQ=255)
                         ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGc                                                                               <  1:222686‑M1/7‑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)
                                          ttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                    >  2:569331‑M1/45‑46 (MQ=255)
                                          ttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                    <  1:569331‑M1/2‑1 (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                         <  2:549928‑M1/90‑53 (MQ=255)
                                                         aGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgtt                         <  1:461915‑M1/90‑53 (MQ=255)
                                                              aaaaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                    <  2:203988‑M1/90‑58 (MQ=255)
                                                                aaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaa                  <  1:583509‑M1/90‑60 (MQ=255)
                                                                aaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaa                  <  2:405674‑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:43512‑M1/1‑15 (MQ=255)
                                                                                aGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatct  <  1:316880‑M1/90‑76 (MQ=255)
                                                                                      |                                                                                   
ATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTA  >  NZ_CP009273/4532413‑4532582

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