Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,236,456 C→G 10.2% L312F (TTG→TTC ycgB ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,236,4560C→G10.2% 2.1 ‑2.2 ‑0.141936ycgBhypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (16/18):  total (23/18)
Fisher's exact test strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.00e-01

TAGGGGGGCTGGAAGACCACATTGGTGTGGCTGTGCAAAAACTCCAGCATAAAACGTTCCGTTACTTTCCC  >  REL606/1236421‑1236491
                                   |                                   
taGGGGGGCTGGAAGACCACATTGGTGTGGCTGTGc                                     <  1:4336049/36‑1 (MQ=255)
  ggggggCTGGAAGACCACATTGGTGTGGCTGTGCaa                                   <  1:3715088/36‑1 (MQ=255)
    ggggCTGGAAGACCACATTGGTGTGGCTGTGCaaaa                                 <  1:1065808/36‑1 (MQ=255)
    ggggCTGGAAGACCACATTGGTGTGGCTGGGCgaaa                                 >  1:5291673/1‑36 (MQ=37)
      ggCTGGAAGACCACATTGGTGTGGCTGTGGaaaaaa                               >  1:1435438/1‑35 (MQ=255)
      ggCTGGAAGACCACATTGGTGTGGCTGTGCAAAAAc                               <  1:1313197/36‑1 (MQ=255)
       gCTGGAAGACCACATTGGTGTGGCTGTGCAAAAACt                              <  1:4060089/36‑1 (MQ=255)
        cTGGAAGACCACATTGGTGTGGCTGGGAAAAAACtg                             >  1:3602956/1‑35 (MQ=37)
        cTGGAAGACAACATCGGTGGGGCTGTGGAAAAAcaa                             >  1:431521/1‑34 (MQ=16)
         tagAAGCCAACAATGGTGAGGCTGTGCAAAAACTcc                            <  1:2449301/34‑1 (MQ=16)
         tGGAAGACCACATTGGTGTGGCTGTGCAAAAACTcc                            <  1:1701743/36‑1 (MQ=255)
          ggAAGACCACATTGGTGTGGCTGTGTAAAAACTCCa                           >  1:2762200/1‑36 (MQ=255)
          ggAAGACCACATTGGTGTGGCTGGGCAAAAACGccc                           >  1:4231427/1‑35 (MQ=37)
           gAAGACCACATTGGTGTGGCTGTGCAAAAACTCCAg                          <  1:3591982/36‑1 (MQ=255)
            aaGACCACATTGGTGTGGCTGTGGAAAAACTCCAGc                         >  1:3506700/1‑36 (MQ=255)
            aaGACCACATTGGTGTGGCTGTGCAAAAACTCCAGc                         >  1:5750942/1‑36 (MQ=255)
            aaGACCACATTGGTGTGGCTGTGCAAAAACTCCAGc                         <  1:5211517/36‑1 (MQ=255)
            aaGACCAAATTGGTGTGGCTTTTCAAAAACTCCAGc                         <  1:5493516/36‑1 (MQ=25)
              gACCACATTGGTGTGGCTGTGCAAAAACTCCAGGAt                       >  1:1808395/1‑36 (MQ=255)
              gACCACATTGGTGTGGCTGTGCAAAAACTCCAGCAt                       <  1:3067001/36‑1 (MQ=255)
               aCCACATTGGTGTGGCTGGGCAAAAACTCCAGCATa                      >  1:567489/1‑36 (MQ=255)
                  acaTTGGTGTGGCTGTGCAAAAACTCCAGCATaaaa                   <  1:873152/36‑1 (MQ=255)
                  acaTTGGTGTGGCTGTGCAAAAACTCCAGCATaaaa                   <  1:5296674/36‑1 (MQ=255)
                  acaTTGGTGTGGCTGGGCAAAAAGTCCAGCAAaaaa                   >  1:3575835/1‑36 (MQ=37)
                   caTTGGTGTGGCTGTGCAAAAACTCCAGAAGAAAAc                  >  1:3473062/1‑36 (MQ=37)
                    aTTGGTGTGGCTGTGGAAAAAATCCAGCATAAAACg                 >  1:813857/1‑36 (MQ=37)
                    aTTGGTGTGGCTGTGAAAAAACTCCAGCATAAAACg                 >  1:3485810/1‑36 (MQ=255)
                       ggTTTGGCTATGCAAAAACTCCAGCATAAAACGTTc              <  1:1368971/36‑1 (MQ=37)
                        gtgtGGCTGTGCAAAAACTCCAGCATAAAACGtttc             >  1:2381425/1‑34 (MQ=255)
                        gtgtGGCTGTGCAAAAACTCCAGCATAAAACGTTcc             >  1:1308445/1‑36 (MQ=255)
                         tgtgGCTGTGCAAAAACTTCAGCATAAAACGTTCCg            >  1:1647393/1‑36 (MQ=255)
                          gtgGCTGTGCAAAAACTCCAGCATAAAACGTTCCGt           <  1:5867123/36‑1 (MQ=255)
                           tgGCTGTGCAAAAACTCCAGCATAAAACGTTCCGtt          <  1:3767340/36‑1 (MQ=255)
                           tgGCTGTGCAAAAACTCCAGCATAAAACGTTCCGtt          <  1:5444912/36‑1 (MQ=255)
                             gCTGTGCAAAAACTCCAGCATAAAACGTTCCGTTAc        >  1:227186/1‑36 (MQ=255)
                             gCTGTGCAAAAACTCCAGCATAAAACGTTCCGTTAc        >  1:324303/1‑36 (MQ=255)
                               tgtgCAAAAACTCCAGCATAAAACGTTCCGTTACtt      >  1:5004770/1‑36 (MQ=255)
                               tgtgCAAAAACTCCAGCATAAAACGTTCCGTTACtt      >  1:1293741/1‑36 (MQ=255)
                               tgtcCAAAAACTCCAGCATAAAACGTTCCGTTACtt      <  1:3891730/36‑1 (MQ=255)
                                gtgCAAAAACTCCAGCATAAAACGTTCCGTTACttt     <  1:276301/36‑1 (MQ=255)
                                  gCAAAAACTCCAGCATAAAACGTTCCGTTACTTTcc   >  1:4992262/1‑36 (MQ=255)
                                  gCAAAAACTCCAGCATAAAACGTTCCGTTACTTTcc   >  1:5160640/1‑36 (MQ=255)
                                   cAAAAACTCCAGCATAAAACGTTCCGTTACTTTccc  <  1:2427747/36‑1 (MQ=255)
                                   |                                   
TAGGGGGGCTGGAAGACCACATTGGTGTGGCTGTGCAAAAACTCCAGCATAAAACGTTCCGTTACTTTCCC  >  REL606/1236421‑1236491

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 16 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG