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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 3289962 3289977 16 28 [2] [0] 27 gltB glutamate synthase, large subunit

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 32899612 (0.040)26 (0.580) 21/70 NT 96.3% coding (95/4554 nt) gltB glutamate synthase, large subunit
?REL606 3289978 = 0 (0.000)coding (112/4554 nt) gltB glutamate synthase, large subunit

GGGTTCCCGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTG  >  REL606/3289927‑3289996
                                  |                                   
gggTTCCCGCAGAGCCTGGGGGAGGTTCACGATATc                                    <  1:114273‑M1/36‑2 (MQ=255)
   ttCCCGCAGAGCCTGGGGGAGGTTCACGATATcttg                                 >  1:2549687‑M1/1‑32 (MQ=255)
   tgcccGCAGAGCCTGGGAGCGGTTCAAGATAtcttg                                 >  1:1964158/3‑36 (MQ=16)
    tCCCGCAGAGCCTGGGGGAGGTTCACGATATcttga                                <  1:3141426‑M1/36‑6 (MQ=255)
      ccGGAGAGCCTGGGGGGGGTTCACGATAtcttgaga                              <  1:2356011/36‑4 (MQ=19)
      ccGCAGAGCCTGGGGGAGGTTCACGATATcttgaga                              >  1:3537964‑M1/1‑29 (MQ=255)
      ccGCAGAGCCTGGGGGAGGTTCACGATATcttgaga                              >  1:212150‑M1/1‑29 (MQ=255)
        gCAGAGCCTGGGGGAGGTTCACGATATcttgagagg                            <  1:3666502‑M1/36‑10 (MQ=255)
         cAGAGCCTGGGGGAGGTTCACGATATcttgagaggg                           >  1:479698‑M1/1‑26 (MQ=255)
           gagGCTGGGGGAGGTTCACGATATcttgagagggat                         <  1:2718163‑M1/36‑13 (MQ=255)
            agCCTGGGGGAGGTTCACTATATcttgagagggata                        <  1:2284033‑M1/36‑14 (MQ=255)
            agCCTGGGGGAGGTTCACGATATcttgagagggata                        >  1:2769377‑M1/1‑23 (MQ=255)
            agCCTGGGGGAGGTTCACGATATcttgagagggata                        >  1:1091927‑M1/1‑23 (MQ=255)
             gCCTGGGGGAGGTTCACGATATcttgagagggataa                       >  1:725672‑M1/1‑22 (MQ=255)
               cTGGGGGAGGTTCACGATATcttgagagggataact                     >  1:1537984‑M1/1‑20 (MQ=255)
               cTGGGGGAGGTTCACGATATcttgagagggataact                     <  1:734381‑M1/36‑17 (MQ=255)
                 gggggAGGTTCACGATATcttgagagggataactgt                   <  1:2331652‑M1/36‑19 (MQ=255)
                 gggggAGGTTCACGATATcttgagagggataactgt                   <  1:4453292‑M1/36‑19 (MQ=255)
                 gggggAGGTTCACGATATcttgagagggataactgt                   <  1:675965‑M1/36‑19 (MQ=255)
                  ggggAGGTTCACGATATcttgagagggataactgtg                  >  1:4836491‑M1/1‑17 (MQ=255)
                        gTTCACGATATcttgagagggataactgtggtttcg            >  1:2439790‑M1/1‑11 (MQ=255)
                         ttCACGATATcttgagagggataactgtggtttcgg           >  1:2284574‑M1/1‑10 (MQ=255)
                             cGATATcttgagagggataactgtggtttcggcctg       <  1:3836818‑M1/36‑31 (MQ=255)
                              gATATcttgagagggataactgtggtttcggcctga      >  1:2561014‑M1/1‑5 (MQ=255)
                                tatcttgagagggataactgtggtttcggcctgatc    >  1:4012450‑M1/1‑3 (MQ=255)
                                 atcttgagagggataactgtggtttcggcctgatcg   >  1:2794408‑M1/1‑2 (MQ=255)
                                 atcttgagagggataactgtggtttcggcctgatcg   >  1:2071972‑M1/1‑2 (MQ=255)
                                  tcttgagagggataactgtggtttcggcctgatcgc  >  1:3576220‑M1/1‑1 (MQ=255)
                                  |                                   
GGGTTCCCGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTG  >  REL606/3289927‑3289996

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG