Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 612,409 +G 7.5% coding (402/414 nt) ybdB → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606612,4091.→G 7.5% 4.5 ‑3.1 0.040coding (402/414 nt)ybdBhypothetical protein
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (12/25):  total (15/25)
Fisher's exact test strand distribution p-value = 4.61e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCA.GTTTTGGGATGAGCCCGACAGAATTAAATGAGA  >  REL606/612374‑612442
                                    |                                 
gggCGGCGTTGCTGCACTTGTCGGCTGGGTAGGGCa                                    >  1:2944272/1‑36 (MQ=255)
  gcggcgTTGCTGCACTTGTCGGCTGGGTACGGCA.gg                                 >  1:3140533/1‑35 (MQ=255)
   cggcgTTGCTGCACTTGTCGGCTGGGTACGGCA.ggt                                >  1:738177/1‑34 (MQ=255)
   cggcgTTGCTGCACTTGTCGGCTGGGTACGGCA.ggt                                >  1:2330619/1‑34 (MQ=255)
   cggcgTTGCTGCACTTGTCGGCTGGGTACGGCA.Gtt                                <  1:590902/36‑1 (MQ=255)
     gcgTTGCTGCACTTGTCGGCTGGGTACGGCA.Gtttt                              <  1:827415/36‑1 (MQ=255)
     gcgTTGCTGCACTTGTCGGCTGGGTACGGCA.Gtttt                              <  1:6221937/36‑1 (MQ=255)
     gcgTTGCTGCACTTGTCGGCTGGGTACGGCA.Gtttt                              <  1:6139127/36‑1 (MQ=255)
     gcgTTGCTGCACTTGTCGGCTGGGTACGGCA.Gtttt                              <  1:2099001/36‑1 (MQ=255)
     gcgTTGCTGCACTTGTCGGCTGGGTACGGCA.Gtttt                              <  1:2348666/36‑1 (MQ=255)
     gcgTTGCTGCACTTGTCGGCTGGGTACGGCA.Gtttt                              <  1:3348053/36‑1 (MQ=255)
     gcgTTGCTGCACTTGTCGGCTGGGTACGGCA.Ggttt                              >  1:1373785/1‑36 (MQ=255)
      cgTTGCTGCACTTGTCGGCTGGGTACGGCA.Ggtttt                             >  1:4036058/1‑35 (MQ=255)
      cgTTGCTGCACTTGTCGGCTGGGTACGGCA.GTTTTg                             <  1:726250/36‑1 (MQ=255)
        ttGCTGCACTTGTCGGCTGGGTACGGCA.GTTTTtgg                           >  1:4436113/1‑36 (MQ=255)
        ttGCTGCACTTGTCGGCTGGGTACGGCA.GGTTTtgg                           >  1:554567/1‑36 (MQ=37)
          gctgcACTTGTCGGCTGGGTACGGCAGGTTTTgggt                          >  1:4872520/1‑35 (MQ=37)
          gctgcACTTGTCGGCTGGGTACGGCAGGTTTTgggt                          >  1:2800714/1‑35 (MQ=37)
          gctgcACTTGTCGGCTGGGTACGGCA.GTTTTtgggt                         >  1:6495130/1‑34 (MQ=37)
          gctgcACTTGTCGGCTGGGTACGGCA.GTTTTGGGAt                         <  1:4273013/36‑1 (MQ=255)
           ctgcACTTGTCGGCTGGGTACGGCA.GTTTTGGGGTg                        >  1:3178339/1‑36 (MQ=255)
           ctgcACTTGTCGGCTGGGTACGGCA.GTTTTGGGATg                        <  1:4187015/36‑1 (MQ=255)
           ctgcACTTGTCGGCTGGGTACGGCA.GTTTTGGGATg                        <  1:4298077/36‑1 (MQ=255)
            tgcACTTGTCGGCTGGGTACGGCA.GTTTTGGGATGa                       <  1:6622182/36‑1 (MQ=255)
             gcACTTGTCGGCTGGGTACGGCA.GTTTTTGGGTggg                      >  1:1615551/1‑34 (MQ=25)
             gcACTTGTCGGCTGGGTACGGCA.GTTTTGGGATGAg                      <  1:3094922/36‑1 (MQ=255)
             gcACTTGTCGGCTGGGTACGGCA.GTTTTGGGATGAg                      <  1:5241339/36‑1 (MQ=255)
              cACTTGTCGGCTGGGTACGGCAGGTTTTGGGTTggg                      >  1:4117220/1‑34 (MQ=21)
               aCTTGTCGGCTGGGTACGGCA.GTTTTTGGGTGAGcc                    >  1:7184939/1‑36 (MQ=37)
                cTTGTCGGCTGGGTACGGCA.GTTTTGGGATGAGccc                   <  1:1790817/36‑1 (MQ=255)
                 ttGTCGGCTGGGTACGGCA.GTTTTGGGATGAGccct                  <  1:6717283/36‑2 (MQ=255)
                 ttGTCGGCTGGGTACGGCA.GTTTTGGGATGAGCCCg                  <  1:4117542/36‑1 (MQ=255)
                   gTCGGCTGGGTACGGCA.GTTTTTGGGTGTGGCCggc                >  1:2598916/1‑34 (MQ=16)
                    tCGGCTGGGTACGGCA.GTTTTTTGATGAGGCCGGCa               >  1:1694049/1‑36 (MQ=16)
                    tCGGCTGGGTACGGCA.GTTTTGGGATGAGCCCGACa               <  1:8220894/36‑1 (MQ=255)
                      ggCTGGGTACGGCA.GGTTTGGGGTGAGGCCGGCAGa             >  1:2449206/1‑36 (MQ=255)
                        cTGGGTACGGCA.GGTTTTGGGTGAGCCCGACAGAAt           >  1:2337999/1‑36 (MQ=25)
                         tGGGTACGGCA.GTTTTGGGATGAGCCCGTCAGAAtt          <  1:6476032/36‑1 (MQ=255)
                           ggTACGGCA.GTTTTGGGATGAGCCCGACAGAATTaa        <  1:4991938/36‑1 (MQ=255)
                             tACGGCA.GTTTTGGGATGAGCCCGACAGAATTAAAt      <  1:5274067/36‑1 (MQ=255)
                             tACGGCA.GTTTTGGGATGAGCCCGACAGAATTAAAt      <  1:7076191/36‑1 (MQ=255)
                               cGGCA.GTTTTGGGATGAGCCCGACAGAATTAAATga    <  1:3681826/36‑1 (MQ=255)
                               cGGCA.GTTTTGGGATGAGCACGACAGAATTAAATga    <  1:2799521/36‑1 (MQ=255)
                                 gCA.GTTTTGGGATGAGCCCGACAGAATTAAATgaga  <  1:5149453/36‑1 (MQ=255)
                                    |                                 
GGGCGGCGTTGCTGCACTTGTCGGCTGGGTACGGCA.GTTTTGGGATGAGCCCGACAGAATTAAATGAGA  >  REL606/612374‑612442

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG