New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 45325074 (0.150)18 (0.740) 14/150 0.4 64.8% intergenic (+57/‑425) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 = 4532803 16 (0.650)intergenic (+353/‑129) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

TAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATAT  >  NZ_CP009273/4532419‑4532594
                                                                                        |                                                                                       
ggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCaa                                                                                        >  1:416423‑M1/81‑90 (MQ=255)
ggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCaa                                                                                        >  2:238767‑M1/81‑90 (MQ=255)
  cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                      >  1:175364‑M1/79‑90 (MQ=255)
  cagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAAc                                                                                      >  1:206485‑M1/79‑90 (MQ=255)
    cGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttg                                                                                    <  1:676817‑M1/90‑6 (MQ=255)
    cGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttg                                                                                    <  1:26136‑M1/90‑6 (MQ=255)
     gAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttga                                                                                   <  1:213111‑M1/90‑7 (MQ=255)
     gAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttga                                                                                   <  1:442032‑M1/90‑7 (MQ=255)
       gttctgtacactttgttttgtaaacgagtttgactgccaaccctgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTcc                                                                                 >  2:232868‑M1/74‑90 (MQ=255)
         tctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCat                                                                               >  2:463403‑M1/72‑90 (MQ=255)
         tctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCat                                                                               >  1:172930‑M1/72‑90 (MQ=255)
          aTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcat                                                                              <  2:668435‑M1/90‑12 (MQ=255)
            gtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATc                                                                            >  1:564706‑M1/69‑90 (MQ=255)
                 ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaa                                                                       <  2:15058‑M1/27‑1 (MQ=255)
                 agagaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaa                                                                       <  2:348985‑M1/90‑19 (MQ=255)
                   ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaata                                                                     >  2:525565‑M1/62‑90 (MQ=255)
                    tgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaataa                                                                    <  2:771498‑M1/30‑1 (MQ=255)
                                 agtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAttttt                                                       >  1:771498‑M1/48‑90 (MQ=255)
                                  gtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAtttttt                                                      <  2:327663‑M1/44‑1 (MQ=255)
                                      gactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCa                                                  >  1:351487‑M1/43‑90 (MQ=255)
                                       actgccaacactgcacagtttccccccaaaagatgaaacatTTGGGGCCAAACTCTCCATATCATAAATAAGTTACGTATTTTTTCTCaa                                                 >  1:67937‑M1/42‑90 (MQ=255)
                                        aTTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttca                                                >  1:759024‑M1/1‑49 (MQ=255)
                                                   aGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgtt                                     <  1:388312‑M1/90‑53 (MQ=255)
                                                    cacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATAtt                                    <  2:4784‑M1/62‑1 (MQ=255)
                                                     ttAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttgg                                   <  1:754850‑M1/90‑55 (MQ=255)
                                                        aaaaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                                <  1:414014‑M1/90‑58 (MQ=255)
                                                        aaaaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt                                <  2:759024‑M1/90‑58 (MQ=255)
                                                          aaaaCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatctcggacaaactttttcagtttatttgttggcttaa                              <  1:184525‑M1/90‑60 (MQ=255)
                                                           aaaCAAATAAATACAAGCCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAAc                             <  1:327689/90‑1 (MQ=255)
                                                             aCAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatat                           <  2:547962‑M1/90‑63 (MQ=255)
                                                               aaataaatACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattc                         <  1:21489‑M1/90‑65 (MQ=255)
                                                                  taaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                      <  1:299006/90‑1 (MQ=255)
                                                                  taaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                      <  1:27349/90‑1 (MQ=255)
                                                                  aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaa                      <  2:195960‑M1/76‑1 (MQ=255)
                                                                         gaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCt               <  2:4033‑M1/83‑1 (MQ=255)
                                                                          aaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa              >  1:527262‑M1/7‑90 (MQ=255)
                                                                          aGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatct              >  2:220840‑M1/1‑15 (MQ=255)
                                                                          aGACAATTGGGGCCAAACTGTCCACATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa              >  1:35268/1‑90 (MQ=255)
                                                                             cAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatcttta           <  1:39768‑M1/90‑79 (MQ=255)
                                                                               actggGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaa                              <  2:432117‑M1/67‑60 (MQ=255)
                                                                               actggGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaa                              >  1:432117‑M1/3‑10 (MQ=255)
                                                                               aTTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatt         >  2:556395‑M1/1‑10 (MQ=255)
                                                                               aTTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatt         >  2:700993‑M1/1‑10 (MQ=255)
                                                                                TTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGtt        >  1:172055‑M1/1‑90 (MQ=255)
                                                                                  ggggCCAttttgactcatagaggaaagaatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttat      <  1:724871‑M1/90‑84 (MQ=255)
                                                                                   GGGCCAAACTGTCCATATCATAAATAAGTTACATATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGa     >  2:492613‑M1/1‑90 (MQ=255)
                                                                                    ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatag    >  1:183261‑M1/1‑5 (MQ=255)
                                                                                    ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatag    <  1:173338‑M1/90‑86 (MQ=255)
                                                                                    ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattc                         <  1:236507‑M1/69‑65 (MQ=255)
                                                                                    ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattc                         >  2:236507‑M1/1‑5 (MQ=255)
                                                                                      CCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGatat  >  1:171363‑M1/1‑90 (MQ=255)
                                                                                        |                                                                                       
TAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATAT  >  NZ_CP009273/4532419‑4532594

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