Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,248,957 A→T 100% 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→T100.0% 223.3 652292infBtranslation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  new base (27/36):  ref base (1/0):  total (29/36)

AATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGC  >  REL606/3248924‑3248989
                                 |                                
aaTCAGGTTATAGATGACGGAGTAGTAACGCAGTTc                                >  1:5545904/1‑36 (MQ=255)
aaTCAGGTTATAGATGACGGAGTAGTAAAGCAGTTc                                >  1:735686/1‑36 (MQ=37)
 aTCAGGTTATAGATGACGGAGTAGTAACGCAGTTcc                               <  1:5095806/36‑1 (MQ=255)
 aTCAGGTTATAGATGACGGAGTAGTAACGCAGTTcc                               <  1:733074/36‑1 (MQ=255)
  tCAGGTTATAGATGACGGAGTAGTAACGCAGTTCCa                              <  1:3322892/36‑1 (MQ=255)
  tCAGGTTATAGATGACGGAGTAGTAACGCAGTTCCa                              <  1:5835985/36‑1 (MQ=255)
  tCAGGTTATAGATGACGGAGTAGTAACGCAGATCCa                              >  1:5893022/1‑36 (MQ=255)
    aGGTTATAGATGACGGAGTAGTAACGCAGTTCCAgg                            <  1:3660613/36‑1 (MQ=255)
    aGGTTATAGATGACGGAGTAGTAACGCAGTTCCAgg                            >  1:5196207/1‑36 (MQ=255)
     ggTTATAGATGACGGAGTAGTAACGCAGTTCCAGGc                           >  1:4397549/1‑36 (MQ=255)
      gTTATAGATGACGGGGTAGTAACGCAGTTCCAGGCt                          <  1:2246302/36‑1 (MQ=37)
      gTTATAGATGACGGAGTAGTAACGCAGTTCCAGGCt                          <  1:3254785/36‑1 (MQ=255)
       ttATAGATGACGGAGTAGTAACGCAGTTCCAGGCtt                         <  1:691411/36‑1 (MQ=255)
       ttATAGATGACGGAGTAGTAACGCAGTTCCAGGCtt                         <  1:3750236/36‑1 (MQ=255)
        tctaGATGACGGAGTAGTACCGCAGTTCCAGGCttt                        <  1:1158810/34‑1 (MQ=25)
        tataGATGACGGAGTAGTAACGCAGTTCCAGGCttt                        <  1:1581948/36‑1 (MQ=255)
         atcgATGACGGAGTAGTAACGCAGTTCCAGGCTTTc                       <  1:3858938/33‑1 (MQ=37)
         ataGATGCCGGAGTAGTAACGCAGTTCCAGGCTTTc                       <  1:781179/36‑1 (MQ=37)
         ataGATGACGGAGTAGTAACGCAGTTCCAGGCTTTc                       >  1:3443534/1‑36 (MQ=255)
         ataGATGACGGAGTAGTAACGCAGTTCCAGGCTTTc                       >  1:4647257/1‑36 (MQ=255)
         ataGATGACGGAGTAGTAACGCAGTTCCAGGCTTTc                       <  1:5641324/36‑1 (MQ=255)
         ataGATGACGGAGTAGTAACGCAGTTCCAGGCTTTc                       <  1:830996/36‑1 (MQ=255)
          taGATGACGTAGTAGTAACGCAGTTCCAGGCTTTcc                      <  1:720759/36‑1 (MQ=37)
          taGATGACGGAGTAGTAACGCAGTTCCAGGCTTTcc                      >  1:638054/1‑36 (MQ=255)
          taGATGACGGAGTAGTAACGCAGTTCCAGGCTTTcc                      >  1:3390141/1‑36 (MQ=255)
           aGATGACGGAGTAGTAACGCATTTCCAGGCTTTCCg                     <  1:3831029/36‑1 (MQ=37)
           aGATGACGGAGTAGTAACGCAGTTCCAGGCTTCCCg                     >  1:2460299/1‑36 (MQ=37)
           aGATGACGGAGTAGTAACGCAGGTCCAGGCTTTCCg                     >  1:5276372/1‑36 (MQ=255)
            gATGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGc                    >  1:2536323/1‑36 (MQ=255)
            gATGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGc                    <  1:1654655/36‑1 (MQ=255)
            gATGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGc                    <  1:108632/36‑1 (MQ=255)
            gATGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGc                    <  1:3989665/36‑1 (MQ=255)
             aTGACGGAGTAGTAACGCAGTTCCAGGCTTTCCgat                   >  1:1468998/1‑34 (MQ=37)
             aTGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCt                   >  1:5067054/1‑36 (MQ=255)
             aTGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCt                   <  1:397073/36‑1 (MQ=255)
              tGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCtt                  <  1:1080416/36‑1 (MQ=255)
               gACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTc                 >  1:3923990/1‑36 (MQ=255)
               gACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTc                 >  1:10500/1‑36 (MQ=255)
               gACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTc                 <  1:3414423/36‑1 (MQ=255)
                aCGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCa                <  1:4378049/36‑1 (MQ=255)
                aCGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCa                <  1:1657613/36‑1 (MQ=255)
                 cGGAGTAGTAACGCATTTCCAGGCTTTCCGCTTCaa               <  1:995849/36‑1 (MQ=37)
                 cGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCaa               >  1:3343180/1‑36 (MQ=255)
                  ggAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCAAt              >  1:2832657/1‑36 (MQ=255)
                  ggAGGAGTAACGCAGTTCCAGGCTTTCCGCTTAAAt              >  1:5607744/1‑36 (MQ=25)
                   gAGTAGTAACGCAGTTCCAGGCTTTCCGCTTTAATc             >  1:1754079/1‑36 (MQ=37)
                   gAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCAATc             <  1:3920070/36‑1 (MQ=255)
                   gAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCAATc             >  1:841980/1‑36 (MQ=255)
                    agtagtAACGCAGTTCCAGGCTTTCCGCTTCAATCa            <  1:5771882/36‑1 (MQ=255)
                    agtagtAACGCAGTTCCAGGCTTTCCGCTTCAATCa            <  1:5347999/36‑1 (MQ=255)
                    agtagtAACGCAGTTCCAGGCTTTCCGCTTCAATCa            <  1:869192/36‑1 (MQ=255)
                     gtagtaACGCAGTTCCAGGCTTTCCGCTTCAATCAc           <  1:5585646/36‑1 (MQ=255)
                     gtagtaACGCAGTTCCAGGCTTTCCGCTTCAATCAc           >  1:274362/1‑36 (MQ=255)
                      tagtaACGCAGTTCCAGGCTTTCCGCTTCAATCACt          <  1:5303374/36‑1 (MQ=255)
                       agtaACGCAGTTCCAGGCTTTCCGCTTCAATCAcct         >  1:3202369/1‑34 (MQ=37)
                       agtaACGCAGTTCCAGGCTTTCCGCTTCAATCACtt         >  1:2808695/1‑36 (MQ=255)
                       agtaACGCAGTTCCAGGCTTTCCGCTTCAATCACtt         <  1:2386582/36‑1 (MQ=255)
                         taACGCAGTTCCAGGCTTTCCGCTTCAATCACTTTa       <  1:5657444/36‑1 (MQ=255)
                         taACGCAGTTCCAGGCTTTCCGCTTCAATCACTTTa       <  1:3693610/36‑1 (MQ=255)
                          aaCGCAGTTCCAGGCTTTCCGCTTCAATCACTTTAc      >  1:503536/1‑36 (MQ=255)
                          aaCGCAGTTCCAGGCTTTCCGCTTCAATCACTTTAc      <  1:3668056/36‑1 (MQ=255)
                           aCGCAGTTCCAGGCTTTCCGCTTCAATCACTTTACg     >  1:779870/1‑36 (MQ=255)
                             gCAGTTCCAGGCTTTCCGCTTTAATCACTTTACGTg   >  1:3154523/1‑36 (MQ=37)
                             gCAGTTCCAGGCTTTCCGCTTCAATCACTTTACGTg   >  1:1502121/1‑36 (MQ=255)
                              cAGTTCCAGGCTTTCCGATTCAATCACTTTACGTgt  >  1:4439152/1‑35 (MQ=37)
                                 |                                
AATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTGC  >  REL606/3248924‑3248989

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