New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	= 4542672	26 (0.710)	3 (0.090)	3/256	NT	10.5%	intergenic (+57/‑425)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	NC_000913_3_hsa_tpiA	= 4542968	27 (0.790)					intergenic (+353/‑129)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

TTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTG  >  NC_000913_3_hsa_tpiA/4542820‑4543114                                                                                                                                                      |                                                                                                                                                   ttCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCaaaa                                                                                                                                                    <  1:345588/149‑1 (MQ=255)  tCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAAT‑AGAAGATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAt                                                                                                                                                   >  1:581170‑M2/1‑148 (MQ=255)               gCAACTTACTCATTCTACTATCCAAAAACCACATCATTCAATATAAACATCTATAAATAAAGAAAACAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGtt                                                                                                                                      <  2:313221/149‑1 (MQ=255)                cATCTTTCTTAGTCGTCTATCCACAAACCACCTCATGCAATAAAAACAGCTAGAACTAAAGATAACAAT‑AAAATATTAAGCCAACACATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGttt                                                                                                                                     <  2:12086/149‑1 (MQ=255)                cATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAAT‑AGAATATTAAGCCAACAAATAAACTGGAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGttt                                                                                                                                     <  1:457373/149‑1 (MQ=255)                          atccttctCTCCAAAAACCACCACATGCAATATAAACATCTATAAATAAAGATAACAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTgg                                                                                                                           <  1:419354/146‑1 (MQ=255)                          aTTTTTCTCTCCAAAAACCACCTCATTCAATATAAACATCTATAAATAAAGTTAACAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTgg                                                                                                                           <  2:359762/149‑1 (MQ=255)                          aTTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAAT‑AGAATATTAAGCCAACAAATAACCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTgg                                                                                                                           <  2:635903/149‑1 (MQ=255)                                ctctcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCACCAAATGTTTCATCTTTTGGGGGaaa                                                                                                                     >  1:613622/1‑149 (MQ=255)                                 tctcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAAAGTTTCATCTTCTGTGGGaaca                                                                                                                    >  1:122828/1‑147 (MQ=255)                                        aaaCCACCTCATGCAATATAAACATCTATAAATAAAGATAACAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCa                                                                                                             <  1:493508/149‑1 (MQ=255)                                           ccaccTCATGCAATATAAACATCTATAAATAAAGATAACAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAgtg                                                                                                          >  1:496308/1‑149 (MQ=255)                                               cTCATGCAATATAAACATCTAAAAAAAAAGAT‑ACACTCATAGCATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTTTGCAGCGTTgg                                                                                                      >  2:293486/1‑149 (MQ=255)                                                               aTCTATAAATAAAGATAACAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAATCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaa                                                                                      <  1:607712‑M2/149‑63 (MQ=255)                                                                tCTATAAATAAAGATAACAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTac                                                                                     >  1:247671/1‑149 (MQ=255)                                                                      aataaaGATAACAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTCTGCGTCAAAATCGCCCCAAACGTGTTAATATGCCGTTTAAAACACTGCAGTCTCAGGCGTCAAACTCGTTTacaaaaca                                                                               >  1:679124/1‑149 (MQ=255)                                                                       ataaaGATAACAAT‑AGAATATTAAGCCAAGAAATAAACTGAAAAAGTTTGTCCACGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTTGCAGTCAAACTCGTTTacaaaacga                                                                              >  2:313534/1‑147 (MQ=255)                                                                        taaaGATAACAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctc                                                                             <  1:635450‑M2/149‑72 (MQ=255)                                                                        taaaGATAACAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaaa                                                                             <  2:357647/149‑1 (MQ=255)                                                                                aCAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCGATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAg                                                                     >  1:451795/1‑149 (MQ=255)                                                                                              aacccAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCAt                                                        <  1:284139/146‑1 (MQ=255)                                                                                              aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCGCAAATGTTTCATCTGTGGGGGGAAAACGGTGCAGTGTTGGCAGTCAGACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCag                                                        >  2:548226/1‑148 (MQ=255)                                                                                                   aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGa                                                   >  1:22603/1‑149 (MQ=255)                                                                                                    acaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTGTACAAAACAAAGTGTACAGAACAACTGCCCATTGCGAt                                                  >  1:234364/1‑149 (MQ=255)                                                                                                    acaaAGAAACTGAAAAAGTTGGTCCGCGATGCTTTCCTCTATGAGTCAAAGTGGCCCCAAATGTTTCATTTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCACACTCGTTTACAAAACGAAGTGTACAGAACGACTGCCCATGTTGAt                                                  <  2:427066/149‑1 (MQ=255)                                                                                                           aaCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGGtttgggaa                                           >  2:543629/1‑144 (MQ=255)                                                                                                             cTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtattgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgatattatttctttccc                                         <  2:425025‑M2/149‑108 (MQ=255)                                                                                                                aaaaaGATAATCCGCGATGCTTTCCACAATGAGTCAAAATGGCCACAAACGAATCATGTCCTGGGAGAAAACTGTGCAGATTTGGCAGTCAAGCTCGTTTACAAAACAAAGTGTACAGAACGAATGCCCAGGTCGATTTAGaagtgcaa                                      >  2:531739/1‑143 (MQ=255)                                                                                                                       ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaa                               <  1:690649/149‑1 (MQ=255)                                                                                                                                    tttCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATGTAGAAATAGATGTTTAACAGGAAAGAAGCAt                  >  1:665387/1‑149 (MQ=255)                                                                                                                                               aGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTaaaa       >  2:230929/1‑149 (MQ=255)                                                                                                                                               aGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTaaaa       <  1:533081/149‑1 (MQ=255)                                                                                                                                                gTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAATGTTTACGGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAAc      >  1:310366/1‑149 (MQ=255)                                                                                                                                                    aaaTGGTCCCAAATGTTTCATCTTTTGGGGTAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTg  <  2:150541/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAAT‑AGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTG  >  NC_000913_3_hsa_tpiA/4542820‑4543114

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.