Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,311,493 C→A 18.4% intergenic (+1/+38) yfaO → / ← ais predicted NUDIX hydrolase/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,311,4930C→A18.4% 3.6 ‑4.4 ‑0.041intergenic (+1/+38)yfaO/aispredicted NUDIX hydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/8):  ref base (17/16):  total (17/24)
Fisher's exact test strand distribution p-value = 1.28e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.99e-01

CGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGA  >  REL606/2311460‑2311526
                                 |                                 
cGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACaa                                 >  1:5107176/1‑36 (MQ=255)
cGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACaa                                 >  1:5387769/1‑36 (MQ=255)
cGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACaa                                 <  1:1545660/36‑1 (MQ=255)
cGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACaa                                 >  1:1966877/1‑36 (MQ=255)
 gAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAat                                <  1:4922143/36‑1 (MQ=255)
 gAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAat                                >  1:3966438/1‑36 (MQ=255)
  aaaaaCGTTACGTTTGAAAGGTCTTCTGTAACAatg                               >  1:4335230/1‑36 (MQ=255)
   aaaaCGTTACGTTTGAAAGGTCTTCTGTAACAatga                              >  1:825562/1‑36 (MQ=255)
   aaaaCGTTACGTTTGAAAGGTCTTCTGTAACAatga                              <  1:5413030/36‑1 (MQ=255)
    aaaCGTTACGTTTGAAAGGTCTTCTGTACCAatgac                             <  1:7083518/36‑2 (MQ=255)
    aaaCGTTACGTTTGAAAGGTCTTCTGTAACAatgat                             <  1:2615094/36‑1 (MQ=255)
     aaCGTTACGTTTGAAAGGTCTTCTGTAACAatgatg                            >  1:1749153/1‑36 (MQ=255)
         ttacgttTGAAAGGTCTTCTGTAACAATGATGACAg                        <  1:1731585/36‑1 (MQ=255)
          tacgttTGAAAGGTCTTCTGTAACAATGATGACAGc                       >  1:8369035/1‑36 (MQ=255)
          tacgttTGAAAGGTCTTCTGTAACAATGATGACAGc                       <  1:2633686/36‑1 (MQ=255)
           acgttTGAAAGGTCTTCTGTAACAATGATGACAGcc                      >  1:923694/1‑36 (MQ=255)
           acgttTGAAAGGTCTTCTGTAACAATGATGACAGcc                      <  1:1343790/36‑1 (MQ=255)
            cgttTGAAAGGTCTTCTGTAACAATGATGACAGCCa                     <  1:4782202/36‑1 (MQ=255)
            cgttTGAAAGGTCTTCTGTAACAATGATGACAGCCa                     >  1:2073012/1‑36 (MQ=255)
            cgttTGAAAGGTCTTCTGTAACAATGATGACAGCCa                     <  1:7220553/36‑1 (MQ=255)
             gttTGAAAGGTCTTCTGTAACAATGATGACAGCCAt                    <  1:6044715/36‑1 (MQ=255)
                agAACGGTCTTCTGTAACAATGATGCCAGCCATTAc                 <  1:2023177/35‑1 (MQ=37)
                 gAAAGGTCTTCTGTAACAATGATGACAGCCATTACt                >  1:2466790/1‑36 (MQ=255)
                   aaGGTCTTCTGTAAAAATGATGACAGCCATTACTgg              <  1:646969/36‑1 (MQ=255)
                   aaGGTCTTCGGTAAAAATGATGACAGCCATTGCTgg              <  1:4817673/36‑1 (MQ=25)
                    aGGTCTTCTGTAACAATGATGACAGCCATTACTGGc             >  1:6392292/1‑36 (MQ=255)
                     ggTTTTCTGAAAAAATGATGACAGCCATTACTGGCt            <  1:2345918/36‑1 (MQ=25)
                     ggTCTTCTGTAACAATGATGACAGCCATTACTGGCt            <  1:1321525/36‑1 (MQ=255)
                     ggTCTTCTGTAAAAATGATGACAGCCATTACTGGCt            <  1:2518733/36‑1 (MQ=255)
                      gTCTTCTGTAACAATGATGACAGCCATTACTGGCTg           >  1:3480776/1‑36 (MQ=255)
                      gTCTTCTGTAACAATGATGACAGCCATTACTGGCTg           >  1:936874/1‑36 (MQ=255)
                       ttttctGTAAAAATAATGACAGCCATTACTGGCTGt          <  1:2173439/34‑1 (MQ=37)
                       tcttctATAAAAATGATTACAGCCATTACTGGCTGt          <  1:3119634/36‑1 (MQ=21)
                        cttctAAAAAAATGATGACAGCCATTACTGGCTGTc         <  1:1453461/36‑1 (MQ=25)
                           ttGTAAAAATGACGACAGCCATTACTGGCTGTCAGa      <  1:7621405/35‑1 (MQ=37)
                           ctGTAACAATGGTGACAGCCATTACTGGATGTCAGa      <  1:3669005/36‑1 (MQ=37)
                           ctGTAACAATGATGACAGCCATTACTGGCTGTCCGa      >  1:3529361/1‑36 (MQ=255)
                             ggaaCAATGATGACAGCCATTACTGGCTGCAAGAtt    >  1:6230423/3‑36 (MQ=37)
                             gTAACAATGATGACAGCCATTACTGGCTGTCAGAtt    <  1:1324682/36‑1 (MQ=255)
                              aaaCAATGATGACAGCCATTACTGGCTGTCAGATTg   <  1:6623187/35‑1 (MQ=255)
                               aaCAATGATGACAGCCATTACTGGCTGGCAGAttta  >  1:4664490/1‑34 (MQ=38)
                                 |                                 
CGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGA  >  REL606/2311460‑2311526

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 17 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG