New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	2065027 =	34 (0.880)	5 (0.130)	5/274	4.5	13.7%	coding (196/744 nt)	cobS	cobalamin synthase
	?	NC_000913_3_hsa_tpiA	2065048 =	30 (0.800)					coding (175/744 nt)	cobS	cobalamin synthase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGAAAAGTAATAATACCGCGAGAATAATGCTCGAAATCCAGTCCCTGGGACCA  >  NC_000913_3_hsa_tpiA/2064880‑2065162                                                                                                                                                    |                                                                                                                                        gACCGCTGTGGGTGTCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGCATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGt                                                                                                                                        <  1:320048/149‑1 (MQ=255) caCCGCCGGGGCTGCCTCAACGGCTATCACGCCTTATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGt                                                                                                                                        <  1:320058/148‑1 (MQ=255)         tGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAg                                                                                                                                <  1:578545/149‑1 (MQ=255)          gggTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAgt                                                                                                                               >  1:441291/1‑149 (MQ=255)           ggTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAgtg                                                                                                                              <  2:538476/149‑1 (MQ=255)                  aaaCGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCACCACCAGCACAAGTACGCTAAAAAGGGCCGCCAGAGGTGcccc                                                                                                                       >  1:647970/1‑149 (MQ=255)                        cTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCa                                                                                                                 <  1:592460/149‑1 (MQ=255)                              cGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTg                                                                                                           <  1:98721/149‑1 (MQ=255)                                  ggATCTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGAACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTgcagc                                                                                                       >  1:643165/2‑149 (MQ=255)                                                 gATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGc                                                                                        <  2:131812/149‑1 (MQ=255)                                                              ggCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGTAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCgcgc                                                                           <  2:326398/148‑1 (MQ=255)                                                                  gAAAATAGGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCCAAACAGGGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAg                                                                       <  2:623123/149‑1 (MQ=255)                                                                     aaTACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGTCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCaa                                                                    <  2:66004/149‑1 (MQ=255)                                                                         cGCCATCACAGGTATCAGCGAGGCCGTCCAGAGGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCaataat                                                                <  1:564565/149‑1 (MQ=255)                                                                          gCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATc                                                               <  2:411045/149‑1 (MQ=255)                                                                           ccATCACGGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATcc                                                              <  1:443334/149‑1 (MQ=255)                                                                                  aGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATCaaa                                                       <  2:441291/149‑1 (MQ=255)                                                                                        cggcgcGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACCCCATGCCTGCAGCACCATGAAGACCAGTCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGaaaa                                                 <  2:608730/147‑1 (MQ=255)                                                                                          gCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCCAGCAATAATCCAATCAAAGGAAAAGt                                               >  1:106292/1‑149 (MQ=255)                                                                                            gAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGAAAAGtaa                                             <  1:37108/149‑1 (MQ=255)                                                                                               gCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGAAAAGtaataa                                          >  2:439500/1‑149 (MQ=255)                                                                                                  gTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGAAAAGTAATAATAc                                       <  1:177434/149‑1 (MQ=255)                                                                                                     cAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGAAAAGTAATAATACcgc                                    <  2:386639/149‑1 (MQ=255)                                                                                                      aGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCACCCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGAAAAGTAATATTACcgcg                                   >  2:329153/1‑149 (MQ=255)                                                                                                       gATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTgca                                                                                                         >  1:599550/1‑78 (MQ=255)                                                                                                       gATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTgca                                                                                                         <  2:599550/78‑1 (MQ=255)                                                                                                             aTCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTACTCGCGCCAAGCAATAATCCAATCAAAGGAAAAGTAATAATACCGCGAGaataa                            >  2:524/1‑149 (MQ=255)                                                                                                                 cccGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGAAAAGTAATAATACCGCGAGAATAATGCt                        >  2:388898/1‑149 (MQ=255)                                                                                                                      tcatcaGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAATACCAGCCCGCTAATCGCGCCAAGCAATAATCCTTTCAAAGGACAAGTAATAATACCGCGAGAATAATGCTCGaaa                   <  2:643165/149‑1 (MQ=255)                                                                                                                        atcaGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGAACAAAAAAGACCAGCCCGCTAATCGCGCCACGCAATAATCCAATCAAAGGAAAAGTAATAATACCGCGAGAATAATGCTCGAAATc                 >  1:301012/1‑149 (MQ=255)                                                                                                                             caccagcacAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGAAAAGTAATAATACCGCGAGAATAATGCTCGAAATCCAGTc            >  2:97702/1‑149 (MQ=255)                                                                                                                                      aaGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGAAAAGTAATAATACCGCGAGAATAATGCTCGAAATCCAGTCCCTGGGAcc   >  1:315710/1‑149 (MQ=255)                                                                                                                                      aaGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATAAAAGGAAAAGTAATAATACCGCTAGAATAATGGTCGAAATCCAGTCCCTgggccg   >  1:107835/1‑146 (MQ=255)                                                                                                                                       aggaCGCTCACCAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGCAGACCAGCCCGCTAATCGCGCCAAGCAAGAATCCAATCAAAGGAAAAGTAATAATACCGCGAGAATAATGCTCGAAATCCAGTCCCTGGGACCa  <  2:218094/146‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                        GACCGCCGTGGGTGCCTAAACGGCTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGAAAAGTAATAATACCGCGAGAATAATGCTCGAAATCCAGTCCCTGGGACCA  >  NC_000913_3_hsa_tpiA/2064880‑2065162

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.