New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4542682 =	21 (0.200)	3 (0.030)	3/452	9.6	8.5%	intergenic (+49/‑433)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913	4542996 =	44 (0.450)					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.
Rejected: Frequency below cutoff threshold.

AGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATC  >  NC_000913/4542451‑4542902                                                                                                                                                                                                                                        |                                                                                                                                                                                                                             aGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGACATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATGCAAGACAATTGGGGCCAAACTGTCCata                                                                                                                                                                                                           >  1:580587/1‑251 (MQ=255)      tGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATCCAGGATTATCTCGAGCGTCGAAACATTCGCCATACTGTGCGACATACCGCTAGTAATGCTGCACGTTTTGCTGTATTATGGGAAAGAAATAATCTCATAAACGAACAATTAAAAAGAGAAGAGGTTCGATTTAACTTATTGATATTAAAGTTAAAAAAACAAATAAATAAAAGCCACTGGGGGCAAAACTCGacataacata                                                                                                                                                                                                      >  2:131029/1‑251 (MQ=255)               ttATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGTGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCAGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTa                                                                                                                                                                                             <  2:422152/251‑1 (MQ=255)                   gAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTa                                                                                                                                                                                         >  1:994614/1‑251 (MQ=255)                                           aCTCGTTTAATTTAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATAtt                                                                                                                                                                 <  1:204664/251‑1 (MQ=255)                                                      tctagATTCTCGCGGGCATCGAAATATTCGCCATACTGTGCGTTAAACCGCCAGTAATGCTGCTCGTTTTGTCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACa                                                                                                                                                      <  2:373005/247‑1 (MQ=255)                                                           aTTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTGTCTCAAGCATAAACATATTAAAAAACGACAAAAAg                                                                                                                                                 >  1:372929/1‑251 (MQ=255)                                                              aTCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATGCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCAt                                                                                                                                              <  2:580703/251‑1 (MQ=255)                                                                            aaTATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCAGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGata                                                                                                                                >  1:46515/1‑251 (MQ=255)                                                                                      aTACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTATCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAattat                                                                                                                      >  2:1062274/1‑251 (MQ=255)                                                                                                ttATACCGCCAGTAATGCTGCTCGATTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGCTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTACATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAATTGTTTGATATATAAATTATTTCTATTGTa                                                                                                            >  1:949951/1‑251 (MQ=255)                                                                                                 tataCCTCCAGGAATTCTGCTATGTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTaa                                                                                                           <  1:871315/251‑1 (MQ=255)                                                                                                     ccGCAAGTAATGCTGCTCGTTTTTCCAGATGATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAAGTTTTATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAAACAATTGGGGCCAAACTGTCCATATCAGAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAaatta                                                                                                       <  1:729199/251‑1 (MQ=255)                                                                                                      cGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAaattaa                                                                                                      <  1:1062112/251‑1 (MQ=255)                                                                                                           gTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAACAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATAGCATGAATAAGTTACGTATTTTTTCTCAAGCATAAAAATACTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCAATTGTAAATTAATTTca                                                                                                 >  2:488359/1‑251 (MQ=255)                                                                                                                                                                  aTTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTc                                          <  2:994776/251‑1 (MQ=255)                                                                                                                                                                                atgTTTGAATTAACTTATTTATAATAAAGTAAAAAAAACAAATAAATACATGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTATCAAGCATAAAAATATTAAAAAACAACTAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCaa                            <  2:46556/249‑1 (MQ=255)                                                                                                                                                                                agGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGGCCATATCATGAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAGAAGCATCAAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCaa                            <  1:488249/251‑1 (MQ=255)                                                                                                                                                                                                          aaattAAAAAAACAAATAAATACAGGACAATATGGACCAAACTGTCCATAACATAAATAAGTTACGTATTTTTTCTCAAGCATACAAGTATTAAAAAACGACAAAAAGCATCTAACTGTTTTATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCACCGCTATATGTAAAGCTAACATTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATc  <  1:130986/248‑1 (MQ=255)                                                                                                                                                                                                                      cAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCt                                           <  1:432952/198‑1 (MQ=255)                                                                                                                                                                                                                                  aGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCaaaa                          <  1:4167/203‑1 (MQ=255)                                                                                                                                                                                                                                        ttGGGGCCAAACTGTCCATATCATCAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTg                                                                                                                                       >  2:1061655/1‑88 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                             AGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATC  >  NC_000913/4542451‑4542902

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