New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4542683 =	31 (1.110)	5 (0.190)	4/266	3.2	15.2%	intergenic (+49/‑433)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913	4542997 =	27 (1.030)					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/above cutoff threshold.

TCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAAT  >  NC_000913/4542862‑4543138                                                                                                                                        |                                                                                                                                              tCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGCGTCAATATGGCCCCAAATGTTTCa                                                                                                                                        >  1:91474/1‑143 (MQ=255) tCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCa                                                                                                                                        >  1:480930/1‑143 (MQ=255)            tcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggggg                                                                                                                             >  1:48668/1‑143 (MQ=255)              cAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTCTCATCTTTTGGGGGaa                                                                                                                           >  2:228508/1‑143 (MQ=255)               aaaaaCCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCACATGTTTCATCTTTTgggggcac                                                                                                                          >  1:206621/1‑140 (MQ=255)                    ccacaTCATGAAATATAATGATCTCTAAATTAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                                                                                                     <  2:98358/143‑1 (MQ=255)                           aTGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtgt                                                                                                              <  1:241490/143‑1 (MQ=255)                            tGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCTGCGATGCGCTCCTCTCTGAGTCAAAATGGCCCCAAATGTTTCATCTATTGGGGGAAAACTGTGCAGTGtt                                                                                                             >  1:158483/1‑143 (MQ=255)                              cAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAATATGGCGCCAAATGGTTCATCTTTTGGGGGAAAACTGTGCAGTGTTgg                                                                                                           <  2:396316/143‑1 (MQ=255)                                   taAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTc                                                                                                      <  1:356901/143‑1 (MQ=255)                                           tataAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGCAAACTGTGCAGTGTTGGCAGTCAAACTCGt                                                                                              >  2:96982/1‑143 (MQ=255)                                           tataAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTATATGAGTCAAAATGGCACCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAAATCGt                                                                                              >  2:453282/1‑143 (MQ=255)                                            ataaataaAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGCCAAACTCGtt                                                                                             >  2:414265/1‑143 (MQ=255)                                                ataaaGATAACAATAGAATATTCAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAACTGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCCCACTCGTTTacc                                                                                         >  1:193179/1‑142 (MQ=255)                                                      aTAAGAAGAGAATATTAAGCCAAAAACTAATCTGAAAAAGTTTGTACGCGATGCCTTCCTCTATGATTCAAAATGGCACCAAATGTTTCATCTTTTTGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTTGTTTacaaaacaa                                                                                   <  2:39720/143‑1 (MQ=255)                                                           aaTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCACAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACACAgtgt                                                                              >  1:241656/1‑143 (MQ=255)                                                            aTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTa                                                                             <  2:318111/143‑1 (MQ=255)                                                                 atatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa                                                                        <  2:273109/143‑1 (MQ=255)                                                                        gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCCAACTCGTTTACAAAACAAAGTGTACAGAACGACTGc                                                                 <  1:313551/143‑1 (MQ=255)                                                                                      tGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAg                                                   >  1:342681/1‑143 (MQ=255)                                                                                        aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAGACTCGTTTACAAAACAAAGTGTACAGAACGTCTGCCCATGTCGATTTAGaa                                                 <  2:253562/143‑1 (MQ=255)                                                                                                  ccctcgATGCTTTCCTCTATGAGTGAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGtttttt                                       <  2:112334/139‑1 (MQ=255)                                                                                                    cgcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACACAGTGTACAGAACGATTGGTCATGACAATATAGAAATAGTTTTTTga                                     >  1:179583/1‑143 (MQ=255)                                                                                                       gATGCTTTCCTCTATTATTCAAAATGGCCCCAGATGTTTCAGCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACGCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaag                                  <  1:7722/143‑1 (MQ=255)                                                                                                                   tATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGCAATAGTTTTTTGCACGGAAAGCAGCATg                      >  1:74149/1‑143 (MQ=255)                                                                                                                        gTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAAt                 >  1:300838/1‑143 (MQ=255)                                                                                                                           aaaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTaa              <  2:480930/143‑1 (MQ=255)                                                                                                                                   tactaaTGTTTCATCTTTTGGGGGAAAACTGTGCAGGGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATTTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTgg      <  1:47713/139‑1 (MQ=255)                                                                                                                                   cccAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCGCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAAATCTga      >  2:223704/1‑142 (MQ=255)                                                                                                                                    ccAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGc     >  1:137825/1‑143 (MQ=255)                                                                                                                                       aaTGTTTCATCTTTTGGGAGAAAACTGTTCAGTGTTGGCAGTCAATCTCGTTTACAAAACAAAGTGTATAGAACGACTGCCCATGTCGATTTAGAAATATTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAAt  <  2:71133/143‑1 (MQ=255)                                                                                                                                        |                                                                                                                                              TCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAAT  >  NC_000913/4542862‑4543138

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