New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 454269022 (0.620)5 (0.150) 5/246 2.6 13.6% intergenic (+57/‑425) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
?NC_000913 = 4542986 43 (1.310)intergenic (+353/‑129) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
Rejected: Frequency below/above cutoff threshold.

ATGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTAT  >  NC_000913/4542560‑4542822
                                                                                                                                     |                                                                                                                                    
aTGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaaa                                                                                                                                    <  2:585952/133‑1 (MQ=255)
   ctgctCGTTTTGCCGGATTATGGCTGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATAGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAtt                                                                                                                                    <  2:350549‑M1/133‑3 (MQ=255)
     gctCGTTTTGCCGGACTAAG‑‑‑GGGAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTc                                                                                                                               <  1:281831/133‑1 (MQ=255)
      ctCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTc                                                                                                                               <  1:556225/132‑1 (MQ=255)
      ctCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTc                                                                                                                               >  2:556225/1‑132 (MQ=255)
        cGTTTTGCTGGATTACG‑‑‑GCAAACAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactca                                                                                                                            <  1:286469‑M1/133‑11 (MQ=255)
           tttGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                                                                                                         <  2:23205/133‑1 (MQ=255)
            ttcctttcaaaaaactatttctaaatcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCata                                                                                                                           >  2:456987‑M1/114‑133 (MQ=255)
                 ggATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaat                                                                                                                   <  2:118772/133‑1 (MQ=255)
                       tG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTa                                                                                                             >  2:521016/1‑133 (MQ=255)
                                       gacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtctc                                                                                                >  1:508101‑M1/87‑133 (MQ=255)
                                       atCTCATAAATGAAATATTAAAAACAGACGAGTTTTGATTTCACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtctc                                                                                                <  2:164552/133‑1 (MQ=255)
                                            aTAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCa                                                                                           >  1:603089/1‑133 (MQ=255)
                                            aTAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCa                                                                                           >  1:5093/1‑133 (MQ=255)
                                                              ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCTATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACg                                                                         <  2:591977‑M1/70‑1 (MQ=255)
                                                                agaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAATGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGAc                                                                       >  1:441344/1‑133 (MQ=255)
                                                                  aagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattct                                                                     >  2:16918‑M1/1‑68 (MQ=255)
                                                                     agGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaaa                                                                  <  2:603089/133‑1 (MQ=255)
                                                                          tGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATc                                                             >  1:394047/1‑133 (MQ=255)
                                                                           gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCt                                                            <  2:5093/133‑1 (MQ=255)
                                                                            aTTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatct                                                           <  1:16918‑M1/133‑76 (MQ=255)
                                                                             tttAACTTACTGATATTAAAGTTAAACAGACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTaa                                                          <  2:278189/133‑1 (MQ=255)
                                                                                 aCTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGTCTCTAACttt                                                      >  2:304899/1‑131 (MQ=255)
                                                                                  cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGtt                                                     >  2:536996/1‑133 (MQ=255)
                                                                                     aTTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttata                                                  <  1:597675‑M1/133‑85 (MQ=255)
                                                                                      ttGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGat                                                 >  1:127898/1‑133 (MQ=255)
                                                                                                   tAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTAttt                                    <  1:521016/133‑1 (MQ=255)
                                                                                                        ttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTAtt                               <  2:531074‑M1/112‑1 (MQ=255)
                                                                                                                   taCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTGTATTGTAaattaatt                    >  2:247576/1‑133 (MQ=255)
                                                                                                                        gACAATTGGGGCCAACCTGTTCATATCAGAAAAAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTccca               <  2:253470/133‑1 (MQ=255)
                                                                                                                        gACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTcaca               >  1:62277/1‑133 (MQ=255)
                                                                                                                             ttGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatgtttatattgcatgaggtggtttttggagag                <  2:603417‑M1/127‑119 (MQ=255)
                                                                                                                             ttGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatgtttatattgcatgaggtggtttttggagag                >  1:603417‑M1/1‑9 (MQ=255)
                                                                                                                             TTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCAcc          <  1:419489‑M1/133‑1 (MQ=255)
                                                                                                                             TTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAGTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATCATTTCTATTGTAAACTAATTTGACATCAcc          <  1:253470‑M1/133‑1 (MQ=255)
                                                                                                                             TTGGGGCCAAACTCTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCAcc          <  2:382185‑M1/133‑1 (MQ=255)
                                                                                                                              TGGGGCCAAACTGTCTATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCt         >  1:591977‑M1/1‑133 (MQ=255)
                                                                                                                              TGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCt         >  1:531074‑M1/1‑133 (MQ=255)
                                                                                                                                GGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGACATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTcc       <  1:307468‑M1/133‑1 (MQ=255)
                                                                                                                                  gCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatgtttatattgcatgaggtggtttttggagagaagaatgagga     >  2:141730‑M1/1‑4 (MQ=255)
                                                                                                                                  GCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGc     <  1:349561‑M1/133‑1 (MQ=255)
                                                                                                                                   CCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCt    <  1:504615‑M1/133‑1 (MQ=255)
                                                                                                                                     AAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTAGTGTAAATAAATGTCGCATCACCTCCGCtat  <  1:169036‑M1/133‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                    
ATGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTAT  >  NC_000913/4542560‑4542822

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