Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 179,945 +G 4.8% coding (522/624 nt) yadS ← conserved inner membrane protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606179,9451.→G 4.8% 3.9 ‑1.4 0.063coding (522/624 nt)yadSconserved inner membrane protein
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (30/30):  total (33/30)
Fisher's exact test strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ATCAATAGCGTCACGACCATGCCCATCATACTGGCT.G.TTTCCAGTGGTACGGAAAATGTGTAATAAGCCG  >  REL606/179910‑179979
                                    |                                   
aTCAATAGCGTCACGACCATGCCCATCATACTGGCt                                      <  1:6561838/36‑1 (MQ=255)
aTCAATAGCGTCACGACCATGCCCATCATACTGGCt                                      >  1:181290/1‑36 (MQ=255)
aTCAATAGCGTCACGACCATGCCCATCATACTGGCt                                      <  1:6813636/36‑1 (MQ=255)
aTCAATAGCGTCACGACCATGCCCATCATACTGGCt                                      <  1:7403600/36‑1 (MQ=255)
aTCAATAGCGTCACGACCATGCCCATCATACTGGCt                                      <  1:7320314/36‑1 (MQ=255)
 tCAATAGCGTCACTACCATGCCCATCATACTGGCT.g                                    <  1:2409191/36‑1 (MQ=255)
 tCAATAGCGTCACGACCATGCCCATCATACTGGCT.g                                    <  1:7968654/36‑1 (MQ=255)
  taaTAGCGTCACGACCATGCCCATCATACTGGCT.G.t                                  <  1:6557178/35‑1 (MQ=255)
  taaTAGCGTCACGACCATGCCCATCATACTGGCT.G.t                                  <  1:7247958/35‑1 (MQ=255)
  cAATAGCGTCACGACCATGCCCATCATACTGGCT.gg                                   >  1:2681691/1‑35 (MQ=255)
  cAATAGCGTCACGACCATGCCCATCATACTGGCT.gg                                   >  1:7461641/1‑35 (MQ=255)
   aaTAGCGTCACGACCATGCCCATCATACTGGCT.G.tt                                 <  1:3470134/36‑1 (MQ=255)
   aaTAGCGTCACGACCATGCCCATCATACTGGCT.G.tt                                 <  1:1324355/36‑1 (MQ=255)
    aTAGCGTCACGACCATGCCCATCATACTGGCT.gggt                                 >  1:4927482/1‑33 (MQ=255)
     tAGCGTCACGACCATGCCCATCATACTGGCT.G.TTTc                               <  1:1694629/36‑1 (MQ=255)
      aGCGTCACGACCATGCCCATCATACTGGCT.G.gtttc                              >  1:1492315/1‑34 (MQ=38)
       gCGTCACGACCATGCCCATCCTACTGGCT.G.ttttcc                             >  1:4310210/1‑33 (MQ=37)
       gCGTCACGACCATGCCCATCATACTGGCT.G.gtttgc                             >  1:7608117/1‑33 (MQ=37)
       gCGTCACGACCATGCCCATCATACTGGCT.G.TTTCCa                             <  1:3598916/36‑1 (MQ=255)
       gCGTCACGACCATGCCCATCATACTGGCT.G.TTTCCa                             <  1:405049/36‑1 (MQ=255)
       gCGTCACGACCATGCCCATCATACTGGCT.G.GTTccc                             >  1:7685222/1‑35 (MQ=255)
       gCGTCACGACCATGCCCATCATACTGGCT.G.GTTTCa                             >  1:2038529/1‑36 (MQ=37)
        cGTCACGACCATGCCCATCATACTGGCTGG.TTTccg                             >  1:6277935/1‑35 (MQ=37)
         gTCACGACCATGCCCATCATACTGGCTGG.TTTccgt                            >  1:4393700/1‑34 (MQ=25)
         gTCACGACCATGCCCATCATACTGGCTGG.TTTccct                            >  1:7187834/1‑34 (MQ=37)
         gTCACGACCATGCCCATCATACTGGCT.G.TTTCCAGt                           <  1:97151/36‑1 (MQ=255)
         gTCACGACCATGCCCATCATACTGGCT.G.TTTCCAGt                           <  1:2615436/36‑1 (MQ=255)
          tCACGACCATGCCCATCATACTGGCT.G.GTTCCCGTg                          >  1:707176/1‑36 (MQ=37)
           cACGACCATGCCCATCATACTGGCT.GTTTTCCAGTg                          >  1:6182437/1‑36 (MQ=37)
           cACGACCATGCCCATCATACTGGCT.G.TTTCCCGTgg                         >  1:1889746/1‑36 (MQ=255)
            aCGACCATGCCCATCATACTGGCT.G.TTTCCAGTGGt                        <  1:7338071/36‑1 (MQ=255)
            aCGACCATGCCCATCATACTGGCT.G.GTTTCCGTggg                        >  1:6484258/1‑35 (MQ=21)
            aCGACCATGCCCATCATACTGGCT.G.GTTCCCGTGGt                        >  1:1916690/1‑36 (MQ=37)
             cGACCATGCCCATCATACTGGCT.G.TTTCCAGTGGTa                       >  1:639958/1‑36 (MQ=255)
             cGACCATGCCCATCATACTGGCT.G.GTTCCCGTGGTa                       >  1:2899752/1‑36 (MQ=37)
              gACCATGCCCATCATACTGGCT.G.TTTCCAGTGGTAc                      <  1:6916260/36‑1 (MQ=255)
              gACCATGCCCATCATACTGGCT.G.TTTCCAGTGGTAc                      <  1:4133366/36‑1 (MQ=255)
              gACCATGCCCATCATACTGGCT.G.TTTCCAGTGGTAc                      <  1:3406305/36‑1 (MQ=255)
               aCCATGCCCATCATACTGGCT.G.TTTCCAGTGGTACg                     >  1:8180735/1‑36 (MQ=255)
               aCCATGCCCATCATACTGGCT.G.GTACCAGTGGTTCg                     >  1:830803/1‑36 (MQ=25)
                ccATGCCCATCATACTGGCT.G.TTTCCCGTGGTACgg                    >  1:2232745/1‑36 (MQ=255)
                 cATGCCCATCATACTGGCT.G.TTTTCCGTGGTACggg                   >  1:4634689/1‑35 (MQ=37)
                 cATGCCCATCATACTGGCT.G.TTTCCAGTGGTACGGa                   <  1:1977095/36‑1 (MQ=255)
                  aTGCCCATCATACTGGCT.G.TTTCCATTGGTACGGaa                  <  1:7304446/36‑1 (MQ=255)
                  aTGCCCATCATACTGGCT.G.TTTCCAGTGGTACGGaa                  <  1:3071262/36‑1 (MQ=255)
                  aTGCCCATCATACTGGCT.G.TTTCCAGTGGTACGGaa                  <  1:4625143/36‑1 (MQ=255)
                    gCCCATCATACTGGCT.G.TTTCCCGTGGTACGGgaaa                >  1:4621020/1‑36 (MQ=37)
                    gCCCATCATACTGGCT.G.TTTCCAGTGGTACGGaaaa                <  1:942097/36‑1 (MQ=255)
                    gCCCATCATACTGGCT.G.TTTCCAGTGGTACGGaaaa                <  1:3593305/36‑1 (MQ=255)
                    gCCCATCATACTGGCT.G.GTTCCCGTGGTTCGGgaaa                >  1:6661600/1‑36 (MQ=16)
                    gCCCATCATACTGGCT.G.GTTCCCGTGGTACGGaaaa                >  1:7518707/1‑36 (MQ=37)
                     cccATCATACTGGCT.G.TTTCCAGTGGTACGGAAAAt               <  1:1675734/36‑1 (MQ=255)
                     cccATCATACTGGCT.G.GTTTCAGTGGTACGGAAAAt               >  1:4375574/1‑36 (MQ=37)
                      ccATCATACTGGCT.G.TTTCCAGTGGTACGGAAAAtt              >  1:7374234/1‑35 (MQ=255)
                      ccATCATACTGGCT.G.GTTCCAGTGGTGCGGAAAAtg              >  1:7709891/1‑36 (MQ=37)
                       catcatACTGGCT.G.TTTCCAGTGGTACGGAAAAtgt             <  1:1749110/36‑1 (MQ=255)
                         tcatACTGGCT.G.TTTCCAGTGGTACGGAAAAtgttt           >  1:4851171/1‑34 (MQ=255)
                         tcatACTGGCT.G.TTTCCAGTGGTACGGAAAATgtgt           <  1:8023068/36‑1 (MQ=255)
                          catACTGGCT.G.TTTCCAGTGGTACGGGAAATGGGta          >  1:8299303/1‑36 (MQ=37)
                             aCTGGCT.G.TTTCCAGTGGTACGGAAAATGTGtaata       <  1:2134382/36‑1 (MQ=255)
                                ggCT.G.TTTCCAGTGGTACGGAAAATGTGTAATAAGc    >  1:5023858/1‑36 (MQ=255)
                                ggCT.G.TTTCCAGTGGTACGGAAAATGTGTAATAAGc    <  1:5794434/36‑1 (MQ=255)
                                ggCT.G.TTTCCAGTGGTACGGAAAATGTGTAATAAGc    >  1:8005240/1‑36 (MQ=255)
                                ggCT.G.TTTCCAGTGGTACGGAAAATGTGTAATAAGc    >  1:7424152/1‑36 (MQ=255)
                                ggCT.G.TTTCCAGTGGTACGGAAAATGTGTAATAAGc    <  1:197869/36‑1 (MQ=255)
                                ggCT.G.TTTCCAGTGGTACGGAAAATGTGTAATAAGc    >  1:4424491/1‑36 (MQ=255)
                                 gCT.G.TTTCCAGTGGTACGTAAAATGTGTAATAAGcc   <  1:6248526/36‑1 (MQ=255)
                                 gCT.G.TTTCCAGTGGTACGGAAAATGTGTAATAAGcc   <  1:2832104/36‑1 (MQ=255)
                                 gCT.G.TTTCCAGTGGTACGGAAAATGTGTAATAAGcc   <  1:1963000/36‑1 (MQ=255)
                                 gCT.G.TTTCCAGTGGTACGGAAAATGTGTAATAAGcc   <  1:2386130/36‑1 (MQ=255)
                                 gCT.G.TTTCCAGTGGTACGGAAAATGTGTAATAAGcc   >  1:5517056/1‑36 (MQ=255)
                                 gCT.G.TTTCCAGTGGTACGGAAAATGTGTAATAAGcc   >  1:7172456/1‑36 (MQ=255)
                                  cT.G.TTTCCAGTGGTACGGAAAATGTGTAATAAGCCg  <  1:6012111/36‑1 (MQ=255)
                                    |                                   
ATCAATAGCGTCACGACCATGCCCATCATACTGGCT.G.TTTCCAGTGGTACGGAAAATGTGTAATAAGCCG  >  REL606/179910‑179979

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