New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 45325070 (0.000)10 (0.390) 8/146 1.1 57.1% intergenic (+57/‑425) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 = 4532803 15 (0.580)intergenic (+353/‑129) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

TAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATG  >  NZ_CP009273/4532419‑4532595
                                                                                        |                                                                                        
ggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCaa                                                                                         >  2:838245‑M1/81‑90 (MQ=255)
ggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCaa                                                                                         >  2:689768‑M1/81‑90 (MQ=255)
  cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                       >  1:816950‑M1/79‑90 (MQ=255)
  cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                       >  1:17858‑M1/79‑90 (MQ=255)
        ttctgtacactttgttttttaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCa                                                                                 >  1:412797‑M1/73‑90 (MQ=255)
          aTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcat                                                                               <  1:403035‑M1/90‑12 (MQ=255)
          aTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcat                                                                               <  2:592682‑M1/90‑12 (MQ=255)
           tgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCatat                                                                              >  1:693470‑M1/70‑90 (MQ=255)
            gtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATc                                                                             >  1:550343‑M1/69‑90 (MQ=255)
                 ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                       <  2:823009‑M1/12‑1 (MQ=255)
                 ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                       >  1:823009‑M1/64‑75 (MQ=255)
                 ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaa                                                                        <  2:747552‑M1/27‑1 (MQ=255)
                        ttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTa                                                                 >  1:111385‑M1/57‑90 (MQ=255)
                        agGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcg                                                                 <  1:871611‑M1/90‑26 (MQ=255)
                         tgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTAc                                                                >  1:704054‑M1/56‑90 (MQ=255)
                               cgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCata                                                                               >  2:433516‑M1/50‑69 (MQ=255)
                               cgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCata                                                                               <  1:433516‑M1/20‑1 (MQ=255)
                                     tgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtctc                                                    >  1:747552‑M1/44‑90 (MQ=255)
                                       tATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttc                                                  <  2:285761‑M1/90‑41 (MQ=255)
                                        aTTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttca                                                 >  1:691206‑M1/1‑49 (MQ=255)
                                              aataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttat                                           >  2:265842‑M1/1‑43 (MQ=255)
                                                        aaaaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                                 <  2:415983‑M1/90‑58 (MQ=255)
                                                             cccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGa                            <  2:629958‑M1/71‑1 (MQ=255)
                                                                  aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                       >  1:798010‑M1/15‑90 (MQ=255)
                                                                    aataCAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattg                     <  1:133635‑M1/90‑70 (MQ=255)
                                                                               aTTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatt          >  1:201254‑M1/1‑10 (MQ=255)
                                                                               aTTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatt          >  1:65783‑M1/1‑10 (MQ=255)
                                                                                       cAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatg  >  2:325780‑M1/1‑2 (MQ=255)
                                                                                       cAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatg  >  2:662393‑M1/1‑2 (MQ=255)
                                                                                        |                                                                                        
TAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATG  >  NZ_CP009273/4532419‑4532595

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