New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 45325078 (0.320)3 (0.130) 3/150 2.2 18.7% intergenic (+57/‑425) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
?NZ_CP009273 = 4532803 19 (0.850)intergenic (+353/‑129) fimE/fimA type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
Rejected: Frequency below/above cutoff threshold.

AAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATG  >  NZ_CP009273/4532420‑4532595
                                                                                       |                                                                                        
aaaCGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaaa                                                                                        <  2:455153/90‑1 (MQ=255)
      aaaaTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgact                                                                                  <  2:38621‑M1/90‑9 (MQ=255)
           tAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                                                             <  2:269075/90‑1 (MQ=255)
                  agaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGGCAATTGGGGCCAAACTGTCCATATCataaata                                                                      <  2:531977/90‑1 (MQ=255)
                  agaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaata                                                                      <  1:112446/90‑1 (MQ=255)
                       agGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcg                                                                 <  2:663203‑M1/90‑26 (MQ=255)
                                                  aGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgtt                                      <  2:553251‑M1/90‑53 (MQ=255)
                                                                         aGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTa               >  1:12165/1‑90 (MQ=255)
                                                                                  gggCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGa      <  1:12916/90‑1 (MQ=255)
                                                                                      cAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATg  >  1:551747/1‑90 (MQ=255)
                                                                                      cAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATg  <  2:551747/90‑1 (MQ=255)
                                                                                       attttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctt              >  1:582893‑M1/1‑1 (MQ=255)
                                                                                       attttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctt              <  2:582893‑M1/77‑77 (MQ=255)
                                                                                       |                                                                                        
AAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATG  >  NZ_CP009273/4532420‑4532595

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