New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4542690	33 (0.680)	9 (0.200)	6/260	5.4	20.5%	intergenic (+57/‑425)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913	= 4542986	39 (0.850)					intergenic (+353/‑129)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
Rejected: Coverage evenness skew score above cutoff.

GCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGT  >  NC_000913/4542553‑4542837                                                                                                                                          |                                                                                                                                                    gCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAACCACGTCAACAAATACAAGACAATTGGGGCCAAACTGTCCata                                                                                                                                          <  1:209699/149‑1 (MQ=255)    aGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATAGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTTGGGCCAAACTGTCCATATCa                                                                                                                                       >  1:392759/1‑149 (MQ=255)      tAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCata                                                                                                                                     >  1:50862/1‑149 (MQ=255)           ctgctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaataa                                                                                                                                >  1:857605/1‑149 (MQ=255)            tgctCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAg                                                                                                                               <  2:50862/149‑1 (MQ=255)                       ccGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATCGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTAGGGCCAAACTGTCCATATAATAAATAAGTTACGTAtttt                                                                                                                    >  2:709587/1‑149 (MQ=255)                                  gaaagaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGATTTTCCATTACGTATTAATATTGAAGGTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCAt                                                                                                         >  1:206197/1‑149 (MQ=255)                                                       gAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaa                                                                                    >  2:563288/1‑149 (MQ=255)                                                         aaaaTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtctc                                                                                                               <  2:606945/120‑1 (MQ=255)                                                         aaaaTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtctc                                                                                                               >  1:606945/1‑120 (MQ=255)                                                         aaaaTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTTTTATTTTTTAAAATATTAAAAAATAACaaac                                                                                  >  2:68180/1‑148 (MQ=255)                                                             ttAAGAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCACATCATAGATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCa                                                                              >  2:658518/1‑149 (MQ=255)                                                                 aaaGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATAATAAAAAACGACAAAATTCATCTa                                                                          >  2:387004/1‑149 (MQ=255)                                                                        agagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGttt                                                                   <  1:509384/149‑1 (MQ=255)                                                                         gagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTCTTTTGTAGCAACCAAAAAAATCTCAAAAAGCGACAAAAAGCATCTAACTGTTTg                                                                  >  1:198962/1‑149 (MQ=255)                                                                           ggTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGat                                                                <  1:708408/149‑1 (MQ=255)                                                                            gTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGata                                                               >  1:477331/1‑149 (MQ=255)                                                                                        ttATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTATATTAttt                                                   >  1:503577/1‑149 (MQ=255)                                                                                        ttATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTAttt                                                   <  2:857605/149‑1 (MQ=255)                                                                                             gATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCCTACACAAACTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTAtt                                              <  1:679729/149‑1 (MQ=255)                                                                                             gATAATAAAGTTAAAAAAACAAATAAATACAAGACAAGTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTAtt                                              <  1:709587/149‑1 (MQ=255)                                                                                                     aGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAaatta                                      <  2:477331/149‑1 (MQ=255)                                                                                                       ttAAAAAAACAAATTAATACAAGAAAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAaattaat                                    >  1:747873/1‑149 (MQ=255)                                                                                                           aaaaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTca                                >  2:260353/1‑149 (MQ=255)                                                                                                           aaaaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTATCAAGCATAAAAATATTAAAAAACGACAAAAAGCATATAACTATTTGATATGTATATTATTTCTATTGTAATTTAATTTct                                >  2:520896/1‑148 (MQ=255)                                                                                                              aaCAAATAAATACAAGACAATTGGGGCCAAACTGTACATACCATAAAAAAGATACGTATTTTTTCTCACGCATAAAAATATTAAAAAACGACAAAAAGCATCTAAATGTTTGATATGTAAATTATTTATATTGTAAATTAATTTcactt                             >  2:658965/1‑147 (MQ=255)                                                                                                                           aaGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACTACACCTCCGCtata                >  2:348410/1‑149 (MQ=255)                                                                                                                           aaGACAAATAGGGCAAAACTGTCCAGAGCATAAAAAAGTTACGCAATTTTGAACAAGCATAAAAATATTAAAAAACGACAAAAAACATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCtata                <  1:68180/149‑1 (MQ=255)                                                                                                                             gACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATg              >  1:392740/1‑149 (MQ=255)                                                                                                                             gACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATg              <  1:260353/149‑1 (MQ=255)                                                                                                                                  ttGGGGCCAAACTGTCCAGATCATAAATAAGTTACGTATGTTTTCTCTAGCATAAAAATATTAAAAAACGACAAAAAGCTTCTAACTGTTTGATATGTAAATTATGTCTATTGTAAATTAATTGCACATCACCTCCGCTATATGTAAAg         <  1:658518/149‑1 (MQ=255)                                                                                                                                   tGGGGCCAAACTGTCAATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGc        <  2:392759/149‑1 (MQ=255)                                                                                                                                         cAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATCTGTAAACCTAACGt  >  2:77740/1‑149 (MQ=255)                                                                                                                                          |                                                                                                                                                    GCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGT  >  NC_000913/4542553‑4542837

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