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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606161,0410T→G90.6% 6.8 ‑0.5 ‑0.453904pcnBpoly(A) polymerase I
Reads supporting (aligned to +/- strand):  new base (28/20):  ref base (4/1):  total (32/21)
Fisher's exact test strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.94e-01

ACCAGAACATGGCGGCAAACAGGAACGCCGGGTTCACGCGCATATCGTTATGGATACGCGTATCG  >  REL606/161008‑161072
                                 |                               
aCCAGAACATGGCGGCAAACAGGAACGCCGGGTGCa                               >  1:8252990/1‑36 (MQ=255)
aCCAGAACATGGCGGCAAACAGGAACGCCGGGTGCa                               >  1:2263740/1‑36 (MQ=255)
 ccAGAACATGGCTGCAAACAGGAACGCCGGGTGCAc                              <  1:8275835/36‑1 (MQ=38)
 ccAGAACATGGCGGCAAACAGGAACGCCGGGTGCAc                              >  1:5927258/1‑36 (MQ=255)
  cAGAACATGGCGGCAAACAGGAACGCCGGGTGCAcg                             >  1:4864853/1‑36 (MQ=255)
  cAGAACATGGCGGCAAACAGGAACGCCGGGTGCAcg                             >  1:6200498/1‑36 (MQ=255)
   aGAACATGGCGGCAAACAGGAACGCCGGGTGCAcgc                            >  1:2349362/1‑36 (MQ=255)
   aGAACATGGCGGCAAACAGGAACGCCGGGTGCAcgc                            >  1:5618804/1‑36 (MQ=255)
    gAACATGGCGGCAAACAGGAACGTCGGGTGCAcgcg                           >  1:2641840/1‑36 (MQ=37)
    gAACATGGCGGCAAACAGGAACGCCGGGTGCAcgcg                           <  1:6521304/36‑1 (MQ=255)
    gAACATGGCGGCAAACAGGAACGCCGGGTGCAcgcg                           <  1:6196700/36‑1 (MQ=255)
    gAACATGGCGGCAAACAGGAACGCCGGGTGCAcgcg                           <  1:5683200/36‑1 (MQ=255)
      aCATGGCGGCAAACAGGAACGCCGGGTGCACGCGCa                         <  1:6202838/36‑1 (MQ=255)
      aCATGGCGGCAAACAGGAACGCCGGGTGCACGCGCa                         >  1:5440861/1‑36 (MQ=255)
       cATGGCGGCAAACAGGAACGCCGGGTGCACGCGCat                        <  1:417629/36‑1 (MQ=255)
       cATGGCGGCAAACAGGAACGCCGGGGGCACGCGCat                        >  1:5224936/1‑36 (MQ=37)
        aTGGCGGCAAACAGGAACGCCGGGTGCACGCGCata                       <  1:1626174/36‑1 (MQ=255)
        aTGGCGGCAAACAGGAACGCCGGGTGCAAGCGAata                       >  1:8413638/1‑36 (MQ=25)
          ggcggcAAACAGGAACGCCGGGTGCACGCGCATATc                     <  1:6331659/36‑1 (MQ=255)
          ggcggcAAACAGGAACGCCGGGTGCACGCGCATATc                     <  1:2293821/36‑1 (MQ=255)
          ggcggcAAACAGGAACGCCGGGTGCACGCGCATATc                     >  1:1020135/1‑36 (MQ=255)
           gcggcAAACAGGAACGCCGGGTTCACGCGCATATCg                    >  1:1352495/1‑36 (MQ=255)
           gcggcAAACAGGAACGCCGGGTGCACGCGCATATCg                    >  1:8360029/1‑36 (MQ=255)
            cggcAAACAGGAACGCCGGGTTCACGCGCATATCGt                   >  1:7793927/1‑36 (MQ=255)
            cggcAAACAGGAACGCCGGGTGCACGCGCATATCGt                   <  1:745444/36‑1 (MQ=255)
            cggcAAACAGGAACGCCGGGTGCACGCGCATATCGt                   >  1:4640352/1‑36 (MQ=255)
            cggcAAACAGGAACGCCGGGTGCACGCGCATATCGt                   >  1:3956664/1‑36 (MQ=255)
             ggcAAACAGGAACGCCGGGTGCACGCGCATATCGtt                  <  1:7653049/36‑1 (MQ=255)
               cAAACAGGAACGCCGGGTTCACGCGCATATCGTTAt                >  1:4647381/1‑36 (MQ=255)
               cAAACAGGAACGCCGGGTGCACGCGCATATCGTTAt                >  1:4295446/1‑36 (MQ=255)
                  aCAGGAACGCCGTGTGCACGCGCATATCGTTATGGa             <  1:4474563/36‑1 (MQ=38)
                  aCAGGAACGCCGGGTGCACGCGCATATCGTTATGGa             >  1:5967611/1‑36 (MQ=255)
                  aCAGGAACGCCGGGTGCACGCGCATATCGTTATGGa             <  1:6007690/36‑1 (MQ=255)
                  aCAGGAACGCCGGGTGCACGCGCATATCGTTATGGa             >  1:866956/1‑36 (MQ=255)
                  aCAGGAACGCCGGGTGCACGCGCATATCGTTATGGa             <  1:2587581/36‑1 (MQ=255)
                   cAGGAACGCCGGGTGCACGCGCATATCGTTATGGAt            >  1:7721583/1‑36 (MQ=255)
                   cAGGAACGCCGGGTGCACGCGCATATCGTTATGGAt            >  1:2406568/1‑36 (MQ=255)
                    aGGAACGCCGGGTGCACGCGCATATCGTTATGGATa           >  1:5273837/1‑36 (MQ=255)
                     ggAACGCCGGGTGCACGCGCATATCGTTATGGATAc          <  1:3504966/36‑1 (MQ=255)
                      gAACGCCGGGTGCACGCGCATATCGTTATGGATAcg         >  1:618252/1‑36 (MQ=255)
                       aaCGCCGGGTTCACGCGCATATCGTTATGGATAcgc        >  1:667577/1‑36 (MQ=255)
                       aaCGCCGGGTGCACGCGCATATCGTTATGGATAcgc        <  1:4888397/36‑1 (MQ=255)
                        aCGCCGGGTGCACGCGCATATCGTTATGGATAcgcg       >  1:925660/1‑36 (MQ=255)
                        aCGCCGGGTGCACGCGCATATCGTTATGGATAcgcg       >  1:1175817/1‑36 (MQ=255)
                        aCGCCGGGTGCACGCGCATATCGTTATGGATAcgcg       >  1:7403472/1‑36 (MQ=255)
                        aCGCCGGGTGCACGCGCATATCGTTATGGATAcgcg       >  1:18307/1‑36 (MQ=255)
                         cGCCGGGTGCACGCGCATATCGTTATGGATACGCGt      <  1:3377774/36‑1 (MQ=255)
                          gCCGGGTGCACGCGCATATCGTTATGGATACGCGTa     <  1:890901/36‑1 (MQ=255)
                           ccGGGTTCACGCGCATATCGTTATGGATACGCGTAt    <  1:5872341/36‑1 (MQ=255)
                            cGGGTGCACGCGCATATCGTTATGGATATGCGTATc   >  1:4823298/1‑36 (MQ=37)
                             gggTGCACGCGCATATCGTTATGGATACGCGTATCg  <  1:6714889/36‑1 (MQ=255)
                             gggTGCACGCGCATATCGTTATGGATACGCGTATCg  >  1:4888561/1‑36 (MQ=255)
                             gggTGCACGCGCATATCGTTATGGATACGCGTATCg  <  1:3339783/36‑1 (MQ=255)
                                 |                               
ACCAGAACATGGCGGCAAACAGGAACGCCGGGTTCACGCGCATATCGTTATGGATACGCGTATCG  >  REL606/161008‑161072

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG