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% 123.7 39576hslUATP‑dependent protease ATP‑binding subunit
Reads supporting (aligned to +/- strand):  new base (28/11):  ref base (0/0):  total (28/11)

GAACATGGACTGCAGCTGGCTGGTCATCTCTTCCATGCCCGGAGGAGCCATAATTTCAACGCCCAT  >  REL606/4100622‑4100687
                                 |                                
gAACATGGACTGCAGCTGGCTGGTCATGTCTTCTAt                                >  2:3096256/1‑36 (MQ=37)
gAACATGGACTGCAGCTGGCTGGTCATCTCTTCTAt                                >  1:1662551/1‑36 (MQ=255)
gAACAGGGGCTGCAGCTGGCTGGTCATCTCTTCTAt                                <  1:3340617/36‑1 (MQ=25)
 aaCATGGACTGCAGCTGGCTGGTCATCTCTTCTATg                               >  1:3380696/1‑36 (MQ=255)
   cATGGACTGCAGCTGGCTGGTCATCTCTTCTATGcc                             >  1:3131550/1‑36 (MQ=255)
   cATGGACTGCAGCTAGCTGGTCATCTCTTCTATGcc                             >  2:1581912/1‑36 (MQ=37)
    aTGGACTGCAGCTGGCTGGTCATCTCTTCTATGccc                            >  1:1324348/1‑36 (MQ=255)
     tGGACTGCAGCTGGCTGGTCATCTCTTCTATGCCCg                           <  2:2011577/36‑1 (MQ=255)
      ggACTGCAGCTGGCTGGTCATCTCTTCTATGCCCgg                          >  2:2425683/1‑36 (MQ=255)
      ggACTGCAGCTGGCTGGTCATCTCTTCTACGCCCgg                          >  2:3713907/1‑36 (MQ=37)
       gACTGCAGCTGGCTGGTCATCTCTTCTATGCCCgga                         >  2:1578292/1‑36 (MQ=255)
        aCTGCAGCTGGCTGGTCCTCTCTTCTATGCCCggag                        >  2:738610/1‑36 (MQ=37)
        aCTGCAGCTGGCTGGTCATCTCTTCTATGCCCggag                        <  2:1316959/36‑1 (MQ=255)
        aCTGCAGCTGGCGGGTCATCACTTCTATGCCCggag                        >  2:289620/1‑36 (MQ=25)
         cTGCAGCTGGCTGGTCATCTCTTCTATGCCCggagg                       >  2:3564844/1‑36 (MQ=255)
         cTGCAGCTGGCTGGTCATCTCTTCTATGCCCggagg                       <  2:2777770/36‑1 (MQ=255)
         cTGCAGCTGGCTGGTCATCTCTTCTATGCCCggaag                       >  1:2735439/1‑34 (MQ=37)
          tGCAGCTGGCTGGTCATCTCTTCTATGCCCggagga                      >  2:3392853/1‑36 (MQ=255)
          tGCAGCTGGCTGGTCATCTCTTCTATGCCCggagga                      <  2:3019063/36‑1 (MQ=255)
           gCAGCTGGCTGGTCATCTCTTCTATGCCCGgaggag                     <  2:106450/36‑1 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   >  2:444456/1‑36 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   <  2:829130/36‑1 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   >  1:1028490/1‑36 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   >  1:2266289/1‑36 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   <  1:3661005/36‑1 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTACGCCCGGAGCAGcc                   >  1:1735116/1‑36 (MQ=25)
              gctggttgGTCATCTCTTCTATGCCCGGAGGAGCCa                  >  1:400876/1‑36 (MQ=37)
              gctggctgGTCATCTCTTCTATGCCCGGAGGAGCCa                  <  2:3300957/36‑1 (MQ=255)
                tggctggTCATCTCTTGTATGCCGGGAGGAGCCATa                <  1:82341/36‑1 (MQ=37)
                   ctggTCATCTCTTCTATGCCCGGAGGTGCCATAAtt             >  2:2862256/1‑36 (MQ=37)
                    ttttCATCTCTTCTATTCCGGGAGGAGCCATAAttt            <  1:428870/33‑1 (MQ=19)
                    tggTCATCTCTTCTATGCCCGGAGGAGCCATAAttt            >  2:2649753/1‑36 (MQ=255)
                     ggTCATCTCTTCTATGCCCGGAGGAGCCATAATTTc           >  2:406998/1‑36 (MQ=255)
                      gTCATCTCTTCTATGCCCGGAGGAGCCATAATTTCa          >  1:3230597/1‑36 (MQ=255)
                       tCATCTCTTCTATGCCCGGAGGAGCCATAATTTCaa         >  1:2030543/1‑36 (MQ=255)
                        cATCTCTTCTATGCCCGGCGGAGCCATACTTTCAAc        >  2:2107654/1‑36 (MQ=37)
                          tctcTTCTATGCCCGGAGGAGCCATAATTTCAACGc      >  2:2284953/1‑36 (MQ=255)
                             ctTCTATGCCCGGAGGAGCCATAATTTCAACGCCCa   >  1:2668308/1‑36 (MQ=255)
                              ttCTATGCCCGGAGGAGCCATAATTTCAACGCCCAt  >  2:375429/1‑36 (MQ=255)
                                 |                                
GAACATGGACTGCAGCTGGCTGGTCATCTCTTCCATGCCCGGAGGAGCCATAATTTCAACGCCCAT  >  REL606/4100622‑4100687

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