Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,100,655 C→T 16.6% 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→T16.6% 10.7 ‑0.9 ‑0.042576hslUATP‑dependent protease ATP‑binding subunit
Reads supporting (aligned to +/- strand):  new base (2/5):  ref base (19/16):  total (21/21)
Fisher's exact test strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.63e-01

GGAACATGGACTGCAGCTGGCTGGTCATCTCTTCCATGCCCGGAGGAGCCATAATTTCAACGCCCATCGG  >  REL606/4100621‑4100690
                                  |                                   
ggAACATGGACTGCAGCTGGCTGGTCATCTCTTCCa                                    >  1:1499729/1‑36 (MQ=255)
ggAACATGGACTGCAGCTGGCTGGTCATCTCTTCCa                                    <  1:6267884/36‑1 (MQ=255)
 gAACATGGACTGCAGCTGGCTGGTCATCTCTTCCAt                                   >  1:3711380/1‑36 (MQ=255)
 gAACATGGACTGCAGCTGGCTGGTCATCTCTTCCAt                                   <  1:175019/36‑1 (MQ=255)
   aCATGGACTGCCGCTGTCTGGTCCTCTCTTCCCTGc                                 <  1:1935198/36‑1 (MQ=16)
   aCATGGACTGCAGCTGGCTGGTCATCTCTTCCATGc                                 >  1:6679590/1‑36 (MQ=255)
    cATGGACTGCAGCTGGCTGGTCATCTCTTCCATGcc                                >  1:4273478/1‑36 (MQ=255)
     aTGGACTGCAGCTGGCTGGTCATCTCTTCCATGccc                               <  1:2231750/36‑1 (MQ=255)
      tGGACTGCAGCTGGCTGGTCATCTCTTCTATGCCCg                              <  1:7056406/36‑1 (MQ=255)
      tGGACTGCAGCTGGCTGGTCATCTCTTCCATGCCCg                              >  1:6501996/1‑36 (MQ=255)
       ggACTGCAGCTGGCTGGTCATCTCTTCCATGCCCgg                             <  1:3698191/36‑1 (MQ=255)
        gACTGCAGCTGGCTGGTCATCTCTTCCATGCCCgga                            <  1:3033757/36‑1 (MQ=255)
           tGCAGCTGGCTGGTCATCTCTTCCATGCCCggagga                         <  1:6817440/36‑1 (MQ=255)
           tGCAGCTGGCTGGTCATCTCTTCCATGCCCggagga                         >  1:7345954/1‑36 (MQ=255)
           tGCAGCTGGCTGGTCATCTCTTCCATGCCCggagga                         >  1:5073226/1‑36 (MQ=255)
            gCAGCTGGCTGGTCATCTCTTCCATGCCCGgaggag                        >  1:6399846/1‑36 (MQ=255)
            gCAGCTGGCTGGTCATCTCTTCCATGCCCGgaggag                        >  1:5286033/1‑36 (MQ=255)
             cAGCTGGCTGGTCATCTCTTCTATGCCCGGAGGAGc                       <  1:1488665/36‑1 (MQ=255)
             cAGCTGGCTGGTCATCTCTTCCATGCCCGGAGGAGc                       <  1:3789132/36‑1 (MQ=255)
              aGCTGGCTGGTCATCTCTTCCATGCCCGGAGGAGcc                      >  1:8569/1‑36 (MQ=255)
               gctggctgGTCATCTCTTCCATGCCCGGAGGAGCCa                     <  1:1891228/36‑1 (MQ=255)
                ctggctggTCATCTCTTCCATGCCCGGAGGAGCCAt                    >  1:4252343/1‑36 (MQ=255)
                ctggctggTCATCTCTTCCATGCCCGGAGGAGCCAt                    <  1:5012123/36‑1 (MQ=255)
                  ggctggTCATCTCTTCCATGCCCGGAGGAGCCATaa                  <  1:3753787/36‑1 (MQ=255)
                   gctgtTCATCTCTTCTAGGCCCGGAGGAGCCATAAt                 <  1:5259283/36‑1 (MQ=25)
                   gctggTCATCTCTTCTATGCCCGGAGGAGCCATAAt                 >  1:6481932/1‑36 (MQ=255)
                   gctggTCATCTCTTCCATGCCCGGAGGAGCCATAAt                 >  1:5767164/1‑36 (MQ=255)
                     tggTCATCTCTTCTATGCCCGGAGGAGCCATAAttt               >  1:1668528/1‑36 (MQ=255)
                      ggTCATCTCTTCCATGCCCGGAGGAGCCATAATTTc              >  1:306833/1‑36 (MQ=255)
                      ggTCATCTCTTCCATGCCCGGAGGAGCCATAATCTc              >  1:3285841/1‑36 (MQ=255)
                          aTCTCTTCCATGCCCGGAGGAGCCATAATTTCAACg          >  1:558551/1‑36 (MQ=255)
                          aTCTCTTCCATGCCCGGAGGAGCCATAATTTCAACg          <  1:3176193/36‑1 (MQ=255)
                          aTCTCTTCCATGCCCGGAGGAGCCATAATTTCAACg          >  1:4765933/1‑36 (MQ=255)
                            ctctTCCATGCCCGGAGGAGCCATAATTTCAACGcc        >  1:4592858/1‑36 (MQ=255)
                             tctTCTAGGCCCGGAGGAGCCATAATTTCAACGccc       <  1:2933779/36‑1 (MQ=37)
                             tctTCCATGCCCGGAGGAGCCATAATTTCAACGccc       <  1:5833660/36‑1 (MQ=255)
                             tctTCCAGGCCCGGAGGAGCCATAATTTCAACGccc       <  1:3218053/36‑1 (MQ=255)
                               ttCTATGCCCGGAGGAGCCTTAATTTCAACGCCCAt     <  1:4443982/36‑1 (MQ=37)
                               ttCCATGCCCGGAGGATCCATAATTTCAACGCCCAt     <  1:7451115/36‑1 (MQ=255)
                               ttCCATGCCCGGAGGAGCCATAATTTCAACGCCCAt     <  1:2164982/36‑1 (MQ=255)
                               ttCCATGCCCGGAGGAGCCATAATTTCAACGCCCAt     >  1:6253872/1‑36 (MQ=255)
                                tCCATGCCCGGAGGTGCCATAATTTCAACGCCCAtt    >  1:7035561/1‑35 (MQ=255)
                                 ccATGCCCGGAGGAGCCATAATTTCAACGCCCATCg   <  1:6670597/36‑1 (MQ=255)
                                 ccATGCCCGGAGGAGCCATAATTTCAACGCCCATCg   >  1:6582365/1‑36 (MQ=255)
                                 ccATGCCCGGAGGAGCCATAATTTCAACGCCCATCg   <  1:2761946/36‑1 (MQ=255)
                                  cATGCCCGGAGGAGCCATAATTTCAACGCCCATCgg  >  1:5772656/1‑36 (MQ=255)
                                  |                                   
GGAACATGGACTGCAGCTGGCTGGTCATCTCTTCCATGCCCGGAGGAGCCATAATTTCAACGCCCATCGG  >  REL606/4100621‑4100690

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