New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4542682 =33 (0.750)7 (0.180) 6/116 1.8 17.2% intergenic (+49/‑433) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
?NC_000913 4542996 = 39 (1.020)intergenic (+363/‑119) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
Rejected: Frequency below/above cutoff threshold.

AAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAA  >  NC_000913/4542615‑4542748
                                                                   |                                                                  
tacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatt                                                                    >  1:3081814‑M1/68‑68 (MQ=255)
 acactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                   <  1:835695‑M1/2‑1 (MQ=255)
 acactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacattt                                                                   <  1:2121911‑M1/2‑1 (MQ=255)
  cactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTg                                                                  <  1:2705502‑M1/3‑1 (MQ=255)
   actttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTtgt                                                                 <  1:3459135‑M1/4‑2 (MQ=255)
   actttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgg                                                                 <  1:3097471‑M1/4‑1 (MQ=255)
   actttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgg                                                                 <  1:3409507‑M1/4‑1 (MQ=255)
   actttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgg                                                                 <  1:1573379‑M1/4‑1 (MQ=255)
   actttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgg                                                                 <  1:2324745‑M1/4‑1 (MQ=255)
   actttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgg                                                                 <  1:229834‑M1/4‑1 (MQ=255)
    agagaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTggg                                                                >  1:679660/1‑68 (MQ=255)
      ttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGc                                                              <  1:3063158‑M1/7‑1 (MQ=255)
      agaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGc                                                              <  1:1224824/68‑1 (MQ=255)
      agaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGc                                                              <  1:1541766/68‑1 (MQ=255)
         ttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCaa                                                           >  1:3107973‑M1/59‑68 (MQ=255)
         agagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATGCAAGACAATTGGGGCCaa                                                           <  1:108726/68‑1 (MQ=255)
         agagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                           >  1:510182/1‑68 (MQ=255)
         agagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaa                                                           <  1:2552722/68‑1 (MQ=255)
           agGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAc                                                         <  1:553505/68‑1 (MQ=255)
            tgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACt                                                        >  1:2902726‑M1/56‑68 (MQ=255)
            ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAAcc                                                        >  1:2993752/1‑67 (MQ=255)
              tttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                      >  1:6710/1‑68 (MQ=255)
              tttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                      <  1:1283287/68‑1 (MQ=255)
              tttGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGt                                                      <  1:921645/68‑1 (MQ=255)
              taaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGt                                                      >  1:1416269‑M1/54‑68 (MQ=255)
              taaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGt                                                      >  1:1319573‑M1/54‑68 (MQ=255)
                   tttAACTTATTGATAATAAAGTTAAAAACACAAATAAATACAAGACAATTGGGGCCAAACTGTCCata                                                 >  1:1725185/1‑68 (MQ=255)
                   tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCata                                                 >  1:659814/1‑68 (MQ=255)
                   tttAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCata                                                 >  1:1516787/1‑68 (MQ=255)
                             tGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaataa                                       >  1:731829/1‑68 (MQ=255)
                               ataataAAGTTAGAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGt                                     <  1:2986374/68‑1 (MQ=255)
                               ataataAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGt                                     >  1:396994/1‑68 (MQ=255)
                                  ataaAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTAc                                  >  1:2417305/1‑68 (MQ=255)
                                   actgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACg                                 <  1:225700‑M1/36‑1 (MQ=255)
                                      aGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAt                              >  1:1387273/1‑68 (MQ=255)
                                        acagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAttt                            <  1:1732295‑M1/41‑1 (MQ=255)
                                             aaaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtct                       >  1:2688468/1‑68 (MQ=255)
                                                   aataaataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCa                 <  1:3453185/68‑1 (MQ=255)
                                                     aaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATa               >  1:195431‑M1/15‑68 (MQ=255)
                                                        ataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaaa            >  1:330594/1‑68 (MQ=255)
                                                        ataCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATaaaa            >  1:3481694/1‑68 (MQ=255)
                                                               aCAATTGGGGCCAAACTGTCCATATCATAAATACGTTACGTATTTTTTCTCAAGCATAAAAATATTaa     >  1:690512/1‑68 (MQ=255)
                                                               aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa     >  1:3332551/1‑68 (MQ=255)
                                                               aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa     >  1:1210513/1‑68 (MQ=255)
                                                               aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa     >  1:895070/1‑68 (MQ=255)
                                                               aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa     >  1:865259/1‑68 (MQ=255)
                                                               aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa     >  1:732208/1‑68 (MQ=255)
                                                               aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa     >  1:1559801/1‑68 (MQ=255)
                                                               aCAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa     >  1:690433/1‑68 (MQ=255)
                                                               aCAATTGGGGCCAAACTGTCCACATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTaa     >  1:1877982/1‑68 (MQ=255)
                                                                  aTTGGGGCCAAACTGTCCATATCATAAATAAGCTACGTATTTTTTCTCAAGCATAAAAATATTaaaaa  >  1:2797094/1‑68 (MQ=255)
                                                                   |                                                                  
AAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAA  >  NC_000913/4542615‑4542748

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