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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,082,6850G→A16.2% 9.7 0.0 ‑0.1371481yegIhypothetical protein
Reads supporting (aligned to +/- strand):  new base (3/3):  ref base (13/18):  total (16/21)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.96e-01

GTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAGCACGTTTTTCTTCCTCCGGTAACCCGAGCATTTCG  >  REL606/2082653‑2082720
                                |                                   
gTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAACAc                                  >  1:6543328/1‑36 (MQ=255)
 ttCCCTTGCGGTGTCGTGAAGTGCTGCCAGAGCACg                                 >  1:662505/1‑36 (MQ=255)
    ccTTGCGGTGTCGTGAAGTGCTGCCAGAGCACGttt                              <  1:6756310/36‑1 (MQ=255)
    ccTTGCGGTGTCGTGAAGTGCTGCCAGAGCACGttt                              <  1:2773525/36‑1 (MQ=255)
      ttGCGGTGTCGTGAAGTTCTGCCAGAGCACGttttt                            >  1:7269327/1‑36 (MQ=255)
      ttGCGGTGTCGTGAAGTGCTGCCAGAGCACGttttt                            <  1:4365815/36‑1 (MQ=255)
      ttGCGGTGTCGTGAAGTGCTGCCAGAGCACGttttt                            >  1:1106252/1‑36 (MQ=255)
       tGCGGTGTCGTGAAGTGCTGCCAGAGCACGTTTttc                           <  1:1420016/36‑1 (MQ=255)
       tGCGGTGTCGTGAAGTGCTGCCAGAGCACGTTTttc                           <  1:3069131/36‑1 (MQ=255)
       tGCGGTGTCGTGAAGTGCTGCCAGAGCACGTTTttc                           <  1:4706792/36‑1 (MQ=255)
       tGCGGTGTCGTGAAGTGCTGCCAGAACACGTTTttc                           <  1:5272493/36‑1 (MQ=255)
        gCGGTGTCGTGAAGTGCTGCCAGAGCGCGTTTttct                          <  1:4295023/36‑1 (MQ=255)
        gCGGTGTCGTGAAGTGCTGCCAGAGCACGTTTttct                          >  1:1420938/1‑36 (MQ=255)
          ggTGTCGTGAAGTGCTGCCAGAGCACGTTTttcttc                        >  1:7001360/1‑36 (MQ=255)
          ggTGTCGTGAAGTGCTGCCAGAGCACGTTTttcttc                        <  1:3992158/36‑1 (MQ=255)
           gtgtCGTGAAGTGCTGCCAGAGCACGTTTTTCTtcc                       >  1:1141839/1‑36 (MQ=255)
            tgtTGTGAAGTGCTGCCAGAGCACGTTTTTCTtcct                      <  1:3075074/36‑1 (MQ=255)
            tgtCGTGCAGTGCTGCCAGAGCACGTTTTTCTtcct                      <  1:1469337/36‑1 (MQ=255)
            tgtCGTGAAGTGCTGCCAGAGCACGTTTTTCTtcct                      >  1:3442349/1‑36 (MQ=255)
            tgtCGTGAAGTGCTGCCAGAGCACGTTTTTCTtcct                      <  1:842258/36‑1 (MQ=255)
             gtCGTGAAGTGCTGCCAGAGCACGTTTTTCTtcctc                     <  1:4224851/36‑1 (MQ=255)
             gtCGTGAAGTGCTGCCAGAACACGTTTTTCTtcctc                     >  1:4110373/1‑36 (MQ=255)
              tCGTGAAGTGCTGCCAGAGCACGTTTTTCTtcctcc                    <  1:3733583/36‑1 (MQ=255)
               cGTGAAGTGCTGCCAGAGCACGTTTTTCTTCCTCCg                   <  1:3934834/36‑1 (MQ=255)
                 tGAAGTGCTGCCAGAGCACGTTTTTCTTCCTCCGGt                 <  1:3057046/36‑1 (MQ=255)
                  gcaGTGCTGCCATAGCACGTTTTTCTTCCTCCGGTa                <  1:4171847/34‑1 (MQ=38)
                   aaGTGCTGCCAGAGCACGTTTTTCTTCCTCCGGTaa               >  1:4938913/1‑36 (MQ=255)
                   aaGTGCTGCCAGAGCACGTTTTTCTTCCTCCGGCaa               >  1:899952/1‑36 (MQ=255)
                   aaGTGCTGCCAGAACACGTTTTTCTTCCTCCGGTaa               >  1:5406806/1‑36 (MQ=255)
                    aGTGCTGCCAGAGCACGTTTTTCTTCCTCCGGTAAc              <  1:5918503/36‑1 (MQ=255)
                    aGTGCTGCCAGAGCACGTTTTTCTTCCTCCGGTAAc              >  1:707818/1‑36 (MQ=255)
                     gTGCTGCCAGAGCACGTTTTTCTTCCTCCGGTAAcc             >  1:4864240/1‑36 (MQ=255)
                     gTGCTGCCAGAGCACGTTTTTCTTCCTCCGGTAAcc             >  1:2938115/1‑36 (MQ=255)
                     gTGCTGCCAGAACACGTTTTTCTTCCTCCGGTAAcc             <  1:5351094/36‑1 (MQ=255)
                      tgctgcCAGAGCACGTTTTTCTTCCTCCGGTAAccc            <  1:927489/36‑1 (MQ=255)
                        ctgcCAGAACACGTTTTTCTTCCTCCGGTAACCCGa          <  1:668526/36‑1 (MQ=255)
                         tgcCAGAGCACGTTTTTCTTCCTCCGGTAACCCGAg         >  1:4004212/1‑36 (MQ=255)
                             agagCACGTTTTTCTTCCTCCGGTAACCCGAGCAtt     >  1:5823479/1‑36 (MQ=255)
                             agagCACGTTTTTCTTCCTCCGGTAACCCGAGCAtt     >  1:5568026/1‑36 (MQ=255)
                             agagCACGTTTTTCTTCCTCCGGTAACCCGAGCAtt     <  1:4696914/36‑1 (MQ=255)
                             agagCACGTTTTTCTTCCTCCGGTAACCCGAGCAtt     <  1:1757870/36‑1 (MQ=255)
                               agCACGTTTTTCTTCCTCCGGTAACCCGAGCATTTc   <  1:1201449/36‑1 (MQ=255)
                               agCACGTTTTTCTTCCTCCGGTAACCCGAGCATTTc   <  1:7485775/36‑1 (MQ=255)
                               agCACGTTTTTCTTCCTCCGGTAACCCGAGCATTTc   <  1:923164/36‑1 (MQ=255)
                               agCACGTTTTTCTTCCTCCGGTAACCCGAGCATTTc   <  1:986746/36‑1 (MQ=255)
                                gCACGTTTTTCTTCCTCCGTTAACCCGCGCATTTCg  >  1:3361229/1‑36 (MQ=37)
                                gCACGTTTTTCTTCCTCCGGTAACCCGAGCATTTCg  >  1:3735891/1‑36 (MQ=255)
                                gCACGTTTTTCTTCCTCCGGTAACCCGAGCATTTCg  >  1:3138300/1‑36 (MQ=255)
                                |                                   
GTTCCCTTGCGGTGTCGTGAAGTGCTGCCAGAGCACGTTTTTCTTCCTCCGGTAACCCGAGCATTTCG  >  REL606/2082653‑2082720

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 37 ≤  ATCG