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
*	REL606	161,041	0	T→G	100.0%	186.0	59	N302H (AAC→CAC)	pcnB	poly(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