Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 9,972 T→G 4.8% N174T (AAC→ACC)  yaaH ← conserved inner membrane protein associated with acetate transport

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6069,9720T→G 4.8% 2.4 ‑0.5 ‑0.263521yaaHconserved inner membrane protein associated with acetate transport
Reads supporting (aligned to +/- strand):  new base (1/2):  ref base (34/26):  total (35/28)
Fisher's exact test strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.64e-01

CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTG  >  REL606/9937‑10007
                                   |                                   
cACCAATCGGCAGAACGGTGCGACCAAACTGCTCGt                                     >  1:5712408/1‑36 (MQ=255)
 aCCAATCGGCAGAACGGTGCGACCAAACTGCTCGtt                                    >  1:3643824/1‑36 (MQ=255)
  ccAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                   >  1:4820465/1‑36 (MQ=255)
  ccAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                   <  1:4869725/36‑1 (MQ=255)
  ccAATCGGCAGAACGGTGCGACCAAACTGCTCGTTc                                   <  1:1960707/36‑1 (MQ=255)
   cAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCa                                  <  1:3561371/36‑1 (MQ=255)
   cAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCa                                  <  1:5323245/36‑1 (MQ=255)
    aaTCGGCAGAACGGTGCGACCAAACTGCTTGTTCAg                                 >  1:944053/1‑36 (MQ=255)
    aaTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAg                                 >  1:6048356/1‑36 (MQ=255)
    aaTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAg                                 >  1:2778409/1‑36 (MQ=255)
    aaTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAg                                 >  1:6208748/1‑36 (MQ=255)
    aaTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAg                                 >  1:6407607/1‑36 (MQ=255)
     aTCGGCAGAACGGTGCGACCAAACTGCTCGTTCatt                                <  1:1844506/36‑3 (MQ=255)
     aTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                                >  1:5019646/1‑36 (MQ=255)
     aTCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGt                                <  1:4752155/36‑1 (MQ=255)
     aTCGGCAGAACGGGGAGCCCAAACGACTCGTTCCGt                                >  1:7250244/1‑36 (MQ=12)
     aTCGGCAGAACGGCGCGACCAAACTGCTAGTTCAGt                                >  1:8315767/1‑36 (MQ=37)
      tCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTa                               >  1:3929556/1‑36 (MQ=255)
      tCGGCAGAACGGTGCGACCAAACTGCTCGGTCAGTa                               <  1:6448051/36‑1 (MQ=255)
       cGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTAc                              >  1:2393815/1‑36 (MQ=255)
       cGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTAc                              >  1:1222186/1‑36 (MQ=255)
        ggCAGAACGGTGCGACCAAACTGCTCGTTCAGTACt                             >  1:3142682/1‑36 (MQ=255)
        ggCAGAACGGTGCGACCAAACTGCTCGTTCAGTACt                             <  1:6635957/36‑1 (MQ=255)
        ggCAGAACGGTGCGACCAAACTGCTCGTTCAGTACt                             <  1:8412138/36‑1 (MQ=255)
         gCAGAACGGTGCGACCAAACTGCTCGTTCAGTACtt                            >  1:8442061/1‑36 (MQ=255)
           aGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCa                          <  1:1995355/36‑1 (MQ=255)
           aGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCa                          >  1:2003071/1‑36 (MQ=255)
           aGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCa                          >  1:6997926/1‑36 (MQ=255)
           aGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCa                          >  1:3229791/1‑36 (MQ=255)
            gAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAc                         >  1:4434961/1‑36 (MQ=255)
            gAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAc                         >  1:5811163/1‑36 (MQ=255)
            gAACGGTGCGACCAAACTGCTCGTTCAGTACTTCAc                         >  1:6576068/1‑36 (MQ=255)
            gAACGGTGCAACCAACCTGCTCCTTCAGTACTTCAc                         <  1:5477431/36‑1 (MQ=25)
              aCGGTGCGACCAAACTGCTCGTTCAGTACTTCAccc                       <  1:2110180/36‑1 (MQ=255)
               cGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCa                      <  1:3150175/36‑1 (MQ=255)
                ggTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                     <  1:2444849/36‑1 (MQ=255)
                ggTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                     <  1:7573488/36‑1 (MQ=255)
                ggTGCGACCAAACTGCTCGTTCAGTACTTCACCCAt                     <  1:2181976/36‑1 (MQ=255)
                 gTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                    >  1:6847658/1‑36 (MQ=255)
                 gTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                    <  1:1361276/36‑1 (MQ=255)
                 gTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                    <  1:7502501/36‑1 (MQ=255)
                 gTGCGACCAAACTGCTCGTTCAGTACTTCACCCATc                    <  1:1086614/36‑1 (MQ=255)
                   gCGACCAAACTGCTCGTTCAGTACTTCACCCATCGc                  >  1:3281443/1‑36 (MQ=255)
                   gCGACCAAACTGCTCGTTCAGTACTTCACCCATCGc                  >  1:3481436/1‑36 (MQ=255)
                    cGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc                 >  1:3178957/1‑36 (MQ=255)
                    cGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc                 <  1:1345453/36‑1 (MQ=255)
                    cGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc                 >  1:3023679/1‑36 (MQ=255)
                     gACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa                >  1:5300175/1‑36 (MQ=255)
                     gACCAAACTGCTCGTTCAGTACTTCACCCATCGCCa                <  1:236558/36‑1 (MQ=255)
                      aCCAAACTGCTCGTTCAGTACTTCACCCATCGCCag               >  1:2221855/1‑36 (MQ=255)
                      aCCAAACTGCTCGTTCAGTACTTCACCCATCGCCag               >  1:5689976/1‑36 (MQ=255)
                        cAAACTGCTCGTTCAGTACTTCACCCATCGCCagat             >  1:4014271/1‑36 (MQ=255)
                         aaaCTGCTCGTTCAGTACTTCACCCATCGCCagata            <  1:7297892/36‑1 (MQ=255)
                         aaaCTGCTCGTTCAGTACTTCACCCATCGCCagata            >  1:694693/1‑36 (MQ=255)
                         aaaCTGCTCGTTCAGTACTTCACCCATCGCCagata            >  1:5708081/1‑36 (MQ=255)
                         aaaCTGCTCGTTCAGTACTTCACACCTCGCCatata            <  1:2477179/36‑1 (MQ=37)
                           aCTGCTCGTTCAGTACTTCACCCATCGCCagataga          <  1:6792277/36‑1 (MQ=255)
                           aCTGCTCGTTCAGTACTTCACCCATCGCCagataga          >  1:662981/1‑36 (MQ=255)
                            cTGCTCGTTCAGTACTTCACCCATCGCCagatagat         <  1:1357613/36‑1 (MQ=255)
                            cTGCTCGGTCAGTACTTCACCCATCGCCagatagat         >  1:4197617/1‑36 (MQ=255)
                               cTCGTTCAGTACTTCACCCATCGCCAGATAGATTGc      <  1:1343183/36‑1 (MQ=255)
                               cTCGGTCAGTACTTCACCCATCGCCAGATAGATTGc      <  1:4190042/36‑1 (MQ=255)
                                tCGTTCATTACTTCACCCATCGCCAGATAGATTGCa     <  1:3504595/36‑1 (MQ=255)
                                tCGTTCAGTACTTCATCCATCGCCAGATAGATTGCa     >  1:8299418/1‑36 (MQ=255)
                                  gTTCAGTACTTCACCCATCGCCAGATAGATTGCACt   >  1:36949/1‑36 (MQ=255)
                                   ttCAGTACTTCACCCATCGCCAGATAGATTGCACTg  <  1:6628827/36‑1 (MQ=255)
                                   ttCAGTACTTCACCCATCGCCAGATAGATTGCACTg  <  1:8054188/36‑1 (MQ=255)
                                   ttCAGTACTTCACCCATCGCCAGATAGATTGCACTg  >  1:4341066/1‑36 (MQ=255)
                                   |                                   
CACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATAGATTGCACTG  >  REL606/9937‑10007

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG