Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 161,041 T→G N302H (AAC→CAC)  pcnB poly(A) polymerase I

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606161,0410T→G100.0% 186.0 59N302H (AAC→CAC) pcnBpoly(A) polymerase I
Reads supporting (aligned to +/- strand):  new base (26/33):  ref base (0/0):  total (26/33)

ACCAGAACATGGCGGCAAACAGGAACGCCGGGTTCACGCGCATATCGTTATGGATACGCGTATCGGTA  >  REL606/161008‑161075
                                 |                                  
aCCAGAACATGGCTGCAAAAAGGAATGCAGGGTGCa                                  >  1:1576226/1‑36 (MQ=12)
aCCAGAACATGGCGGCAAACAGGAACGCCGGGTGCa                                  <  1:88784/36‑1 (MQ=255)
 ccAGAACATGGCGGCAAACAGGAACGCCGGGTGCAc                                 <  1:2974933/36‑1 (MQ=255)
  cAGAACATGGCGGCAAACAGGAACGCCGGGTGCAcg                                >  1:4867991/1‑36 (MQ=255)
  cAGAACATGGCGGCAAACAGGAACGCCCGGCGCAcg                                >  1:3646268/1‑36 (MQ=21)
  cAGAACATGGCGGCAAACAAGAACGCCGGGTGCAcg                                >  1:1184111/1‑36 (MQ=25)
   aGAACCGGGCGGCAAACAGGAACGCCGGGTGCAcgc                               <  1:5560199/36‑1 (MQ=21)
   aGAACATGGCGGCAAACAGGAACGCCGGGTGCAcgc                               >  1:707186/1‑36 (MQ=255)
   aGAACATGGCGGCAAACAGGAACGCCGGGTGCAcgc                               <  1:4943435/36‑1 (MQ=255)
    gAACATGGCGGCAAACAGGAACGCCGGGTGCAcgcg                              <  1:6503247/36‑1 (MQ=255)
    gAACATGGCGGCAAACAGGAACGCCGGGTGCAcgcg                              <  1:2602443/36‑1 (MQ=255)
    gAACATCGCGGCAAACAGGAACGCCGGCTGCAcgcg                              >  1:4292397/1‑36 (MQ=18)
     aaCATGGCGGCAAACAGGAACGCCGGGTGCACgcgc                             <  1:1413358/36‑1 (MQ=255)
      aCATGGCGGCAAACAGGAACGCCGGGTGCACGCGCa                            >  1:5202255/1‑36 (MQ=255)
      aCATGGCGGCAAACAGGAACGCCGGGTGCACGCGCa                            >  1:1755262/1‑36 (MQ=255)
       cATGGCGGCAAACAGGAACGCCGGGTGCACGCGCat                           >  1:5362701/1‑36 (MQ=255)
       cATGGCGGCAAACAGGAACGCCGGGGGCACgcgcct                           >  1:7175196/1‑34 (MQ=21)
          ggcggcAAACAGGAACGCCGGGTGGACGCGCATATc                        >  1:4337762/1‑36 (MQ=25)
          ggcggcAAACAGGAACGCCGGGTGCACGCGCATATc                        <  1:7424582/36‑1 (MQ=255)
           gcggcAGACAGGAGCGGCGGGTGCACGCGCATATCg                       <  1:1129391/36‑1 (MQ=18)
           gcggcAAACAGGAACGCCGGGTGCACGCGCATATCg                       >  1:6220732/1‑36 (MQ=255)
           gcggcAAACAGGAACGCCGGGTGCACGCGCATATCg                       >  1:4081437/1‑36 (MQ=255)
           acggcAAACAGGAACGCCGGGTGCACGCGCATATCg                       <  1:1401653/35‑1 (MQ=25)
            cggcAAGCAGGAACGCCGGGTGCACGCGCATATCGt                      <  1:3361544/36‑1 (MQ=25)
            cggcAAACAGGAACGCCGGGTGCACGCGCATATCGt                      <  1:3518270/36‑1 (MQ=255)
            cggcAAACAGGAACGCCGGGTGCACGCGCACATCGt                      >  1:3851736/1‑36 (MQ=25)
             ggcAAACAGGAACGCCGGGTGCACGCGCATATCGtt                     <  1:3484875/36‑1 (MQ=255)
             ggcAAACAGGAACGCCGGGTGCACGCGCATATCGtt                     <  1:4597152/36‑1 (MQ=255)
             ggcAAACAGGAACGCCGGGTGCACGCGCATATCGtt                     >  1:5574343/1‑36 (MQ=255)
             ggcAAACAGGAACGCCGGGTGCAAGCGCACATCGtt                     >  1:6679867/1‑36 (MQ=21)
              gcAAACAGGAACGCCGGGTGCACGCGCATATCGTTa                    >  1:5993522/1‑36 (MQ=255)
              gcAAACAGGAACGCCGGGTGCACGCGCATATCGTTa                    <  1:4798805/36‑1 (MQ=255)
               cAAACAGGAACGCCGGGTGCACGCGCATATCGTTAt                   <  1:3812773/36‑1 (MQ=255)
               cAAACAGGAACGCCGGGTGCACGCGCATATCGTTAt                   >  1:1547665/1‑36 (MQ=255)
               cAAACAGGAACGCCGGGTGCACGCGCATATCGTTAt                   <  1:2381744/36‑1 (MQ=255)
                aaaCAGGAACGCCGGGTGCACGCGCATATCGTTATg                  >  1:7214297/1‑36 (MQ=255)
                 aaCAGGAACGCCGGGTGCACGCGCATATCGTTATgg                 <  1:3973592/36‑1 (MQ=255)
                 aaCAGGAACGCCGGGTGCACGCGCATATCGTTATgg                 <  1:414799/36‑1 (MQ=255)
                 aaCAGGAACGCCGGGTGCACGCGCATATCGTTACgg                 >  1:6608063/1‑36 (MQ=25)
                 aaCAGCAACGCCGGGTGCACGCGCATATCGTTAttg                 <  1:6526862/36‑3 (MQ=25)
                  aCAGGAACGCCGGGTGCACGCGCATATCGTTATGGa                <  1:1114715/36‑1 (MQ=255)
                  aCAGGAACGCCGGGTGCACGCGCATATCGTTATGGa                >  1:5019627/1‑36 (MQ=255)
                  aCAGGAACGCCGGGTGCACGCGCATATCGTTATGGa                <  1:3000300/36‑1 (MQ=255)
                  aCAGGAACGCCGGGTGCACGCGCACATCGTTATGGa                <  1:6177465/36‑1 (MQ=25)
                   cAGGAGCGCCGGGTGCACGCGCATATCGTTATGGAt               <  1:2940633/36‑1 (MQ=25)
                   cAGGAACGCCGGGTGCACGCGCATATCGTTATGGAt               >  1:5796108/1‑36 (MQ=255)
                   cAGGAACGCCGGGTGCACGCGCATATCGTTATGGAt               <  1:7130508/36‑1 (MQ=255)
                     ggAACGCCGGGTGCACGCGCATATCGTTATGGATAc             <  1:1447968/36‑1 (MQ=255)
                       aaCGCCGGGTGCACGCGCATATCGTTATGGATAcgc           <  1:2183580/36‑1 (MQ=255)
                       aaCGCCGGGTGCACGCGCATATCGTTATGGATAccc           >  1:1575721/1‑34 (MQ=37)
                        aggCCGGGTGCACGCGCATATCGTTATGGATAcgcg          <  1:224693/34‑1 (MQ=25)
                        aCGCCGGGTGCACGCGCATATCGTTATGGATAggcg          >  1:4329936/1‑36 (MQ=25)
                        aCGCCGGGTGCACGCGCATATCGTTATGGATAcgcg          >  1:2407913/1‑36 (MQ=255)
                         cGCCGGGTGCACGCGCATATCGTTATGGATACGCGt         <  1:4692262/36‑1 (MQ=255)
                          gCCGGGTGCACGCGCATATCGTTATGGATACGCGTa        <  1:5014364/36‑1 (MQ=255)
                          gCCGGGTGCACGCGCATATCGTTATGGATACGCGTa        <  1:3285470/36‑1 (MQ=255)
                           ccGGGTGCACGCGCATATCGTTATGGATACGCGTAt       <  1:6546850/36‑1 (MQ=255)
                              ggTGCACGCGCATATCGTTATGGATACGCGTATCgg    >  1:2974942/1‑36 (MQ=255)
                              ggTGCACGCGCATATCGTTATGGATACGCGTATCgg    <  1:2307038/36‑1 (MQ=255)
                                atCACTCGCATATCGTTATGGATGCGCGTATCGGTa  <  1:7257490/35‑1 (MQ=21)
                                 |                                  
ACCAGAACATGGCGGCAAACAGGAACGCCGGGTTCACGCGCATATCGTTATGGATACGCGTATCGGTA  >  REL606/161008‑161075

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 16 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG