New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4532499 =7 (0.300)16 (0.770) 14/152 0.3 62.4% intergenic (+49/‑433) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 4532813 = 13 (0.630)intergenic (+363/‑119) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

ATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTT  >  NZ_CP009273/4532410‑4532588
                                                                                         |                                                                                         
ataatCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGTTTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT                                                                                           <  1:47378‑M1/90‑1 (MQ=255)
 tcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                                          >  1:137284‑M1/89‑90 (MQ=255)
 taatCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATT                                                                                          <  2:450751‑M1/90‑1 (MQ=255)
   atCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGG                                                                                        <  2:31251‑M1/90‑1 (MQ=255)
      tcATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGC                                                                                     <  2:499714‑M1/90‑1 (MQ=255)
       cATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCC                                                                                    >  2:52043‑M1/1‑90 (MQ=255)
             gtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTg                                                                              <  1:633820‑M1/14‑1 (MQ=255)
             gtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTg                                                                              <  1:183388‑M1/14‑1 (MQ=255)
                 ttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCa                                                                          >  1:468661‑M1/73‑90 (MQ=255)
                   ctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCata                                                                        <  1:527034‑M1/20‑1 (MQ=255)
                     gtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATc                                                                      >  1:410837‑M1/69‑90 (MQ=255)
                       acactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCat                                                                    <  2:563735‑M1/24‑1 (MQ=255)
                          ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaa                                                                 >  2:415555‑M1/64‑90 (MQ=255)
                          ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaa                                                                 <  2:274843‑M1/27‑1 (MQ=255)
                           tttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaat                                                                <  1:197426‑M1/28‑1 (MQ=255)
                            ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaata                                                               >  1:274843‑M1/62‑90 (MQ=255)
                                    taaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGt                                                       <  2:370172‑M1/37‑1 (MQ=255)
                                     ttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcgga                                                      >  2:621395‑M1/1‑61 (MQ=255)
                                     ttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTa                                                      >  2:472816/1‑90 (MQ=255)
                                        aTTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaa                                                   <  2:592976‑M1/90‑33 (MQ=255)
                                        aTTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaa                                                   <  2:501005‑M1/90‑33 (MQ=255)
                                            aaCTTATTGGTTGTAACGGTAAAAAAACAAATAAAAACAAGACAGTTGGGGACAAAATTTCCATATCATAAGTAAGCTACGTATTTGTtc                                               >  1:51391/1‑90 (MQ=255)
                                                       aataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttat                                    >  2:354272‑M1/1‑43 (MQ=255)
                                                           aaGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAata                                >  1:506992/1‑90 (MQ=255)
                                                            aGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgtt                               <  1:621395‑M1/90‑53 (MQ=255)
                                                            aGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgtt                               <  2:412938‑M1/90‑53 (MQ=255)
                                                             cacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATAtt                              >  1:159283‑M1/29‑90 (MQ=255)
                                                                 aaaaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                          <  2:386709‑M1/90‑58 (MQ=255)
                                                                   aaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcagacaaactttttcagtttatttgttggcttaa                        <  1:52043‑M1/90‑60 (MQ=255)
                                                                       cAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGAc                    >  1:493144/1‑90 (MQ=255)
                                                                        aaataaatACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattc                   <  1:483998‑M1/90‑65 (MQ=255)
                                                                        aaataaatACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattc                   <  1:223273‑M1/90‑65 (MQ=255)
                                                                          ataaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaa                 <  2:493144/90‑1 (MQ=255)
                                                                           aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                >  2:311054‑M1/15‑90 (MQ=255)
                                                                           aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                >  2:183388‑M1/15‑90 (MQ=255)
                                                                             aagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAg              >  1:370172‑M1/13‑90 (MQ=255)
                                                                                 cAAGACAATTGGGGCCAttttgactcataggggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttat          >  2:285746‑M1/1‑17 (MQ=255)
                                                                                 cAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATc          <  1:540089/90‑1 (MQ=255)
                                                                                     acatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAAc      >  2:298724‑M1/5‑90 (MQ=255)
                                                                                        aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGt   >  1:189191/1‑90 (MQ=255)
                                                                                         TTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGtt  <  1:548549‑M1/90‑1 (MQ=255)
                                                                                         |                                                                                         
ATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTT  >  NZ_CP009273/4532410‑4532588

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