New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 250217336 (0.650)4 (0.080) 4/146 3.5 10.4% coding (733/966 nt) glk glucokinase
?NZ_CP009273 = 2502211 37 (0.760)coding (695/966 nt) glk glucokinase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTT  >  NZ_CP009273/2502122‑2502296
                                                                                         |                                                                                     
ccgcAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATgcgc                                                                                       <  2:967793/90‑1 (MQ=255)
   cAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCgcggc                                                                                    >  1:1046296/1‑90 (MQ=255)
   cAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCgcggc                                                                                    >  2:149078/1‑90 (MQ=255)
   cAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCgcggc                                                                                    >  1:1146428/1‑90 (MQ=255)
        aaCACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCaa                                                                               <  1:230575/90‑1 (MQ=255)
            cgccgcCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCgg                                                                           >  2:255978/1‑90 (MQ=255)
             gccgccAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGt                                                                          >  1:434093/1‑90 (MQ=255)
             gccgccAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGt                                                                          >  1:443413/1‑90 (MQ=255)
                    aaTGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTg                                                                   <  2:90074/90‑1 (MQ=255)
                        tCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTcag                                                               >  1:535524/1‑90 (MQ=255)
                          ccGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTcagcc                                                             >  1:950998/1‑90 (MQ=255)
                          ccGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTcagcc                                                             >  2:1474821/1‑90 (MQ=255)
                                 tGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAgcgcgc                                                      >  2:1311254/1‑90 (MQ=255)
                                  gAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGcgcgcg                                                     >  1:473174/1‑90 (MQ=255)
                                   aGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGt                                                    <  1:200670/90‑1 (MQ=255)
                                   aGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGt                                                    <  1:346265/90‑1 (MQ=255)
                                    gcgcCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGtt                                                   <  2:1146428/90‑1 (MQ=255)
                                    gcgcCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGtt                                                   <  2:1497016/90‑1 (MQ=255)
                                    gcgcCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGtt                                                   <  1:219608/90‑1 (MQ=255)
                                    gcgcCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGtt                                                   <  1:1041683/90‑1 (MQ=255)
                                     cgcCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTc                                                  <  1:361352/90‑1 (MQ=255)
                                     cgcCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTc                                                  <  1:99770/90‑1 (MQ=255)
                                            ttGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAata                                           <  2:10053/90‑1 (MQ=255)
                                                 gccAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCttt                                      <  1:1387571/90‑1 (MQ=255)
                                                    aaaaCGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCa                                                                       <  2:755167/54‑1 (MQ=255)
                                                    aaaaCGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCa                                                                       >  1:755167/1‑54 (MQ=255)
                                                    aaaaCGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTgg                                   >  2:1455841/1‑90 (MQ=255)
                                                          gCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTgag                             <  2:443413/90‑1 (MQ=255)
                                                                    aCGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGc                   >  2:1106241/1‑90 (MQ=255)
                                                                    aCGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGc                   >  1:253368/1‑90 (MQ=255)
                                                                    aCGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGc                   >  1:389320/1‑90 (MQ=255)
                                                                    aCGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGc                   >  1:1227567/1‑90 (MQ=255)
                                                                      gCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAg                 >  1:677061/1‑90 (MQ=255)
                                                                             aGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAat                                            <  2:1553919/56‑1 (MQ=255)
                                                                             aGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAat                                            >  1:1553919/1‑56 (MQ=255)
                                                                             aGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTATGGCTTGAGATTTTCTGGCAGGCGGGTg          >  2:1296045/1‑90 (MQ=255)
                                                                                gACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCa       <  1:701827/90‑1 (MQ=255)
                                                                                 aCAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAg      >  2:1087618/1‑90 (MQ=255)
                                                                                   aaTGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCt    <  2:16685/90‑1 (MQ=255)
                                                                                   aaTGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCt    >  2:1208807/1‑90 (MQ=255)
                                                                                   aaTGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCt    <  2:959063/90‑1 (MQ=255)
                                                                                     tGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCttt  >  2:1408108/1‑90 (MQ=255)
                                                                                         |                                                                                     
CCGCAATAAACACGCCGCCAAATGTCCCGAGATTGAGCGCCAGATTGCCGCCAAAACGGCCCATAATGACGCAAAACAGCGACAATGCGCGGCGGCAATCGGTGCAGCTGTCAGCCAGCGCGCGTTCGGTAATATCTTTTGGCTTGAGATTTTCTGGCAGGCGGTTGTCAGCTTT  >  NZ_CP009273/2502122‑2502296

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