New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	= 4542672	30 (0.780)	3 (0.080)	3/254	NT	9.8%	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.750)					intergenic (+353/‑129)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

TCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGG  >  NC_000913_3_hsa_tpiA/4542821‑4543115                                                                                                                                                    |                                                                                                                                                    tCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAt                                                                                                                                                    <  1:65035‑M2/149‑2 (MQ=255)   tgtgGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTAAAAAtag                                                                                                                                                  <  1:226057‑M2/149‑4 (MQ=255)    gtgGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggc                                                                                                                                                 >  1:154409‑M2/1‑145 (MQ=255)            aCGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATACAAATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATTGCCCCaaacg                                                                                                                                         >  2:329005/1‑147 (MQ=255)               cATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATATACGTCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGGCCGCGATGCTTTCCTCTCTGAGTCAAAATGGCCCCAAATGttt                                                                                                                                      <  1:553270/149‑1 (MQ=255)                aTCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTc                                                                                                                                     >  1:284465/1‑149 (MQ=255)                        cATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTg                                                                                                                             <  2:642106/149‑1 (MQ=255)                        cATGCTTCACTCCAAAAACCACCTCATGCAATATAAACATGTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGTCCCCAAATGTTTCATCTTTTg                                                                                                                             <  2:642009/149‑1 (MQ=255)                         aTTCTTCTGTCATAAAACCTCCTGATGAAATATACATATCTATAAATAAAGATAAAAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTgg                                                                                                                            <  1:332351/149‑1 (MQ=255)                                  tcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACt                                                                                                                   >  1:6821/1‑149 (MQ=255)                                            cccTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAAAAAACTGAAAAAGGTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGtt                                                                                                         <  1:28497/148‑1 (MQ=255)                                              cTCAAGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTgg                                                                                                       <  1:251612/149‑1 (MQ=255)                                                 aTGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCACCTTTTGGGGAAAAACGATGCAGTGGTGGCAg                                                                                                    >  2:414871/1‑149 (MQ=255)                                                                   taaataaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa                                                                                  >  2:668769/1‑149 (MQ=255)                                                                            aTAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAATGTTGGCAGTCAAAATCGTTTACAAAACACAGTGTa                                                                         >  2:686843/1‑149 (MQ=255)                                                                              aaCAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACa                                                                       <  1:527632/149‑1 (MQ=255)                                                                                 aaTAGAATATTAAGCCAACAACTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa                                                                    <  1:197414/149‑1 (MQ=255)                                                                                 aaTAGAATATTAAGCCAACAAATAAAATGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATAGCGCCAAATGTGTCATCTTTTGCGGCAAAACTGTACAGGGTTGGCAGTCAGCAGCGGTTACACAACACAGTGTACAGac                                                                    >  2:130819/1‑148 (MQ=255)                                                                                   tAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgttttttttatttttttatcaataagtgaaatcacacctcttctttttttaattttt                                                                  <  2:599734‑M2/149‑84 (MQ=255)                                                                                   tAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACg                                                                  >  1:403177/1‑149 (MQ=255)                                                                                             aGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgag                                                        <  2:2519‑M2/149‑94 (MQ=255)                                                                                                 aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCCGTGCAGCCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGa                                                    >  1:568320/1‑149 (MQ=255)                                                                                                      aTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAACTGTTTCATCTTCTGGGGGAAAACTAAGCCGTGTTAGCCGTCAAACTCGTTTACAAAACAAAGCGTACAGAACGACTGCCCATGTCGATCTAg                                               >  1:695691/1‑149 (MQ=255)                                                                                                         aaCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGaaa                                            >  2:696162/1‑149 (MQ=255)                                                                                                                  aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGtttttt                                   <  1:592913/149‑1 (MQ=255)                                                                                                                     ttGGCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTGCATCTATTGGTGGAAAAGTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaa                                <  2:327617/149‑1 (MQ=255)                                                                                                                         acgcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAAATGTGCAGTGTTGTCAGTCAAACTCGTTTACAAAAAAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaagga                            <  1:269930/148‑1 (MQ=255)                                                                                                                           gcgATGTTTTCCTCTCTGAGTCAAAAGGGCCCCAAATGTTTTATCTTTTGGGGGAAAAATGTGCCGTTTTGGCAGTGCAACTCGTTTACAAATTAAAGTGTAAAGTACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaaggaaa                          <  2:223742/149‑1 (MQ=255)                                                                                                                                gctctgctcgATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAagcagc                     <  1:195417/143‑1 (MQ=255)                                                                                                                                       tctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaaggaaa                          <  1:704419/137‑1 (MQ=255)                                                                                                                                       tctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaaggaaa                          >  2:704419/1‑137 (MQ=255)                                                                                                                                       tctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGaaaa              >  1:33981/1‑149 (MQ=255)                                                                                                                                              gTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagcagcattactggcggta       <  1:607519‑M2/149‑143 (MQ=255)                                                                                                                                              gTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAGTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAAc       >  1:147230/1‑149 (MQ=255)                                                                                                                                                   aaTGGCCCCAAATGTGTCATATTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTgg  <  1:129022/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGG  >  NC_000913_3_hsa_tpiA/4542821‑4543115

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.