New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	= 3614829	32 (0.940)	5 (0.160)	5/266	3.7	14.6%	coding (1163/1575 nt)	nikA	nickel/heme ABC transporter periplasmic binding protein
	?	NC_000913_3_hsa_tpiA	= 3614832	28 (0.870)					coding (1166/1575 nt)	nikA	nickel/heme ABC transporter periplasmic binding protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTGCTGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGC  >  NC_000913_3_hsa_tpiA/3614684‑3614966                                                                                                                                                     |                                                                                                                                       ttGCTGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGAt                                                                                                                                        <  1:241864/149‑1 (MQ=255)   gcaggAAAAAGCTGGTCGGACGTTGCCTTCGGGCAAAGACAGCCGCCAGTAAAAGGGTCAGCCGCTGCGCATTGTACTTTCGTTCATCGTCACCGATGCTTTAAGCAAATCTATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTg                                                                                                                                      <  2:469119/146‑1 (MQ=255)     tGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGgc                                                                                                                                    >  2:281262/1‑149 (MQ=255)        aaaaaGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCa                                                                                                                                 <  2:446642/149‑1 (MQ=255)             gCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCCTTCAGGCTGATATGCGCCAGGTTGGCGCAGATGt                                                                                                                            >  1:109267/1‑149 (MQ=255)                     gACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGa                                                                                                                    <  2:297860/149‑1 (MQ=255)                          tgcctgcGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGc                                                                                                               <  1:281262/149‑1 (MQ=255)                          ggccggcGGGCAAAGACATCCGCGAGAAAAATGGTAAGTCTCAGCGGATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAACTCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGc                                                                                                               <  2:214382/148‑1 (MQ=255)                             ctgcGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAATCTGGCGCatat                                                                                                                              >  1:124782/1‑131 (MQ=255)                             ctgcGGGCAAAGACATCCGCGAGAAAAAGGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCgaa                                                                                                            <  1:300869/149‑1 (MQ=255)                                   gCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATGCAGGCTGATATGCGCCAGATTGGCGCCGATGCCTCGCTGATTGGCGAagaagaa                                                                                                      >  2:496229/1‑148 (MQ=255)                                     aaaGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAagagag                                                                                                    >  1:124464/1‑149 (MQ=255)                                     aaaGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTCTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAagagag                                                                                                    >  1:328953/1‑149 (MQ=255)                                       aGACATCCTCGAGAAAAATGGTCGGCGGCTGCTCATTGAACTTGCGTTCATCGGCCCCGATGCGTTAAGAAAATCGATGGCTGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAAACATGTCGATGACTGGCGAAGAAGAGAGCa                                                                                                  <  2:284347/149‑1 (MQ=255)                                          cATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTa                                                                                               <  1:207750/149‑1 (MQ=255)                                            tCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTCTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATc                                                                                             <  2:328953/149‑1 (MQ=255)                                               gggagaAAAATGGTCGGCCGATGCGGATGGAACTTTCGTTCATGGGCACCGATGCGTTAAGCAAATCGATGGCGGACATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTAt                                                                                          <  2:46960/147‑1 (MQ=255)                                                   gaAAAATGGTCAGCCGCTGCGTATTGAACTTTTGTTCATGGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGAAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTc                                                                                      <  2:74841/149‑1 (MQ=255)                                                    aaaaaTGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCg                                                                                     >  1:224487/1‑149 (MQ=255)                                                      aaaTGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGCTTGGCGAAGAAGAGAGCAGTATCTATGCTCGTc                                                                                   >  1:558836/1‑149 (MQ=255)                                                         tGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATCGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAgc                                                                                >  2:92193/1‑149 (MQ=255)                                                               gCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACg                                                                          >  2:528875/1‑149 (MQ=255)                                                                          ttGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCt                                                                                                                      <  1:513091/94‑1 (MQ=255)                                                                          ttGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCt                                                                                                                      >  2:513091/1‑94 (MQ=255)                                                                                    gTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGTGTTTct                                                     >  1:329616/1‑148 (MQ=255)                                                                                      tCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCAc                                                   >  1:40409/1‑149 (MQ=255)                                                                                                    gCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGgcgc                                     <  2:362793/149‑1 (MQ=255)                                                                                                      gTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGTCGCCGCTTGTGGCATGATTTTCCACCGCACCTGGGgcgcgc                                   >  1:466480/1‑149 (MQ=255)                                                                                                             taaTCGATGGCTGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGAACTGTTTGCAGCCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGa                            <  1:291239/148‑1 (MQ=255)                                                                                                             aaaTCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGa                            <  1:528875/149‑1 (MQ=255)                                                                                                                                ttCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCCCCTGGGGCGCGCCATATGATCAACACGCCTTCCTCAGt         >  2:79281/1‑149 (MQ=255)                                                                                                                                    gggtGATATGCGGCAGATTGGGTCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTCTCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACCCCCCTGTCTCAGTTCaa     <  2:521117/146‑1 (MQ=255)                                                                                                                                       tGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATgc  >  1:558368/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                       TTGCTGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGC  >  NC_000913_3_hsa_tpiA/3614684‑3614966

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