New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 4542691	41 (0.860)	9 (0.200)	9/266	4.2	19.9%	intergenic (+57/‑425)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913_3_bme_pgi	= 4542987	34 (0.760)					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.
Rejected: Frequency below/above cutoff threshold.

GACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTG  >  NC_000913_3_bme_pgi/4542850‑4543133                                                                                                                                          |                                                                                                                                                   gACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGAAATATAAACATCTATTAATAAAGATAACAATATAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaatt                                                                                                                                         <  2:542213/149‑3 (MQ=255)     cATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGttt                                                                                                                                     <  1:752410/149‑1 (MQ=255)      aTCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAATTGtctt                                                                                                                                    <  1:63524/149‑4 (MQ=255)      aTCTTCCTCATTCTTATGTCCACAAAACAGTGCACACAAAACAAATATTTATAAAGAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTc                                                                                                                                    <  2:345849/149‑1 (MQ=255)      aTCTTCCTAATTCTCCTCTCCAAAAACCACCTAATGCAATATAAACATCTCTAAATCAAGATCACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCAGCGATGCTTTTCTCTATGTGGTGAAGTTGCCCCAATTGtctt                                                                                                                                    >  2:590851/1‑146 (MQ=255)             tCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCtttt                                                                                                                             <  1:74372/149‑1 (MQ=255)                   ttctCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTTGGGGa                                                                                                                       >  2:714434/1‑149 (MQ=255)                        tcCAAACACCACCTCATGCAATATAAACATCTATAACTAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAACAAGTTTGTCCGCGATGCTTTCCTCCATGAGTCAAAATGCCCCCAAATGTTTCAACTCTTGGAGGTAAACt                                                                                                                  >  2:111447/1‑149 (MQ=255)                        tcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACt                                                                                                                  >  1:612207/1‑149 (MQ=255)                           aaaaaCCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgtg                                                                                                               <  1:146684/149‑1 (MQ=255)                                   ccTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTg                                                                                                       <  1:355212/149‑1 (MQ=255)                                    cTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTgg                                                                                                      <  1:732382/149‑1 (MQ=255)                                           aaTATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCaa                                                                                               >  1:669322/1‑149 (MQ=255)                                                    aTCTAAAAATAAAGATAAAAAAAAAATAATAAGCCAACAAAAAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTAAGAATCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTa                                                                                      <  2:329841/149‑1 (MQ=255)                                                            ataaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTCGTCGCCGTCAAACTTCTTTacaaaacaa                                                                              >  2:836814/1‑149 (MQ=255)                                                             taaaGATAACAAAAGATTATTAATCCAACAAATAAACTTAAAACGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaaa                                                                             <  1:132055/149‑1 (MQ=255)                                                               aaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgt                                                                           <  1:833891/149‑1 (MQ=255)                                                                  ctAACAGTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTCCAGTGTTGGGCGTCACACTCGTGTACAAAACCAACTGTa                                                                        >  2:286867/2‑149 (MQ=255)                                                                    aaCAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACa                                                                      <  1:328244/149‑1 (MQ=255)                                                                            aaTATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGtt                                                                                                        <  2:331200/107‑1 (MQ=255)                                                                            aaTATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTAAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGtt                                                                                                        >  1:331200/1‑107 (MQ=255)                                                                                    gCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATgc                                                      <  1:585139/149‑2 (MQ=255)                                                                                      caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCCAATGTTTCATCTTTTGGGGGAAAACTGTTCAGTGTTGGCAGTCAACCTAGTTTACAAAACAAAGTGTCCAGAACGACTGCGCATGTCg                                                    >  1:223922/1‑149 (MQ=255)                                                                                      caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCg                                                    >  2:241824/1‑149 (MQ=255)                                                                                       aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGa                                                   >  2:319185/1‑149 (MQ=255)                                                                                       aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACGGTGCAGTGTTGGCAGTCAAACTCGTTTAGAAAACAAAGTGTACAGAACGACTGCCCAGGTCGa                                                   >  1:744751/1‑149 (MQ=255)                                                                                        aacactAAACTTAAAAAGCTTGTCCGCGATGGTTCCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGCCGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAt                                                  <  1:207430/144‑1 (MQ=255)                                                                                          aaaTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAttt                                                <  1:836814/149‑1 (MQ=255)                                                                                               aaCTGAAAACGTTTGTCCGCGATGCTTTCCTCTATGAGCCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTGGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGaaa                                           <  1:700054/149‑1 (MQ=255)                                                                                                  tGAAAAAGTTTGTCCGCGATGCTTTACTCTATGAGTCAAAATGGCCGCAAATGTTTCATCTTTTGGCAGAAAACTGTGCAGTGTAGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACTACTGCCCATGTCGATTTACAAATAg                                        >  1:757000/1‑149 (MQ=255)                                                                                                    aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGtt                                      >  2:446590/1‑149 (MQ=255)                                                                                                        aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGtttttt                                  <  1:369153/149‑1 (MQ=255)                                                                                                                        tttCCTCTATGAGTCAAAATTGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGc                                                            <  1:86572/107‑1 (MQ=255)                                                                                                                        tttCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGc                                                            >  2:86572/1‑107 (MQ=255)                                                                                                                                   aGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTaaaa       <  1:287741/149‑1 (MQ=255)                                                                                                                                     tCAAAATGGCCACAGAAGTTTCATCTTTTAGGGGAAAACTGTGCAATGTCGGCAGTCAAACTTGTTTACAAAACAAAGTGTATAGTACGACTGCCCATGTAGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAAct     >  2:755998/1‑149 (MQ=255)                                                                                                                                        aaaCGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTg  <  1:82212/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                   GACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTG  >  NC_000913_3_bme_pgi/4542850‑4543133

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.