Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,524,736 +G 4.3% coding (126/885 nt) 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,7361.→G 4.3% 2.9 ‑1.7 ‑0.246coding (126/885 nt)fdnHformate dehydrogenase‑N, Fe‑S (beta) subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (18/26):  total (20/26)
Fisher's exact test strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.41e-01

GCAAAACTTATCGACGTTTCCACCTGTATCGGCTGT.AAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGA  >  REL606/1524701‑1524771
                                    |                                   
gCAAAACTTATCGACGTTTCCACCTGTATCGGCTGt                                      >  1:3359224/1‑36 (MQ=255)
gCAAAACTTATCGACGTTTCCACCTGTATCGGCTGt                                      <  1:1521721/36‑1 (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:61850/36‑1 (MQ=255)
    aaCTTATCGACGTTTCCACCTGTATCGGCTGT.AAAg                                 <  1:2512765/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:1683190/36‑1 (MQ=255)
      cTTATCGACGTTTCCACCTGTATCGGCTGT.AAAGcc                               <  1:6732450/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:2761321/1‑36 (MQ=35)
                 ttCCACCTGTATCGGCTGT.AAAGCCTGTCAGGTGGc                    <  1:3600150/36‑1 (MQ=35)
                  tCCACCTGTATCGGCTGT.AAAACCTGTCAGGTgggg                   >  1:5035737/1‑34 (MQ=11)
                        tGTATTGGCTGT.AAATCCTGCGTGGTGGCCTTTTct             >  1:6769911/1‑35 (MQ=37)
                          tATCGGCTGT.AAAGCCTGTCAGGTGGCGTGTTCGGa           <  1:7910732/36‑1 (MQ=38)
                          tATCGGCTGT.AAAGCCTGTCAGGTGGCGTGTTCGGa           <  1:2262081/36‑1 (MQ=38)
                           aTCGGCTGT.AAAGCCTGTCAGGTGGCGTGTTCGGAg          <  1:1792758/36‑1 (MQ=38)
                           aTCGGCTGT.AAAGCCTGTCAGGTGGCGTGTTCGGAg          <  1:212006/36‑1 (MQ=38)
                            tCGGCTGT.AAAGCCTGTCAGGTGGCGTGTTCGGAGt         <  1:7706261/36‑1 (MQ=38)
                            tCGGCTGT.AAAGCCTGTCAGGTGGCGTGTTCGGAGt         <  1:1750570/36‑1 (MQ=38)
                            tCGGCTGT.AAAGCCTGTCAGGTGGCGTGTTCGGAGt         <  1:1749406/36‑1 (MQ=38)
                            tCGGCTGT.AAAGCCTGTCAGGTGGCGTGTTCGGAGt         >  1:4728363/1‑36 (MQ=38)
                             cGGCTGT.AAAGCCTGTCAGGTGGCGTGTTCGGAGTg        <  1:7694211/36‑1 (MQ=38)
                             cGGCTGT.AAAGCCTGTCAGGTGGCGTGTTCGGAGTg        <  1:7837006/36‑1 (MQ=38)
                                 tGT.AAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAAc    <  1:6072162/36‑1 (MQ=38)
                                 tGT.AAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAAc    <  1:6881447/36‑1 (MQ=38)
                                  gT.AAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACg   <  1:1801276/36‑1 (MQ=38)
                                  gT.AAAGCCTGTCAGGCGGCGTGTTCGGAGTGGAACg   >  1:6324778/1‑36 (MQ=255)
                                  gT.AAAGCCTGTCAGGCGGCGTGGTCGGAGTGGAACg   >  1:7289424/1‑36 (MQ=14)
                                   t.AAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGa  >  1:3701783/1‑36 (MQ=38)
                                   t.AAAGCCTGTCAGGTGGCGTGTTCGGAGGGGAACgc  >  1:2988443/1‑35 (MQ=255)
                                    |                                   
GCAAAACTTATCGACGTTTCCACCTGTATCGGCTGT.AAAGCCTGTCAGGTGGCGTGTTCGGAGTGGAACGA  >  REL606/1524701‑1524771

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