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.

TGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGAC  >  NZ_CP009273/4331843‑4332017
                                                                                               |                                                                                       
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:1387373/1‑90 (MQ=255)
        gACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCaa                                                                                >  1:529208/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                                                            >  2:227729/1‑90 (MQ=255)
                            gTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGa                                                            >  1:85648/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                                           >  1:1307503/1‑90 (MQ=255)
                                             gTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCgtgg                                           >  2:125464/1‑90 (MQ=255)
                                              tAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCgtggt                                          <  2:581231/90‑1 (MQ=255)
                                                aGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAATCCTTAGAGGAt                                                           >  1:1252901/1‑71 (MQ=255)
                                                aGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAATCCTTAGAGGAt                                                           <  2:1252901/71‑1 (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                       >  1:149624/1‑90 (MQ=255)
                                                                 ggggCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGt                       >  2:1039614/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)
                                                                                               |                                                                                       
TGGACATTGACCCCACCCGCCTGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCA‑‑‑‑CCTG‑‑‑CCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGAC  >  NZ_CP009273/4331843‑4332017

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