New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4532499 =5 (0.090)25 (0.530) 17/148 1.0 50.8% intergenic (+49/‑433) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 4532813 = 44 (0.930)intergenic (+363/‑119) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

TAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGT  >  NZ_CP009273/4532411‑4532587
                                                                                        |                                                                                        
tcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                         >  1:6637‑M1/89‑90 (MQ=255)
   acatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTggg                                                                                      <  2:1172400‑M1/5‑1 (MQ=255)
     tcATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                                                    <  2:988251‑M1/90‑1 (MQ=255)
     tcATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                                                    <  2:737006‑M1/90‑1 (MQ=255)
     atgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGc                                                                                    >  2:1164079‑M1/84‑90 (MQ=255)
      cATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCC                                                                                   >  2:192367‑M1/1‑90 (MQ=255)
          cagtcgttctgtacactttgttttgtgaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                               >  2:27559‑M1/79‑90 (MQ=255)
          cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                               >  1:775111‑M1/79‑90 (MQ=255)
                aaaTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCa                                                                         >  2:43518/1‑90 (MQ=255)
                 aaTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCat                                                                        <  1:1568717/90‑1 (MQ=255)
                  aTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcat                                                                       <  2:116788‑M1/90‑12 (MQ=255)
                   tgtaccctttgttttgtaaacgagtttgactgccaacactgcacagtttgcccccaaaagatgaaacatTTGGGGCCAAACTGTCCatat                                                                      >  2:893008‑M1/70‑90 (MQ=255)
                   tgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCatat                                                                      >  1:455867‑M1/70‑90 (MQ=255)
                         ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaa                                                                >  1:1270631‑M1/64‑90 (MQ=255)
                              ttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                         <  1:1248029‑M1/2‑1 (MQ=255)
                              ttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                         >  2:1248029‑M1/59‑60 (MQ=255)
                              ttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCaa                                                                                 <  1:1534815‑M1/10‑1 (MQ=255)
                              ttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCaa                                                                                 >  2:1534815‑M1/59‑68 (MQ=255)
                              ttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGt                                                           <  1:248670‑M1/32‑1 (MQ=255)
                                  gTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACg                                                       <  2:285538/90‑1 (MQ=255)
                                       cgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAttt                                                  <  2:1258086‑M1/41‑1 (MQ=255)
                                          gtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                               >  1:444014‑M1/47‑58 (MQ=255)
                                          gtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                               <  2:444014‑M1/12‑1 (MQ=255)
                                              gactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa                                           >  1:1310661‑M1/43‑90 (MQ=255)
                                              gactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa                                           >  1:630204‑M1/43‑90 (MQ=255)
                                               tATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcgtacaaactttttc                                          <  2:1245696‑M1/90‑41 (MQ=255)
                                               tATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttc                                          <  2:81977‑M1/90‑41 (MQ=255)
                                               tATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttc                                          <  2:1111469‑M1/90‑41 (MQ=255)
                                                  gccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCa                                       <  2:667023‑M1/52‑1 (MQ=255)
                                                      acactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaa                                   >  2:861330‑M1/35‑90 (MQ=255)
                                                           gcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat                              <  2:679610‑M1/61‑1 (MQ=255)
                                                            cacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa                                           >  2:421492‑M1/29‑76 (MQ=255)
                                                            cacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa                                           <  1:421492‑M1/48‑1 (MQ=255)
                                                             ttAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttgg                            <  1:260280‑M1/90‑55 (MQ=255)
                                                             ttAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTa                            >  1:1526957/1‑90 (MQ=255)
                                                                aaaaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                         <  2:1354583‑M1/90‑58 (MQ=255)
                                                                aaaaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                         <  2:928857‑M1/90‑58 (MQ=255)
                                                                  aaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaa                       <  2:595110‑M1/90‑60 (MQ=255)
                                                                          aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa               >  2:1236572‑M1/15‑90 (MQ=255)
                                                                          aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa               <  2:1393513‑M1/76‑1 (MQ=255)
                                                                          aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa               <  1:1042019‑M1/76‑1 (MQ=255)
                                                                          aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa               >  1:1169913‑M1/15‑90 (MQ=255)
                                                                            aagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAg             >  1:1393513‑M1/13‑90 (MQ=255)
                                                                                   gACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctt      <  1:543267‑M1/90‑77 (MQ=255)
                                                                                    acatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAAc     >  1:481030‑M1/5‑90 (MQ=255)
                                                                                     catTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACt    >  2:472702‑M1/4‑90 (MQ=255)
                                                                                     cAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatcttta    >  1:822325‑M1/1‑12 (MQ=255)
                                                                                       aTTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatt  >  2:152501‑M1/1‑10 (MQ=255)
                                                                                       aTTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatt  >  1:1054335‑M1/1‑10 (MQ=255)
                                                                                       aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGt  >  2:1296759/1‑90 (MQ=255)
                                                                                        |                                                                                        
TAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGT  >  NZ_CP009273/4532411‑4532587

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