Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,248,957 A→T 93.4% D764E (GAT→GAA infB ← translation initiation factor IF‑2

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,248,9570A→T93.4% 5.3 ‑0.5 ‑0.0612292infBtranslation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  new base (30/27):  ref base (3/1):  total (33/28)
Fisher's exact test strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.54e-01

CAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGC  >  REL606/3248923‑3248989
                                  |                                
cAATCAGGTTATAGATGACGGAGTAGTAACGCAGAt                                 >  1:4819749/1‑36 (MQ=255)
 aaTCAGGTTATAGATGACGGAGTAGTAACGCAGTTc                                >  1:1592010/1‑36 (MQ=255)
 aaTCAGGTTATAGATGACGGAGTAGTAACGCAGTTc                                >  1:6385111/1‑36 (MQ=255)
  aTCAGGTTATAGATGACGGAGTAGTAACGCAGTTcc                               <  1:6069279/36‑1 (MQ=255)
   tCAGGTTATAGATGACGGAGTAGTAACGCAGTTcaa                              >  1:814068/1‑34 (MQ=38)
   tCAGGTTATAGATGACGGAGTAGTAACGCAGTTCCa                              >  1:5622569/1‑36 (MQ=255)
    cAGGTTATAGATGACGGAGTAGTAACGCAGTTCCAg                             <  1:489058/36‑1 (MQ=255)
    cAGGTTATAGATGACGGAGTAGTAACGCAGTTCCAg                             <  1:8390786/36‑1 (MQ=255)
    cAGGTTATAGATGACGGAGTAGTAACGCAGATCCAg                             >  1:2217421/1‑36 (MQ=255)
      ggTTATAGATGACGGAGTAGTATCGCAGTTCCAGGc                           >  1:5800617/1‑36 (MQ=37)
      ggTTATAGATGACGGAGTAGTAACGCAGTTCCAGGc                           <  1:60049/36‑1 (MQ=255)
      ggTTATAGATGACGGAGTAGTAACGCAGTTCCAGGc                           <  1:1093973/36‑1 (MQ=255)
      ggTTATAGATGACGGAGTAGTAACGCAGTTCCAGGc                           >  1:2340402/1‑36 (MQ=255)
      ggTTATAGATGACGGAGTAGTAACGCAGTTCCAGGc                           >  1:2019034/1‑36 (MQ=255)
        ttATAGATGACGGAGTAGTAACGCAGTTCCAGGCtt                         >  1:2912757/1‑36 (MQ=255)
        ttATAGATGACGGAGTAGTAACGCAGTTCCAGGCtt                         >  1:3053068/1‑36 (MQ=255)
        ttATAGATGACGGAGTAGTAACGCAGTTCCAGGCtt                         >  1:375727/1‑36 (MQ=255)
         tataGATGACGGAGTAGTAACGCAGTTCCAGGCttt                        >  1:3828831/1‑36 (MQ=255)
          ataGATGACGGAGTAGTAACGCAGTTCCAGGCTTTc                       <  1:7344302/36‑1 (MQ=255)
           taGATGACGGAGTAGTAACGCAGTTCCAGGCTTTcc                      >  1:5756357/1‑36 (MQ=255)
           taGATGACGGAGTAGTAACGCAGTTCCAGGCTTTcc                      <  1:6236143/36‑1 (MQ=255)
            aGATGACGGAGTAGTAACGCAGTTCCAGGCTTTCCg                     >  1:5434361/1‑36 (MQ=255)
            aGATGACGGAGTAGTAACGCAGTTCCAGGCTTTCCg                     <  1:1975395/36‑1 (MQ=255)
             gATGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGc                    >  1:600454/1‑36 (MQ=255)
             gATGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGc                    >  1:1850846/1‑36 (MQ=255)
             gATGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGc                    <  1:1951682/36‑1 (MQ=255)
              aTGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCt                   <  1:4840192/36‑1 (MQ=255)
               tGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCtt                  <  1:4067702/36‑1 (MQ=255)
               tGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCtt                  >  1:6611201/1‑36 (MQ=255)
               tGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCtt                  <  1:6807180/36‑1 (MQ=255)
               tGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCtt                  <  1:2357009/36‑1 (MQ=255)
                 aCGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCa                >  1:914480/1‑36 (MQ=255)
                 aCGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCa                >  1:8080544/1‑36 (MQ=255)
                 aCGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCa                >  1:4920265/1‑36 (MQ=255)
                  cGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCaa               >  1:2423516/1‑36 (MQ=255)
                   ggAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCAAt              <  1:8324339/36‑1 (MQ=255)
                   ggAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCAAt              <  1:4957945/36‑1 (MQ=255)
                    gAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCAATc             >  1:7917982/1‑36 (MQ=255)
                    gAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCAATc             >  1:7016523/1‑36 (MQ=255)
                    gAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCAATc             >  1:1427826/1‑36 (MQ=255)
                     agtagtAACGCAGTTCCAGGCTTTCCGCTTCAATCa            >  1:7421828/1‑36 (MQ=255)
                       tagtaACGCAGTTCCAGGCTTTCCGCTTCAATCACt          <  1:2359648/36‑1 (MQ=255)
                         gtaACGCAGTTCCAGGCTTTCCGCTTCAATCATttt        >  1:3875464/1‑36 (MQ=37)
                         gtaACGCAGTTCCAGGCTTTCCGCTTCAATCACttt        <  1:250047/36‑1 (MQ=255)
                         gtaACGCAGTTCCAGGCTTTCCGCTTCAATCACttt        <  1:5450271/36‑1 (MQ=255)
                           aaCGCAGTTCCAGGCTTTCCGCTTCAATCACTTTAc      <  1:8364691/36‑1 (MQ=255)
                           aaCGCAGTTCCAGGCTTTCCGCTTCAATCACTTTAc      >  1:2694136/1‑36 (MQ=255)
                            aCGCAGTTCCAGGCTTTCCGCTTCAATCACTTTACg     <  1:5328027/36‑1 (MQ=255)
                            aCGCAGTTCCAGGCTTTCCGCTTCAATCACTTTACg     <  1:3589826/36‑1 (MQ=255)
                             cGCAGTTCCAGGCTTTCCGCTTCAATCACTTTACGt    <  1:7657515/36‑1 (MQ=255)
                             cGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGt    >  1:5404642/1‑36 (MQ=255)
                              gCAGTTCCAGGCTTTCCGCTTCAATCACTTTACGTg   >  1:2203265/1‑36 (MQ=255)
                              gCAGTTCCAGGCTTTCCGCTTCAATCACTTTACGTg   <  1:8128774/36‑1 (MQ=255)
                              gCAGTTCCAGGCTTTCCGCTTCAATCACTTTACGTg   <  1:717597/36‑1 (MQ=255)
                              gCAGTTCCAGGCTTTCCGCTTCAATCACTTTACGTg   >  1:5516408/1‑36 (MQ=255)
                               cAGTTCCAGGCTTTCCGCTTCAATCACTTTACGTGc  <  1:7665600/36‑1 (MQ=255)
                               cAGTTCCAGGCTTTCCGCTTCAATCACTTTACGTGc  >  1:6619723/1‑36 (MQ=255)
                               cAGTTCCAGGCTTTCCGCTTCAATCACTTTACGTGc  >  1:6609298/1‑36 (MQ=255)
                               cAGTTCCAGGCTTTCCGCTTCAATCACTTTACGTGc  <  1:4052658/36‑1 (MQ=255)
                               cAGTTCCAGGCTTTCCGCTTCAATCACTTTACGTGc  <  1:4390442/36‑1 (MQ=255)
                               cAGTTCCAGGCTTTCCGCTTCAATCACTTTACGTGc  <  1:8420600/36‑1 (MQ=255)
                                  |                                
CAATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGC  >  REL606/3248923‑3248989

Base quality scores:  ATCG  < 3 ≤  ATCG  < 22 ≤  ATCG  < 28 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG