New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4532499 =2 (0.030)7 (0.120) 6/148 2.6 29.9% intergenic (+49/‑433) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 4532813 = 31 (0.540)intergenic (+363/‑119) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

AGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAT  >  NZ_CP009273/4532729‑4532901
                                                                                    |                                                                                        
aGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGttt                                                                                     >  1:978437/1‑90 (MQ=255)
          aGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTgg                                                                           <  2:978437/90‑1 (MQ=255)
           gAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                                          >  2:71967/1‑90 (MQ=255)
             atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggggg                                                                        >  2:530613/1‑90 (MQ=255)
                ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAGAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                                                     >  1:1445732/1‑90 (MQ=255)
                ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                                                     >  1:1285265/1‑90 (MQ=255)
                ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                                                     >  2:1501123/1‑90 (MQ=255)
                ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                                                     >  1:1364670/1‑90 (MQ=255)
                  aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                                                   >  2:264425/1‑90 (MQ=255)
                  aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAc                                                                   <  2:1285265/90‑1 (MQ=255)
                      caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgtg                                                               <  2:448364/90‑1 (MQ=255)
                       aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                                                              >  1:1708117/1‑90 (MQ=255)
                       aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                                                              >  1:547312/1‑90 (MQ=255)
                                aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGc                                                     >  1:1776821/1‑90 (MQ=255)
                                 gcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCa                                                    >  2:593251‑M2/52‑90 (MQ=255)
                                   gAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                                                  >  1:644673/1‑90 (MQ=255)
                                   gAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                                                  >  2:160636/1‑90 (MQ=255)
                                   gAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                                                  >  1:1449235/1‑90 (MQ=255)
                                        aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAAc                                             <  2:679121/90‑1 (MQ=255)
                                                 gcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTaca                                    <  2:1160150/90‑1 (MQ=255)
                                                       tatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa                              <  1:593251‑M2/61‑1 (MQ=255)
                                                         tttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                            <  1:580388‑M2/63‑1 (MQ=255)
                                                         ttCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                            <  1:1063902/90‑1 (MQ=255)
                                                         ttCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                            >  2:1063902/1‑90 (MQ=255)
                                                         ttCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                            <  1:530613/90‑1 (MQ=255)
                                                              gatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTAc                       >  1:865283‑M2/23‑90 (MQ=255)
                                                               atatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACa                      >  1:1337698‑M2/22‑90 (MQ=255)
                                                                aTGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa                              <  1:241530/81‑1 (MQ=255)
                                                                aTGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa                              >  2:241530/1‑81 (MQ=255)
                                                                aTGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAg                     >  1:860784/1‑90 (MQ=255)
                                                                   aGTCAAAGTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                            <  1:928316/80‑1 (MQ=255)
                                                                   aGTCAAAGTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                            >  2:928316/1‑80 (MQ=255)
                                                                      cAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCa                                                             >  2:721592/1‑44 (MQ=255)
                                                                      cAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCa                                                             <  1:721592/44‑1 (MQ=255)
                                                                        gtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTg             >  2:1581854‑M2/13‑90 (MQ=255)
                                                                        gtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTg             >  2:258761‑M2/13‑90 (MQ=255)
                                                                        aaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTa                                      >  1:1698106/1‑65 (MQ=255)
                                                                        aaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTa                                      <  2:1698106/65‑1 (MQ=255)
                                                                             gccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCAt        >  2:209545‑M2/8‑90 (MQ=255)
                                                                             gccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCAt        >  1:427219‑M2/8‑90 (MQ=255)
                                                                             gccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCAt        >  1:164829‑M2/8‑90 (MQ=255)
                                                                                ccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTc     >  2:1346937‑M2/5‑90 (MQ=255)
                                                                                  aaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                            <  2:149893‑M2/63‑1 (MQ=255)
                                                                                  aaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                            >  1:149893‑M2/3‑65 (MQ=255)
                                                                                  aaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGa   >  2:499346‑M2/3‑90 (MQ=255)
                                                                                   aATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAt  <  2:164829‑M2/89‑1 (MQ=255)
                                                                                    |                                                                                        
AGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAT  >  NZ_CP009273/4532729‑4532901

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