New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 704375 =18 (0.810)3 (0.150) 3/116 1.8 16.3% coding (432/1947 nt) nagE N‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components
?NC_000913 704405 = 15 (0.770)coding (462/1947 nt) nagE N‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components
Rejected: Frequency below/above cutoff threshold.

GTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATC  >  NC_000913/704340‑704471
                                                                 |                                                                  
gtgCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGt                                                                  >  2:179868/1‑68 (MQ=255)
gtgCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGt                                                                  <  2:204594/68‑1 (MQ=255)
 tgCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTa                                                                 >  1:669747/1‑68 (MQ=255)
      aTTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCa                                                            <  2:537660/68‑1 (MQ=255)
       ttGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCAc                                                           <  2:441531/68‑1 (MQ=255)
        tGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACg                                                          <  2:788318/68‑1 (MQ=255)
          gcctgcatggatagcgtgctgtaccggcggccagacgtaaccaaaaatggccgccGGTACAGCACGCt                                                        <  1:232697‑M2/13‑1 (MQ=255)
           cctgcatggatagcgtgctgtaccggcggccagacgtaaccaaaaatggccgccGGTACAGCACGCTa                                                       <  1:502090‑M2/14‑1 (MQ=255)
            aCCGGTTTCTTCTGCCTGGTGCTGGCGGCCATTTTAGGTTACGTCTGGCCGCCGGTACAGCACGCTAt                                                      <  2:221048/68‑1 (MQ=255)
            aCCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTAt                                                      <  2:99590/68‑1 (MQ=255)
            aCCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTAt                                                      <  1:750866/68‑1 (MQ=255)
                     tagcgtgctgtaccGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAgg                                             <  1:502090‑M1/54‑1 (MQ=255)
                      agcgtgctgtaccGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAggc                                            <  1:232697‑M1/55‑1 (MQ=255)
                          ccTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGgcggcg                                        >  1:625212/1‑68 (MQ=255)
                           cTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGa                                       <  1:671475/68‑1 (MQ=255)
                                gCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATc                                                     >  1:219877/1‑49 (MQ=255)
                                gCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATc                                                     <  2:219876/49‑1 (MQ=255)
                                gCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTggg                                  >  1:21162/1‑67 (MQ=255)
                                     cggccagacgtaaccaaaaatggccgccGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGtt                             >  1:31006‑M2/29‑68 (MQ=255)
                                     CGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGtt                             <  1:792683‑M1/68‑1 (MQ=255)
                                                ggTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGgcgc                  <  2:228878/68‑1 (MQ=255)
                                                          ggCCGCCGGTACAGCTCGCTTTCCATTCAGGCGGCGAGTGTATCGTTTCTGCGGTCGCGCTTGGTTcc        <  1:551996/68‑1 (MQ=255)
                                                            ccgccGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCgg      >  2:786591‑M2/6‑68 (MQ=255)
                                                               ccGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTAt   <  2:792685‑M2/66‑1 (MQ=255)
                                                                cGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATc  <  1:94478‑M2/67‑1 (MQ=255)
                                                                 |                                                                  
GTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATC  >  NC_000913/704340‑704471

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 37 ≤ 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.