New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4532499 =29 (0.660)9 (0.230) 8/150 1.6 18.2% intergenic (+49/‑433) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 4532813 = 55 (1.400)intergenic (+363/‑119) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
Rejected: Frequency below/above cutoff threshold.

AAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAT  >  NZ_CP009273/4532727‑4532901
                                                                                      |                                                                                        
gttagatgctttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGt                                                                                       >  1:125352‑M2/87‑90 (MQ=255)
  aGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAATTGtct                                                                                     <  1:802430/90‑3 (MQ=255)
  aGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGttt                                                                                     >  1:1096848/1‑90 (MQ=255)
  aGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGttt                                                                                     <  1:684847/90‑1 (MQ=255)
    aTAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCa                                                                                   >  1:419883/1‑90 (MQ=255)
          aTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCtttt                                                                             <  1:1257530/90‑1 (MQ=255)
          aTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCtttt                                                                             <  2:880824/90‑1 (MQ=255)
          aTAGAATATTAAGCCAACAAATAAACTGAAAAAGCTTGTCGGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATTTTTCATCtttt                                                                             >  1:742684/1‑90 (MQ=255)
             gAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                                          >  2:523099/1‑90 (MQ=255)
             gAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                                          >  1:549692/1‑90 (MQ=255)
             gAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                                          >  1:459022/1‑90 (MQ=255)
               tcgttttttaatatttttatgcttgagaaaaaatacgtaacttgtttatgatatggacagtttggccccaaATGTTTCATCTTTTggggg                                                                        <  2:1205135‑M2/19‑1 (MQ=255)
               tcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTggggg                                                                        <  1:116722‑M2/19‑1 (MQ=255)
               atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggggg                                                                        <  1:285178/90‑1 (MQ=255)
               atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggggg                                                                        >  2:640177/1‑90 (MQ=255)
               atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggggg                                                                        >  1:355277/1‑90 (MQ=255)
               atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggggg                                                                        <  1:523099/90‑1 (MQ=255)
               atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggggg                                                                        <  2:1239511/90‑1 (MQ=255)
               atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggggg                                                                        >  1:181276/1‑90 (MQ=255)
                  ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                                                     >  1:22543/1‑90 (MQ=255)
                  ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                                                     >  2:433756/1‑90 (MQ=255)
                  ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                                                     >  2:434017/1‑90 (MQ=255)
                    ttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAAc                                                                   <  1:578102‑M2/24‑1 (MQ=255)
                    aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggggg                                                                        <  2:368660/85‑1 (MQ=255)
                    aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggggg                                                                        >  1:368660/1‑85 (MQ=255)
                    aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                                                   >  1:39982/1‑90 (MQ=255)
                    aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                                                   >  2:930762/1‑90 (MQ=255)
                     aGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACt                                                                  >  1:392871/1‑90 (MQ=255)
                        caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgtg                                                               <  2:1138334/90‑1 (MQ=255)
                        caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgtg                                                               <  2:1096848/90‑1 (MQ=255)
                        caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgtg                                                               >  2:454953/1‑90 (MQ=255)
                        caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgtg                                                               <  2:894763/90‑1 (MQ=255)
                        caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgtg                                                               <  1:433756/90‑1 (MQ=255)
                         aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                                                              >  2:557220/1‑90 (MQ=255)
                         aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                                                              >  1:402425/1‑90 (MQ=255)
                         aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAATGTGc                                                              >  2:335590/1‑90 (MQ=255)
                          acaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCa                                                             <  1:70841/90‑1 (MQ=255)
                          acaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCa                                                             <  1:316961/90‑1 (MQ=255)
                          acaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCa                                                             <  1:166435/90‑1 (MQ=255)
                             aaTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                                          >  2:1144303/1‑90 (MQ=255)
                                    tGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAg                                                   >  1:822284/1‑90 (MQ=255)
                                     ttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                                                  >  1:429990‑M2/50‑90 (MQ=255)
                                     gAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                                                  >  2:846325/1‑90 (MQ=255)
                                     gAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                                                  >  2:872369/1‑90 (MQ=255)
                                                      aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa                                 <  2:918501/90‑1 (MQ=255)
                                                         tatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTttacaaaacga                              >  1:509611‑M2/30‑88 (MQ=255)
                                                           tttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                            <  2:412142‑M2/63‑1 (MQ=255)
                                                           ttCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                            <  1:323166/90‑1 (MQ=255)
                                                                ctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTAc                       >  2:637751/1‑90 (MQ=255)
                                                                ctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTAc                       >  2:637087/1‑90 (MQ=255)
                                                                ctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTAc                       <  2:503742/90‑1 (MQ=255)
                                                                ctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTAc                       >  1:558162/1‑90 (MQ=255)
                                                                  aTGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAg                     >  2:490786/1‑90 (MQ=255)
                                                                   tGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGa                    <  2:770388/90‑1 (MQ=255)
                                                                    tggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa                   >  1:416625‑M2/19‑90 (MQ=255)
                                                                    gAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa                   >  1:1112316/1‑90 (MQ=255)
                                                                     aGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAAc                  >  1:669427/1‑90 (MQ=255)
                                                                       tCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                                           >  2:1161579/1‑47 (MQ=255)
                                                                       tCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgt                                                           <  1:1161579/47‑1 (MQ=255)
                                                                       tCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGa                >  2:4493/1‑90 (MQ=255)
                                                                       tCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGa                >  1:847044/1‑90 (MQ=255)
                                                                       tCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGa                >  1:217274/1‑90 (MQ=255)
                                                                        cAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGAc               <  2:558162/90‑1 (MQ=255)
                                                                          gtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTg             >  2:156636‑M2/13‑90 (MQ=255)
                                                                               gccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTGCAAAACAAAGTGTACAGAACGACTGCCCAt        >  2:1213356‑M2/8‑90 (MQ=255)
                                                                                    aaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGa   >  2:531621‑M2/3‑90 (MQ=255)
                                                                                     aATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAt  <  2:967615‑M2/89‑1 (MQ=255)
                                                                                      |                                                                                        
AAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAT  >  NZ_CP009273/4532727‑4532901

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