New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	4542683 =	40 (0.840)	10 (0.220)	9/266	4.2	21.4%	intergenic (+49/‑433)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913_3_bme_pgi	4542997 =	36 (0.810)					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.

CATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTT  >  NC_000913_3_bme_pgi/4542854‑4543145                                                                                                                                                |                                                                                                                                                     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)                                                                                                                                                aTGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCgttgtt  <  1:250923/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     CATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTT  >  NC_000913_3_bme_pgi/4542854‑4543145

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.