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% 165.9 49576hslUATP‑dependent protease ATP‑binding subunit
Reads supporting (aligned to +/- strand):  new base (32/16):  ref base (0/1):  total (32/17)

GAACATGGACTGCAGCTGGCTGGTCATCTCTTCCATGCCCGGAGGAGCCATAATTTCAACGCCCAT  >  REL606/4100622‑4100687
                                 |                                
gAACATGGACTGCAGTTGGCTGGTCATCTCTTCTAt                                <  1:2617771/36‑1 (MQ=37)
gAACATGGACTGCAGCTGGCTGGTCATCTCTTCTAt                                >  1:3567849/1‑36 (MQ=255)
 aaCATGGACTGCAGCTGGCTGGTCATCTCTTCTATg                               >  1:3723329/1‑36 (MQ=255)
  aCATGGACTGCAGCTGGCTGGTCATCTCTTCTATGc                              >  1:5004391/1‑36 (MQ=255)
   cATGGACTGCAGCTGGCTGGTCATCTCTTCTATGcc                             >  1:1229915/1‑36 (MQ=255)
    aTGGACTGCAGCTGGCTGGTCATCTCTTCTATGccc                            >  1:7503824/1‑36 (MQ=255)
    aTGGACTGCAGCTGGCTGGAATTGTCTTCTATGcac                            >  1:1574154/1‑34 (MQ=12)
     tGGACTGCAGCTGGCTGGTCATCTCTTCTATGCCCg                           >  1:2645529/1‑36 (MQ=255)
      ggACTGCAGCTGGCTGGTCATCTCTTCTATGCCCgg                          >  1:5032911/1‑36 (MQ=255)
       gccTGCAGCTGGCTGGTCATCTCTTCTATGCCCgga                         <  1:4128253/34‑1 (MQ=37)
       gACTGCAGCTGGCTGGTCATCTCTTCTATTCCCgga                         >  1:1193887/1‑36 (MQ=37)
        aCTGCAGCTGGCTGGTCATCTCTTATATGCCCggag                        >  1:7747901/1‑36 (MQ=37)
         cTGCAGCTGGCTGGTCATCTCTTCTATGCCCggagg                       <  1:3527625/36‑1 (MQ=255)
          tGCAGCTGGCTGGTCATCTCTTCTATGCCCggagga                      >  1:1553857/1‑36 (MQ=255)
            cAGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGc                    >  1:2577206/1‑36 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   >  1:1337765/1‑36 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   <  1:456060/36‑1 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   >  1:7490552/1‑36 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   >  1:3292361/1‑36 (MQ=255)
             aGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGcc                   >  1:6205748/1‑36 (MQ=255)
              gctggctgGTCATCTCTTCTATGCCCGGAGGAGCCa                  >  1:5822343/1‑36 (MQ=255)
              gctggctgGTCATCTCTTCTATGCCCGGAGGAGCCa                  >  1:7008856/1‑36 (MQ=255)
              gctggctgGTCATCTCTTCTATGCCCGGAGGAGCCa                  >  1:1348845/1‑36 (MQ=255)
               ctggctggTCATCTCTTCTATGCCCGGAGGAGCCAt                 <  1:4685856/36‑1 (MQ=255)
                tggctggTCATCTCTTCTATGCCCGGAGGAGCCATa                >  1:5060902/1‑36 (MQ=255)
                  gctggTCATCTCTTCTATGCCCGGAGGAGCCATAAt              >  1:6234200/1‑36 (MQ=255)
                  gctggTCATCTCTTCTATGCCCGGAGGAGCCATAAt              <  1:6719552/36‑1 (MQ=255)
                  gctggTCATCTCTTCTATGCCCGGAGGAGCCATAAt              <  1:248926/36‑1 (MQ=255)
                   ctggTCATCTCTTCTATGCCCGGAGGAGCCATAAtt             >  1:4179646/1‑36 (MQ=255)
                   ctggTCATCTCTTCTATGCCCGGAGGAGCCATAAtt             <  1:3418425/36‑1 (MQ=255)
                   ctggCCATCTCTTCTACGCCAGGAGGAGCCATAAtt             <  1:2962344/36‑1 (MQ=19)
                     ggTCATCTCTTCTATGCCCGGGGGAGCCATAATTTc           >  1:6759582/1‑36 (MQ=38)
                      ggcATCTCTTCTATGCCCGGAGGAGCCATAATTTCa          >  1:7094846/3‑36 (MQ=255)
                        cATCTCTTCTATGCCCGGAGGAGCCATAATTTCAAc        >  1:5988300/1‑36 (MQ=255)
                        cATCTCTTCTATGCCCGGAGGAGCCATAATTTCAAc        >  1:6953980/1‑36 (MQ=255)
                        cATCTCTTCTATGCCCGGAGGAGCCATAATTTCAAc        >  1:2837640/1‑36 (MQ=255)
                        cATCTCTTCTATGCCCGGAGGAGCCATAATTTCAAc        <  1:2130829/36‑1 (MQ=255)
                         aTCTCTTCTATGCCCGGAGGAGCCATAATTTCAACg       >  1:6978514/1‑36 (MQ=255)
                         aTCTCTTCTATGCCCGGAGGAGCCATAATTTCAACg       >  1:1069720/1‑36 (MQ=255)
                          tctcTTCTATGCCCGGAGGAGCCATAATTTCAACGc      >  1:403472/1‑36 (MQ=255)
                           ctctTCTATGCCCGGAGGAGCCATAATTTCAACGcc     <  1:700124/36‑1 (MQ=255)
                           ctctTCTATGCCCGGAGGAGCCATAATTTCAACGcc     <  1:8196272/36‑1 (MQ=255)
                            ttttCTATGCCCGTGGGAGCCATAATTTCAACGccc    <  1:8397614/34‑1 (MQ=21)
                            tctTCTATGCCCGGAGGAGCCATAATTTCAACGccc    <  1:1949012/36‑1 (MQ=255)
                              ttCTATGCCCGGAGGAGCCATAATTTCAACGCCCAt  >  1:6168273/1‑36 (MQ=255)
                              ttCTATGCCCGGAGGAGCCATAATTTCAACGCCCAt  >  1:5892211/1‑36 (MQ=255)
                              ttCTATGCCCGGAGGAGCCATAATTTCAACGCCCAt  <  1:4009289/36‑1 (MQ=255)
                              ttCTATGCCCAGAGGAGCCATAATTTCAACGCCCAt  <  1:1758705/36‑1 (MQ=37)
                              ttCCATGCCCGGAGGAGCCATAATTTCAACGCCCAt  <  1:8247527/36‑1 (MQ=255)
                                 |                                
GAACATGGACTGCAGCTGGCTGGTCATCTCTTCCATGCCCGGAGGAGCCATAATTTCAACGCCCAT  >  REL606/4100622‑4100687

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