New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 335560612 (0.650)5 (0.310) 4/116 1.3 24.7% coding (882/4461 nt) gltB glutamate synthase, large subunit
?NC_000913 = 3355642 20 (1.230)coding (918/4461 nt) gltB glutamate synthase, large subunit

GGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCAT  >  NC_000913/3355585‑3355707
                                                         |                                                                 
ggAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGc                                                         >  1:250776/1‑68 (MQ=255)
   aCTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGa                                                      >  2:121193/1‑68 (MQ=255)
        tgctgGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGaacaac                                                 >  1:318346/1‑68 (MQ=255)
         gctgGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAAcc                                                <  1:122065/68‑1 (MQ=255)
             gCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGa                                            <  2:42897/68‑1 (MQ=255)
             gCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGa                                            <  1:497357/68‑1 (MQ=255)
             gCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACAACCCGCCTGGCAGAACAACCCGGa                                            >  1:290179/1‑68 (MQ=255)
             gCAGGCGGGATGGATATCATCCGCGCCATGCGGCTATTAGTACCACCCGCCTGGCAGAACAACCCGGa                                            <  1:452649/68‑1 (MQ=255)
              cAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTGCCACCCGCCTGGCAGAACAACCCGGat                                           >  2:500411/1‑68 (MQ=255)
               aGGCGGGATGGACATCATCCGCGCCATGCGTCTATTAGTACCAcccgcctgccagcagcagttccagc                                          >  1:441086‑M2/1‑43 (MQ=255)
                    ggatggatATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGAc                                     <  2:1493/68‑1 (MQ=255)
                      atggatATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGAccc                                   <  2:281754/68‑1 (MQ=255)
                           tatCATCCGCGCCATGCGTCTATTAGTACCAcccgcctgccagcagcagttccagcatgttatccatc                              <  2:412264‑M2/68‑38 (MQ=255)
                           tatCATCCGCGCCATGCGTCTATTAGTACCAcccgcctgccagcagcagttccagcatgttatccatc                              <  1:25165‑M2/68‑38 (MQ=255)
                               atcCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTgag                          >  1:626009/1‑66 (MQ=255)
                                tcCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAGCCCGGATATGGACACGGAACTGCGt                         >  2:409487/1‑68 (MQ=255)
                                tcCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGt                         <  1:273627/68‑1 (MQ=255)
                                 ccGCGCCATGCGTCTATTAGTACCAcccgcctgccagcagcagttccagcatgttatccatcgaactg                        >  1:462890‑M2/1‑25 (MQ=255)
                                     gcCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCctt                    <  2:250775/68‑1 (MQ=255)
                                          gCGTCTATTAGTACCAcccgcctggcagcagcagttccagcatgttagccatcgaactggagtcagag               >  2:44111‑M2/1‑16 (MQ=255)
                                                   aGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACt      <  1:659067/68‑1 (MQ=255)
                                                   aGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACt      <  1:659787/68‑1 (MQ=255)
                                                     tACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTcc    >  1:538776/1‑68 (MQ=255)
                                                       ccaccCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCAt  >  1:337179/1‑68 (MQ=255)
                                                       ccaccCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCAt  >  2:403727/1‑68 (MQ=255)
                                                         |                                                                 
GGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGTCTATTAGTACCACCCGCCTGGCAGAACAACCCGGATATGGACCCGGAACTGCGTGCCTTCTTTGACTTTAACTCCAT  >  NC_000913/3355585‑3355707

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