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.

CTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGAAAAGTAATAATACCGCGAGAATAATGCTCGAAATCCAGTCCCTGGGACCAGCGACGCGGTACAGGCAGACGCGTAATAAATGA  >  NC_000913_3_hsa_tpiA/2064903‑2065195
                                                                                                                                                 |                                                                                                                                                   
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)
                                                                                                                              gccgccAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGAAAAGTAATAATAcc                                                                       >  2:388041/1‑98 (MQ=255)
                                                                                                                              gccgccAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGAAAAGTAATAATAcc                                                                       <  1:388041/98‑1 (MQ=255)
                                                                                                                                                cacCATCCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGAAAAGTAATAATACCGCGAGAATAATGCTCGAAATCCAGTCCCTGGGACCAGCGACGCGGTACAGGCAGACGCGTAATAAATga  <  1:389594/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTATCACGCATGATTTCCAACATGCGATCGCGGCTACGTGCAGAAAATACGCCATCACAGGTATCAGCGAGGCCGTCCAGATGGAATCCCCCGGTCATCAGCACCAGCACAAGTACGCTAAACAGTGCCGCCAGTGGTGCGCCACACCATGCCTGCAGCACCATGAAGACCAGCCCGCTAATCGCGCCAAGCAATAATCCAATCAAAGGAAAAGTAATAATACCGCGAGAATAATGCTCGAAATCCAGTCCCTGGGACCAGCGACGCGGTACAGGCAGACGCGTAATAAATGA  >  NC_000913_3_hsa_tpiA/2064903‑2065195

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.