Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,045,269 C→A 18.8% L67F (TTG→TTT yihP ← predicted transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,045,2690C→A18.8% 2.3 ‑2.5 0.033201yihPpredicted transporter
Reads supporting (aligned to +/- strand):  new base (1/6):  ref base (15/11):  total (16/17)
Fisher's exact test strand distribution p-value = 8.54e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ATCGGTAAACGCAGTAAAGAATTTTGAAATCAAAAAGATAATGCCGCCATAGGTGCCTGGCAGCC  >  REL606/4045239‑4045303
                              |                                  
atCGGTAAACGCAGTAAAGAATTTTGAAATCaaaaa                               <  1:6015603/36‑1 (MQ=255)
atCGGTAAACGCAGTAAAGAATTTTGAAATCaaaaa                               <  1:3999885/36‑1 (MQ=255)
atCGGTAAACGCAGTAAAGAATTTTGAAATCaaaaa                               >  1:2416478/1‑36 (MQ=255)
   ggTAAACGCAGTAAAGAATTTTGAAATCAAAAAGAt                            >  1:3868219/1‑36 (MQ=255)
   ggTAAACGCAGTAAAGAATTTTGAAATCAAAAAGAt                            <  1:5654416/36‑1 (MQ=255)
   ggTAAACGCAGTAAAGAATTTTGAAATCAAAAAGAt                            >  1:2720588/1‑36 (MQ=255)
      aaaCGCAGTAAAGAATTTTGAAATCAAAAAGATAAt                         <  1:5569426/36‑1 (MQ=255)
       aaCGCAGTAAAGAATTTTGAAATCAAAAAGATAATg                        >  1:2488836/1‑36 (MQ=255)
         cGCAGTAAAGAATTTTGAAATCAAAAAGATAATgcc                      <  1:4142718/36‑1 (MQ=255)
         cGCAGTAAAGAATTTTCAAATCAAAAAAATAATgcc                      <  1:2978146/36‑1 (MQ=37)
          gCAGTAAAGAATTTTGAAATCAAAAAGATAATgccg                     <  1:3685266/36‑1 (MQ=255)
           cAGTAAAGAATTTTGGAATCAAAACGATAATgccgc                    >  1:3936937/1‑36 (MQ=37)
            aGTAAAGAATTTTGAAATCGAAAAGATAATgccgcc                   >  1:4694124/1‑36 (MQ=255)
            aGTAAAGAATTTTGAAATCAAAAAGATAATgccgcc                   >  1:2920764/1‑36 (MQ=255)
            aGTAAAGAATTTTGAAATCAAAAAGATAATgccgcc                   <  1:2450965/36‑1 (MQ=255)
             gTAAAGAATTTTGAAATCAAAAAGATAATGCCGCCa                  >  1:2719777/1‑36 (MQ=255)
             gTAAAGAATTTTGAAATCAAAAAGATAATGCCGCCa                  >  1:5591557/1‑36 (MQ=255)
             gTAAAGAATTTTGAAATCAAAAAGATAATGCCGCCa                  >  1:4947638/1‑36 (MQ=255)
              tAAAGAATTTTGAAATAAAAAAGATAATGCCGCCAt                 >  1:2613734/1‑36 (MQ=255)
              tAAAGAATTTTGAAATAAAAAAGATAATGCCGCCAt                 <  1:2197273/36‑1 (MQ=255)
              tAAAGAATGTTGAAATAAAAAAGATAATGCCGCCAt                 <  1:5406672/36‑1 (MQ=37)
               aaaGAATTTAGAAAAAAAAAAGATAAGGCCGCCATa                <  1:468409/36‑1 (MQ=16)
                 aGAATTTTGAAATCAAAAAGATAATGCCGCCATAgg              <  1:1890418/36‑1 (MQ=255)
                   aaTTTTGAAATCAAAAAGATAATGCTGCCATATGTg            >  1:1546598/1‑36 (MQ=37)
                   aaTTTTGAAATCAAAAAGATAATGCCGCCATAGGTg            >  1:1404868/1‑36 (MQ=255)
                   aaTTTTGAAATAAAAAAGATAAGGCCGCCATAGGTg            <  1:2747370/36‑1 (MQ=37)
                      tttGAAATCAAAAAGATAATGCCGCCATAGGTGCCt         >  1:4821530/1‑36 (MQ=255)
                      tttAAAATAAAAAAGATAATGCCGCCATAGGTGCCt         <  1:6841083/36‑1 (MQ=37)
                         gAAATAAAAAAGATAATGCCGCCATAGGTGCCTGGc      <  1:6782477/36‑1 (MQ=255)
                            aTCAAAAAGATAATGCCGCCATAGGTGCCTGGCAGc   <  1:3185025/36‑1 (MQ=255)
                            aTCAAAAAGATAATGCCGCCATAGGTGCCTGGCAGc   >  1:6453986/1‑36 (MQ=255)
                            aTCAAAAAGACAATGCCGCCAAAAGCGCCGGGAAGc   >  1:6127451/1‑36 (MQ=9)
                             tCAAAAAGATAGTTCCGCCATAGGTGCCTGGCAGcc  <  1:6471934/36‑1 (MQ=38)
                              |                                  
ATCGGTAAACGCAGTAAAGAATTTTGAAATCAAAAAGATAATGCCGCCATAGGTGCCTGGCAGCC  >  REL606/4045239‑4045303

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 21 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG