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% 166.5 55G290S (GGC→AGC) nadRnicotinamide‑nucleotide adenylyltransferase
Reads supporting (aligned to +/- strand):  new base (34/19):  ref base (1/1):  total (35/20)

TTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAA  >  REL606/4616363‑4616428
                                 |                                
ttGCAGTATTCCGATTACGATAAAATCGCGCTGAGc                                >  1:4053025/1‑36 (MQ=255)
  gCAGTATTCCGATTACGATAAAATCGCGCTGAGcaa                              >  1:7000973/1‑34 (MQ=25)
  gCAGTATTCCGATTACGATAAAATCGCGCTGAGCCa                              <  1:1728794/36‑1 (MQ=255)
  gCAGTATTCCGATTACGATAAAATCGCGCTGAGCCa                              >  1:1735257/1‑36 (MQ=255)
  gCAGTATTCCGATTACGATAAAATCGCGCTGAGCCa                              >  1:7001514/1‑36 (MQ=255)
  gCAGTATTCCGATTACGATAAAATCGCGCTGAGCCa                              <  1:698032/36‑1 (MQ=255)
    aGTATTCCGATTACGATAAAATCGCGCTGAGCCACg                            >  1:5071174/1‑36 (MQ=255)
     gTATTCCGATTACGATAAAATCGCGCTGAGCCACGc                           <  1:3760760/36‑1 (MQ=255)
     gTATTCCGATTACGATAAAATCGCGCTGAGCCACGc                           <  1:1561291/36‑1 (MQ=255)
       aTTCCGATTACGATAAAATCGCGCTGAGCCACGcac                         >  1:5750874/1‑36 (MQ=255)
         tCCGATTACGATAAAATCGCGCTGAGCCACGCACaa                       <  1:1924497/36‑1 (MQ=255)
          ccGATTACGATAAAATCGCGCTGAGCCACGCACAAt                      <  1:1448935/36‑1 (MQ=255)
          ccGATTACGATAAAATCGCGCTGAGCCACCCACAAt                      >  1:4407146/1‑36 (MQ=25)
           cGATTACGATAAAATCGCGCTGAGCCACGCACAATa                     <  1:2083970/36‑1 (MQ=255)
            gATTACGATAAAATCGCGCTGAGCCACGCACAATAc                    <  1:3696316/36‑1 (MQ=255)
             aTTACGATAAAATCGCGCTGAGCCACGCACAATACa                   <  1:2265582/36‑1 (MQ=255)
             aTTACGATAAAATCGCGCTGAGCCACGCACAATACa                   >  1:140752/1‑36 (MQ=255)
             aTTACGATAAAATCGCGCTGAGCCACGCACAATACa                   <  1:1167370/36‑1 (MQ=255)
              ttACGATAAAATCGCGCTGAGCCACGCATAATACAt                  >  1:4530314/1‑36 (MQ=25)
              ttACGATAAAATCGCGCTGAGCCACGCACAATACAt                  <  1:6632169/36‑1 (MQ=255)
              ttACGATAAAATCGCGCTGAGCCACGCACAATACAt                  >  1:3720805/1‑36 (MQ=255)
               tACGATAAAATCGCGCTGGGCCACGCACAATACAtt                 >  1:7251599/1‑36 (MQ=255)
               tACGATAAAATCGCGCTGAGCCACGCACAATACAtt                 <  1:186046/36‑1 (MQ=255)
                aCGATAAAATCGCGCTGAGCTACGCACAATACATTg                >  1:987676/1‑36 (MQ=25)
                aCGATAAAATCGCGCTGAGCCACGCACAATACATTg                <  1:6552078/36‑1 (MQ=255)
                aCGATAAAATCGCGCTGAGCCACGCACAATACATTg                >  1:5541506/1‑36 (MQ=255)
                aCGATAAAATCGCGCTGAGCCACGCACAATACATTg                <  1:990717/36‑1 (MQ=255)
                aCGATAAAATCGCGCTGAGCCACGCACAATACATTg                >  1:1947677/1‑36 (MQ=255)
                 cGATAAAATCGCGCTGAGCCACGCACAATACATTGa               <  1:6523659/36‑1 (MQ=255)
                  gATAAAATCGCGCTGAGCCACGCACAATACATTGAt              <  1:6795178/36‑1 (MQ=255)
                  gATAAAATCGCGCTGAGCCACGCACAATACATTGAt              >  1:229018/1‑36 (MQ=255)
                  gATAAAATCGCGCTGAGCCACGCACAATACATTGAt              >  1:3342035/1‑36 (MQ=255)
                   aTAAAATCGCGCTGAGCCACGCACAATACATTGAtt             >  1:115898/1‑36 (MQ=255)
                   aTAAAATCGCGCTGAGCCACGCACAATACATTGAtt             >  1:516929/1‑36 (MQ=255)
                    tAAAATCGCGCTGAGCCACGCACCATACATTGAttt            >  1:6712816/1‑36 (MQ=25)
                    tAAAATCGCGCTGAGCCACGCACAATACATTGAttt            >  1:6286351/1‑36 (MQ=255)
                    tAAAATCGCGCTGAGCCACGCACAATACATTGAttt            >  1:6087213/1‑36 (MQ=255)
                     aaaaTCGCGCTGAGCCACGCACAATACATTGAtttt           >  1:4134416/1‑36 (MQ=255)
                     aaaaTCGCGCTGAGCCACGCACAATACATTGAtttt           >  1:210263/1‑36 (MQ=255)
                      aaaTCGCGCTGAGCCACGCACAATACATTGATTTTg          >  1:4823627/1‑36 (MQ=255)
                       aaTCGCGCTGAGCCACGCACAATACATTGTTTTTgt         >  1:736659/1‑35 (MQ=25)
                        aTCGCGCTGAGCCACGCACAATACATTGATTTTGCa        <  1:1777168/36‑1 (MQ=255)
                        aTCGCGCTGAGCCACGCACAATACATTGATTTTGCa        >  1:1623730/1‑36 (MQ=255)
                         tCGCGCTAAGCCACGCACAATACATTGATTTTGCAg       <  1:3197348/36‑1 (MQ=25)
                          cgcgCTGGGCCACGCGCAATACATTGATTTTGCAGt      <  1:7349442/36‑1 (MQ=255)
                          cgcgCTGAGCCACGCACAATACATTGATTTTGCAGt      >  1:3847110/1‑36 (MQ=255)
                          cgcgCGGAGCCACGCACAATACATTGATTTTGCAGt      <  1:3847385/36‑1 (MQ=25)
                           gcgcTGAGCCACGCACAATACATTGATTTTGCTGTg     >  1:5278341/1‑36 (MQ=25)
                           gcgcTGAGCCACGCACAATACATTGATTTTGCTGTg     >  1:4003210/1‑36 (MQ=25)
                            cgcTGAGCCACGCACAATACATTGATTTTGCAGTGa    >  1:3196539/1‑36 (MQ=255)
                            cgcTGAGCCACGCACAATACATTGATTTTGAAGTGa    >  1:5693924/1‑36 (MQ=25)
                             gcTGAGCCACGCCCAATACATTGATTTTGCAGTGaa   >  1:4651127/1‑36 (MQ=25)
                             gcTGAGCCACGCACAATACATTGATTTTGCAGTGaa   >  1:972411/1‑36 (MQ=255)
                             gcTGAGCCACGCACAATACATTCATTCTGCAGTGaa   >  1:1734716/1‑36 (MQ=21)
                              cTGAGCCACGCACAATACATTGATTTTGCAGCGaaa  >  1:2734989/1‑36 (MQ=25)
                                 |                                
TTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAA  >  REL606/4616363‑4616428

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