New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA = 455746027 (1.140)4 (0.180) 4/262 2.6 14.3% coding (812/912 nt) hypT hypochlorite‑responsive transcription factor
?NC_000913_3_hsa_tpiA = 4557471 23 (1.050)coding (801/912 nt) hypT hypochlorite‑responsive transcription factor
Rejected: Frequency below/above cutoff threshold.

GCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGAAGAGACAAAAAAGGTGCTG  >  NC_000913_3_hsa_tpiA/4557339‑4557606
                                                                                                                                    |                                                                                                                                       
gcAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGGCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACaa                                                                                                                         >  1:291248/1‑149 (MQ=255)
              gagaAGGTCAGCTAGGCACAATCTCCGGCTCGCGCAGTTCACGGCAGAAGCGTTCGGTAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTg                                                                                                           <  2:11414/149‑1 (MQ=255)
                           aaGCACAATCTGCCGCTGGCGCAGTTCACGCCAGAATCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCAGAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCttgtt                                                                                              <  2:224512/149‑1 (MQ=255)
                              cacaATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTggc                                                                                           >  2:306184/1‑148 (MQ=255)
                              cacaATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGAATCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGa                                                                                           >  2:69039/1‑149 (MQ=255)
                                 aaTCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGAGTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTGCTGGTGGATgg                                                                                        >  1:159614/1‑149 (MQ=255)
                                         gcgcgcgCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACTAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTcc                                                                                <  2:346957/146‑1 (MQ=255)
                                                  gTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGGTCACCAGTTCGTCCCGTTTCAGTACAACGAGCTTCCCGCTGCGAGTTTCTTGTTGGATGGCGTACTCCGGCAGCCAg                                                                       <  2:121358/149‑1 (MQ=255)
                                                  gTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAg                                                                       <  2:121350/149‑1 (MQ=255)
                                                             aaGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCcaca                                                            <  1:83824/149‑1 (MQ=255)
                                                                     ggCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGACTTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGa                                                    <  1:123194/149‑1 (MQ=255)
                                                                      gCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAg                                                   <  1:81196/149‑1 (MQ=255)
                                                                        aaCGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGCggg                                                 >  2:294968/1‑149 (MQ=255)
                                                                        aaCGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAggg                                                 >  2:219923/1‑149 (MQ=255)
                                                                          cGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCa                                               >  2:213596/1‑149 (MQ=255)
                                                                                  aTGCGGGTACTCCTCCGGTATGCATCGGCCTGAACCGGGATCACCAGTTCGTCCCGATTCAGTAACACGAGCTTCCCGCTGCGAATGTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCtt                                       <  1:412697/149‑1 (MQ=255)
                                                                                       ggTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTaaaa                                  <  2:291248/149‑1 (MQ=255)
                                                                                                               cTGAATCGGGATCACCAGTTCGtcccgattcaggcctatgcataccggatgaatacccgcatgaatcccgttgccgaac                                                                                <  2:193075‑M2/79‑58 (MQ=255)
                                                                                                               cTGAATCGGGATCACCAGTTCGtcccgattcaggcctatgcataccggatgaatacccgcatgaatcccgttgccgaac                                                                                >  1:193075‑M2/1‑22 (MQ=255)
                                                                                                               cTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGTTAGCTTACTCCGGCAGCCAGGCAATCCAACAGCCGTCGAGAGCAACATGCTTGCAAAGCTAGCTCATCGTAGAGACaaaaa          >  1:183726/1‑149 (MQ=255)
                                                                                                               cTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGAAGAGAcaacaa          >  2:67146/1‑149 (MQ=255)
                                                                                                               cTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGAAGAGACaaaaa          >  1:190276/1‑149 (MQ=255)
                                                                                                                tGAATCGGGATCACCAGTTCGtcccgattcaggcctatgcataccggatgaatacccgcatgaatcccgttgccgaacgcttctggcgtgaactgcgc                                                            <  2:91507‑M2/98‑78 (MQ=255)
                                                                                                                tGAATCGGGATCACCAGTTCGtcccgattcaggcctatgcataccggatgaatacccgcatgaatcccgttgccgaacgcttctggcgtgaactgcgc                                                            >  1:91507‑M2/1‑21 (MQ=255)
                                                                                                                       ggATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCATCCTGCTTTAAAAGCTCGCTCATCGAAGAGACAAAAAAGGTGCTg  >  2:312381/1‑149 (MQ=255)
                                                                                                                              cAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGt                                                       <  2:352734/89‑1 (MQ=255)
                                                                                                                              cAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGt                                                       >  1:352734/1‑89 (MQ=255)
                                                                                                                                    |                                                                                                                                       
GCAGTGCGACGCTGGAGAATGTCAGCTAAGCACAATCTCCAGCTCGCGCAGTTCACGCCAGAAGCGTTCGGCAACGGGATTCATGCGGGTATTCATCCGGTATGCATAGGCCTGAATCGGGATCACCAGTTCGTCCCGATTCAGTACAACGAGCTTCCCGCTGCGAATTTCTTGTTGGATGGCGTACTCCGGCAGCCAGGCAATCCCACAGCCGTCGAGGGCAACCTGCTTTAAAAGCTCGCTCATCGAAGAGACAAAAAAGGTGCTG  >  NC_000913_3_hsa_tpiA/4557339‑4557606

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.