New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	3186096 =	NA (NA)	59 (1.410)	48/284	0.0	100%	intergenic (+6/‑91)	yqiG/insC1	pseudogene; fimbrial export usher family;putative membrane; Not classified; putative membrane protein/IS2 repressor TnpA
	?	NC_000913_3_hsa_tpiA	= 3725928	0 (0.000)					coding (42/996 nt)	wecH	O‑acetyltransferase for enterobacterial common antigen (ECA)

GATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGA  >  NC_000913_3_hsa_tpiA/3185950‑3186236                                                                                                                                                   |                                                                                                                                             gATTGCTTACGAGTAATCGCATTTGGCCGGGGTCGCCCGTCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGAATAAAAATAATTATTTTTACTTATCAGTTATTAg                                                                                                                                            >  2:220718/1‑149 (MQ=255)       ttACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGc                                                                                                                                      >  1:666108/1‑149 (MQ=255)                 tCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAAtc                                                                                                                            >  2:326135/1‑149 (MQ=255)                    cATTTGGCCGGGTTCGCCCGGCCAAATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTcc                                                                                                                         <  2:592751/149‑1 (MQ=255)                          gCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACaa                                                                                                                   <  2:387623/149‑1 (MQ=255)                              ggTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCaa                                                                                                               <  1:393525/149‑1 (MQ=255)                                           acaTCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCAGTTATTAGACTGGCCCCCTGAATCTCCGGACAACCAATATCACTTAAATa                                                                                                  >  2:11528/1‑149 (MQ=255)                                             atcatTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAg                                                                                                >  2:145964/1‑149 (MQ=255)                                                atTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGa                                                                                             <  1:326135/149‑1 (MQ=255)                                                  tAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATa                                                                                           <  2:10295/149‑1 (MQ=255)                                                   aaCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAg                                                                                          >  2:485585/1‑149 (MQ=255)                                                    aCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGt                                                                                         >  1:590008/1‑149 (MQ=255)                                                    aCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGt                                                                                         >  2:278528/1‑149 (MQ=255)                                                    aCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACATAAATAAGTGATAGt                                                                                         <  2:590008/149‑1 (MQ=255)                                                     cGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTc                                                                                        <  2:299050/149‑1 (MQ=255)                                                         ttAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTaa                                                                                    <  2:534288/149‑1 (MQ=255)                                                          tAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                                                                                   >  1:381408/1‑149 (MQ=255)                                                          tAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                                                                                   <  2:229402/149‑1 (MQ=255)                                                           aaTGATAAGTGGATCAGATGTATAAAAAATTAAAACTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATa                                                                                  <  1:15802/149‑1 (MQ=255)                                                             tGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTTTTAATACt                                                                                >  1:160947/1‑149 (MQ=255)                                                                tAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGt                                                                             <  1:485585/149‑1 (MQ=255)                                                                 aaGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                                                                            <  2:666108/149‑1 (MQ=255)                                                                   gTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTAAAATAAGTGATAGTCTTAATACTAGtttt                                                                          >  1:439227/1‑149 (MQ=255)                                                                       aTCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTTATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGa                                                                      <  1:493958/149‑1 (MQ=255)                                                                                    aaaattaaaattAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgag                                                         <  2:320493/149‑1 (MQ=255)                                                                                    aaaattaaaattAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgag                                                         >  1:116874/1‑149 (MQ=255)                                                                                              ttaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATAAGACTGGCCCCCTGAATCTCCAGACAACCAATATCCCTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGCTGTTGATGAACTgatgat                                               >  1:476774/1‑149 (MQ=21)                                                                                                  ccACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgattgat                                           <  2:381408/149‑1 (MQ=255)                                                                                                               aTTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCTGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTTACTTTTCTTTGGAGAACAGCTGGTTTATGTCTTAGGGTCgg                              >  2:503971/1‑149 (MQ=18)                                                                                                                    ttAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGaaa                         >  2:617376/1‑149 (MQ=255)                                                                                                                       aaaaTATTTATTGTTCATTATCCGTTGTTAGACTGGCCCCCTGAATCTCCCGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGc                      <  1:635816/149‑1 (MQ=255)                                                                                                                                       actATCCGTTATTAGACTGGCCCCCTGAATCTCCTGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGCTGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACcacac      <  1:503971/147‑1 (MQ=255)                                                                                                                                        ttATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACcacaca     >  2:316162/1‑149 (MQ=38)                                                                                                                                         tATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCTCTTAAATAAGTGATAGTCTTAATACTATTTTGTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAg    <  1:11528/149‑1 (MQ=255)                                                                                                                                           tCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGa  >  1:303621/1‑149 (MQ=35)                                                                                                                                              gTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGaa                                                       <  2:205924/93‑1 (MQ=34)                                                                                                                                              gTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGaa                                                       >  1:205924/1‑93 (MQ=34)                                                                                                                                                   |                                                                                                                                             GATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCACATCATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGA  >  NC_000913_3_hsa_tpiA/3185950‑3186236

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ 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.