New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4532499 =4 (0.130)21 (0.790) 16/148 0.4 83.1% intergenic (+49/‑433) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 4532813 = 5 (0.190)intergenic (+363/‑119) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

ATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGT  >  NZ_CP009273/4532410‑4532587
                                                                                         |                                                                                        
atcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatt                                                                                          <  1:102261‑M1/1‑1 (MQ=255)
atcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatt                                                                                          <  1:407372‑M1/1‑1 (MQ=255)
atcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatt                                                                                          <  2:375879‑M1/1‑1 (MQ=255)
 tcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                         >  2:102261‑M1/89‑90 (MQ=255)
 tcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                         >  2:85817‑M1/89‑90 (MQ=255)
 tcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                         >  2:482326‑M1/89‑90 (MQ=255)
   atCTCATAAAAGAAAAATTAAAAAGAGAAAAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAAACAATTgg                                                                                       >  1:845309/1‑90 (MQ=255)
     catgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgggg                                                                                     >  2:629937‑M1/85‑90 (MQ=255)
      tcATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGc                                                                                    <  2:350875/90‑1 (MQ=255)
      atgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGc                                                                                    >  2:64620‑M1/84‑90 (MQ=255)
          gcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCaaa                                                                                >  2:185088‑M1/80‑90 (MQ=255)
            agtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACt                                                                              <  2:38572‑M1/13‑1 (MQ=255)
             gtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCaaa                                                                                <  1:504115‑M1/11‑1 (MQ=255)
             gtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCaaa                                                                                >  2:504115‑M1/77‑87 (MQ=255)
             gtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTg                                                                             >  1:47839‑M1/77‑90 (MQ=255)
             gtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTg                                                                             >  1:199947‑M1/77‑90 (MQ=255)
                    ttAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCatat                                                                      <  1:754864/90‑1 (MQ=255)
                     gtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATc                                                                     <  1:24975‑M1/22‑1 (MQ=255)
                       acactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGc                                                                                    >  2:515630‑M1/67‑73 (MQ=255)
                       acactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGc                                                                                    <  1:515630‑M1/7‑1 (MQ=255)
                            ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaata                                                              >  2:530932‑M1/62‑90 (MQ=255)
                            ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaata                                                              <  2:63331‑M1/29‑1 (MQ=255)
                               ttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGt                                                           <  1:575842‑M1/32‑1 (MQ=255)
                               ttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGt                                                           >  1:524634‑M1/59‑90 (MQ=255)
                               ttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGt                                                           <  2:855069‑M1/32‑1 (MQ=255)
                                  tgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTAc                                                        >  2:24975‑M1/56‑90 (MQ=255)
                                                             cacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATAtt                             >  1:14872‑M1/29‑90 (MQ=255)
                                                                 gttttcccccaaaagatgaaacatTTGGGGCCAAACTg                                                                             >  1:726646‑M1/25‑38 (MQ=255)
                                                                 gttttcccccaaaagatgaaacatTTGGGGCCAAACTg                                                                             <  2:726646‑M1/14‑1 (MQ=255)
                                                                 gttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaaa                         <  1:273725‑M1/66‑1 (MQ=255)
                                                                       ccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGAc                   >  1:320215‑M1/19‑90 (MQ=255)
                                                                        cccaaaagatgaaacatTTGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaa                 >  1:37168‑M1/18‑90 (MQ=255)
                                                                           aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa               >  2:92861‑M1/15‑90 (MQ=255)
                                                                                atgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCAt          >  1:521872‑M1/10‑90 (MQ=255)
                                                                                        aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGt  >  2:36156/1‑90 (MQ=255)
                                                                                         |                                                                                        
ATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGT  >  NZ_CP009273/4532410‑4532587

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