Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,535,004 G→T 11.8% C84F (TGC→TTC)  focB → predicted formate transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,535,0040G→T11.8% 4.8 ‑4.6 ‑0.042251focBpredicted formate transporter
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (13/24):  total (18/24)
Fisher's exact test strand distribution p-value = 1.01e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.94e-01

TTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACG  >  REL606/2534969‑2535038
                                   |                                  
ttgctttgcACTCGGCTTTATTTTGCTGGCGGTTTg                                    <  1:425698/36‑1 (MQ=255)
ttgctttgcACTCGGCTTTATTTTGCTGGCGGTTTg                                    <  1:2046454/36‑1 (MQ=255)
 tgctttgcACTCGGCTTTATTTTGCTGGCGGTTTGc                                   <  1:6012737/36‑1 (MQ=255)
   ctttgcACTCGGCTTTATTTTGCTGGCGGTTTGCgg                                 <  1:4308127/36‑1 (MQ=255)
    tttgcACTCGGCTTTATTTTGCTGGCGGTTTGCGGc                                <  1:7201571/36‑1 (MQ=255)
    tttgcACTCGGCTTTATTTTGCTGGCGGTTTGCGGc                                <  1:3444840/36‑1 (MQ=255)
       gcACTCGGCTTTATTTTGCTGGCGGTTTGCGGCAcc                             <  1:6258640/36‑1 (MQ=255)
         aCTCGGCTTTATTTTGCTGGCGGTTTTCGGGAcccg                           >  1:4252514/1‑34 (MQ=18)
           tCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGcc                         <  1:7673072/36‑1 (MQ=255)
           tCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGcc                         <  1:5844039/36‑1 (MQ=255)
            cGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCt                        <  1:3587206/36‑1 (MQ=255)
            cGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCt                        <  1:7485327/36‑1 (MQ=255)
              gCTTTATTTTGCTGGCGGTTTTTTGCACTAGccata                      >  1:7359571/1‑33 (MQ=9)
              gCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGt                      <  1:517871/36‑1 (MQ=255)
              gCTTTATTTTGCCGGCGGTTTGCGGCACCAGCCTGt                      <  1:2855350/36‑1 (MQ=255)
                tttATTTTGCTGGCGGTTTTTGTCACCAGCCTTttt                    >  1:6903336/1‑35 (MQ=12)
                tttATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTc                    <  1:7422921/36‑1 (MQ=255)
                  tATTTTGCTGGCGGTTTTTGGGACCAGCCTTTTTAc                  >  1:6433864/1‑36 (MQ=9)
                   aTTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCAcc                 <  1:3214022/36‑1 (MQ=255)
                     tttGCTGGCGGTTTTCGGCACCAGCCTGTTTACCtc               >  1:4454855/1‑36 (MQ=37)
                       tGCGGGCGGTTTGCGTCCCCAGCCTGTTCCCCtcgt             <  1:7810895/36‑1 (MQ=19)
                        gCTGGGGTTTTGCGGCACCAGACTGTTCACCtcgtc            <  1:560910/36‑1 (MQ=21)
                        gCTGGCGGTTTGCGGCACCAGCCTGTTCACCtcgtc            <  1:2433685/36‑1 (MQ=255)
                        gCTGGCGGTTTGCGGCACCAGCCTGTTCACCtcgtc            <  1:2565204/36‑1 (MQ=255)
                         cTGGCGGTTTGCGGCACCAGCCTGTTCACCtcgtcg           <  1:5628230/36‑1 (MQ=255)
                         cTGGCGGTTTGCGGCACCAGCCTGTTCACCtcgtcg           >  1:6123060/1‑36 (MQ=255)
                          tGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCgg          <  1:1802043/36‑1 (MQ=255)
                          tGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCgg          <  1:425405/36‑1 (MQ=255)
                           ggCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGt         <  1:2217039/36‑1 (MQ=255)
                           ggCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGt         <  1:2856074/36‑1 (MQ=255)
                           ggCGGTTTGCGGCACCAGCCGGTTTACCTCGTCTGt         >  1:389259/1‑36 (MQ=25)
                            gCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTa        >  1:5439213/1‑36 (MQ=255)
                            gCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTa        >  1:5183248/1‑36 (MQ=255)
                             cGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTaa       >  1:582157/1‑36 (MQ=255)
                             cGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTaa       >  1:134449/1‑36 (MQ=255)
                              ggTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAAt      >  1:5095765/1‑36 (MQ=255)
                              ggTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAAt      >  1:2407271/1‑36 (MQ=255)
                               gTTTGCGGCACCAGCCTGTTCACCTCGTCGGTGATg     >  1:6810212/1‑36 (MQ=255)
                               gTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATg     >  1:2412316/1‑36 (MQ=255)
                               gTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATg     >  1:4271114/1‑36 (MQ=255)
                                 ttGCGGCACCAGCCTGTTCACCTCGTCGGTAATGAc   <  1:835661/36‑1 (MQ=255)
                                 ttGCGGCACCAGCCTGTTCACCTCGTCGGTAATGAc   >  1:1263222/1‑36 (MQ=255)
                                  tGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACg  <  1:6927567/36‑1 (MQ=255)
                                  tGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACg  >  1:4711790/1‑36 (MQ=255)
                                   |                                  
TTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACG  >  REL606/2534969‑2535038

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 19 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG