New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 250869436 (0.930)4 (0.120) 4/110 2.9 11.7% coding (733/966 nt) glk glucokinase
?NC_000913 = 2508732 30 (0.910)coding (695/966 nt) glk glucokinase
Rejected: Frequency below/above cutoff threshold.

GCGGCACGATACCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGC  >  NC_000913/2508628‑2508761
                                                                  |                                                                   
gcgGCACGATACCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGCTTGCCGCCa                                                                    >  1:35479/1‑68 (MQ=255)
       gATACCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGc                                                             <  2:1068042/68‑1 (MQ=255)
       gATACCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGc                                                             <  2:262241/68‑1 (MQ=255)
       gATACCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGc                                                             <  2:262242/68‑1 (MQ=255)
       gATACCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGc                                                             <  1:846118/68‑1 (MQ=255)
       gATACCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGc                                                             <  1:563044/68‑1 (MQ=255)
       gATACCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGc                                                             <  2:62068/68‑1 (MQ=255)
       gATACCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGc                                                             <  1:136756/68‑1 (MQ=255)
       gATACCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGc                                                             <  1:1137692/68‑1 (MQ=255)
       gATACCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGc                                                             <  1:1168659/68‑1 (MQ=255)
         tACCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGccc                                                           <  1:216061/68‑1 (MQ=255)
          aCCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCa                                                          >  1:532669/1‑68 (MQ=255)
            cgcccgcAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATa                                                        <  1:805755/68‑1 (MQ=255)
              cccgcAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAAt                                                      >  1:31738/1‑68 (MQ=255)
              cccgcAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAAt                                                      <  1:1097792/68‑1 (MQ=255)
                 gcAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGAc                                                   <  2:279546/68‑1 (MQ=255)
                       aaCACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCaaaa                                             <  2:836211/68‑1 (MQ=255)
                         cacGCCGCCAAATGTCCCGAGATTGCGCGCCAGATTGCCGCCAAAACGGCCCATAATGCCGCAAACCa                                           >  2:1371548/1‑68 (MQ=255)
                         cacGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGAc                                                   <  2:1167051/60‑1 (MQ=255)
                         cacGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGAc                                                   >  1:1167055/1‑60 (MQ=255)
                         cacGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACa                                           >  1:410641/1‑68 (MQ=255)
                         cacGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACa                                           >  2:55974/1‑68 (MQ=255)
                         cacGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACa                                           >  1:217434/1‑68 (MQ=255)
                         cacGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACa                                           >  1:1379584/1‑68 (MQ=255)
                         cacGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACa                                           >  2:824071/1‑68 (MQ=255)
                         cacGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGCCCCATAATGACGCAAAACa                                           >  1:831058/1‑68 (MQ=255)
                           cgccgcCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGc                                         <  2:431692/68‑1 (MQ=255)
                              cgccAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCgcgcattgtcgctgttttgcgtcattatggg                                      <  1:1405605‑M1/68‑32 (MQ=255)
                              cgccAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGAc                                      <  1:1394183/68‑1 (MQ=255)
                              cgccAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGAc                                      <  1:740939/68‑1 (MQ=255)
                                      gTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGCCGCAAAACAGCGACAATGCGCg                              <  1:1181453‑M2/68‑2 (MQ=255)
                                        cccGAGATTGAGCGCCAGATTGCCGCCAAAACGTCCCATAATGACGCAAAACAGCGACAATGCGCggc                            >  1:1356565‑M2/1‑65 (MQ=255)
                                        cccGAGATAGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGCCGCAAAACAGCGACAATGCGCggc                            <  1:444433‑M2/68‑4 (MQ=255)
                                        ccagagaTTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGCCAATGCgcggc                            >  1:442541/4‑68 (MQ=255)
                                         ccGAGATTGAGCGCCAGATTGCCGCCgcgcattgtcgctgttttgcgtcattatgggccg                                   >  2:1383064‑M1/1‑26 (MQ=255)
                                         ccGAGATTGAGCGCCAGATTGCCGCCgcgcattgtcgctgttttgcgtcattatgggccg                                   <  1:1383069‑M1/60‑35 (MQ=255)
                                            agaTTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCggcggca                        <  1:1087192‑M2/68‑8 (MQ=255)
                                             gaTTGCGCGCCAGATTGCCGCCAAAGCGGCCCATAATGACGCAAAACAGCGACAATGGGCggcggcaa                       <  2:903803‑M2/68‑9 (MQ=255)
                                             gaTTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAGAACAGCGACAATGCGCggcggcgg                       <  2:908710‑M2/68‑9 (MQ=255)
                                             gaTTGAGCGCCAGATAGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCggcggcaa                       >  1:670572‑M2/1‑60 (MQ=255)
                                              aTTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCggcggcagt                      <  2:106394‑M2/68‑10 (MQ=255)
                                                  aGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGt                  <  1:968035/68‑1 (MQ=255)
                                                  aGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGt                  <  1:949231/68‑1 (MQ=255)
                                                   gcgcCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTg                 >  1:1235877/1‑68 (MQ=255)
                                                   gcgcCAGATTGCCGCAgcgcattgtcgctgttttgagtcattatgggccgttttggcggcaatctgga                 >  1:998098‑M1/1‑16 (MQ=255)
                                                     tccaGATTGCCGCCAAAACGGCCCATAATGACTCAAAACAGCGACAATGCGCtgcggcaatctggcgc               >  1:998098‑M2/2‑52 (MQ=255)
                                                        aGATTGCCGCCAAAACGGCCCATAAAGACGCAAAACAGCGACAATGCGCGGCTGCAATCGGTGCAGCt            <  2:1344130/68‑1 (MQ=255)
                                                          actgcCGCCgcgcattgtcgctgttttgcgtcattatgggccgttttggcggcaatctggcgctcaat          <  2:106394‑M1/66‑60 (MQ=255)
                                                          aTTGCCGCCAAAACGCCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGt          <  2:831053/68‑1 (MQ=255)
                                                           ttGCCGCCgcgcattgtcgctgttttgcgtcattatgggccgttttggcggctatctggcgctcaatc         >  1:670572‑M1/1‑8 (MQ=255)
                                                           ttGCCGCCgcccattgtcgctgttttgcgtcattatgggccgctttggcggcaatctggcgcgcaatc         <  2:903803‑M1/68‑61 (MQ=255)
                                                           ccgccgccgcgcattgtcgctgttctgcgtcattatgggccgttttggcggcaatctggcgctcaatc         <  2:908710‑M1/66‑61 (MQ=255)
                                                            tGCCGCCgcgcattgtcgctgttttgcgtcattatgggccgttttggcggcaatctggcgctcaatct        <  1:1087192‑M1/68‑62 (MQ=255)
                                                                gccgcgcattgtcgctgttttgcgtcattatgggacgttttggcggcaatctggcgctcaatctcggg    >  1:1356565‑M1/1‑3 (MQ=255)
                                                                gccgcgcattgtcgctgttttgcggcattatgggccgttttggcggcaatctggcgctctatctcggg    <  1:444433‑M1/68‑66 (MQ=255)
                                                                 ccAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTcagccag   >  1:640298/1‑68 (MQ=255)
                                                                  cgcgcattgtcgctgttttgcggcattatgggccgttttggcggcaatctggcgctcaatctcgggac  <  1:1181453‑M1/68‑68 (MQ=255)
                                                                  |                                                                   
GCGGCACGATACCGCCCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGC  >  NC_000913/2508628‑2508761

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.