New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2416377 =12 (0.440)13 (0.550) 11/84 0.6 53.6% coding (1631/2145 nt) pta phosphate acetyltransferase
?NC_000913 = 2416373 12 (0.510)coding (1627/2145 nt) pta phosphate acetyltransferase

CAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCG  >  NC_000913/2416330‑2416423
                                               |                                              
gatcaaaactgcaccgggcagctccctggtatcttccgtgttcttcatg                                               <  1:1023910‑M1/2‑1 (MQ=255)
gatcaaaactgcaccgggcagctccctggtatcttccgtgttcttcatg                                               <  1:371269‑M1/2‑1 (MQ=255)
 atcaaaactgcaccgggcagctccctggtatcttccgtgttcttcatgc                                              >  1:946869‑M1/47‑49 (MQ=255)
 atcaaaactgcaccgggcagctccctggtatcttccgtgttcttcatgc                                              >  1:1907350‑M1/47‑49 (MQ=255)
 atcaaaactgcaccgggcagctccctggtatcttccgtgttcttcat                                                >  1:1153667‑M1/47‑47 (MQ=255)
  aaaCTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgctgct                                             <  1:1098341/49‑1 (MQ=255)
  aaaCTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCAtgctgct                                             >  1:179499/1‑49 (MQ=255)
       actgcaccgggcagctccctggtatcttccgtgttcttcaTGCTGCCgg                                        >  1:1098944‑M1/41‑49 (MQ=255)
           caccgggcagctccctggtatcttccgtgttcttcaTGCTGCCGGaa                                      <  1:2079671‑M1/11‑1 (MQ=255)
                ggcagctccctggtatcttccgtgttcttcaTGCTGCCGGAACAGgttt                               >  1:1023182‑M1/32‑49 (MQ=255)
                  cTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgtttac                             <  1:445390/49‑1 (MQ=255)
                   agctccctggtatcttccgtgttcttcaTGCTGCCGGAACAGgtttacg                            <  1:777622‑M1/21‑1 (MQ=255)
                    cccTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgtttacgt                           >  1:1092707/1‑49 (MQ=255)
                        cctggtatcttccgtgttcttcaTGCTGCCGGAACAGgtttacgttt                         >  1:1120719‑M1/24‑47 (MQ=255)
                        cctggtatcttccgtgttcttcaTGCTGCCGGAACAGgtttacgttt                         >  1:417248‑M1/24‑47 (MQ=255)
                         gTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgtttacgtttac                       <  1:1092578/48‑1 (MQ=255)
                         gTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgtttacgtttac                       <  1:1236691/48‑1 (MQ=255)
                         gTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGgtttacgttta                        <  1:190145/47‑1 (MQ=255)
                          tATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACgg                     >  1:1403267/1‑49 (MQ=255)
                             tatcttccgtgttcttcaTGCTGCCGGAACAGGTTTACGTTTACGGTGa                  <  1:1256243‑M1/31‑1 (MQ=255)
                                   tgtTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGc            <  1:139013/49‑1 (MQ=255)
                                   tgtTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGc            >  1:600508/1‑49 (MQ=255)
                                   ccgtgttcttcaTGCTGCCGGAACAGGTTTACGTTTACGGTGACtgt              >  1:1481099‑M1/13‑47 (MQ=255)
                                   ccgtgttcttcaTGCTGCCGGAACAGGTTTACGTTTACGGTGACtgt              >  1:1269868‑M1/13‑47 (MQ=255)
                                    cgtgttcttcaTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGc            <  1:501493‑M1/37‑1 (MQ=255)
                                      tgttcttcaTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGAt         >  1:1353221‑M1/10‑49 (MQ=255)
                                          cATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATc        <  1:217365/46‑1 (MQ=255)
                                          cATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATc        <  1:498699/46‑1 (MQ=255)
                                           ttcaTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAAcc    <  1:111077‑M1/45‑1 (MQ=255)
                                              aTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCg  <  1:109822‑M1/47‑1 (MQ=255)
                                               |                                              
CAAAACTGCACCGGGCAGCTCCCTGGTATCTTCCGTGTTCTTCATGCTGCTGCCGGAACAGGTTTACGTTTACGGTGACTGTGCGATCAACCCG  >  NC_000913/2416330‑2416423

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