Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,082,685 G→A 100% A494V (GCT→GTT)  yegI ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,082,6850G→A100.0% 127.0 351481yegIhypothetical protein
Reads supporting (aligned to +/- strand):  new base (20/14):  ref base (1/0):  total (21/14)

CGTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAGCACGTTTTTCTTCCTCCGGTAACCCGAGCATT  >  REL606/2082652‑2082717
                                 |                                
cGTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACa                                <  1:1093216/36‑1 (MQ=255)
 gTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACAc                               >  1:4059332/1‑36 (MQ=255)
   tCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACACGt                             >  1:4495920/1‑36 (MQ=255)
   tCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACACGt                             <  1:2075272/36‑1 (MQ=255)
     ccTTGCGGTGTCGTGAAGTGCTGCCAGAACACGttt                           <  1:3152982/36‑1 (MQ=255)
      cTTGCGGTGTCGTGAAGTGCTGCCAGAACACGtttt                          <  1:6287093/36‑1 (MQ=255)
      cTTGCGGTGTCGTGAAGTGCTGCCAGAACACGtttt                          <  1:5280550/36‑1 (MQ=255)
      cTTGCGGTGTCGTGAAGCGCTGCCAGAACACGtttt                          >  1:1103839/1‑36 (MQ=37)
       ttGCGGTGTCGTGAAGTGCTGCCAGAACACGttttt                         >  1:1480743/1‑36 (MQ=255)
       ttGCGGTGTCGTGAAGTGCTGCCAGAACACGttttt                         <  1:17543/36‑1 (MQ=255)
       ttGCGGTGTCGTGAAGTGCTGCCAGAACACGttttt                         <  1:4451173/36‑1 (MQ=255)
         gCGGTGTCGTGAAGTGCTGCCAGAACACGTTTttct                       >  1:2630062/1‑36 (MQ=255)
          cGGTGTCGTGAAGTGCTGCCAGAACACGTTTttctt                      <  1:1583224/36‑1 (MQ=255)
           ggTGTCGTGAAGTGCTGCCAGAACACGTTTttcttc                     >  1:907123/1‑36 (MQ=255)
            gtgtCGTGAAGTGCTGCCAGAACACGTTTTTCTtcc                    <  1:4478838/36‑1 (MQ=255)
             tgtCGTGAAGTGCTGCCAGAACACGTTTTTCTtcct                   >  1:1249844/1‑36 (MQ=255)
             tgtCGTGAAGTGCTGCCAGAACACGTTTTTCTtcct                   >  1:6738851/1‑36 (MQ=255)
             tgtCGTGAAGTGCTGCCAGAACACGTTTTTCTtcct                   <  1:4786512/36‑1 (MQ=255)
              gtCGTGAAGTGCTGCCAGAACACGTTTTTCTtcctc                  <  1:1846128/36‑1 (MQ=255)
                cGTGAAGTGCTGCCAGAACACGTTTTTCTTCCTCCg                >  1:3288282/1‑36 (MQ=255)
                 gTGAAGTGCTGCCAGAGCACGTTTTTCTTCCTCCgg               >  1:760878/1‑36 (MQ=255)
                 gTGAAGTGCTGCCAGAACACGTTTTTCTTCCTCCgg               <  1:2908401/36‑1 (MQ=255)
                 gTGAAGTGCTGCCAGAACACGTTTTTCTTCCTCCgg               >  1:7038909/1‑36 (MQ=255)
                 gTGAAGTGCTGCCAGAACACGTTTTTCTTCCTACgg               >  1:4798640/1‑36 (MQ=37)
                  tGAAGTGCTGCCAGAACACGTTTTTCTTCCTCCGGt              >  1:6752775/1‑36 (MQ=255)
                   gAAGTGCTGCCAGAACACGTTTTTCTTCCTCCGGTa             <  1:5768542/36‑1 (MQ=255)
                    aaGTGCTGCCAGAACACGTTTTTCTTCCTCCGGTaa            >  1:3092154/1‑36 (MQ=255)
                         ctgcCAGAACACGTTTTTCTTCATCCGGTAACCCGa       >  1:3532122/1‑36 (MQ=37)
                          tgcCAGAACACGTTTTTCTTCCTCCGGTAACCCGAg      >  1:7378504/1‑36 (MQ=255)
                          tgcCAGAACACGTTTTTCTTCCTCCGGTAACCCGAg      >  1:3844572/1‑36 (MQ=255)
                           gcCAGAACACGTTTTTCTTCCTCCGGTAACCCGAGc     <  1:3341914/36‑1 (MQ=255)
                           gcCAGAACACGTTTTTCTTCCTCCGGTAACCCGAGc     >  1:3585004/1‑36 (MQ=255)
                           gcCAGAACACGTTTTTCTTCCTCCGGTAACCCGAGc     >  1:2500758/1‑36 (MQ=255)
                            ccAGAACACGTTTTTCTTCCTCCGGTAACCCGAGCa    >  1:1578419/1‑36 (MQ=255)
                             cAGAACACGTTTTTCTTCCTCCGGTAACCCGAGCAt   >  1:5679881/1‑36 (MQ=255)
                              agaaCACGTTTTTCTTCCTCCGGTAACCCGAGCAtt  <  1:5624726/36‑1 (MQ=255)
                              agaaCACGTTTTTCTTCCTCCGGTAACCCGAGCAtt  <  1:2067248/36‑1 (MQ=255)
                              agaaCACGTTTTTCTTCCTCCGGTAACCCGAGCAtt  <  1:7806097/36‑1 (MQ=255)
                              agaaCACGTTTTTCTTCCTCCGGTAACCCGAGCAtt  <  1:886132/36‑1 (MQ=255)
                                 |                                
CGTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAGCACGTTTTTCTTCCTCCGGTAACCCGAGCATT  >  REL606/2082652‑2082717

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG