Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,100,655 C→T M192I (ATG→ATA hslU ATP‑dependent protease ATP‑binding subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,100,6550C→T100.0% 114.1 37M192I (ATG→ATAhslUATP‑dependent protease ATP‑binding subunit
Reads supporting (aligned to +/- strand):  new base (26/11):  ref base (0/0):  total (26/11)

GAACATGGACTGCAGCTGGCTGGTCATCTCTTCCATGCCCGGAGGAGCCATAATTTCAACGCCCAT  >  REL606/4100622‑4100687
                                 |                                
gAACATGGACTGCAGCTGGCTGGTCATGTCTTCTAt                                >  1:6835466/1‑36 (MQ=25)
gAACATGGACTGCAGCTGGCTGGTCATCTCTTCTAt                                >  1:1662551/1‑36 (MQ=255)
gAACAGGGGCTGCAGCTGGCTGGTCATCTCTTCTAt                                <  1:3340617/36‑1 (MQ=21)
 aaCATGGACTGCAGCTGGCTGGTCATCTCTTCTATg                               >  1:3380696/1‑36 (MQ=255)
   cATGGACTGCAGCTGGCTGGTCATCTCTTCTATGcc                             >  1:3131550/1‑36 (MQ=255)
   cATGGACTGCAGCTAGCTGGTCATCTCTTCTATGcc                             >  1:5321122/1‑36 (MQ=25)
    aTGGACTGCAGCTGGCTGGTCATCTCTTCTATGccc                            >  1:1324348/1‑36 (MQ=255)
     tGGACTGCAGCTGGCTGGTCATCTCTTCTATGCCCg                           <  1:5750787/36‑1 (MQ=255)
      ggACTGCAGCTGGCTGGTCATCTCTTCTATGCCCgg                          >  1:6164893/1‑36 (MQ=255)
      ggACTGCAGCTGGCTGGTCATCTCTTCTACGCCCgg                          >  1:7453117/1‑36 (MQ=25)
       gACTGCAGCTGGCTGGTCATCTCTTCTATGCCCgga                         >  1:5317502/1‑36 (MQ=255)
        aCTGCAGCTGGCTGGTCCTCTCTTCTATGCCCggag                        >  1:4477820/1‑36 (MQ=25)
        aCTGCAGCTGGCTGGTCATCTCTTCTATGCCCggag                        <  1:5056169/36‑1 (MQ=255)
         cTGCAGCTGGCTGGTCATCTCTTCTATGCCCggagg                       >  1:7304054/1‑36 (MQ=255)
         cTGCAGCTGGCTGGTCATCTCTTCTATGCCCggagg                       <  1:6516980/36‑1 (MQ=255)
         cTGCAGCTGGCTGGTCATCTCTTCTATGCCCggaag                       >  1:2735439/1‑34 (MQ=25)
          tGCAGCTGGCTGGTCATCTCTTCTATGCCCggagga                      >  1:7132063/1‑36 (MQ=255)
          tGCAGCTGGCTGGTCATCTCTTCTATGCCCggagga                      <  1:6758273/36‑1 (MQ=255)
           gCAGCTGGCTGGTCATCTCTTCTATGCCCGgaggag                     <  1:3845660/36‑1 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   <  1:4568340/36‑1 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   >  1:4183666/1‑36 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   <  1:3661005/36‑1 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   >  1:2266289/1‑36 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   >  1:1028490/1‑36 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTACGCCCGGAGCAGcc                   >  1:1735116/1‑36 (MQ=21)
              gctggttgGTCATCTCTTCTATGCCCGGAGGAGCCa                  >  1:400876/1‑36 (MQ=25)
              gctggctgGTCATCTCTTCTATGCCCGGAGGAGCCa                  <  1:7040167/36‑1 (MQ=255)
                tggctggTCATCTCTTGTATGCCGGGAGGAGCCATa                <  1:82341/36‑1 (MQ=21)
                   ctggTCATCTCTTCTATGCCCGGAGGTGCCATAAtt             >  1:6601466/1‑36 (MQ=25)
                    ttttCATCTCTTCTATTCCGGGAGGAGCCATAAttt            <  1:428870/33‑1 (MQ=18)
                    tggTCATCTCTTCTATGCCCGGAGGAGCCATAAttt            >  1:6388963/1‑36 (MQ=255)
                     ggTCATCTCTTCTATGCCCGGAGGAGCCATAATTTc           >  1:4146208/1‑36 (MQ=255)
                      gTCATCTCTTCTATGCCCGGAGGAGCCATAATTTCa          >  1:3230597/1‑36 (MQ=255)
                       tCATCTCTTCTATGCCCGGAGGAGCCATAATTTCaa         >  1:2030543/1‑36 (MQ=255)
                          tctcTTCTATGCCCGGAGGAGCCATAATTTCAACGc      >  1:6024163/1‑36 (MQ=255)
                             ctTCTATGCCCGGAGGAGCCATAATTTCAACGCCCa   >  1:2668308/1‑36 (MQ=255)
                              ttCTATGCCCGGAGGAGCCATAATTTCAACGCCCAt  >  1:4114639/1‑36 (MQ=255)
                                 |                                
GAACATGGACTGCAGCTGGCTGGTCATCTCTTCCATGCCCGGAGGAGCCATAATTTCAACGCCCAT  >  REL606/4100622‑4100687

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 16 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG