Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,616,396 G→A 9.9% G290S (GGC→AGC)  nadR → nicotinamide‑nucleotide adenylyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,616,3960G→A 9.9% 7.2 0.0 ‑0.051868nadRnicotinamide‑nucleotide adenylyltransferase
Reads supporting (aligned to +/- strand):  new base (3/2):  ref base (25/20):  total (28/23)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.95e-01

CATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATA  >  REL606/4616361‑4616430
                                   |                                  
cATTGCAGTATTCCGATTACGATAAAATCGCGCTgg                                    <  1:6688738/36‑1 (MQ=255)
  ttGCAGTATTCCGATTACGATAAAATCGCGCTGGGc                                  <  1:4278100/36‑1 (MQ=255)
  ttGCAGTATTCCGATTACGATAAAATCGCGCTGGGc                                  >  1:4369115/1‑36 (MQ=255)
   tGCAGTATTCCGATTACGATAAAATAGCGCTGGGcc                                 >  1:1486014/1‑36 (MQ=255)
    gCAGTATTCCGATTACGATAAAATCGCGCTGGGCCa                                <  1:1648051/36‑1 (MQ=255)
     cAGTATTCCGATTACGATAAAATCGCGCTGGGCCAc                               <  1:5051677/36‑1 (MQ=255)
      aGTATTCCGATTACGATAAAATCGCGCTGGGCCACg                              >  1:4655648/1‑36 (MQ=255)
       gTATTCCGATTACGATAAAATCGCGCTGGGCCACGc                             <  1:6336030/36‑1 (MQ=255)
        tATTCCGATTACGCTAAAATCGCGCTGAGCCACGca                            <  1:7089274/36‑1 (MQ=25)
        tATTCCGATTACGATAAAATCGCGCTGAGCCACGca                            <  1:1419482/36‑1 (MQ=255)
         aTTCCGATTACGATAAAATCGCGCTGGGCCACGcac                           <  1:2108897/36‑1 (MQ=255)
         aTTCCGATTACGATAAAATCGCGCTGGGCCACGcac                           >  1:5980604/1‑36 (MQ=255)
         aTTCCGATTACGATAAAATCGCGCTGGGCCACGcac                           <  1:4867480/36‑1 (MQ=255)
         aTTCCGATTACGATAAAATCGCGCTGGGCCACGcac                           <  1:3454044/36‑1 (MQ=255)
          ttCCGATTACGATAAAATCGCGCTGGGCCACGcaca                          >  1:797262/1‑36 (MQ=255)
           tCCGATTACGATAAAATCGCGCTGGGCCACGCACaa                         <  1:4463556/36‑1 (MQ=255)
             cGATTACGATAAAATCGCGCTGGGCCACGCACAATa                       <  1:863036/36‑1 (MQ=255)
             cGATTACGATAAAATCGCGCTGGGCCACGCACAATa                       <  1:3541531/36‑1 (MQ=255)
              gATTACGATAAAATCGCGCTGGGCCACGCACAATAc                      >  1:7497035/1‑36 (MQ=255)
              gATTACGATAAAATCGCGCTGGGCCACGCACAATAc                      <  1:3031930/36‑1 (MQ=255)
              gATTACGATAAAATCGCGCTGGGCCACGCACAATAc                      <  1:549074/36‑1 (MQ=255)
               aTTACGATAAAATCGCGCTGGGCCACGCACAATACa                     >  1:6293868/1‑36 (MQ=255)
                 tACGATAAAATCGCGCTGGGCCACGCACAATACAtt                   >  1:7609552/1‑36 (MQ=255)
                  aCGATAAAATCGCGCTGGGCCACGCACAATACATTg                  >  1:4744292/1‑36 (MQ=255)
                   cGATAAAATCGCGCTGGGCCACGCACAATACATTGa                 <  1:6187468/36‑1 (MQ=255)
                    gATAAAATCGCGCTGGGCCACGCACAATACATTGAt                >  1:7809857/1‑36 (MQ=255)
                    gATAAAATCGCGCTGAGCCACGCACAATACATTGAt                >  1:4508332/1‑36 (MQ=255)
                    gATAAAATCGCGCGGGGAAATGCACAATACATTGAt                >  1:5017429/1‑36 (MQ=16)
                     aTAAAATCGCGCTGGGCCACGCACAATACATTGAtt               >  1:5001624/1‑36 (MQ=255)
                      tAAAATCGCGCTGGGCCACGCACAATACATTGAttt              <  1:6864459/36‑1 (MQ=255)
                      tAAAATCGCGCTGGGCCACGCACAATACATTGAttt              <  1:4580653/36‑1 (MQ=255)
                       aaaaTCGCGCTGGGCCACGCACAATACATTGAtttt             >  1:6053553/1‑36 (MQ=255)
                       aaaaTCGCGCTGGGCCACGCACAATACATTGAtttt             >  1:5833745/1‑36 (MQ=255)
                       aaaaTCGCGCTGGGCCACGCACAATACATTGAtttt             <  1:3840116/36‑1 (MQ=255)
                       aaaaTCGCGCTGGGCCACGCACAATACATTGAtttt             >  1:1909986/1‑36 (MQ=255)
                       aaaaTCGCGCTGGGCCAAGCACAATACATTGAtttt             >  1:7095555/1‑36 (MQ=255)
                         aaTCGCGCTGGGCCACGCACAATACATTGATTTTGc           >  1:3008891/1‑36 (MQ=255)
                           tCGCGCTGGGCCACGCACAATACATTGATTTTGCAg         <  1:3004346/36‑1 (MQ=255)
                            cgcgCTGGGCCACGCACAATACATTGATTTTGCAGt        <  1:6204465/36‑1 (MQ=255)
                             tcgcTGTGCCACGCCCAATACATTGATTTTGCAGTg       <  1:2356794/35‑1 (MQ=37)
                             gcgcTGGGCCACGCACAATGCATTGATTTTGCAGTg       >  1:2330088/1‑36 (MQ=255)
                             gcgcTGGGCCACGCACAATACATTGATTTTGCAGTg       >  1:5015567/1‑36 (MQ=255)
                             gcgcTGGGCCACGCACAATACATTGATTTTGCAGTg       >  1:7872994/1‑36 (MQ=255)
                             gcgcTGGGCCACGCACAATACATTGATTTTGCAGTg       <  1:3646965/36‑1 (MQ=255)
                               gcTGGGCCACGCACAATACATTGATTTTGCAGTGaa     >  1:3448612/1‑36 (MQ=255)
                               gcTGGGCCACGCACAATACATTGATTTTGCAGTGaa     <  1:2030570/36‑1 (MQ=255)
                                cTGAGCCACGCACAATACATTGATTTTGCAGTGaaa    >  1:3086259/1‑36 (MQ=255)
                                cTGAGCCACGCACAATACATTGATTTTGCAGTGaaa    >  1:3627210/1‑36 (MQ=255)
                                 tGGGCCACGCACAATACATTGATTTTGCAGTGAAat   >  1:7217221/1‑36 (MQ=255)
                                 tGGGCCACGCACAATACATTGATTTTGCAGTGAAat   >  1:7827371/1‑36 (MQ=255)
                                 tGGGCCACGCACAATACATTGATTTTGCAGTGAAat   >  1:3617284/1‑36 (MQ=255)
                                 tGGGCCACGCACAATACATTGATTTTGCAGTGAAat   >  1:1480542/1‑36 (MQ=255)
                                  gggCCACGCAAAATACATTGATTTTGCAGTGAAata  <  1:4625418/36‑1 (MQ=255)
                                  gggCCAAGCACAATACANNGNNTTTGCAGTGAAata  <  1:6994154/36‑1 (MQ=37)
                                   |                                  
CATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATA  >  REL606/4616361‑4616430

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG