New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4542682 =	17 (0.380)	8 (0.190)	8/264	3.9	28.6%	intergenic (+49/‑433)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913	4542996 =	24 (0.570)					intergenic (+363/‑119)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
Rejected: Coverage evenness skew score above cutoff.

TACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTA  >  NC_000913/4542549‑4542827                                                                                                                                      |                                                                                                                                                  taCCGCCAGTAATGCTGCTAGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTc                                                                                                                                    <  1:154403/149‑1 (MQ=255)   tgcagaGTAATGCTGCTCGTTTTGCCTGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTTATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCa                                                                                                                                  <  2:119442/144‑1 (MQ=255)     gcaaGTAATGCCGGTCGTTTTGCGGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCata                                                                                                                                <  1:605530/146‑1 (MQ=255)      ccAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCatat                                                                                                                               >  1:612135/1‑149 (MQ=255)               ctgctCGTATAGTCGGAATAAGGGAAAAAAATAATCTCATAAACGAAAAAATAAAAAAAGAAGAGGTTTGATAAAACTGATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaataa                                                                                                                      <  2:12570/149‑1 (MQ=255)                    cGTTTTGACGGATTATGGGAAAGAAATAATCTCATAAACGAAAAAGTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTAc                                                                                                                 <  2:159688/149‑1 (MQ=255)                                      gaaagaaaTAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCCCAAGCAt                                                                                               >  2:228097/1‑149 (MQ=255)                                                                    aaaaGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCt                                                                 >  2:287962/1‑149 (MQ=255)                                                                           aggcgggttgattgaGCGTATTGGTAAAAAAGGTAAAAAAACAAAGCGATACAAGACGATTGGGGCCGAACTGTGCATATCTTAAATAAGGTACGTATTTTTTGTCAAGCATAAAAGTATTAAAAAACGACAAAAAGCGTCTAACTGtt                                                          <  2:578341/142‑1 (MQ=255)                                                                                gTTTAATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCCCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGata                                                     <  2:162496/149‑1 (MQ=255)                                                                                            ttATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATGAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTAttt                                         <  1:511658/149‑1 (MQ=255)                                                                                                tGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTAt                                     <  1:193888/149‑1 (MQ=255)                                                                                                  ataataAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTAttt                                   >  1:258547/1‑148 (MQ=255)                                                                                                             aaaaaaaCAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAACATAATTGCAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAAttt                        >  1:764797/1‑149 (MQ=255)                                                                                                                     aaataaatACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGAATAAAAATATTAAAAAACGACAAAAAGCAACTAACTGTTTGATATGTAAATTATTTCTATTGGAAATTAATTTcaccgaac                >  2:87610/1‑144 (MQ=255)                                                                                                                               aaGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCtata      >  2:249222/1‑149 (MQ=255)                                                                                                                                  aCATTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTATATTGTAAATTAATTTCACATCACCTCCGCTatattt   >  2:433640/1‑147 (MQ=255)                                                                                                                                   cAAGTGGGGCAAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTa  <  2:585861/149‑1 (MQ=255)                                                                                                                                      |                                                                                                                                                  TACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTA  >  NC_000913/4542549‑4542827

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