Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,524,719 T→A 7.4% S37T (TCC→ACC)  fdnH → formate dehydrogenase‑N, Fe‑S (beta) subunit, nitrate‑inducible

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,524,7190T→A 7.4% 5.8 ‑3.4 ‑0.054109fdnHformate dehydrogenase‑N, Fe‑S (beta) subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (20/30):  total (24/30)
Fisher's exact test strand distribution p-value = 3.36e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.54e-01

ATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGT.AAAGCCTGTCAGGTGGCG  >  REL606/1524684‑1524754
                                   |                                    
aTTACAAAGCAGAAGTCGCAAAACTTATCGACGttt                                      >  1:8337/1‑36 (MQ=255)
aTTACAAAGCAGAAGTCGCAAAACTTATCGACGttt                                      >  1:1358155/1‑36 (MQ=255)
  taaaaaGCAGAAGTCGCTAAACTTATCGACGTTTcc                                    <  1:5911187/33‑1 (MQ=38)
  tACAAAGCAGAAGTCGCAAAACTTATCGACGTTTcc                                    <  1:7652710/36‑1 (MQ=255)
  tACAAAGCAGAAGTCGCAAAACTTATCGACGTTTcc                                    >  1:5008019/1‑36 (MQ=255)
  tACAAAGCAGAAGGCGCAAAACTTATCGACGtttac                                    >  1:1397292/1‑34 (MQ=38)
   aCAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCa                                   <  1:2343752/36‑1 (MQ=255)
    cAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCAc                                  <  1:3979957/36‑1 (MQ=255)
     aaaGCAGAAGTCGCAAAACTTATCGACGTTTCCAcc                                 <  1:4430199/36‑1 (MQ=255)
      aaGCAGAAGTCGCAAAACTTATCGACGTTTCCACCt                                <  1:2342097/36‑1 (MQ=255)
      aaGCAGAAGTCGCAAAACTTATCGACGTTTCCACCt                                >  1:1271256/1‑36 (MQ=255)
       aGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTg                               <  1:6515724/36‑1 (MQ=255)
       aGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTg                               >  1:5886344/1‑36 (MQ=255)
        gCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGt                              >  1:7498719/1‑36 (MQ=255)
         cAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTa                             <  1:6743106/36‑1 (MQ=255)
          aGAAGTCGCAAAACTTATCGACGTTTCCACCTGTAt                            <  1:7100347/36‑1 (MQ=255)
           gAAGTCGCAAAACTTATCGACGTTTCCACCTGTATc                           >  1:4315616/1‑36 (MQ=255)
             aGTCGCAAAACTTATCGACGTTTCCACCTGTATCgg                         <  1:5208451/36‑1 (MQ=255)
              gTCGCAAAACTTATCGACGTTTCCACCTGTATCGGc                        <  1:5667992/36‑1 (MQ=255)
              gTCGCAAAACTTATCGACGTTTCCACCTGTATCGGc                        <  1:548678/36‑1 (MQ=255)
               tCGCAAAACTTATCGACGTTTCCACCTGTATCGGCt                       <  1:2402946/36‑1 (MQ=255)
               tCGAAAAACTTATCGACGTTTCCACCTGTATCGGCt                       <  1:4489421/36‑1 (MQ=255)
               acGCGAAGCTTATCGACGTTTCCACCTATATCGGCt                       <  1:5902956/35‑1 (MQ=14)
                cGCAAAACTTATCGACGTTTCCACCTGTATCGGCTg                      <  1:4049262/36‑1 (MQ=255)
                 gCAAAACTTATCGACGTTTCCACCTGTATCGGCTgg                     >  1:100201/1‑35 (MQ=255)
                 gCAAAACTTATCGACGTTTCCACCTGTATCGGCTGt                     <  1:1521721/36‑1 (MQ=255)
                 gCAAAACTTATCGACGTTTCCACCTGTATCGGCTGt                     >  1:3359224/1‑36 (MQ=255)
                  cAAAACTTATCGACGTTTCCACCTGTATCGGCTGtt                    >  1:6915303/1‑35 (MQ=255)
                   aaaaCTTATCGACGTTTCCACCTGTATCGGCTGT.aa                  <  1:1397198/36‑1 (MQ=255)
                    aaaCTTATCGACGTTTCCACCTGTATCGGCTGT.taa                 >  1:7484211/1‑36 (MQ=255)
                    aaaCTTATCGACGTTTCCACCTGTATCGGCTGT.aaa                 <  1:1478245/36‑1 (MQ=255)
                    aaaCTTATCGACGTTTCCACCTGTATCGGCTGT.aaa                 <  1:1837302/36‑1 (MQ=255)
                     aaCTTATCGACGTTTCCACCTGTATCGGCTGT.TAAg                >  1:6623378/1‑36 (MQ=255)
                     aaCTTATCGACGTTTCCACCTGTATCGGCTGT.AAAg                <  1:2512765/36‑1 (MQ=255)
                     aaCTTATCGACGTTTCCACCTGTATCGGCTGT.AAAg                <  1:61850/36‑1 (MQ=255)
                      aCTTATCGACGTTTCCACCTGTATCGGCTGT.AAAGc               <  1:3758079/36‑1 (MQ=255)
                      aCTTATCGACGTTTCCACCTGTATCGGCTGG.TAAGc               >  1:2142534/1‑36 (MQ=17)
                      aCTCATCGACGTTACCACCTGTATCGGCTGT.taaac               >  1:2713/1‑34 (MQ=17)
                       cTTATCGACGTTTCCACCTGTATCGGCTGT.AAAGcc              <  1:6732450/36‑1 (MQ=255)
                       cTTATCGACGTTTCCACCTGTATCGGCTGT.AAAGcc              <  1:1683190/36‑1 (MQ=255)
                        tcaTCGACGTTACCACCTGTATCGGCTGTGAAAGct              >  1:5418475/3‑35 (MQ=17)
                        tcaTCGACGTTACCACCTGTATCGGCTGTGAAAGct              >  1:3784773/3‑35 (MQ=17)
                         tATCGACGTTTCCACCTGTATCGGCTGT.AAAGCCTg            <  1:7067727/36‑1 (MQ=40)
                          aTCGACGTTTCCACCTGTATCGGCTGT.AAAGCCTGt           <  1:7462499/36‑1 (MQ=38)
                           tCGACGTTTCCACCTGTATCGGCTGT.TAAGCCTGTc          >  1:3870989/1‑36 (MQ=255)
                           tCGACGTTTCCACCTGTATCGGCTGT.AAAGCCTGTc          <  1:1225770/36‑1 (MQ=38)
                           tCGACGTTTCCACCTGTATCGGCTGG.GAAGCCTTTc          >  1:5763801/1‑36 (MQ=14)
                            cggcGTTTCCACCTGTATCTGCCGT.AAAGCCTGTCa         <  1:1572591/33‑1 (MQ=14)
                            cGACGTTTCCACCTGTATCGGCTGT.AAAGCCTTTct         >  1:4968419/1‑35 (MQ=255)
                             gACGTTTCCACCTGTATCGGCTGT.AAAAGCTGTCAg        >  1:3812951/1‑36 (MQ=14)
                             gACGTTACCACCTGTATCGGCTGT.CAATGCTGTGAg        >  1:7817899/1‑36 (MQ=9)
                              aCGTTTCCACCTGTATCGGCTGT.AAAACCTTTCCgg       >  1:3415054/1‑36 (MQ=14)
                               cGTTTCCACCTGTATCGGCTGT.AAAGCCTGTCAGGt      <  1:3733389/36‑1 (MQ=38)
                                gTTTCCACCTGTATCGGCTGT.TAAGACTGTCAgggg     >  1:6251438/1‑34 (MQ=14)
                                gTTTCCACCTGTATCGGCTGT.AAAGCCTTTCAgggg     >  1:169607/1‑34 (MQ=14)
                                gTTTCCACCTGTATCGGCTGT.AAAACCTGTCAgggg     >  1:6249114/1‑34 (MQ=14)
                                gTTTCCACCTGTATCGGCTGG.TAAGCCTGTTAgggg     >  1:6177761/1‑34 (MQ=14)
                                 tttCCTCCTGTATCGGCTGT.AAAGCCTGTCAGGTgg    <  1:3721980/36‑1 (MQ=255)
                                 tttCCACCTGTATCGGCTGT.AAAGCCTGTCAGGTgg    <  1:2435315/36‑1 (MQ=38)
                                 tttCCACCTGTATCGGCTGT.AAAGCCTGTCAGGTgg    <  1:3467384/36‑1 (MQ=38)
                                  ttCCACCTGTATCGGCTGT.AAAGCCTGTCAGGTGGc   <  1:3600150/36‑1 (MQ=35)
                                  ttCCACCTGTATCGGCTGT.AAAGCCTGTCAGGTGGc   >  1:2761321/1‑36 (MQ=35)
                                   tCCACCTGTATCGGCTGT.AAAACCTGTCAGGTgggg  >  1:5035737/1‑34 (MQ=11)
                                   |                                    
ATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGT.AAAGCCTGTCAGGTGGCG  >  REL606/1524684‑1524754

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 18 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG