New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 454269042 (0.760)10 (0.210) 7/110 1.9 18.3% intergenic (+57/‑425) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
?NC_000913 = 4542986 53 (1.120)intergenic (+353/‑129) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
Rejected: Frequency below/above cutoff threshold.

AGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGA  >  NC_000913/4542626‑4542752
                                                                |                                                              
agGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCA                                                                >  1:2066455‑M1/1‑65 (MQ=255)
 ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                               >  1:651865/1‑65 (MQ=255)
 ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                               >  1:3598676/1‑65 (MQ=255)
 ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                               <  1:2022620/65‑1 (MQ=255)
 ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                               <  1:1976878/65‑1 (MQ=255)
 ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                               >  1:1765737/1‑65 (MQ=255)
 ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAt                                                               <  1:609895‑M1/65‑2 (MQ=255)
   tttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                             >  1:3461997/1‑65 (MQ=255)
   tttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                             >  1:190849/1‑65 (MQ=255)
   tttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATAGGGGCCAAAc                                                             <  1:4276524/65‑1 (MQ=255)
   tttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAAAAAATACAAGACAATTGGGGCCAAAc                                                             <  1:82632/65‑1 (MQ=255)
    ttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACt                                                            <  1:2256732/65‑1 (MQ=255)
      gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                          <  1:1594416/65‑1 (MQ=255)
      gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                          >  1:1730254/1‑65 (MQ=255)
      gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                          <  1:2634357/65‑1 (MQ=255)
      gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                          <  1:1365577/65‑1 (MQ=255)
      gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                          <  1:3737646/65‑1 (MQ=255)
      acgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGt                                                          >  1:806915‑M1/51‑65 (MQ=255)
        tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgact                                                        >  1:1762820‑M1/1‑57 (MQ=255)
        tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgact                                                        >  1:738927‑M1/1‑57 (MQ=255)
        tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgact                                                        >  1:485840‑M1/1‑57 (MQ=255)
        tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTcc                                                        >  1:2771383/1‑65 (MQ=255)
        tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTcc                                                        >  1:3487280/1‑65 (MQ=255)
        tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTcc                                                        >  1:2081061/1‑65 (MQ=255)
        tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTcc                                                        >  1:1576406/1‑65 (MQ=255)
            aCTTATTGATAATAAAGTTAAAAAAGCAAATAAGTACAAGACAATTGGGGCCAAACTGTCCatat                                                    <  1:4218365/65‑1 (MQ=255)
             cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatag                                                   >  1:598295‑M1/1‑52 (MQ=255)
             cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                                   <  1:903781/65‑1 (MQ=255)
             cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGCACAAACTGTCCATATc                                                   >  1:70140/1‑65 (MQ=255)
              ttATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCa                                                  >  1:928751/1‑65 (MQ=255)
               tATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCat                                                 <  1:4093712/65‑1 (MQ=255)
                  tGATAATAAAGTTAAAAAAACAAATAAATACAAGACAACTCGCGCCACACTGCCCATATCataaa                                              >  1:2725938/1‑65 (MQ=255)
                    ataataAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaata                                            >  1:3046867/1‑65 (MQ=255)
                     aacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCataaataa                                           >  1:3964235‑M1/36‑65 (MQ=255)
                       ataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcat                                         >  1:3015839‑M1/1‑42 (MQ=255)
                       ataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGt                                         <  1:1225875/65‑1 (MQ=255)
                         ctgcacagttttaccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTa                                       <  1:2687233‑M1/34‑1 (MQ=255)
                         aaaGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTa                                       <  1:1889924/65‑1 (MQ=255)
                            gTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGt                                    >  1:97170/1‑65 (MQ=255)
                             ttAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTa                                   <  1:3677577/65‑1 (MQ=255)
                             ttAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTa                                   <  1:481972/65‑1 (MQ=255)
                             acagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTa                                   >  1:2206769‑M1/28‑65 (MQ=255)
                                  tttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAttttt                              <  1:726194‑M1/43‑1 (MQ=255)
                                          taaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGc                      <  1:1884989/65‑1 (MQ=255)
                                          taaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGc                      >  1:3169116/1‑65 (MQ=255)
                                           aaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCa                     <  1:3017368/65‑1 (MQ=255)
                                               aCAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttat                 >  1:725162‑M1/1‑18 (MQ=255)
                                               aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaa                 >  1:1338566/1‑65 (MQ=255)
                                               aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaa                 <  1:1449187/65‑1 (MQ=255)
                                               aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATCTTTTCTCAAGCATaaa                 >  1:2243365/1‑65 (MQ=255)
                                                cAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaaa                <  1:273176/65‑1 (MQ=255)
                                                    aCAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgtt            >  1:3306156‑M1/1‑13 (MQ=255)
                                                    aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat            >  1:1212046/1‑65 (MQ=255)
                                                    aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat            >  1:1938563/1‑65 (MQ=255)
                                                    aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat            <  1:4209879/65‑1 (MQ=255)
                                                    aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat            >  1:2468757/1‑65 (MQ=255)
                                                       ttgggggccattttgactcatagaggaaagcatcgcgggcaaactttttcagtttatttgttggc         <  1:4215688‑M1/62‑56 (MQ=255)
                                                       aTTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggc         <  1:3162879‑M1/65‑56 (MQ=255)
                                                       aTTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggc         <  1:2797648‑M1/65‑56 (MQ=255)
                                                         tGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt       >  1:61819‑M1/1‑8 (MQ=255)
                                                         tGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctt       <  1:747715‑M1/65‑58 (MQ=255)
                                                         TGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaaaa       <  1:3459220‑M1/65‑1 (MQ=255)
                                                          ggggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggctta      >  1:700462‑M1/1‑7 (MQ=255)
                                                            ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaat    <  1:4242026‑M1/65‑61 (MQ=255)
                                                            ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaat    <  1:1639028‑M1/65‑61 (MQ=255)
                                                            ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaat    <  1:1221630‑M1/65‑61 (MQ=255)
                                                            ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaat    <  1:1115975‑M1/65‑61 (MQ=255)
                                                            ggCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaat    <  1:662542‑M1/65‑61 (MQ=255)
                                                            ggCCAttttgacacatagaggaaagcatcgcggacaaactttttcagttgatttgttggcttaat    <  1:122481‑M1/65‑61 (MQ=255)
                                                              ccAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatat  >  1:1149638‑M1/1‑3 (MQ=255)
                                                              ccAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatat  <  1:2069234‑M1/65‑63 (MQ=255)
                                                              ccAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatat  <  1:2424467‑M1/65‑63 (MQ=255)
                                                              ccAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatat  >  1:2740541‑M1/1‑3 (MQ=255)
                                                              ccAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatat  <  1:3906431‑M1/65‑63 (MQ=255)
                                                              ccAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatat  <  1:3865714‑M1/65‑63 (MQ=255)
                                                              ccAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatat  <  1:3047849‑M1/65‑63 (MQ=255)
                                                              ccAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatat  <  1:3288880‑M1/65‑63 (MQ=255)
                                                              CCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGa  <  1:2160150‑M1/65‑1 (MQ=255)
                                                              CCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGa  <  1:492006‑M1/65‑1 (MQ=255)
                                                              CCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGa  <  1:2781512‑M1/65‑1 (MQ=255)
                                                              CCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGa  <  1:2827041‑M1/65‑1 (MQ=255)
                                                                |                                                              
AGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGA  >  NC_000913/4542626‑4542752

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