Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,266,652 T→A 7.7% intergenic (‑348/‑45) acs ← / → nrfA acetyl‑coenzyme A synthetase/nitrite reductase, formate‑dependent, cytochrome

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,266,6520T→A 7.7% 2.4 ‑3.7 ‑0.039intergenic (‑348/‑45)acs/nrfAacetyl‑coenzyme A synthetase/nitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  new base (0/3):  ref base (27/9):  total (27/12)
Fisher's exact test strand distribution p-value = 2.41e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.62e-01

GACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGAT  >  REL606/4266617‑4266687
                                   |                                   
gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                     >  1:4766895/1‑36 (MQ=255)
 aCAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATa                                    <  1:936185/36‑1 (MQ=255)
 aCAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATa                                    >  1:1092498/1‑36 (MQ=255)
  cAGTGTAGGTGCGGTTACTTTCGTGAAAAACAATaa                                   <  1:3511098/36‑1 (MQ=255)
  cAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa                                   <  1:4485273/36‑1 (MQ=255)
  cAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa                                   <  1:3836543/36‑1 (MQ=255)
   aGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaaa                                  >  1:3516446/1‑36 (MQ=255)
   aGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaaa                                  >  1:1514912/1‑36 (MQ=255)
   aGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaaa                                  <  1:150647/36‑1 (MQ=255)
    gtgtAGGTGCGGTTACTGTCGTGAAAAACAATacac                                 >  1:1031598/1‑33 (MQ=255)
    gtgtAGGTGCGGTTACTGTCGTGAAAAACAATaaaa                                 <  1:2082630/36‑1 (MQ=255)
    gtgtAGGGGCGGTTACTGTCGTGAAAAACAATaaaa                                 >  1:3300189/1‑36 (MQ=255)
      gtAGGTGCGGTTACTGTCGTTAAAAACAATAAAAAc                               <  1:1351959/36‑1 (MQ=255)
       tAGGTGCGGTTACTGTCGTGAAAAACAATAAAAAcc                              <  1:4924012/36‑1 (MQ=255)
       tAGGTGCGGTTACTGTCGTGAAAAACAATAAAAAcc                              <  1:4514062/36‑1 (MQ=255)
        aGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCg                             >  1:1387587/1‑36 (MQ=255)
        aGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCg                             >  1:802081/1‑36 (MQ=255)
          gTGCGGTTACTGTCGTGAAAAACAATAAAAACCGcc                           >  1:1230607/1‑36 (MQ=255)
           tGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCa                          >  1:5237752/1‑36 (MQ=255)
               gTTACTGTCGTGAAAAACAATAAAAACCGCCATTGc                      >  1:499179/1‑36 (MQ=255)
                  aCTGTCGTGAAAAACAATAAAAACCGCCATTGcaac                   >  1:5565534/1‑36 (MQ=255)
                    tGTCGTGAAAAACAATAAAAACCGCCATTGcaacaa                 >  1:1708651/1‑36 (MQ=255)
                      tCGTGAAAAACAATAAAAACCGCCATTGCAACAATg               >  1:2144126/1‑36 (MQ=255)
                       cGTGAAAAACAATAAAAACCGCCATTGCAACAATgg              >  1:1368058/1‑36 (MQ=255)
                        gTGAAAAACAATAAAAACCGCCATTGCAACAATGgc             >  1:1085601/1‑36 (MQ=255)
                        aaaaaaaaAAAAAAAACCCGCCATTGCAACAATGgc             <  1:4577481/33‑1 (MQ=16)
                          gAAAAACAATAAAAACCGCCATGGCAACAATGgcgc           >  1:3607579/1‑36 (MQ=255)
                           aaaaaCAATAAAAACCGCCATTGCAACAATGGCGCa          >  1:15115/1‑36 (MQ=255)
                           aaaaaCAATAAAAACCGCCATTGCAACAATGGCGCa          >  1:2846415/1‑36 (MQ=255)
                           aaaaaAAAAAAAAACCGCCATTGCAACAATGGCGCa          <  1:243449/36‑1 (MQ=37)
                            aaaaCAATAAAAACCGCCATTGCAACAATGGCGCaa         >  1:3685319/1‑36 (MQ=255)
                            aaaaAAAAAAAAACCGCCATTGCAACAATGGCGCaa         <  1:1115865/36‑1 (MQ=37)
                              aaCAATAAAAACCGCCATTGCAACAATGGCGCAAtt       >  1:4731755/1‑36 (MQ=255)
                              aaCAATAAAAACCGCCATTGCAACAATGGCGCAAtt       >  1:2460680/1‑36 (MQ=255)
                               aCAATAAAAACCGCCATTGCAACAATGGCGCAATTc      >  1:4654501/1‑36 (MQ=255)
                               aCAATAAAAACCGCCATTGCAACAATGGCGCAATTc      >  1:4370186/1‑36 (MQ=255)
                                cAATAAAAACCGCCATCGCAACAATGGCGCAATTCg     >  1:2218423/1‑36 (MQ=255)
                                 aaTAAAAACCGCCATTGCAACAATGGCGCAATTCgg    >  1:2692569/1‑36 (MQ=255)
                                  aTAAAAACCGCCATTGCAACAATGGCGCAATTCGGa   >  1:1452443/1‑36 (MQ=255)
                                  aTAAAAACCGCCATTGCAACAAGGGCGCAATACgg  >  1:802492/1‑35 (MQ=37)
                                   tAAAAACCGCCATTGCAACAATGGCGCAATTCGGAt  >  1:4733469/1‑36 (MQ=255)
                                   tAAAAACCGCCATTGCAACAATGGCGCAATTCGGAt  >  1:4968207/1‑36 (MQ=255)
                                   tAAAAACCGCCATTGCAACAATGGCGCAACTCGGAt  >  1:1292099/1‑36 (MQ=255)
                                   |                                   
GACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGAT  >  REL606/4266617‑4266687

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