New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4532499 =21 (0.340)13 (0.230) 9/144 2.1 23.7% intergenic (+49/‑433) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 4532813 = 65 (1.170)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:227990‑M1/89‑90 (MQ=255)
taatCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATT                                                                                         <  1:1495987‑M1/90‑1 (MQ=255)
 aatCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTG                                                                                        >  2:832567‑M1/1‑90 (MQ=255)
  atCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGG                                                                                       >  2:979068‑M1/1‑90 (MQ=255)
  atCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGG                                                                                       >  2:88351‑M1/1‑90 (MQ=255)
     tcATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                                                    <  2:395781‑M1/90‑1 (MQ=255)
      cATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCC                                                                                   <  1:1792103‑M1/90‑1 (MQ=255)
       gggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCa                                                                                  >  1:640302‑M1/82‑90 (MQ=255)
          cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                               >  1:1629971‑M1/79‑90 (MQ=255)
          aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                                               <  2:1280349/90‑1 (MQ=255)
          aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                                               <  2:479818/90‑1 (MQ=255)
          aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                                               <  2:532216/90‑1 (MQ=255)
          aaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                                               <  2:140266/90‑1 (MQ=255)
            cGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttg                                                                             <  1:1309490‑M1/90‑6 (MQ=255)
               aaaaTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAATACAATTGGGGCCAAACTGTcc                                                                          <  1:44874/90‑1 (MQ=255)
                ttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCa                                                                         >  1:1816017‑M1/73‑90 (MQ=255)
                ttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCa                                                                         >  1:791925‑M1/73‑90 (MQ=255)
                   tgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCatat                                                                      >  2:1240832‑M1/70‑90 (MQ=255)
                         ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGc                                                                                    >  2:175982‑M1/64‑70 (MQ=255)
                         ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGc                                                                                    <  1:175982‑M1/7‑1 (MQ=255)
                           ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCa                                                                         >  1:1094539‑M1/62‑79 (MQ=255)
                           ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCa                                                                         <  2:1094540‑M1/18‑1 (MQ=255)
                           ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaata                                                              <  2:401884‑M1/29‑1 (MQ=255)
                           agaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaata                                                              <  1:852135/90‑1 (MQ=255)
                            tgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaataa                                                             <  1:273844‑M1/30‑1 (MQ=255)
                            tgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaataa                                                             >  2:259298‑M1/61‑90 (MQ=255)
                                    ttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTa                                                     >  2:1719653/1‑90 (MQ=255)
                                      gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggaca                                                   >  2:940534‑M1/1‑59 (MQ=255)
                                           aaCTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtc                                              >  2:1154914/1‑90 (MQ=255)
                                            ttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgggg                                                                                     >  2:767938‑M1/45‑50 (MQ=255)
                                            ttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgggg                                                                                     <  1:767938‑M1/6‑1 (MQ=255)
                                              ttATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTCACGTATTTTTTCTCa                                           <  2:921248/90‑1 (MQ=255)
                                                    ataataAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagttt                                     <  1:889925‑M1/90‑46 (MQ=255)
                                                    ataataAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagttt                                     <  1:1375146‑M1/90‑46 (MQ=255)
                                                      aataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaa                                   <  2:984564/90‑1 (MQ=255)
                                                            cacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATAtt                             <  1:996394‑M1/62‑1 (MQ=255)
                                                             ttAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTa                            >  2:1227071/1‑90 (MQ=255)
                                                                aaaaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                         <  1:1496270‑M1/90‑58 (MQ=255)
                                                                aaaaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                         <  1:940533‑M1/90‑58 (MQ=255)
                                                                       aaataaatAAAAGACAATTGGGGACAAACCCTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACa                  <  1:661369/90‑1 (MQ=255)
                                                                          aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa               >  2:1869243‑M1/15‑90 (MQ=255)
                                                                                  aaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa       >  1:1044747‑M1/7‑90 (MQ=255)
                                                                                  aaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa       >  1:1325573‑M1/7‑90 (MQ=255)
                                                                                  aGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatct       >  2:1444079‑M1/1‑15 (MQ=255)
                                                                                  aGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa       >  1:737829/1‑90 (MQ=255)
                                                                                  aGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa       >  1:1665245/1‑90 (MQ=255)
                                                                                  aGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa       >  2:83414/1‑90 (MQ=255)
                                                                                  aGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa       >  2:1716660/1‑90 (MQ=255)
                                                                                       aTTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatt  >  1:735267‑M1/1‑10 (MQ=255)
                                                                                       aTTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatt  >  1:1040752‑M1/1‑10 (MQ=255)
                                                                                       aTTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatt  >  1:639728‑M1/1‑10 (MQ=255)
                                                                                       aTTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatt  >  2:106938‑M1/1‑10 (MQ=255)
                                                                                       aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAg             <  1:1715955/79‑1 (MQ=255)
                                                                                       aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAg             >  2:1715956/1‑79 (MQ=255)
                                                                                       aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGt  >  2:1109982/1‑90 (MQ=255)
                                                                                       aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGt  >  1:1390131/1‑90 (MQ=255)
                                                                                       aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGt  >  2:1429769/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.