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% 206.3 571481yegIhypothetical protein
Reads supporting (aligned to +/- strand):  new base (28/29):  ref base (0/0):  total (28/29)

CGTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAGCACGTTTTTCTTCCTCCGGTAACCCGAGCATT  >  REL606/2082652‑2082717
                                 |                                
cGTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACa                                >  1:4204713/1‑36 (MQ=255)
cGTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACa                                <  1:6590630/36‑1 (MQ=255)
cGTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACa                                <  1:2233831/36‑1 (MQ=255)
cGTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACa                                >  1:2509915/1‑36 (MQ=255)
cGTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACa                                >  1:2982942/1‑36 (MQ=255)
 gTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACAc                               <  1:618819/36‑1 (MQ=255)
 gTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACAc                               >  1:827723/1‑36 (MQ=255)
 gTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACAc                               >  1:5348921/1‑36 (MQ=255)
 gTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACAc                               >  1:242193/1‑36 (MQ=255)
   tCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACACGt                             <  1:7482514/36‑1 (MQ=255)
   tCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACACGt                             <  1:4991737/36‑1 (MQ=255)
   tCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACACGt                             >  1:3605308/1‑36 (MQ=255)
   tCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACACGt                             <  1:5973349/36‑1 (MQ=255)
    cccTTGCGGTGTCGTGAAGTGCTGCCAGAACACGtt                            >  1:7779345/1‑36 (MQ=255)
    cccTTGCGGTGTCGTGAAGTGCTGCCAGAACACGtt                            >  1:5950282/1‑36 (MQ=255)
     ccTTGCGTTGTCGTGAAGTGCTGCCAGAACACGttt                           <  1:3184112/36‑1 (MQ=37)
     ccTTGCGGTGTCGTGAAGTGCTGCCAGAACACGttt                           <  1:3257956/36‑1 (MQ=255)
      cTTGCGGTGTCGTGAAGTGCTGCCAGAACACGtttt                          <  1:2945933/36‑1 (MQ=255)
      cTTGCGGTGTCGTGAAGTGCTGCCAGAACACGtttt                          <  1:6950322/36‑1 (MQ=255)
       ttGCGGTGTCGTGAAGTGCTGCCATAACACGttttt                         >  1:1864934/1‑36 (MQ=38)
       ttGCGGTGTCGTGAAGTGCTGCCAGAACACGttttt                         >  1:5343691/1‑36 (MQ=255)
       ttGCGGTGTCGTGAAGTGCTGCCAGAACACGttttt                         <  1:6635085/36‑1 (MQ=255)
        tGCGGTGTCGTGAAGTGCTGCCAGAACACGTTTttc                        >  1:2351947/1‑36 (MQ=255)
        tGCGGTGTCGTGAAGTGCTGCCAGAACACGTTTttc                        <  1:2322102/36‑1 (MQ=255)
        tGCGGTGTCGTGAAGTGCTGCCAGAACACGTTTttc                        >  1:4121831/1‑36 (MQ=255)
         gCGGTGTCGTGAAGTGCTGCCAGAACACGTTTttct                       >  1:2650147/1‑36 (MQ=255)
         gCGGTGTCGTGAAGTGCTGCCAGAACACGTTTttct                       >  1:5584904/1‑36 (MQ=255)
          cGGTGTCTTGAAGTGCTGCCAGAACACGTTTttctt                      <  1:6626907/36‑1 (MQ=37)
          cGGTGTCGTTAAGTGCTGCCAGAACACGTTTttctt                      <  1:7023703/36‑1 (MQ=37)
          cGGTGTCGTGCAGTGCTGGCAGAACACGTTTttctt                      <  1:6948690/36‑1 (MQ=25)
          cGGTGTCGTGAAGTGCTGCCAGAACACGTTTttctt                      <  1:4237517/36‑1 (MQ=255)
          cGGTGTCGTGAAGTGCTGCCAGAACACGTTTttctt                      >  1:5534084/1‑36 (MQ=255)
            ttgtCGTGAAGTGCTGCCAGAACACGTTTTTCTtcc                    <  1:7587485/35‑1 (MQ=255)
            gtgtCGTGAAGTGCTGCCAGAACACGTTTTTCTtcc                    >  1:177307/1‑36 (MQ=255)
             tgtCGTGAAGTGCTGCCAGAACGCGTTTTTCTtcct                   <  1:309244/36‑1 (MQ=37)
              gtCGTGAAGTGCTGCCAGAACACGTTTTTCTtcctc                  >  1:2957359/1‑36 (MQ=255)
               tCGTGAAGTGCTGCCAGAACACGTTTTTCTtcctcc                 <  1:3512838/36‑1 (MQ=255)
               tCGTGAAGTGCTGCCAGAACACGTTTTTCTtcctcc                 <  1:5586632/36‑1 (MQ=255)
               tCGTGAAGTGCTGCCAGAACACGTTTTTCTtcctcc                 <  1:832278/36‑1 (MQ=255)
               tCGTGAAGAGCTGCCAGAACACGTTTTTCTtcctcc                 <  1:3510898/36‑1 (MQ=37)
                cttGAAGTGCTGCCAGAACACGTTTTTCTTCCTCCg                <  1:6527301/34‑1 (MQ=37)
                 gTGAAGTGCTGCCAGAACACGTTTTTCTTCCTCCgg               <  1:3405284/36‑1 (MQ=255)
                 gTGAAGTGCTGCCAGAACACGTTTTTCTTCCTCCgg               >  1:5767898/1‑36 (MQ=255)
                 gTGAAGTGCTGCCAGAACACGTTTTTCTTCCTCCgg               >  1:3043822/1‑36 (MQ=255)
                 gTGAAGTGCTGCCAGAACACGTTTTTCTTCCTCCgg               <  1:7167051/36‑1 (MQ=255)
                 gTGAAGTGCTGCCAGAACACGTTTTTCTTCCTCCgg               >  1:1289753/1‑36 (MQ=255)
                 gTGAAGTGCTGCCAGAACACGTGTTTCTTCCTCCgg               >  1:923617/1‑36 (MQ=37)
                   gAAGTGCTGCCAGAACACGTTTTTCTTCCTCCGGTa             >  1:3845917/1‑36 (MQ=255)
                   gAAGTGCTGCCAGAACACGTTTTTCTTCCACCGGTa             >  1:4873615/1‑36 (MQ=37)
                    aaGTGCTGCCAGAACACGTTTTTCTTCCTCCGGTaa            >  1:6131971/1‑36 (MQ=255)
                    aaGTGCTGCCAGAACACGTTTTTCTTCCTCCGGTaa            <  1:4949766/36‑1 (MQ=255)
                    aaGTGCTGCCAGAACACGTTTTTCTTCCTCCGGTaa            >  1:2141163/1‑36 (MQ=255)
                        gctgcCAGAACACGTTTTTCTTCCTCCGGTAACCCg        <  1:5221704/36‑1 (MQ=255)
                         ctgcCAGAACACGTTTTTCTTCCTCCGGTAACCCGa       >  1:2519195/1‑36 (MQ=255)
                          tgcCAGAACACGTTTTTCTTCCTCCGGTAACCCGAg      >  1:1322078/1‑36 (MQ=255)
                            ccAGAACACGTTTTTCTTCCTCCGGTAACCCGAGCa    <  1:5373245/36‑1 (MQ=255)
                             cAGAACACGTTTTTCTTCCTCCGGTAACCCGAGCAt   <  1:2067934/36‑1 (MQ=255)
                              agaaCACGTTTTTCTTCCTCCGGTAACCCGAGCAtt  <  1:3742443/36‑1 (MQ=255)
                                 |                                
CGTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAGCACGTTTTTCTTCCTCCGGTAACCCGAGCATT  >  REL606/2082652‑2082717

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