Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,616,396 G→A 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% 170.7 55868nadRnicotinamide‑nucleotide adenylyltransferase
Reads supporting (aligned to +/- strand):  new base (34/19):  ref base (1/1):  total (35/20)

TTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAA  >  REL606/4616363‑4616428
                                 |                                
ttGCAGTATTCCGATTACGATAAAATCGCGCTGAGc                                >  2:313815/1‑36 (MQ=255)
  gCAGTATTCCGATTACGATAAAATCGCGCTGAGcaa                              >  2:3261763/1‑34 (MQ=37)
  gCAGTATTCCGATTACGATAAAATCGCGCTGAGCCa                              <  1:698032/36‑1 (MQ=255)
  gCAGTATTCCGATTACGATAAAATCGCGCTGAGCCa                              >  2:3262304/1‑36 (MQ=255)
  gCAGTATTCCGATTACGATAAAATCGCGCTGAGCCa                              >  1:1735257/1‑36 (MQ=255)
  gCAGTATTCCGATTACGATAAAATCGCGCTGAGCCa                              <  1:1728794/36‑1 (MQ=255)
    aGTATTCCGATTACGATAAAATCGCGCTGAGCCACg                            >  2:1331964/1‑36 (MQ=255)
     gTATTCCGATTACGATAAAATCGCGCTGAGCCACGc                           <  1:1561291/36‑1 (MQ=255)
     gTATTCCGATTACGATAAAATCGCGCTGAGCCACGc                           <  2:21550/36‑1 (MQ=255)
       aTTCCGATTACGATAAAATCGCGCTGAGCCACGcac                         >  2:2011664/1‑36 (MQ=255)
         tCCGATTACGATAAAATCGCGCTGAGCCACGCACaa                       <  1:1924497/36‑1 (MQ=255)
          ccGATTACGATAAAATCGCGCTGAGCCACGCACAAt                      <  1:1448935/36‑1 (MQ=255)
          ccGATTACGATAAAATCGCGCTGAGCCACCCACAAt                      >  2:667936/1‑36 (MQ=37)
           cGATTACGATAAAATCGCGCTGAGCCACGCACAATa                     <  1:2083970/36‑1 (MQ=255)
            gATTACGATAAAATCGCGCTGAGCCACGCACAATAc                    <  1:3696316/36‑1 (MQ=255)
             aTTACGATAAAATCGCGCTGAGCCACGCACAATACa                   >  1:140752/1‑36 (MQ=255)
             aTTACGATAAAATCGCGCTGAGCCACGCACAATACa                   <  1:1167370/36‑1 (MQ=255)
             aTTACGATAAAATCGCGCTGAGCCACGCACAATACa                   <  1:2265582/36‑1 (MQ=255)
              ttACGATAAAATCGCGCTGAGCCACGCATAATACAt                  >  2:791104/1‑36 (MQ=37)
              ttACGATAAAATCGCGCTGAGCCACGCACAATACAt                  <  2:2892959/36‑1 (MQ=255)
              ttACGATAAAATCGCGCTGAGCCACGCACAATACAt                  >  1:3720805/1‑36 (MQ=255)
               tACGATAAAATCGCGCTGGGCCACGCACAATACAtt                 >  2:3512389/1‑36 (MQ=255)
               tACGATAAAATCGCGCTGAGCCACGCACAATACAtt                 <  1:186046/36‑1 (MQ=255)
                aCGATAAAATCGCGCTGAGCTACGCACAATACATTg                >  1:987676/1‑36 (MQ=37)
                aCGATAAAATCGCGCTGAGCCACGCACAATACATTg                <  2:2812868/36‑1 (MQ=255)
                aCGATAAAATCGCGCTGAGCCACGCACAATACATTg                >  1:1947677/1‑36 (MQ=255)
                aCGATAAAATCGCGCTGAGCCACGCACAATACATTg                <  1:990717/36‑1 (MQ=255)
                aCGATAAAATCGCGCTGAGCCACGCACAATACATTg                >  2:1802296/1‑36 (MQ=255)
                 cGATAAAATCGCGCTGAGCCACGCACAATACATTGa               <  2:2784449/36‑1 (MQ=255)
                  gATAAAATCGCGCTGAGCCACGCACAATACATTGAt              <  2:3055968/36‑1 (MQ=255)
                  gATAAAATCGCGCTGAGCCACGCACAATACATTGAt              >  1:3342035/1‑36 (MQ=255)
                  gATAAAATCGCGCTGAGCCACGCACAATACATTGAt              >  1:229018/1‑36 (MQ=255)
                   aTAAAATCGCGCTGAGCCACGCACAATACATTGAtt             >  1:115898/1‑36 (MQ=255)
                   aTAAAATCGCGCTGAGCCACGCACAATACATTGAtt             >  1:516929/1‑36 (MQ=255)
                    tAAAATCGCGCTGAGCCACGCACCATACATTGAttt            >  2:2973606/1‑36 (MQ=37)
                    tAAAATCGCGCTGAGCCACGCACAATACATTGAttt            >  2:2348003/1‑36 (MQ=255)
                    tAAAATCGCGCTGAGCCACGCACAATACATTGAttt            >  2:2547141/1‑36 (MQ=255)
                     aaaaTCGCGCTGAGCCACGCACAATACATTGAtttt           >  2:395206/1‑36 (MQ=255)
                     aaaaTCGCGCTGAGCCACGCACAATACATTGAtttt           >  1:210263/1‑36 (MQ=255)
                      aaaTCGCGCTGAGCCACGCACAATACATTGATTTTg          >  2:1084417/1‑36 (MQ=255)
                       aaTCGCGCTGAGCCACGCACAATACATTGTTTTTgt         >  1:736659/1‑35 (MQ=37)
                        aTCGCGCTGAGCCACGCACAATACATTGATTTTGCa        >  1:1623730/1‑36 (MQ=255)
                        aTCGCGCTGAGCCACGCACAATACATTGATTTTGCa        <  1:1777168/36‑1 (MQ=255)
                         tCGCGCTAAGCCACGCACAATACATTGATTTTGCAg       <  1:3197348/36‑1 (MQ=37)
                          cgcgCTGGGCCACGCGCAATACATTGATTTTGCAGt      <  2:3610232/36‑1 (MQ=255)
                          cgcgCTGAGCCACGCACAATACATTGATTTTGCAGt      >  2:107900/1‑36 (MQ=255)
                          cgcgCGGAGCCACGCACAATACATTGATTTTGCAGt      <  2:108175/36‑1 (MQ=37)
                           gcgcTGAGCCACGCACAATACATTGATTTTGCTGTg     >  2:264000/1‑36 (MQ=37)
                           gcgcTGAGCCACGCACAATACATTGATTTTGCTGTg     >  2:1539131/1‑36 (MQ=37)
                            cgcTGAGCCACGCACAATACATTGATTTTGCAGTGa    >  1:3196539/1‑36 (MQ=255)
                            cgcTGAGCCACGCACAATACATTGATTTTGAAGTGa    >  2:1954714/1‑36 (MQ=38)
                             gcTGAGCCACGCCCAATACATTGATTTTGCAGTGaa   >  2:911917/1‑36 (MQ=37)
                             gcTGAGCCACGCACAATACATTGATTTTGCAGTGaa   >  1:972411/1‑36 (MQ=255)
                             gcTGAGCCACGCACAATACATTCATTCTGCAGTGaa   >  1:1734716/1‑36 (MQ=25)
                              cTGAGCCACGCACAATACATTGATTTTGCAGCGaaa  >  1:2734989/1‑36 (MQ=37)
                                 |                                
TTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAA  >  REL606/4616363‑4616428

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG