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.

CGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTT  >  NC_000913/4542851‑4543144                                                                                                                                                  |                                                                                                                                                      cGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGATTCAAAATGGCCCGAATTGt                                                                                                                                                    >  2:398424/1‑149 (MQ=255)   cATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGttt                                                                                                                                                  <  2:683515/149‑1 (MQ=255)      cTTCCTCATTCTTCTCTCCAAGAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCACGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCAt                                                                                                                                               >  2:688758/1‑149 (MQ=255)       ttCCTCATTGTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATc                                                                                                                                              <  2:258547/149‑1 (MQ=255)          cTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCttt                                                                                                                                           >  1:183818/1‑149 (MQ=255)            cATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCACCAAATGTTTCATCtttgg                                                                                                                                         >  2:251993/1‑147 (MQ=255)              ttcttcTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                                                                                                       <  1:475395/149‑1 (MQ=255)              ttcttcTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                                                                                                       <  2:183818/149‑1 (MQ=255)                      tcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGCGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACt                                                                                                                              >  1:557640/1‑149 (MQ=255)                                accTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGtt                                                                                                                    <  2:764797/149‑1 (MQ=255)                                          atatAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCaaa                                                                                                          <  1:586361/149‑1 (MQ=255)                                           tataAACATCTATAAATAAAGATAACAAAAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAAc                                                                                                         >  2:441629/1‑149 (MQ=255)                                                 cATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGttt                                                                                                   <  1:750473/149‑1 (MQ=255)                                                           taaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaaa                                                                                         <  1:209357/149‑1 (MQ=255)                                                                             atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGcc                                                                       <  1:526800/149‑1 (MQ=255)                                                                                     aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGa                                                               <  1:445976/149‑1 (MQ=255)                                                                                       caaATAAACTGAAATAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGCGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTCCCCATGTCGact                                                             >  1:676184/1‑147 (MQ=255)                                                                                        aaaTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTAACGAACGACTGCCCATGTCGAttt                                                            >  2:597585/1‑149 (MQ=255)                                                                                                         ttGTGCGGGATGCTTTCCTCTATGAGTCAAAATGG‑CCCATCTCAATTACCAGCTGGGGAAGGAGCCAGGCAGTGTTGGCAGTAAATGTCCTCGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaa                                           <  2:557640/149‑1 (MQ=255)                                                                                                          tGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaa                                          <  1:548870/149‑1 (MQ=255)                                                                                                                          ctctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaag                                         >  1:454464/1‑134 (MQ=255)                                                                                                                          ctctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaag                                         <  2:454464/134‑1 (MQ=255)                                                                                                                              aTGAGGCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCGGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTa                      <  1:251993/149‑1 (MQ=255)                                                                                                                                       aaTGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTCTCCTCGTCATAGCAGCATGAAAATTAAAACTCTgg             <  2:561856/149‑1 (MQ=255)                                                                                                                                             cccAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCg       >  2:464383/1‑149 (MQ=255)                                                                                                                                                 aaTGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCgttgt   >  2:482073/1‑149 (MQ=255)                                                                                                                                                 aaTGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCgttgt   <  1:688758/149‑1 (MQ=255)                                                                                                                                                  aTGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCgttgtt  <  2:670330/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                      CGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGA‑AAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTT  >  NC_000913/4542851‑4543144

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