New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4542682 =62 (1.400)4 (0.100) 4/110 2.1 6.9% intergenic (+49/‑433) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
?NC_000913 4542996 = 54 (1.420)intergenic (+363/‑119) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
Rejected: Frequency below/above cutoff threshold.

AAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAA  >  NC_000913/4542618‑4542745
                                                                |                                                               
actttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatt                                                                 <  1:2716468‑M1/1‑1 (MQ=255)
actttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatt                                                                 <  1:2988647‑M1/1‑1 (MQ=255)
aagagaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAt                                                                 <  1:1954590/65‑1 (MQ=255)
 ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                <  1:977074‑M1/2‑1 (MQ=255)
 ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                <  1:1578615‑M1/2‑1 (MQ=255)
 ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                <  1:876704‑M1/2‑1 (MQ=255)
 ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                <  1:3192207‑M1/2‑1 (MQ=255)
 agagaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAtt                                                                >  1:1799972/1‑65 (MQ=255)
   ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgg                                                              <  1:1757183‑M1/4‑1 (MQ=255)
   ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgg                                                              <  1:582861‑M1/4‑1 (MQ=255)
    tgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTggg                                                             >  1:1362238‑M1/61‑65 (MQ=255)
    gaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAGTTggg                                                             <  1:2330299/65‑1 (MQ=255)
    gaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTggg                                                             <  1:26985/65‑1 (MQ=255)
    gaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTggg                                                             <  1:991413/65‑1 (MQ=255)
       tttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGcc                                                          <  1:2345350‑M1/8‑1 (MQ=255)
        agGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCa                                                         >  1:513621/1‑65 (MQ=255)
         ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCat                                                        <  1:1747590/65‑2 (MQ=255)
         ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                        >  1:318098/1‑65 (MQ=255)
         ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                        <  1:2880776/65‑1 (MQ=255)
         ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                        <  1:8986/65‑1 (MQ=255)
         ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                        <  1:565988/65‑1 (MQ=255)
           tttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                      >  1:1475132/1‑65 (MQ=255)
            ttGATTTAACTTATTGATAATAAAGTTAAAGAAACAAATAAATACAAGACAATTGGGGCCAAACt                                                     <  1:2639063/65‑1 (MQ=255)
            ttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACt                                                     <  1:3401158/65‑1 (MQ=255)
            ttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACt                                                     <  1:3403012/65‑1 (MQ=255)
            ttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACt                                                     <  1:3084535/65‑1 (MQ=255)
            ttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACt                                                     <  1:361859/65‑1 (MQ=255)
            aaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACt                                                     >  1:3436308‑M1/53‑65 (MQ=255)
              gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                   <  1:1530871/65‑1 (MQ=255)
              gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                   <  1:1391848/65‑1 (MQ=255)
              acgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGt                                                   >  1:2947063‑M1/51‑65 (MQ=255)
                tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTcc                                                 >  1:2704555/1‑65 (MQ=255)
                tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTcc                                                 >  1:1033550/1‑65 (MQ=255)
                tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTcc                                                 >  1:2476966/1‑65 (MQ=255)
                tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTcc                                                 >  1:239951/1‑65 (MQ=255)
                tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTcc                                                 >  1:2162304/1‑65 (MQ=255)
                   aaCTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCata                                              <  1:691289/65‑1 (MQ=255)
                   aaCTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCata                                              <  1:1415506/65‑1 (MQ=255)
                     cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                            >  1:1015247/1‑65 (MQ=255)
                     cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                            >  1:2517551/1‑65 (MQ=255)
                     cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                            >  1:957588/1‑65 (MQ=255)
                     cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                            >  1:1095841/1‑65 (MQ=255)
                     cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                            >  1:748783/1‑65 (MQ=255)
                     cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                            <  1:564701/65‑1 (MQ=255)
                     cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                            <  1:3453745/65‑1 (MQ=255)
                     cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                            >  1:1515951/1‑65 (MQ=255)
                        aTTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCata                                         <  1:105630/65‑1 (MQ=255)
                        aTTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAGTGGGGCCAAACTGTCCATATCata                                         <  1:951559/65‑1 (MQ=255)
                          tGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaa                                       >  1:244518/1‑65 (MQ=255)
                          tGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaa                                       >  1:397284/1‑65 (MQ=255)
                          tGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaa                                       >  1:2546785/1‑65 (MQ=255)
                            ataataAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaata                                     >  1:2258116/1‑65 (MQ=255)
                            ataataAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaata                                     >  1:1358947/1‑65 (MQ=255)
                               ataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGt                                  >  1:1618300/1‑65 (MQ=255)
                               ataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGt                                  >  1:1159698/1‑65 (MQ=255)
                                actgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGtt                                 <  1:2150530‑M1/33‑1 (MQ=255)
                                  aaGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTAc                               <  1:1451887/65‑1 (MQ=255)
                                  aaGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTAc                               <  1:1875826/65‑1 (MQ=255)
                                   aGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACg                              >  1:332946/1‑65 (MQ=255)
                                   aGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACg                              >  1:1102660/1‑65 (MQ=255)
                                    gTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGt                             >  1:1814660/1‑65 (MQ=255)
                                        gttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAACTAAGTTACGTAttt                         <  1:2597021‑M1/41‑1 (MQ=255)
                                                   aaataCAAGTCAATTGTGGCCAAACTGTCCATATCATAAATAAGTTACGCATTTTTTCTCAAGCa              <  1:2908982/65‑1 (MQ=255)
                                                   aaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCa              <  1:1823937/65‑1 (MQ=255)
                                                     ataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATa            <  1:83090/65‑1 (MQ=255)
                                                       aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaa          >  1:530880/1‑65 (MQ=255)
                                                       aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaa          >  1:403525/1‑65 (MQ=255)
                                                       aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaa          >  1:894177/1‑65 (MQ=255)
                                                       aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaa          >  1:2172047/1‑65 (MQ=255)
                                                       aCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaa          >  1:2458771/1‑65 (MQ=255)
                                                        cAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaaa         >  1:595045/1‑65 (MQ=255)
                                                          aGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAat       <  1:1829603/65‑1 (MQ=255)
                                                            aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat     >  1:1444254/1‑65 (MQ=255)
                                                            aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat     >  1:2848996/1‑65 (MQ=255)
                                                            aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAatat     >  1:2165256/1‑65 (MQ=255)
                                                               aTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa  >  1:313018/1‑65 (MQ=255)
                                                                |                                                               
AAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAA  >  NC_000913/4542618‑4542745

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