Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,100,655 C→T 100% 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% 171.3 54576hslUATP‑dependent protease ATP‑binding subunit
Reads supporting (aligned to +/- strand):  new base (25/28):  ref base (0/0):  total (25/29)

GAACATGGACTGCAGCTGGCTGGTCATCTCTTCCATGCCCGGAGGAGCCATAATTTCAACGCCCAT  >  REL606/4100622‑4100687
                                 |                                
gAACATGGACTGCAGCTGGCTGGTCATCTCTTCTAt                                <  1:588384/36‑1 (MQ=255)
gAACATGGACTGCAGCTGGCTGGTCATCTCTTCTAt                                >  1:4509658/1‑36 (MQ=255)
 aaCATGGACTGCAGCTGGCTGGTCATCTCTTCTATg                               <  1:3407802/36‑1 (MQ=255)
 aaCATGGACTGCAGCTGGCTGGTCATCTCTTCTATg                               <  1:1877806/36‑1 (MQ=255)
 aaCATGGACTGCAGCTGGCTGGTCATCTCTTCTATg                               <  1:2122919/36‑1 (MQ=255)
 aaCATGGACTGCAGCTGGCTGGTCATCTCTTCTATg                               <  1:4641000/36‑1 (MQ=255)
   aaTGGACTGCAGCTGGCTGGTCATCTCTTCTATgtc                             >  1:731913/2‑34 (MQ=37)
    aTGGACTGCAGCTGGCTGGTCATCTCTTCTATGccc                            <  1:6895741/36‑1 (MQ=255)
    aTGGACTGCAGCTGGCTGGTCATCTCTTCTATGccc                            >  1:3108669/1‑36 (MQ=255)
    aTGGACTGCAGCTGGCTGGTCATCTCTTCTATGccc                            <  1:5337880/36‑1 (MQ=255)
     tGGACTGCAGCTGGCTGGTCATCTCTTCTATGCCCg                           >  1:3152178/1‑36 (MQ=255)
     tGGACTGCAGCTGGCTGGTCATCTCTTCTATGCCCg                           >  1:7155419/1‑36 (MQ=255)
      ggACTGCAGCTGGCTGGTCATCTCTTCTATGCCCgg                          <  1:4044038/36‑1 (MQ=255)
      ggACTGCAGCTGGCTGGTCATCTCTTCTATGCCCgg                          <  1:3002183/36‑1 (MQ=255)
      ggACTGCAGCTGGCTGGTCATCTCTTCTATGCCCgg                          <  1:2989129/36‑1 (MQ=255)
      ggACTGCAGCTGGCTGGCTATCTCTTATATGCCCgg                          >  1:7054857/1‑36 (MQ=19)
       gACTGCAGCTGGCTGGTCATCTCTTCTATGCCCgga                         >  1:2735043/1‑36 (MQ=255)
       gACTGCAGCTGGCTGGTCATCTCTTCTATGCCCgga                         >  1:5927524/1‑36 (MQ=255)
       gACTGCAGCTGGCTGGTCATCTCTTCTATGCCCgga                         >  1:6972915/1‑36 (MQ=255)
        aCTGCAGCTGGCTGGTCATCTCTTCTATGCCCggag                        <  1:4310351/36‑1 (MQ=255)
        aCTGCAGCTGGCTGGTCATCTCTTCTATGCCCggag                        <  1:725804/36‑1 (MQ=255)
        aCTGCAGCTGGCTGGTCATCTCTTCTATGCCCggag                        <  1:5459676/36‑1 (MQ=255)
        aCTGCAGCTGGCTGGTCATCTCTTCTATGCCCggag                        <  1:7087324/36‑1 (MQ=255)
         cTGCAGCTGGCTGGTGATCTCTTCTATGCACggagg                       >  1:2681910/1‑36 (MQ=25)
         cTGCAGCTGGCTGGTCATCTCTTCTATGCCCggagg                       <  1:7545042/36‑1 (MQ=255)
          tGCAGCTGGCTGGTCATCTCTTCTATGCCCggagga                      <  1:6101119/36‑1 (MQ=255)
           gCAGCTGGCTGGTCATCTCTTCTATGCCCGgaggag                     >  1:178648/1‑36 (MQ=255)
           gCAGCTGGCTGGTCATCTCTTCTATGCCCGgaggag                     >  1:2410065/1‑36 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   <  1:2759016/36‑1 (MQ=255)
               ctggctggTCATCTCTTCTATGCCCGGAGGCGCCAt                 >  1:430157/1‑36 (MQ=37)
               ctggctggTCATCTCTTCTATGCCCGGAGGAGCCAt                 >  1:7387400/1‑36 (MQ=255)
                tggctggTCATCTCTTCTATGCCCGGAGGAGCCATa                >  1:6180770/1‑36 (MQ=255)
                  gctggTCATCTCTTCTATGCCCGGAGGAGCCATAAt              >  1:5432138/1‑36 (MQ=255)
                   ctggTCATCTCTTCTATGCCCGGAGGAGCCATAAtt             >  1:4632308/1‑36 (MQ=255)
                    tggTCATCTCTTCAATGCGCGGGGGAGCCATAAttt            <  1:375769/36‑1 (MQ=37)
                     ggTCATCTCTCCTATGCCCGGAGGAGCCATAATTTc           <  1:4146912/36‑1 (MQ=38)
                     ggTCAGCTCTTCTATGCCCGGAGGAGCCATAATTTc           <  1:2331916/36‑1 (MQ=255)
                      gTCATCTCTTCTATGCCCGGAGGAGCCATAATTTCa          >  1:2414124/1‑36 (MQ=255)
                      gTCATCTCTTCTATGCCCGGAGGAGCCATAATTTCa          <  1:2347500/36‑1 (MQ=255)
                       tCATCTCTTCTATGCCCGGAGGAGCCATAATTTCaa         >  1:2184870/1‑36 (MQ=255)
                        cATCTCTTCTATGCCCGGAGGAGCCATAATTTCAAc        <  1:5513521/36‑1 (MQ=255)
                        cATCTCTTCTATGCCCGGAGGAGCCATAATTTCAAc        <  1:6892386/36‑1 (MQ=255)
                        cATCTCTTCTATGCCCGGAGGAGCCATAATTTCAAc        <  1:4270321/36‑1 (MQ=255)
                         aTCTCTTCTATGCCCGGAGGAGCCATAATTTCAACg       >  1:5837521/1‑36 (MQ=255)
                         aTCTCTTCTATGCCCGGAGGAGCCATAATTTCAACg       >  1:6713804/1‑36 (MQ=255)
                           ctctTCTATGCCCGGAGGAGCCATAATTTCAACGcc     >  1:1347694/1‑36 (MQ=255)
                            tctTCTATGCCCGGAGGAGCCATAATTTCAACGccc    <  1:2451465/36‑1 (MQ=255)
                            tctTCTATGCCCGGAGGAGCCATAATTTCAACGccc    <  1:7816773/36‑1 (MQ=255)
                            tctTCTACGCCAGGAGGAGCCATAATTTCAACGccc    <  1:7349093/36‑1 (MQ=25)
                              ttCTATGCCCGGCGCAGCCATAATTTCAACGCCCaa  >  1:3305371/1‑35 (MQ=25)
                              ttCTATGCCCGGAGGAGCCATAATTTCAACGCCCAt  <  1:7104536/36‑1 (MQ=255)
                              ttCTATGCCCGGAGGAGCCATAATTTCAACGCCCAt  <  1:5636922/36‑1 (MQ=255)
                              ttCTATGCCCGGAGGAGCCATAATTTCAACGCCCAt  >  1:5286437/1‑36 (MQ=255)
                              ttCTATGCCCGGAGGAGCCATAATTTCAACGCCCAt  >  1:7686072/1‑36 (MQ=255)
                                 |                                
GAACATGGACTGCAGCTGGCTGGTCATCTCTTCCATGCCCGGAGGAGCCATAATTTCAACGCCCAT  >  REL606/4100622‑4100687

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG