New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA = 166522519 (0.800)4 (0.180) 4/268 2.6 17.5% intergenic (+105/‑90) dmsD/clcB twin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter
?NC_000913_3_hsa_tpiA = 1665236 20 (0.900)intergenic (+116/‑79) dmsD/clcB twin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter
Rejected: Frequency below/above cutoff threshold.

TATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTC  >  NC_000913_3_hsa_tpiA/1665132‑1665384
                                                                                                        |                                                                                                                                                    
tatTGCCGGAGGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCAt                                                                                                          <  2:392421/149‑1 (MQ=255)
tatTGCCGGAGGGGGACGATTACGCGCGACCAGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACGGGTGTTGCGGCTGGCGTCTTTTGCCGGAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCAt                                                                                                          <  2:143470/149‑1 (MQ=255)
        gATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCCCATGCatcatc                                                                                                  >  1:22637/1‑149 (MQ=255)
            ggtaCGACTACGTGCAATCCGGTAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTaca                                                                                              <  2:171697/146‑1 (MQ=255)
                         gcAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGAt                                                                                 >  1:26742/1‑149 (MQ=255)
                                          ggTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTATTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATGTACACATCTATCCGGATCTGCGCACTATGTTCCa                                                                <  2:333913/149‑1 (MQ=255)
                                           gTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCAc                                                               <  2:258286/149‑1 (MQ=255)
                                                    tCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCGCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCAACCGTCTGCtt                                                      >  1:321055/1‑149 (MQ=255)
                                                      aCGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCAACACATCTCTCCGGATCTGCGCACTATGTTCCACCGTCTGCTTAt                                                    >  1:267331/1‑149 (MQ=255)
                                                           tttttATCCCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTGGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCAGGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCaa                                               <  2:18514/149‑1 (MQ=255)
                                                                aTCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTgccgcaaaagcagtgcgtgctttaagagataaaaaaacgtgacactgtaacctattattgc                                                                                         >  1:167340‑M2/1‑41 (MQ=255)
                                                                aTCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTgccgcaaaagcagtgcgtgctttaagagataaaaaaacgtgacactgtaacctattattgc                                                                                         <  2:167340‑M2/102‑62 (MQ=255)
                                                                 tctcTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTAt                                                    <  2:116360/138‑1 (MQ=255)
                                                                 tctcTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTAt                                                    >  1:116360/1‑138 (MQ=255)
                                                                    ctTAAATCCCGCACTGCTTTGGCGGCTGGCCTCTTTTGTCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAGtcgtcg                                      <  1:392652/149‑1 (MQ=255)
                                                                     ttAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCCCACTATGTTCCACCGTCTGCTTATGCCAACAGCCGTCgg                                     >  2:392652/1‑149 (MQ=255)
                                                                           cacgcacTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAGTCGTCGGTATtct                               >  1:191273/1‑149 (MQ=255)
                                                                                 cTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGGGCACTATGTTCCACCGTCTGCTTATCGCAACCGTCGTCGGTATTTTCGCGGc                         >  2:116719/1‑149 (MQ=255)
                                                                                           ggCTGGCCTCCTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTttgccgtt               >  2:216470/1‑149 (MQ=255)
                                                                                                gCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACACTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCgg          >  2:288018/1‑149 (MQ=255)
                                                                                                  ctctTTTgccgcaaaagcagtgcgtgctttaagagataaaaaaacgtgacactgtaacctattattgccggattgcgcgtaatcgtcaccatcc                                                               <  2:250003‑M2/94‑88 (MQ=255)
                                                                                                  ctctTTTgccgcaaaagcagtgcgtgctttaagagataaaaaaacgtgacactgtaacctattattgccggattgcgcgtaatcgtcaccatcc                                                               >  1:250003‑M2/1‑7 (MQ=255)
                                                                                                  ctctTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAGTCTTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGt        >  2:232047/1‑149 (MQ=255)
                                                                                                        tGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTc  <  2:26742/149‑1 (MQ=255)
                                                                                                        |                                                                                                                                                    
TATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTTGCCGGGTTTCGTC  >  NC_000913_3_hsa_tpiA/1665132‑1665384

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.