Predicted mutation
evidence seq id position mutation freq annotation gene description
JC REL606 3,289,962 Δ16 bp 91.9% 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 = 32899613 (0.050)57 (0.920) 41/70 NT 91.9% coding (95/4554 nt) gltB glutamate synthase, large subunit
?REL606 3289978 = 7 (0.110)coding (112/4554 nt) gltB glutamate synthase, large subunit

GGGGTTCCCGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTG  >  REL606/3289926‑3289996
                                   |                                   
ggggTTCCCGCAGAGCCTGGGGGAGGTTCACGatat                                     >  1:6404611/1‑36 (MQ=255)
 gggTTCCCGCAGAGCCTGGGGGAGGTTCACGATATc                                    >  1:2535021‑M1/1‑35 (MQ=255)
  ggTTCCCGCAGAGCCTTGGGGAGGTTCACGATATct                                   >  1:1784363‑M1/1‑34 (MQ=255)
   gTTCCCGCAGAGCCTGGGGGAGGTTCACGATATctt                                  >  1:1935873‑M1/1‑33 (MQ=255)
   gTTCCCGCAGAGCCTGGGGGAGGTTCACGATATctt                                  >  1:4769754‑M1/1‑33 (MQ=255)
    ttCCCGCAGAGCCTGGGGGAGGTTCACGATATcttg                                 <  1:239180‑M1/36‑5 (MQ=255)
    ttCCCGCAGAGCCTGGGGGAGGTTCACGATATcttg                                 >  1:2748387‑M1/1‑32 (MQ=255)
      cccGCAGAGCCTGGGGGAGGTTCACGATATcttgag                               >  1:1939056‑M1/1‑30 (MQ=255)
      cccGCAGAGCCTGGGGGAGGTTCACGATATcttgag                               <  1:7470444‑M1/36‑7 (MQ=255)
      cccGCAGAGCCTGGGGGAGGTTCACGATATcttgag                               >  1:721690‑M1/1‑30 (MQ=255)
      cccGCAGAGCCTGGGGGAGGTTCACGATATcttgag                               <  1:6030321‑M1/36‑7 (MQ=255)
       ccGCAGAGCCTGGGGGAGGTTCACGATATGTTGTag                              >  1:2581656/1‑35 (MQ=255)
        cGCAGAGCCTGGGGGAGGTTCACGATATcttgagag                             <  1:2252758‑M1/36‑9 (MQ=255)
         gCAGAGCCTGGGGGAGGTTCACGATATcttgagagg                            >  1:3754236‑M1/1‑27 (MQ=255)
         gCAGAGCCTGGGGGAGGTTCACGATATcttgagagg                            <  1:3398774‑M1/36‑10 (MQ=255)
          cAGAGCCTGGGGGAGGTTCACGATATcttgagaggg                           <  1:821992‑M1/36‑11 (MQ=255)
          cAGAGCCTGGGGGAGGTTCACGATATcttgagaggg                           >  1:5238448‑M1/1‑26 (MQ=255)
           agagCCTGGGGGAGGTTCACGATATcttgagaggga                          >  1:6940953‑M1/1‑25 (MQ=255)
           agagCCTGGGGGAGGTTCACGATATcttgagaggga                          >  1:1494100‑M1/1‑25 (MQ=255)
            gagCCTGGGGGAGGTTCACGATATcttgagagggat                         <  1:2703217‑M1/36‑13 (MQ=255)
               ccTGGGGGAGGTTCACGATATcttgagagggataac                      >  1:6878744‑M1/1‑21 (MQ=255)
               ccTGGGGGAGGTTCACGATATcttgagagggataac                      <  1:6854547‑M1/36‑16 (MQ=255)
               ccTGGGGGAGGTTCACGATATcttgagagggataac                      <  1:5724242‑M1/36‑16 (MQ=255)
               ccTGGGGGAGGTTCACGATATcttgagagggataac                      <  1:1691600‑M1/36‑16 (MQ=255)
                cTGGGGGAGGTTCACGATATcttgagagggataact                     <  1:7170948‑M1/36‑17 (MQ=255)
                 tGGGGGAGGTTCACGATATcttgagagggataactg                    >  1:465866‑M1/1‑19 (MQ=255)
                  gggggAGGTTCACGATATcttgagagggataactgt                   >  1:2110004‑M1/1‑18 (MQ=255)
                  gggggAGGTTCACGATATcttgagagggataactgt                   <  1:4425170‑M1/36‑19 (MQ=255)
                  gggggAGGTTCAAGATATcttgagagggataactgt                   <  1:7492848‑M1/36‑19 (MQ=255)
                   ggggAGGTTCACGATATcttgagagggataactgtg                  >  1:4479537‑M1/1‑17 (MQ=255)
                    gggAGGTTCACGATATcttgagagggataactgtgg                 >  1:5476304‑M1/1‑16 (MQ=255)
                    gggAGGTTCACGATATcttgagagggataactgtgg                 >  1:3397284‑M1/1‑16 (MQ=255)
                     ggAGGTTCACGATATcttgagagggataactgtggt                >  1:4092632‑M1/1‑15 (MQ=255)
                     ggAGGTTCACGATATcttgagagggataactgtggt                <  1:4073875‑M1/36‑22 (MQ=255)
                      gAGGTTCACGATATcttgagagggataactgtggtt               <  1:2437678‑M1/36‑23 (MQ=255)
                        ggTTCACGATATcttgagagggataactgtggtttc             <  1:6755974‑M1/36‑25 (MQ=255)
                        ggTTCACGATATcttgagagggataactgtggtttc             <  1:3366475‑M1/36‑25 (MQ=255)
                        ggTTCACGATATcttgagagggataactgtggtttc             >  1:7591963‑M1/1‑12 (MQ=255)
                         gTTCACGATATcttgagagggataactgtggtttcg            <  1:7061531‑M1/36‑26 (MQ=255)
                         gTTCACGATATcttgagagggataactgtggtttcg            >  1:7434715‑M1/1‑11 (MQ=255)
                         gTTCACGATATcttgagagggataactgtggtttcg            >  1:8011210‑M1/1‑11 (MQ=255)
                          ttCACGATATcttgagagggataactgtggtttcgg           <  1:5844015‑M1/36‑27 (MQ=255)
                          ttCACGATATcttgagagggataactgtggtttcgg           <  1:6284943‑M1/36‑27 (MQ=255)
                          ttCACGATATGTTGTACGATAAATCCCTTGAGAggg           >  1:5162283/1‑36 (MQ=255)
                           tCACGATATcttgagagggataactgtggtttcggc          >  1:5505034‑M1/1‑9 (MQ=255)
                           tCACGATATcttgagagggataactgtggtttcggc          >  1:238133‑M1/1‑9 (MQ=255)
                            cACGATATcttgagagggataactgtggtttcggcc         <  1:2135234‑M1/36‑29 (MQ=255)
                            cACGATATcttgagagggataactgtggtttcggcc         >  1:478695‑M1/1‑8 (MQ=255)
                             aCGATATcttgagagggatanctgtggtttcggcct        >  1:2025370‑M1/1‑7 (MQ=255)
                              cGATATcttgagagggataactatggtttcggcctg       >  1:6099526‑M1/1‑6 (MQ=255)
                               gATATcttgagagggataactgtggtttcggcctga      <  1:4172767‑M1/36‑32 (MQ=255)
                               gATATcttgagagggataactgtggtttcggcctga      <  1:4189255‑M1/36‑32 (MQ=255)
                                atatcttgagagggataactgtggtttcggcctgat     >  1:862925‑M1/1‑4 (MQ=255)
                                atatGTTGTACGATAAATCCCTTGAGAGGGATAACt     <  1:2871307/36‑1 (MQ=255)
                                 tatcttgagagggataactgtggtttcggcctgatc    <  1:2969510‑M1/36‑34 (MQ=255)
                                 tatcttgagagggataactgtggtttcggcctgatc    <  1:7927168‑M1/36‑34 (MQ=255)
                                  atcttgagagggataactgtggtttcggcctgatcg   >  1:3213557‑M1/1‑2 (MQ=255)
                                  atcttgagagggataactgtggtttcggcctgatcg   <  1:8127944‑M1/36‑35 (MQ=255)
                                   tcttgagagggataactgtggtttcggcctgatcgc  >  1:3005160‑M1/1‑1 (MQ=255)
                                   tcttgagagggataactgtggtttcggcctgatcgc  >  1:6092144‑M1/1‑1 (MQ=255)
                                   tcttgagagggataactgtggtttcggcctgatcgc  >  1:3863835‑M1/1‑1 (MQ=255)
                                   |                                   
GGGGTTCCCGCAGAGCCTGGGGGAGGTTCACGATATGTTGTACGATAAATCCCTTGAGAGGGATAACTGTG  >  REL606/3289926‑3289996

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