New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA = 67155836 (1.040)5 (0.160) 5/262 3.8 13.0% coding (48/1032 nt) holA DNA polymerase III, delta subunit
?NC_000913_3_hsa_tpiA = 671568 34 (1.070)coding (38/1032 nt) holA DNA polymerase III, delta subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTCCAGTCAGTGTTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCTTGG  >  NC_000913_3_hsa_tpiA/671421‑671709
                                                                                                                                                    |                                                                                                                                             
ttCCAGTCAGTGTTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAg                                                                                                                                              <  2:397684/149‑1 (MQ=255)
              gggATCAATGGAAAAAGTGTGTTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACACTGTATTGTTACCAAGTAGAATATACGCCGCGCGCAGCCTTTCATTGAGCTGCGCGCGGAgtt                                                                                                                                <  1:448034/149‑1 (MQ=255)
                gATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAgttgt                                                                                                                              >  1:527999/1‑149 (MQ=255)
                caTCAACGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCCGCTAACAGACTAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAgttgt                                                                                                                              <  1:115097/148‑1 (MQ=255)
                    aaTGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGa                                                                                                                                               >  2:307949/1‑128 (MQ=255)
                    aaTGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGa                                                                                                                                               <  1:307949/128‑1 (MQ=255)
                            aaGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTCACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGCGCTGCGCGCGGAGTTGTTCCGGGTACAAc                                                                                                                  <  2:497787/149‑1 (MQ=255)
                                              aTCCTTGTGCCGCAGCTACCTGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCa                                                                                                <  1:44786/149‑1 (MQ=255)
                                                     ggacgtaGCTACCAGACGAACATGCG‑CCTGGCTTTCCTGCAATCACCGACGATCGTCACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGTGCGGAGTTGTGCCGGGTATATTCGAATCACCAGGTACCCTGCgtggt                                                                                         <  2:265715/143‑1 (MQ=255)
                                                       ccGCAGCTACCTGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGtggtgg                                                                                       >  2:415579/1‑149 (MQ=255)
                                                       ccGCAGCTACCTGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGtggtgg                                                                                       <  1:29137/149‑1 (MQ=255)
                                                          cAGCTACCTGACGAACA‑GCGTCCTGGCTTTCCTGCTATAACATAGGATCGTTACCAAGTATAAGATACGCCGCGCGCAGCCCTTCATTGAGCTACGCGCGGAGTTGTTCCGGGTTCAACCGAATCATCAGCTACCCAGCGTGGTGgaga                                                                                    <  2:479246/149‑1 (MQ=255)
                                                             cTACCTGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGAcgc                                                                                 <  2:527999/149‑1 (MQ=255)
                                                               aCCTGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGtggtgg                                                                                       >  2:446008/1‑141 (MQ=255)
                                                               aCCTGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGtggtgg                                                                                       <  1:446008/141‑1 (MQ=255)
                                                                  tGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAAtcatca                                                                                                        <  2:577775/121‑1 (MQ=255)
                                                                  tGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAAtcatca                                                                                                        >  1:577775/1‑121 (MQ=255)
                                                                   gatgcaca‑cCGGCCTGGCTTGCCTGCAATAACAGCGGCTCGTTATCTAGTAAAATATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAg                                                                           <  1:73431/144‑1 (MQ=255)
                                                                                 cTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGCAGt                                                              >  2:218130/1‑149 (MQ=255)
                                                                                   ggCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTAt                                                            <  1:154696/149‑1 (MQ=255)
                                                                                   ggCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTAt                                                            <  1:560742/149‑1 (MQ=255)
                                                                                           tGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCgtggt                                                                                         <  1:533841/112‑1 (MQ=255)
                                                                                           tGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCgtggt                                                                                         >  2:533841/1‑112 (MQ=255)
                                                                                              aaTAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTGATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGTTc                                                 >  2:559880/1‑149 (MQ=255)
                                                                                              aaTAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTCCCCAGCGAGGTGGAGACGCGTGCAGGCGTTGCCGGAGTACTCGTTCTgggc                                                 >  2:500467/1‑147 (MQ=255)
                                                                                                 aaCAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCgtggtg                                                                                        <  2:347489/107‑1 (MQ=255)
                                                                                                 aaCAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCgtggtg                                                                                        >  1:347489/1‑107 (MQ=255)
                                                                                                            ggtACCAAGTAAAAGATACGCCGGGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTCCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGttcttc                                   <  1:60718/147‑1 (MQ=255)
                                                                                                                ccAAGTAAAAGATACGCCTCGCGCAGCCCTTCATTGAGCCGCGCGCGTAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGTCGTTGCCGGAGTATCCGGTGTGGTCGACGTCTGTTCTtcgtcg                               <  1:461204/149‑1 (MQ=255)
                                                                                                                ccAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTCCCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGCTGTGGTCGACGTCTGTTCTtcgtcg                               <  1:460599/149‑1 (MQ=255)
                                                                                                                        aaGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTGGGAACGAat                       >  1:38595/1‑149 (MQ=255)
                                                                                                                                           ccTTCTTTGTGTTGGGCGCGGTTTTGTGTCTGCTACTACTGAATGGTCTGTTATGGAGCGTGGTGGTGACGCGTTCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCtt    <  2:571158/149‑1 (MQ=255)
                                                                                                                                           ccTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCtt    <  1:415579/149‑1 (MQ=255)
                                                                                                                                            cTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCTTg   >  2:487952/1‑149 (MQ=255)
                                                                                                                                             ttCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCAACAGTTACCCAGCGTGATGGAGCCGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGCCGGAACGAATATCCGCAGCACGGATGCTTgg  >  1:380745/1‑149 (MQ=255)
                                                                                                                                                    aGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAAc                           <  2:31786/117‑1 (MQ=255)
                                                                                                                                                    aGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAAc                           >  1:31786/1‑117 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
TTCCAGTCAGTGTTGGGATCAATGGAAAAAGTGTGGTGTTCTTCGAATCCTTGTGCCGCAGCTACCTGACGAACA‑GCGTCCTGGCTTTCCTGCAATAACAGAGGATCGTTACCAAGTAAAAGATACGCCGCGCGCAGCCCTTCATTGAGCTGCGCGCGGAGTTGTTCCGGGTACAACCGAATCATCAGTTACCCAGCGTGGTGGAGACGCGTGCAGGCGTTGCCGGAGTATCCGTTGTGGTCGACGTCTGTTCTTCGTCGGAACGAATATCCGCAGCACGGATGCTTGG  >  NC_000913_3_hsa_tpiA/671421‑671709

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.