Predicted mutation
evidence seq id position mutation freq annotation gene description
JC REL606 3,289,962 Δ16 bp 100% coding (96‑111/4554 nt) gltB → glutamate synthase, large subunit

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 328996121 (0.350)14 (0.230) 13/70 NT 38.9% coding (95/4554 nt) gltB glutamate synthase, large subunit
?REL606 3289978 = 23 (0.380)coding (112/4554 nt) gltB glutamate synthase, large subunit

GGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATCGCC  >  REL606/3289944‑3290013
                                  |                                   
gggggAGGTTCACGATATGTTGTACGATAAATCCCt                                    >  1:830839/1‑36 (MQ=255)
 ggggAGGTTCACGATATGTTGTACGATAAATCCCtt                                   >  1:989366/1‑36 (MQ=255)
  gggAGGTTCACGATATGTTGTACGATAAATCCCTTg                                  >  1:5579321/1‑36 (MQ=255)
     ccgcagagcctgggggaggttcacgatatCTTgaga                               <  1:5613345‑M2/7‑1 (MQ=255)
     aGGTTCACGATATGTTGTACGATAAATCCCTTgaga                               >  1:1205430/1‑36 (MQ=255)
       gcagagcctgggggaggttcacgatatCTTGAGAgg                             >  1:848703‑M2/28‑36 (MQ=255)
       gTTCACGATATGTTGTACGATAAATCCCTTGAGAgg                             >  1:4759292/1‑36 (MQ=255)
        cagagcctgggggaggttaacgatatCTTGAGAggg                            >  1:359200‑M2/27‑36 (MQ=255)
         tCACGATATGTTGTACGATAAATCCCTTGAGAGGGa                           >  1:938930/1‑36 (MQ=255)
          gagcctgggggaggttcacgatatCTTGAGAGGGAt                          >  1:844267‑M2/25‑36 (MQ=255)
           aCGATATGTTGTACGATAAATCCCTTGAGAGGGATa                         >  1:4412783/1‑36 (MQ=255)
             cctgggggaggttcacgatatCTTGAGAGGGATAAc                       >  1:2599119‑M2/22‑36 (MQ=255)
               tatGTTGTACGATAAATCCCTTGAGAGGGATAACtg                     >  1:7088477/1‑36 (MQ=255)
                atGTTGTACGATAAATCCCTTGAGAGGGATAACtgt                    <  1:5558920/36‑1 (MQ=255)
                 tgttgtACGATAAATCCCTTGAGAGGGATAACtgtg                   <  1:2036259/36‑1 (MQ=255)
                  gttgtACGATAAATCCCTTGAGAGGGATAACTGTgg                  >  1:5833667/1‑36 (MQ=255)
                  gttgtACGATAAATCCCTTGAGAGGGATAACTGTgg                  >  1:1764113/1‑36 (MQ=255)
                  gggaggttcacgatatCTTGAGAGGGATAACTGTgg                  <  1:3878194‑M2/20‑1 (MQ=255)
                  gggaggttcacgatatCTTGAGAGGGATAACTGTgg                  >  1:4842882‑M2/17‑36 (MQ=255)
                   ttgtACGATAAATCCCTTGAGAGGGATAACTGTGGt                 <  1:3249762/36‑1 (MQ=255)
                    tgtACGATAAATCCCTTGAGAGGGATAACTGTGGtt                <  1:1908655/36‑1 (MQ=255)
                    tgtACGATAAATCCCTTGAGAGGGATAACTGTGGtt                <  1:3990760/36‑1 (MQ=255)
                     aggttcacgatatCTTGAGAGGGATAACTGTGGttt               <  1:3571257‑M2/23‑1 (MQ=255)
                      tACGATAAAGCCCTTGAGAGGGATAACTGTGGTTTc              <  1:1937046/36‑1 (MQ=255)
                        ttcacgatatCTTGAGAGGGATAACTGTGGTTTcgt            >  1:1089366‑M2/11‑35 (MQ=255)
                        ttcacgatatCTTGAGAGGGATAACTGTGGTTTCgg            >  1:2437415‑M2/11‑36 (MQ=255)
                          cacgatatCTTGAGAGGGATAACTGTGGTTTCGGcc          <  1:4097779‑M2/28‑1 (MQ=255)
                          aTAAATCCCTTGAGAGGGATAACTGTGGTTTCGGcc          >  1:6554410/1‑36 (MQ=255)
                           tAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCt         <  1:7324801/36‑1 (MQ=255)
                           acgatatCTTGAGAGGGATAACTGTGGTTTCGGCCt         >  1:221285‑M2/8‑36 (MQ=255)
                            aaaTCCCTTGAGAGGGATAACTGTGGTTTCGGCCTg        >  1:6785370/1‑36 (MQ=255)
                             gatatCTTGAGAGGGATAACTGTGGTTTCGGCCTGa       <  1:7078695‑M2/31‑1 (MQ=255)
                             gatatCTTGAGAGGGATAACTGTGGTTTCGGCCTGa       >  1:71720‑M2/6‑36 (MQ=255)
                               tCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATc     <  1:7438362/36‑1 (MQ=255)
                               tCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATc     <  1:613393/36‑1 (MQ=255)
                                cccTTGAGAGGGATAACTGTGGTTTCGGCCTGATCg    <  1:7155875/36‑1 (MQ=255)
                                cccTTGAGAGGGATAACTGTGGTTTCGGCCTGATCg    >  1:3068194/1‑36 (MQ=255)
                                  cTTGAGAGGGATAACTGTGGTTTCGGCCTGATCGcc  >  1:6370210/1‑36 (MQ=255)
                                  cTTGAGAGGGATAACTGTGGTTTCGGCCTGATCGcc  <  1:7268732/36‑1 (MQ=255)
                                  cTTGAGAGGGATAACTGTGGTTTCGGCCTGATCGcc  <  1:692685/36‑1 (MQ=255)
                                  cTTGAGAGGGATAACTGTGGTTTCGGCCTGATCGcc  >  1:6604156/1‑36 (MQ=255)
                                  cTTGAGAGGGATAACTGTGGTTTCGGCCTGATCGcc  >  1:656302/1‑36 (MQ=255)
                                  |                                   
GGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTGGTTTCGGCCTGATCGCC  >  REL606/3289944‑3290013

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG