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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,616,3960G→A100.0% 193.8 54868nadRnicotinamide‑nucleotide adenylyltransferase
Reads supporting (aligned to +/- strand):  new base (22/32):  ref base (0/0):  total (22/32)

TTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAA  >  REL606/4616363‑4616427
                                 |                               
ttGCAGTATTCCGATTACTATAAAATCGCGCTGAGc                               <  1:106687/36‑1 (MQ=37)
ttGCAGTATTCCGATTACGATAAAATCGCGCTGAGc                               >  1:5189374/1‑36 (MQ=255)
ttGCAGTATTCCGATTACGATAAAATCGCGCTGAGc                               >  1:1539276/1‑36 (MQ=255)
 tGCAGTATTCCGATTACGATAAAATCGCGCTGAGcc                              >  1:1292080/1‑36 (MQ=255)
 tGCAGTATTCCGATTACGATAAAATCGCGCTGAGcc                              >  1:152183/1‑36 (MQ=255)
 tGCAGTATTCCGATTACGATAAAATCGCGCTGAGcc                              >  1:3776188/1‑36 (MQ=255)
  gCAGTATTCCGATTACGATAAAATCGCGCTGAGCCa                             >  1:5313585/1‑36 (MQ=255)
  gCAGTATTCCGATTACGATAAAATCGCGCTGAGCCa                             >  1:4782626/1‑36 (MQ=255)
  gCAGTATTCCGATTACGATAAAATCGCGCTGAGCCa                             >  1:4220118/1‑36 (MQ=255)
   cAGTATTCCGATTACGATAAAATCGCGCTGAGCCAc                            <  1:4994876/36‑1 (MQ=255)
    aGTATTCCGATTACGATAAAATCGCGCTGAGCCACg                           >  1:801873/1‑36 (MQ=255)
    aGTATTCCGATTACGATAAAATCGCGCTGAGCCACg                           >  1:450248/1‑36 (MQ=255)
    aGTATTCCGATTACGATAAAATCGCGCTGAGCCACg                           >  1:4363337/1‑36 (MQ=255)
      tATTCCGATTACGATAAAATCGCGCTGAGCCACGca                         <  1:5302156/36‑1 (MQ=255)
      tATTCCGATTACGATAAAATCGCGCTGAGCCACGca                         <  1:3148110/36‑1 (MQ=255)
       aTTCCGATTACGATAAAATCGCGCTGAGCCACGcac                        >  1:5806638/1‑36 (MQ=255)
       aTTCCGATTACGATAAAATCGCGCTGAGCCACGcac                        <  1:1401400/36‑1 (MQ=255)
       aTTCCGATTACGATAAAATCGCGCTGAGCCACGcac                        <  1:4377705/36‑1 (MQ=255)
       aTTCCGATTACGATAAAATCGCGCTGAGCCACGcac                        >  1:1947978/1‑36 (MQ=255)
        ttCCGATTACGATAAAATCGCGGTGATCCACGccca                       <  1:5457904/36‑1 (MQ=16)
        ttCCGATTACGATAAAATCGCGCTGAGCCACGcccc                       <  1:484273/36‑4 (MQ=37)
        ttCCGATTACGATAAAATCGCGCTGAGCCACGcaca                       <  1:4253906/36‑1 (MQ=255)
         tCCGATTACGATAAAATCGCGCTGAGCCACGCACaa                      <  1:3531651/36‑1 (MQ=255)
         tCCGATTACGATAAAATCGCGCTGAGCCACGCACaa                      <  1:4850317/36‑1 (MQ=255)
          ccGATTACGATAAAATCGCGCTGAGCCACGCACAAt                     <  1:4673823/36‑1 (MQ=255)
           cGATTACGATAAAATCGCGCTGAGCCACGCACAATa                    >  1:4027944/1‑36 (MQ=255)
            gATTACGATAAAATCGCGCTGAGCCACGCACAATAc                   <  1:381360/36‑1 (MQ=255)
            gATTACGATAAAATCGCGCTGAGCCACGCACAATAc                   >  1:4353968/1‑36 (MQ=255)
             aTTACGATAAAATCGCGCTGAGCCACGCACAATACa                  <  1:816508/36‑1 (MQ=255)
              ttACGATAAAATCGCGCTGAGCCACGCACAATACAt                 <  1:4510093/36‑1 (MQ=255)
              ttACGATAAAATCGCGCTGAGCCACGCACAATACAt                 <  1:396849/36‑1 (MQ=255)
               tACGATAAAATCGCGCTGAGCCACGCACAATACAtt                <  1:4542917/36‑1 (MQ=255)
               tACGATAAAATCGCGCTGAAACACGCACAATACAtt                >  1:4954154/1‑36 (MQ=25)
                aCGATAAAATCGCGCGGAGCCACGCACAATACATTg               >  1:1382376/1‑36 (MQ=37)
                 cGATAAAATCGCGCTGAGCCACGCACAATACATTGa              <  1:1980341/36‑1 (MQ=255)
                  gATAAAATCGCGCTGAGCCACGCACAATACATTGAt             <  1:2329306/36‑1 (MQ=255)
                  gATAAAATCGCGCTGAGCCACGCACAATACATTGAt             <  1:585469/36‑1 (MQ=255)
                     taaaACGCGATGAGCCCAGCACAATGCATTGAtttt          <  1:3506596/35‑1 (MQ=9)
                      aaaTCGCGCTGAGCCACGCACAATACATTGATTTTg         <  1:2283846/36‑1 (MQ=255)
                      aaaTCGCGCTGAGCCACGCACAATACATTGATTTTg         >  1:3927848/1‑36 (MQ=255)
                      aaaTCGCGCTGAGCCACGCACAATACATTGATTTTg         >  1:1505349/1‑36 (MQ=255)
                      aaaTCGCGCTGAGCCACGCACAATACATTGATTTTg         <  1:2612462/36‑1 (MQ=255)
                       aaTCGCGCTGAGCCACGCACACTACATTGATTTTGc        >  1:5863696/1‑36 (MQ=38)
                       aaTCGCGCTGAGCCACGCACAATACATTGATTTTGc        <  1:3028141/36‑1 (MQ=255)
                       aaTCGCGCTGAGCCACGCACAATACATTGATTTTGc        <  1:720417/36‑1 (MQ=255)
                       aaTCGCGCTGAGCCAAGCACAATACATTGATTTTGc        <  1:3090549/36‑1 (MQ=37)
                        aTCGCGCTGAGCCACGCACAATACATTGATTTTGCa       >  1:450474/1‑36 (MQ=255)
                        aTCGCGCTGAGCCACGCACAATACATTGATTTTGCa       <  1:1605359/36‑1 (MQ=255)
                         tCGCGCTGAGCCACGCACAATACATTGATTTTGCAg      <  1:537249/36‑1 (MQ=255)
                          cgcgCTGAGCCACGCACAATACATTGATTTTGCAGt     <  1:1689688/36‑1 (MQ=255)
                          cgcgCTGAGCCACGCACAATACATTGATTTTGCAGt     <  1:332161/36‑1 (MQ=255)
                           gcgcTGAGCCACGCACAATACATTGATTTTGCAGTg    <  1:5896592/36‑1 (MQ=255)
                             gcTGAGCCACGCACAATACATTGATTTTGCAGTGaa  >  1:2208673/1‑36 (MQ=255)
                             gcTGAGCCACGCACAATACATTGATTTTGCAGTGaa  <  1:3132901/36‑1 (MQ=255)
                                 |                               
TTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAA  >  REL606/4616363‑4616427

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