New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 433193062 (1.380)6 (0.130)
+CCA
3/146 3.2 8.8% coding (203/1356 nt) melA alpha‑galactosidase
?NZ_CP009273 4331931 = 64 (1.350)coding (204/1356 nt) melA alpha‑galactosidase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGAC  >  NZ_CP009273/4331841‑4332017
                                                                                             |                                                                                           
gATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTg                                                                                             >  2:971114/1‑90 (MQ=255)
gATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTg                                                                                             >  1:584300/1‑90 (MQ=255)
  gggccATTGCCCCCCCCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑cc                                                                                        <  2:995538/86‑1 (MQ=255)
   ggACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑Cca                                                                                       >  2:1469380/1‑90 (MQ=255)
   ggACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑Cca                                                                                       <  2:524232/90‑1 (MQ=255)
     aCATTGACCCCACCCGCCTGGAATAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑Ccaca                                                                                     >  1:959939/1‑90 (MQ=255)
     aCATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑Ccaca                                                                                     >  1:782696/1‑90 (MQ=255)
      ctttGCCCCCCCCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCacac                                                                                    <  2:193499/88‑1 (MQ=255)
        ttGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGACAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACAccc                                                                                  >  2:703642/1‑90 (MQ=255)
          gACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCaa                                                                                >  1:529208/1‑90 (MQ=255)
          gACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCaa                                                                                >  1:1387373/1‑90 (MQ=255)
            ccccccccGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCC‑‑‑‑CCTG‑‑‑CCACACCCAACa                                                                              <  1:1403723/90‑1 (MQ=255)
            ccccccccGCCGGGAAGAGTCGCATTTTGTGGTGCGAAAGCTTATGGATTCAGCGGGGGCCAGCGGCAAAATCC‑‑‑‑CCTG‑‑‑CCACACCCAACa                                                                              <  1:1469380/90‑1 (MQ=255)
             cccACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACag                                                                             >  1:1412990/1‑90 (MQ=255)
              ccACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACaga                                                                            >  2:583446/1‑90 (MQ=255)
               cACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACagaa                                                                           >  1:626306/1‑90 (MQ=255)
               cACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACagaa                                                                           >  2:1077552/1‑90 (MQ=255)
                          aaGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTag                                                                <  2:1387373/90‑1 (MQ=255)
                           agagTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCC‑‑‑‑CCTG‑‑‑CCACACCCAACagaaaga                                                                        <  1:561321/81‑1 (MQ=255)
                           agagTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCC‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTaga                                                               <  2:856035/90‑1 (MQ=255)
                           agagTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACagaaaga                                                                        >  2:561321/1‑81 (MQ=255)
                            gagTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTagag                                                              <  2:372974/90‑1 (MQ=255)
                              gTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCC‑‑‑‑CCTG‑‑‑CCACACCCAACagaaagaa                                                                       <  1:955699/79‑1 (MQ=255)
                              gTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTGCCACCACACCCAACAGAAAGAAGCCTTaga                                                               >  1:822698/1‑90 (MQ=255)
                              gTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACagaaagaa                                                                       >  2:955699/1‑79 (MQ=255)
                              gTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAaagaag                                                                      >  2:854974/1‑80 (MQ=255)
                              gTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAaagaag                                                                      <  1:854974/80‑1 (MQ=255)
                              gTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCt                                                                   <  2:1162766/83‑1 (MQ=255)
                              gTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCt                                                                   >  1:1162766/1‑83 (MQ=255)
                              gTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGa                                                            >  1:85648/1‑90 (MQ=255)
                              gTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGa                                                            >  2:227729/1‑90 (MQ=255)
                               tCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGAt                                                           >  1:1294463/1‑90 (MQ=255)
                                        gtggtgCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGAtttt                                                  >  2:31959/1‑90 (MQ=255)
                                           gtgCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTGCCACCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGAtttt                                                  <  2:822698/90‑1 (MQ=255)
                                           gtgCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTc                                               <  2:1294463/90‑1 (MQ=255)
                                              cGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCC‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCgtg                                            <  1:971114/90‑1 (MQ=255)
                                               gaaaGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCgtgg                                           >  2:323393/3‑90 (MQ=255)
                                               gTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCgtgg                                           >  2:125464/1‑90 (MQ=255)
                                               gTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCgtgg                                           >  1:1307503/1‑90 (MQ=255)
                                                tAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCgtggt                                          <  2:581231/90‑1 (MQ=255)
                                                  aGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAATCCTTAGAGGAt                                                           <  2:1252901/71‑1 (MQ=255)
                                                  aGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAATCCTTAGAGGAt                                                           >  1:1252901/1‑71 (MQ=255)
                                                  aGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGtggtgg                                        >  1:563198/1‑90 (MQ=255)
                                                        tGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTc                                  <  2:1307503/90‑1 (MQ=255)
                                                        tGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTc                                  <  1:323393/90‑1 (MQ=255)
                                                         ggATTCAGCAGGGGCCGGCGGAAAAATCC‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCa                                 <  1:148878/90‑1 (MQ=255)
                                                             tCAGCAGGGGCCAGCGGCAAAATCACCTGCCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGCCGTGGGGGCATTTCa                                 >  1:622294/1‑90 (MQ=255)
                                                              cagcagGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTg                            >  2:675018/1‑90 (MQ=255)
                                                                gcgggggCCAGCGGCAAAATCC‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGc                          <  1:141854/87‑1 (MQ=255)
                                                                gcagGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGc                          <  2:959939/90‑1 (MQ=255)
                                                                gcagGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGc                          <  2:858113/90‑1 (MQ=255)
                                                                gcagGGGCAAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCAAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGc                          <  2:107961/90‑1 (MQ=255)
                                                                   ggggCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGt                       >  2:1039614/1‑90 (MQ=255)
                                                                   ggggCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGt                       >  1:149624/1‑90 (MQ=255)
                                                                     ggCCAGCGGCAAAATCC‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTc                                  <  2:476914/77‑1 (MQ=255)
                                                                     ggCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTc                                  >  1:476914/1‑77 (MQ=255)
                                                                     ggCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTa                     >  2:1396073/1‑90 (MQ=255)
                                                                     ggCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTa                     >  1:440068/1‑90 (MQ=255)
                                                                       ccAGGGGAAAAATCC‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATg                   <  1:227729/90‑1 (MQ=255)
                                                                        cAGCGGCAAAATCC‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGc                          <  1:1177127/82‑1 (MQ=255)
                                                                        cAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGc                          >  2:1177127/1‑82 (MQ=255)
                                                                         aGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGaa                <  1:1039614/90‑1 (MQ=255)
                                                                         aGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGAATTTCAGATTGGCGGTTATGaa                <  2:239160/90‑1 (MQ=255)
                                                                              cAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTg           >  2:828349/1‑90 (MQ=255)
                                                                              cAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACttt           >  1:1407221/1‑89 (MQ=255)
                                                                               aaaaTCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGGTTTTAACCTTGc          >  1:381186/1‑90 (MQ=255)
                                                                                 aaTCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCAc        <  2:440068/90‑1 (MQ=255)
                                                                                 aaTCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCAc        <  1:675018/90‑1 (MQ=255)
                                                                                  aTCC‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACg       <  2:584300/90‑1 (MQ=255)
                                                                                    cA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTGTGAAACCTTGccggg      >  2:48376/1‑86 (MQ=255)
                                                                                         cccTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTg    <  1:147748/89‑1 (MQ=255)
                                                                                           cTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGGGAc  >  2:450328/1‑90 (MQ=255)
                                                                                             |                                                                                           
GATGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGAC  >  NZ_CP009273/4331841‑4332017

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