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.

AACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTT  >  NC_000913_3_hsa_tpiA/1665080‑1665371
                                                                                                                                                 |                                                                                                                                                  
aaCTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGAACGCACTGCTTTTGCGGCTgg                                                                                                                                                 >  2:256980/1‑149 (MQ=255)
     ttAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctct                                                                                                                                            >  1:267056/1‑149 (MQ=255)
                                                    tatTGCCGGAGGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCAt                                                                                             <  2:392421/149‑1 (MQ=255)
                                                    tatTGCCGGAGGGGGACGATTACGCGCGACCAGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACGGGTGTTGCGGCTGGCGTCTTTTGCCGGAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCAt                                                                                             <  2:143470/149‑1 (MQ=255)
                                                           ggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCaaaagag                                                                                                                                             <  2:250003‑M1/94‑8 (MQ=255)
                                                           ggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCaaaagag                                                                                                                                             >  1:250003‑M1/1‑87 (MQ=255)
                                                            gATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCCCATGCatcatc                                                                                     >  1:22637/1‑149 (MQ=255)
                                                                ggtaCGACTACGTGCAATCCGGTAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTaca                                                                                 <  2:171697/146‑1 (MQ=255)
                                                                             gcAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGAt                                                                    >  1:26742/1‑149 (MQ=255)
                                                                                     gCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCaaaagaggccagccgcaaaagcagtgcgtgctttaagagat                                                                                                           >  1:167340‑M1/1‑61 (MQ=255)
                                                                                     gCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCaaaagaggccagccgcaaaagcagtgcgtgctttaagagat                                                                                                           <  2:167340‑M1/102‑42 (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                                                                            <  2:167340‑M2/102‑62 (MQ=255)
                                                                                                                    aTCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTgccgcaaaagcagtgcgtgctttaagagataaaaaaacgtgacactgtaacctattattgc                                                                            >  1:167340‑M2/1‑41 (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)
                                                                                                                                                 |                                                                                                                                                  
AACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTTATCGCAACAGTCGTCGGTATTCTCGCGGCCTTTGCCGTT  >  NC_000913_3_hsa_tpiA/1665080‑1665371

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.