Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,560,010 C→A 13.7% L31F (TTG→TTT yfgM ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,560,0100C→A13.7% 2.2 ‑3.0 0.03593yfgMhypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (16/14):  total (16/19)
Fisher's exact test strand distribution p-value = 4.93e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCC  >  REL606/2559975‑2560039
                                   |                             
ttCCAGTAGCGCCAGCCAATCAGTGCGCCAACGccc                               >  1:5485368/1‑36 (MQ=255)
 tCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCa                              >  1:4890714/1‑36 (MQ=255)
 tCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCa                              <  1:140755/36‑1 (MQ=255)
 tCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCa                              >  1:8126034/1‑36 (MQ=255)
 gccAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCa                              <  1:8213017/35‑1 (MQ=255)
  ccAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCaa                             >  1:3758763/1‑36 (MQ=255)
   cAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCaaa                            <  1:4542069/36‑1 (MQ=255)
   cAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCaaa                            <  1:8469361/36‑1 (MQ=255)
     gTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAAt                          <  1:7914530/36‑1 (MQ=255)
      tAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATc                         <  1:3872894/36‑1 (MQ=255)
      tAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATc                         <  1:4393845/36‑1 (MQ=255)
       aGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCa                        >  1:1787448/1‑36 (MQ=255)
       aGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCa                        <  1:1200534/36‑1 (MQ=255)
        gcgcCAGCCAATCAGTGCGCCAACGCCCAAAATCAc                       <  1:3708111/36‑1 (MQ=255)
          gccagccaATCAGTGCGCCAACGCCCAAAATCAccc                     <  1:8259140/36‑1 (MQ=255)
          gccagccaATCAGTGCGCCAACGCCCAAAATCAccc                     >  1:4954463/1‑36 (MQ=255)
          gccagccaATCAGTGCGCCAACGCCCAAAATCAccc                     >  1:5824165/1‑36 (MQ=255)
           ccagccaATCAGTGCGCCAACGCCCAAAATCAcccc                    >  1:6029075/1‑36 (MQ=255)
           ccagccaATCAGTGCGCCAACGCCCAAAATCAcccc                    <  1:5737570/36‑1 (MQ=255)
             agccaATCAGTGCGCCAACGCCCAAAATCACCCCaa                  >  1:6567867/1‑36 (MQ=255)
             agccaATCAGTGCGCCAACGCCCAAAATCACCCCaa                  >  1:6292580/1‑36 (MQ=255)
              gccaATCAGTGCGCCAACGCCCAAAATCACCCCAAc                 <  1:5466509/36‑1 (MQ=255)
               ccaATCAGTGCGCCAACGCCCAAAATCACCCCAACa                <  1:7775666/36‑1 (MQ=255)
               ccaATCAGTGCGCCAACGCCCAAAATCACCCCAACa                <  1:1045669/36‑1 (MQ=255)
                 aaTCAGTGCGCCAACGCCCAAAATCACCCCAACAGc              >  1:5180082/1‑36 (MQ=255)
                 aaTCAGTGCGCCAACGCCCAAAATCACCCCAACAGc              >  1:5096500/1‑36 (MQ=255)
                  aTCAGTGCGCCAACGCCCAAAATCACCCCAACAGcc             >  1:8151745/1‑36 (MQ=255)
                    cAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAg           >  1:5673532/1‑36 (MQ=255)
                     aGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGt          >  1:7070450/1‑36 (MQ=255)
                       tGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGc        >  1:5412058/1‑36 (MQ=255)
                         cgcCAACGCCAAAAATCCCCCCAACAGCCAGTGCtt      <  1:95038/36‑1 (MQ=25)
                         cgcAAACCCCAAAAATCACCCCACCAGCCAGTGCtt      <  1:2960765/36‑1 (MQ=16)
                           ccAACGCCAAAAACCACCCCAACAGCCAGTGCTTTg    <  1:3686849/36‑1 (MQ=37)
                             aaCGCCCAAAATCACCCCAACAGCCAGTGCTTTGcc  >  1:3381070/1‑36 (MQ=255)
                             aaCGCAAAAAACCCCCCCAACAGCCAGTGCTTTGcc  <  1:6585905/36‑1 (MQ=16)
                             aaCCCCAAAAACCCCCCCACCAGCCAGTGCTTTGcc  <  1:5718994/36‑1 (MQ=16)
                                   |                             
TTCCAGTAGCGCCAGCCAATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCC  >  REL606/2559975‑2560039

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG