Predicted mutation
evidence seq id position mutation freq annotation gene description
JC REL606 3,289,962 Δ16 bp 39.4% 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 = 328996120 (0.330)14 (0.230) 13/70 NT 39.4% coding (95/4554 nt) gltB glutamate synthase, large subunit
?REL606 3289978 = 23 (0.380)coding (112/4554 nt) gltB glutamate synthase, large subunit

GGGGTTCCCGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTG  >  REL606/3289926‑3289996
                                   |                                   
ggggTTCCCGCAGAGCCTGGGGGAGGTTCACGatat                                     <  1:7744543/36‑1 (MQ=255)
 gggTTCCCGCAGAGCCTGGGGGAGGTTCACGATAtg                                    >  1:6963722/1‑36 (MQ=255)
   gTTCCAGCAGAGCCTGGGGGAGGTTCACGATAtgtt                                  <  1:2571053/36‑1 (MQ=255)
      cccGCATAACCTGGGGGATGTTGCCGATATGTTGTa                               <  1:6922387/36‑1 (MQ=12)
       ccGCAGAGCCTGGGGGAGGTTCACGATATcttgaga                              <  1:5613345‑M1/36‑8 (MQ=255)
       ccGCAGAGCCTGGGGGAGGTTCACGATATGTTGTAc                              <  1:4483416/36‑1 (MQ=255)
       ccGCAGAGCCTGGGGGAGGTTCACGATATGTTGTAc                              >  1:7621355/1‑36 (MQ=255)
        cGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACg                             <  1:1967300/36‑1 (MQ=255)
         gCAGAGCCTGGGGGAGGTTCACGATATcttgagagg                            >  1:848703‑M1/1‑27 (MQ=255)
          cAGAGCCTGGGGGAGGTTAACGATATcttgagaggg                           >  1:359200‑M1/1‑26 (MQ=255)
           agagCCTGGGGGAGGTTCACGATATGTTGTACGATa                          >  1:6221978/1‑36 (MQ=255)
            gagCCTGGGGGAGGTTCACGATATcttgagagggat                         >  1:844267‑M1/1‑24 (MQ=255)
            gagCCTGGGGGAGGTTCACGATATGTTGTACGATaa                         >  1:6604058/1‑36 (MQ=255)
             agCCTGGGGGAGGTTCACGATATGTTGTACGATaaa                        >  1:3907001/1‑36 (MQ=255)
              gCCTGGGGGAGGTTCACGATATGTTGTACGATAAAt                       <  1:3495515/36‑1 (MQ=255)
               ccTGGGGGAGGTTCACGATATcttgagagggataac                      >  1:2599119‑M1/1‑21 (MQ=255)
               ccTGGGGGAGGTTCACGATATGTTGTACGATAAATc                      >  1:6929097/1‑36 (MQ=255)
                  gggggAGGTTCACGATATGTTGTACGATAAATCCCt                   >  1:830839/1‑36 (MQ=255)
                   ggggAGGTTCACGATATGTTGTACGATAAATCCCtt                  >  1:989366/1‑36 (MQ=255)
                    gggAGGTTCACGATATcttgagagggataactgtgg                 >  1:4842882‑M1/1‑16 (MQ=255)
                    gggAGGTTCACGATATcttgagagggataactgtgg                 <  1:3878194‑M1/36‑21 (MQ=255)
                    gggAGGTTCACGATATGTTGTACGATAAATCCCTTg                 >  1:5579321/1‑36 (MQ=255)
                       aGGTTCACGATATcttgagagggataactgtggttt              <  1:3571257‑M1/36‑24 (MQ=255)
                       aGGTTCACGATATGTTGTACGATAAATCCCTTgaga              >  1:1205430/1‑36 (MQ=255)
                         gTTCACGATATGTTGTACGATAAATCCCTTGAGAgg            >  1:4759292/1‑36 (MQ=255)
                          ttCACGATATcttgagagggataactgtggtttcgt           >  1:1089366‑M1/1‑10 (MQ=255)
                          ttCACGATATcttgagagggataactgtggtttcgg           >  1:2437415‑M1/1‑10 (MQ=255)
                           tCACGATATGTTGTACGATAAATCCCTTGAGAGGGa          >  1:938930/1‑36 (MQ=255)
                            cACGATATcttgagagggataactgtggtttcggcc         <  1:4097779‑M1/36‑29 (MQ=255)
                             aCGATATcttgagagggataactgtggtttcggcct        >  1:221285‑M1/1‑7 (MQ=255)
                             aCGATATGTTGTACGATAAATCCCTTGAGAGGGATa        >  1:4412783/1‑36 (MQ=255)
                               gATATcttgagagggataactgtggtttcggcctga      <  1:7078695‑M1/36‑32 (MQ=255)
                               gATATcttgagagggataactgtggtttcggcctga      >  1:71720‑M1/1‑5 (MQ=255)
                                 tatGTTGTACGATAAATCCCTTGAGAGGGATAACtg    >  1:7088477/1‑36 (MQ=255)
                                  atGTTGTACGATAAATCCCTTGAGAGGGATAACtgt   <  1:5558920/36‑1 (MQ=255)
                                   tgttgtACGATAAATCCCTTGAGAGGGATAACtgtg  <  1:2036259/36‑1 (MQ=255)
                                   |                                   
GGGGTTCCCGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTG  >  REL606/3289926‑3289996

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG