New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_pgi	4542710 =	43 (0.810)	5 (0.100)	5/252	NT	10.9%	intergenic (+49/‑433)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913_3_hsa_pgi	4543024 =	42 (0.840)					intergenic (+363/‑119)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

GACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTT  >  NC_000913_3_hsa_pgi/4542877‑4543169                                                                                                                                                    |                                                                                                                                                  gACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAAt                                                                                                                                                  >  1:41843/1‑149 (MQ=255)   cGCATCTTCCTCATTCTTCTTTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGt                                                                                                                                                <  1:157406/149‑1 (MQ=255)    gCATGTTCCTCAAGCCTATCTCCAAAAACCAACTCATGCAATATAAACATCTATAAATAAAGATAACAAAAGCATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCACGATGCTTTCCTCTCTAGGTTAAAATGGCCCCAAATGtt                                                                                                                                               >  2:723824/1‑149 (MQ=255)      aTCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGtttt                                                                                                                                             >  2:493603/1‑148 (MQ=255)      aTCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAATTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTc                                                                                                                                             >  1:318352/1‑149 (MQ=255)               aTTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAATTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTgg                                                                                                                                    <  1:72910/149‑1 (MQ=255)               aTTCTTCTCTCCAAAAACCACACCATGCAATATAACCCTCTACAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCtttttg                                                                                                                                    >  2:171242/1‑147 (MQ=255)                 tcttctCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTCTGAGTCAAACTAGCCCCAAATGTTTCTTCTTCTcggg                                                                                                                                  >  1:959997/1‑149 (MQ=255)                   ttctCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCGCTATGAGTCAAAATGGCCCCAAACGTTTCATCTTTTGGGGGa                                                                                                                                >  2:436954/1‑149 (MQ=255)                      tctcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaaa                                                                                                                             >  1:843617/1‑149 (MQ=255)                      tctcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAAGAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGATTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaaa                                                                                                                             <  2:881236/149‑1 (MQ=255)                       ctcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCCAATGTTTCTTCTTTCGGGGGAAAAc                                                                                                                            >  2:614768/1‑149 (MQ=255)                        tcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAAAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATCGACCCAAATGTTTCATCTTTTGGGGGaaaatt                                                                                                                           >  2:45989/1‑147 (MQ=255)                          cAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACACTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCAAAATGTTTGATATTTTGGGGGAAAACtgt                                                                                                                         >  1:37234/1‑149 (MQ=255)                          cAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgt                                                                                                                         >  1:646725/1‑149 (MQ=255)                             aaaCCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCa                                                                                                                      <  2:843617/149‑1 (MQ=255)                                      cATGAAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCa                                                                                                             <  1:214537/149‑1 (MQ=255)                                      cATGAAAAATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCa                                                                                                             <  2:142103/149‑1 (MQ=255)                                       aTGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCCTCTTTTGGGGGAAAACTCTGCCGTGTTGGCat                                                                                                            >  2:605834/1‑148 (MQ=255)                                       aTGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAg                                                                                                            <  1:931293/149‑1 (MQ=255)                                             taaaaaCAACTATAAAAAAAAATAACAAAATAATATTAAACCAACAAAAAAACTGAAAAAGTTTGTCCGCGATTCTTTCCTGTAAGAGTCAAAATTGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAAc                                                                                                      <  2:367313/146‑1 (MQ=255)                                                 tacatatcaaacagttagatgctttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGt                                                                                                  >  1:168965‑M2/100‑149 (MQ=255)                                                  aCATCGATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGtt                                                                                                 <  1:945650/149‑1 (MQ=255)                                                    aTCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTa                                                                                               <  2:41843/149‑1 (MQ=255)                                                       tataAAGAAAGATAACAATAGAATAATAAGCCAACAAATGAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaa                                                                                            <  1:493619/149‑1 (MQ=255)                                                         aaacagttagatgctttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa                                                                                          >  1:592588‑M2/92‑149 (MQ=255)                                                               aaGATAACATTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTTCAGTTTTGGCAGTCAAACTCGTTTACAAAACAAAgt                                                                                    <  2:816604/149‑1 (MQ=255)                                                                  aTAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTa                                                                                 >  1:802498/1‑149 (MQ=255)                                                                  aTAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTa                                                                                 <  1:556850/149‑1 (MQ=255)                                                                  aTAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTATTGGCAGTCAAACTCGTTTACAAAACAAAGTGTa                                                                                 >  1:331441/1‑149 (MQ=255)                                                                       cattgAATATTAAGCTAACAAATAAACTGAAAAAGTTTGGTCGCGATGCTTTCCTCTATGAGTCAAAATGGCCTCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa                                                                            <  1:306741/145‑1 (MQ=255)                                                                       aaTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa                                                                            >  2:579076/1‑149 (MQ=255)                                                                       aaTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa                                                                            <  1:529696/149‑1 (MQ=255)                                                                            aaTATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACt                                                                       >  2:816942/1‑149 (MQ=255)                                                                                ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGccc                                                                   <  2:646725/149‑1 (MQ=255)                                                                                 tAAGCCAAGAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGTGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCa                                                                  <  2:436972/149‑1 (MQ=255)                                                                                  ttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatgtgtacagtttgggcccaaATGTTTCATCTTTTGGTGGAAAACTGTGAAGTGTTGGCAGTCAAACTCGTCTACACAACAAGGTTTACAGAACGACtgcccgt                                                                 >  2:166437‑M2/67‑147 (MQ=255)                                                                                  aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGTTTTCCTCTAGGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCAt                                                                 <  1:80183/149‑1 (MQ=255)                                                                                     taatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTc                                                              <  2:168965‑M2/86‑1 (MQ=255)                                                                                      caacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCATCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAATCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTACCCATGTCg                                                             >  1:144826/1‑149 (MQ=255)                                                                                       aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCATACTCGTTTACAAAACAAAGTGTACAGAACGACTGCACCTGTcta                                                            >  2:382209/1‑147 (MQ=255)                                                                                             tAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGa                                                      >  1:546090/1‑149 (MQ=255)                                                                                                       aaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGttttt                                            >  1:644626/1‑149 (MQ=255)                                                                                                            tGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCAGCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTAAAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaa                                       <  2:318352/149‑1 (MQ=255)                                                                                                                  taacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaaggaaag                                 >  1:569911‑M2/35‑149 (MQ=255)                                                                                                                       ctgtCCTCTATGAGTAAAGATGGCCCCAAATTTTTCATCTCTTGGGAGTGAACTGTGCAGTGTTGTCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAgcagca                            <  1:636848/146‑1 (MQ=255)                                                                                                                              ctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAAt                     >  2:294908/1‑149 (MQ=255)                                                                                                                               tATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAAtt                    >  2:704345/1‑149 (MQ=255)                                                                                                                                 tGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGATAAACTGTTGAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTaa                  <  1:689282/149‑1 (MQ=255)                                                                                                                                  gAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTaaa                 <  1:294908/149‑1 (MQ=255)                                                                                                                                      cAAAATGGCCCCAAATGTTTTATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAActc             <  2:427711/149‑1 (MQ=255)                                                                                                                                        aaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTg           <  1:171242/149‑1 (MQ=255)                                                                                                                                         tttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTgg          <  2:786999‑M2/138‑1 (MQ=255)                                                                                                                                                ccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCCGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTAAAAGGTAAACAGCATGAAACTTAAAACTCtcgcaacagt   >  2:430526‑M2/5‑145 (MQ=255)                                                                                                                                                ccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCAAGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATAAAAAACCTGGCAATCgt   >  1:511814‑M2/5‑149 (MQ=255)                                                                                                                                                 caaATGTTTCATCTTTTGGGGGCACCATGTGCGATGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCgtt  >  2:972129‑M2/4‑149 (MQ=255)                                                                                                                                                 caaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAGGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCgtt  >  1:4013‑M2/4‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTT  >  NC_000913_3_hsa_pgi/4542877‑4543169

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.