New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA = 315942118 (0.730)4 (0.170) 4/272 2.8 17.7% coding (1726/2220 nt) ygiQ Radical SAM superfamily protein
?NC_000913_3_hsa_tpiA = 3159442 20 (0.850)coding (1705/2220 nt) ygiQ Radical SAM superfamily protein
Rejected: Frequency below/above cutoff threshold.

GGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTT  >  NC_000913_3_hsa_tpiA/3159298‑3159586
                                                                                                                                                |                                                                                                                                                
ggTTTTTCCCGGTGTAATACATGGTGGTTGAGTTGGCCAGCGGCGACGGATAGAAGTGCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAgatg                                                                                                                                              >  1:376411‑M2/1‑145 (MQ=255)
   ttttCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAgatgaga                                                                                                                                           <  1:143400‑M2/149‑8 (MQ=255)
        ccGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAggg                                                                                                                                      <  1:242122/149‑1 (MQ=255)
        ccGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAggg                                                                                                                                      <  2:46454/149‑1 (MQ=255)
                          gctGAGGTCGGCAGCGGCGACGTATCGAAGTGCTGCCCCTGGGCGAGGCGGAAGCGGTGCTTTTTCAGCCATAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCttt                                                                                                                    <  2:189360/147‑1 (MQ=255)
                                  cGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAAATACTGCTCTTTAcctgcct                                                                                                            >  1:246116/1‑149 (MQ=255)
                                        cggcgACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTg                                                                                                      <  1:411974/149‑1 (MQ=255)
                                               ggATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGagtaagt                                                                                               >  1:368379/1‑149 (MQ=255)
                                                                cTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCAAACAGCTCTTTaa                                                                              >  2:81687/1‑149 (MQ=255)
                                                                 tGGTCGAGGCGGCAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTaaa                                                                             <  1:285981/149‑1 (MQ=255)
                                                                      gAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGt                                                                        >  1:289035/1‑149 (MQ=255)
                                                                           ggAAGCTGTGCTTTTACAGCCACAGGGCCAGAGCCACCACATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGTGATCAGATACTGTTCTTTACTTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAg                                                                   <  2:193157/149‑1 (MQ=255)
                                                                                  gTGCTTTTTCAGCCACAGCGCCAGATTCACCTTATCTTCATCACGCGTACCGGGGTGCGCGGAGCTGAAATCGGGGATCAGATACTGCTTTTTACCGGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCa                                                            >  1:16954/1‑149 (MQ=255)
                                                                                   tGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCAt                                                           >  2:207112/1‑149 (MQ=255)
                                                                                     cTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAAAAGCTCTTTAAAGCGGTCATGGCTGCCCATGc                                                         >  1:154995/1‑149 (MQ=255)
                                                                                            aGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCtt                                                  >  1:231676/1‑149 (MQ=255)
                                                                                                     tccAGACTCACCATATCTTCATCACGCGTACCGGGGTGCGCTGAGATGAAATACGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTTCCCATGCCCGGCTTCATCATCtt                                         <  1:398701/148‑1 (MQ=255)
                                                                                                                 atatCTTCATCACGCGTACCGGGGTGCGCGGAgatgaagatatggtgaatctggcgctgtggctgaaaaagcaccgcttccgcctcgaccaggtgcagaacttctatccgtcgccgctggcgaact                                                    <  2:148038‑M2/126‑95 (MQ=255)
                                                                                                                 atatCTTCATCACGCGTACCGGGGTGCGCGGAgatgaagatatggtgaatctggcgctgtggctgaaaaagcaccgcttccgcctcgaccaggtgcagaacttctatccgtcgccgctggcgaact                                                    >  1:148038‑M2/1‑32 (MQ=255)
                                                                                                                      ttCATCGCGCGTACCGGGGTGCGCGGAGATGAAATAGGGTATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTCAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGTCCGGCTTCATCATCTTCGATAACGGCCcttctt                        <  1:61072/149‑1 (MQ=255)
                                                                                                                          tcaCGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCCTCTTCGATAACCGCCCTTCATCGGt                    >  1:97190/1‑149 (MQ=255)
                                                                                                                            aCGCGTACCGGGGTGCGCGTatgaagatatggtgaatctggcgctgtggctgaaaaagcaccgcttccgcctcgaccaggtgcagaacttctatccgtcgccgctggcgaac                                                       <  2:294117‑M2/113‑93 (MQ=255)
                                                                                                                            aCGCGTACCGGGGTGCGCGTatgaagatatggtgaatctggcgctgtggctgaaaaagcaccgcttccgcctcgaccaggtgcagaacttctatccgtcgccgctggcgaac                                                       >  1:294117‑M2/1‑21 (MQ=255)
                                                                                                                                    cGGGGTGCGGGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCggg          <  2:16954/149‑1 (MQ=255)
                                                                                                                                            gcgGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCtt  <  2:143400/149‑1 (MQ=255)
                                                                                                                                            gcgGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCtt  >  1:306597/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GGTTTTTCCCGGTGTAATACATGGTGGTTGAGTTCGCCAGCGGCGACGGATAGAAGTTCTGCACCTGGTCGAGGCGGAAGCGGTGCTTTTTCAGCCACAGCGCCAGATTCACCATATCTTCATCACGCGTACCGGGGTGCGCGGAGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTT  >  NC_000913_3_hsa_tpiA/3159298‑3159586

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.