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% 196.0 61904pcnBpoly(A) polymerase I
Reads supporting (aligned to +/- strand):  new base (28/33):  ref base (0/0):  total (28/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                                >  2:1128781/1‑36 (MQ=255)
  cAGAACATGGCGGCAAACAGGAACGCCCGGCGCAcg                                >  1:3646268/1‑36 (MQ=25)
  cAGAACATGGCGGCAAACAAGAACGCCGGGTGCAcg                                >  1:1184111/1‑36 (MQ=37)
   aGAACCGGGCGGCAAACAGGAACGCCGGGTGCAcgc                               <  2:1820989/36‑1 (MQ=25)
   aGAACATGGCGGCAAACAGGAACGCCGGGTGCAcgc                               <  2:1204225/36‑1 (MQ=255)
   aGAACATGGCGGCAAACAGGAACGCCGGGTGCAcgc                               >  1:707186/1‑36 (MQ=255)
    gAACATGGCGGCAAACAGGAACGCCGGGTGCAcgcg                              <  2:2764037/36‑1 (MQ=255)
    gAACATGGCGGCAAACAGGAACGCCGGGTGCAcgcg                              <  1:2602443/36‑1 (MQ=255)
    gAACATCGCGGCAAACAGGAACGCCGGCTGCAcgcg                              >  2:553187/1‑36 (MQ=21)
     aaCATGGCGGCAAACAGGAACGCCGGGTGCACgcgc                             <  1:1413358/36‑1 (MQ=255)
      aCATGGCGGCAAACAGGAACGCCGGGTGCACGCGCa                            >  2:1463045/1‑36 (MQ=255)
      aCATGGCGGCAAACAGGAACGCCGGGTGCACGCGCa                            >  1:1755262/1‑36 (MQ=255)
       cATGGCGGCAAACAGGAACGCCGGGTGCACGCGCat                           >  2:1623491/1‑36 (MQ=255)
       cATGGCGGCAAACAGGAACGCCGGGGGCACgcgcct                           >  2:3435986/1‑34 (MQ=25)
          ggcggcAAACAGGAACGCCGGGTGGACGCGCATATc                        >  2:598552/1‑36 (MQ=37)
          ggcggcAAACAGGAACGCCGGGTGCACGCGCATATc                        <  2:3685372/36‑1 (MQ=255)
           gcggcAGACAGGAGCGGCGGGTGCACGCGCATATCg                       <  1:1129391/36‑1 (MQ=19)
           gcggcAAACAGGAACGCCGGGTGCACGCGCATATCg                       >  2:342227/1‑36 (MQ=255)
           gcggcAAACAGGAACGCCGGGTGCACGCGCATATCg                       >  2:2481522/1‑36 (MQ=255)
           acggcAAACAGGAACGCCGGGTGCACGCGCATATCg                       <  1:1401653/35‑1 (MQ=255)
            cggcAAGCAGGAACGCCGGGTGCACGCGCATATCGt                      <  1:3361544/36‑1 (MQ=38)
            cggcAAACAGGAACGCCGGGTGCACGCGCATATCGt                      <  1:3518270/36‑1 (MQ=255)
            cggcAAACAGGAACGCCGGGTGCACGCGCACATCGt                      >  2:112526/1‑36 (MQ=37)
             ggcAAACAGGAACGCCGGGTGCACGCGCATATCGtt                     <  1:3484875/36‑1 (MQ=255)
             ggcAAACAGGAACGCCGGGTGCACGCGCATATCGtt                     >  2:1835133/1‑36 (MQ=255)
             ggcAAACAGGAACGCCGGGTGCACGCGCATATCGtt                     <  2:857942/36‑1 (MQ=255)
             ggcAAACAGGAACGCCGGGTGCAAGCGCACATCGtt                     >  2:2940657/1‑36 (MQ=25)
              gcAAACAGGAACGCCGGGTGCACGCGCATATCGTTa                    >  2:2254312/1‑36 (MQ=255)
              gcAAACAGGAACGCCGGGTGCACGCGCATATCGTTa                    <  2:1059595/36‑1 (MQ=255)
               cAAACAGGAACGCCGGGTGCACGCGCATATCGTTAt                   <  1:2381744/36‑1 (MQ=255)
               cAAACAGGAACGCCGGGTGCACGCGCATATCGTTAt                   >  1:1547665/1‑36 (MQ=255)
               cAAACAGGAACGCCGGGTGCACGCGCATATCGTTAt                   <  2:73563/36‑1 (MQ=255)
                aaaCAGGAACGCCGGGTGCACGCGCATATCGTTATg                  >  2:3475087/1‑36 (MQ=255)
                 aaCAGGAACGCCGGGTGCACGCGCATATCGTTATgg                 <  1:414799/36‑1 (MQ=255)
                 aaCAGGAACGCCGGGTGCACGCGCATATCGTTATgg                 <  2:234382/36‑1 (MQ=255)
                 aaCAGGAACGCCGGGTGCACGCGCATATCGTTACgg                 >  2:2868853/1‑36 (MQ=38)
                 aaCAGCAACGCCGGGTGCACGCGCATATCGTTAttg                 <  2:2787652/36‑3 (MQ=37)
                  aCAGGAACGCCGGGTGCACGCGCATATCGTTATGGa                <  1:1114715/36‑1 (MQ=255)
                  aCAGGAACGCCGGGTGCACGCGCATATCGTTATGGa                >  2:1280417/1‑36 (MQ=255)
                  aCAGGAACGCCGGGTGCACGCGCATATCGTTATGGa                <  1:3000300/36‑1 (MQ=255)
                  aCAGGAACGCCGGGTGCACGCGCACATCGTTATGGa                <  2:2438255/36‑1 (MQ=37)
                   cAGGAGCGCCGGGTGCACGCGCATATCGTTATGGAt               <  1:2940633/36‑1 (MQ=38)
                   cAGGAACGCCGGGTGCACGCGCATATCGTTATGGAt               <  2:3391298/36‑1 (MQ=255)
                   cAGGAACGCCGGGTGCACGCGCATATCGTTATGGAt               >  2:2056898/1‑36 (MQ=255)
                     ggAACGCCGGGTGCACGCGCATATCGTTATGGATAc             <  1:1447968/36‑1 (MQ=255)
                      gAACGCCGGGTGCACGCGCTTATCGCTATGGATaag            >  1:3641925/1‑34 (MQ=21)
                       aaCGCCGGGTGCACGCGCATATCGTTATGGATAcgc           <  1:2183580/36‑1 (MQ=255)
                       aaCGCCGGGTGCACGCGCATATCGTTATGGATAccc           >  1:1575721/1‑34 (MQ=38)
                        aggCCGGGTGCACGCGCATATCGTTATGGATAcgcg          <  1:224693/34‑1 (MQ=37)
                        aCGCCGGGTGCACGCGCATATCGTTATGGATAggcg          >  2:590726/1‑36 (MQ=37)
                        aCGCCGGGTGCACGCGCATATCGTTATGGATAcgcg          >  1:2407913/1‑36 (MQ=255)
                         cGCCGGGTGCACGCGCATATCGTTATGGATACGCGt         <  2:953052/36‑1 (MQ=255)
                          gCCGGGTGCACGTGCCTCTCGTTATGGATAGGCGtg        >  2:3295466/1‑35 (MQ=12)
                          gCCGGGTGCACGCGCATATCGTTATGGATACGCGTa        <  1:3285470/36‑1 (MQ=255)
                          gCCGGGTGCACGCGCATATCGTTATGGATACGCGTa        <  2:1275154/36‑1 (MQ=255)
                           ccGGGTGCACGCGCATATCGTTATGGATACGCGTAt       <  2:2807640/36‑1 (MQ=255)
                              ggTGCACGCGCATATCGTTATGGATACGCGTATCgg    >  1:2974942/1‑36 (MQ=255)
                              ggTGCACGCGCATATCGTTATGGATACGCGTATCgg    <  1:2307038/36‑1 (MQ=255)
                                atCACTCGCATATCGTTATGGATGCGCGTATCGGTa  <  2:3518280/35‑1 (MQ=37)
                                 |                                  
ACCAGAACATGGCGGCAAACAGGAACGCCGGGTTCACGCGCATATCGTTATGGATACGCGTATCGGTA  >  REL606/161008‑161075

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 15 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG