New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 45325075 (0.300)5 (0.340) 5/148 0.8 48.9% intergenic (+57/‑425) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 = 4532803 6 (0.410)intergenic (+353/‑129) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

AACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATG  >  NZ_CP009273/4532421‑4532595
                                                                                      |                                                                                        
cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                       >  1:237481‑M1/79‑90 (MQ=255)
       tctgtacactttgttttgtaaacgagtttgactgccaacactgcccagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCat                                                                                >  2:482185‑M1/72‑90 (MQ=255)
         ttAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCatat                                                                              <  2:52130/90‑1 (MQ=255)
               ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTcc                                                                                  >  2:415289‑M1/64‑80 (MQ=255)
               ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTcc                                                                                  <  1:415289‑M1/17‑1 (MQ=255)
               ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaa                                                                        <  1:245443‑M1/27‑1 (MQ=255)
                                    gactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa                                                   >  2:136361‑M1/43‑90 (MQ=255)
                                     tATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttc                                                  <  1:65907‑M1/90‑41 (MQ=255)
                                             ataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                                                             >  1:409285/1‑55 (MQ=255)
                                             ataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                                                             <  2:409285/55‑1 (MQ=255)
                                                aaGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgt                                       <  1:302250‑M1/90‑52 (MQ=255)
                                                 aGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgtt                                      <  2:13806‑M1/90‑53 (MQ=255)
                                                                aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                       >  1:12302‑M1/15‑90 (MQ=255)
                                                                aaaagatgaaacatTTGGGGCCAAACTGTCCATATCACAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                       >  2:245443‑M1/15‑90 (MQ=255)
                                                                        aGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa               >  2:270826/1‑90 (MQ=255)
                                                                             aTTGGGGCCAttttgactcatagaggaaagcatcgcggac                                                            <  1:484448‑M1/40‑31 (MQ=255)
                                                                             aTTGGGGCCAttttgactcatagaggaaagcatcgcggac                                                            >  2:484448‑M1/1‑10 (MQ=255)
                                                                             aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGt          >  2:360914/1‑90 (MQ=255)
                                                                              ttGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttattt         >  1:360914‑M1/1‑9 (MQ=255)
                                                                                  ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattc                          >  2:335429‑M1/1‑5 (MQ=255)
                                                                                  ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattc                          <  1:335429‑M1/69‑65 (MQ=255)
                                                                                  GGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGat     <  1:335300‑M1/90‑1 (MQ=255)
                                                                                     cAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatg  >  2:83036‑M1/1‑2 (MQ=255)
                                                                                      |                                                                                        
AACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATG  >  NZ_CP009273/4532421‑4532595

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