Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,266,652 T→A 11.3% 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→A11.3% 5.5 ‑4.1 ‑0.044intergenic (‑348/‑45)acs/nrfAacetyl‑coenzyme A synthetase/nitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  new base (1/4):  ref base (30/8):  total (32/12)
Fisher's exact test strand distribution p-value = 1.68e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.94e-01

GACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGA  >  REL606/4266617‑4266686
                                   |                                  
gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                    >  1:7596362/1‑36 (MQ=255)
gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                    <  1:3959198/36‑1 (MQ=255)
gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                    >  1:5463755/1‑36 (MQ=255)
gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                    >  1:3732881/1‑36 (MQ=255)
  cAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa                                  >  1:7925611/1‑36 (MQ=255)
  cAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa                                  <  1:4828138/36‑1 (MQ=255)
  cAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa                                  >  1:2285525/1‑36 (MQ=255)
  cAGTGTAAGTGCGGTTACTGTCGTGAAAAACAATaa                                  <  1:5407866/36‑1 (MQ=255)
   agggtAGGTGCGGTTACTGTCGTGTAAAAAAATaaa                                 >  1:3269903/4‑36 (MQ=25)
   aGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAa                                 >  1:1397330‑M2/1‑35 (MQ=255)
    gtgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaa                                <  1:7278335‑M2/36‑3 (MQ=255)
    gtgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaa                                >  1:2262089‑M2/1‑34 (MQ=255)
     tgtAGGTGCGGTTACTGTCGTGAAAATCAATataaa                               >  1:3913464/1‑36 (MQ=37)
     tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaaa                               <  1:2729121‑M2/36‑4 (MQ=255)
     tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaaa                               >  1:3354270‑M2/1‑33 (MQ=255)
     tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaaa                               <  1:2036716‑M2/36‑4 (MQ=255)
      gtAGGTGCGGTTACTGTTGTGAAAAACAATAAaaac                              >  1:6458527‑M2/1‑32 (MQ=255)
      gtAGGCGCGCTTACTGTAGTGAAAAACAATAAAAAc                              <  1:4256303/36‑1 (MQ=25)
       tAGGTGCGGTTACTGTCGTGAAAAACAATAAaaacc                             >  1:5475810‑M2/1‑31 (MQ=255)
       tAGGTGCGGTTACTGTCGTGAAAAACAATAAaaacc                             <  1:1786231‑M2/36‑6 (MQ=255)
          gggCGGTTACTGCAGTGAAAAACAATAAAAACCGcc                          <  1:6051002/34‑1 (MQ=37)
           tGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCa                         >  1:4912596/1‑36 (MQ=255)
             cGGTTACTGTCGTGAAAAACAATAAAAACCGCCAtt                       >  1:868948/1‑36 (MQ=255)
             cGGTTACTGTCGTGAAAAAAAATAAAAACCGCCAtt                       >  1:61833/1‑36 (MQ=255)
              ggTTACTGTCGTGAAAAACAATAAAAACCGCCATTg                      >  1:6793296/1‑36 (MQ=255)
              ggTTACTGTCGTGAAAAACAATAAAAACCGCCATTg                      >  1:1215817/1‑36 (MQ=255)
               gTTACTGTCGTAAAAAACAATAAAAACCGCCATTGc                     >  1:3563058/1‑36 (MQ=255)
                ttACTGTCGTGAAAAACAATAAAAACCGCCATTGca                    >  1:1563161/1‑36 (MQ=255)
                      tCGTGAAAAACAAGAAAAACCGCCATTGCAACAATg              >  1:2285148/1‑36 (MQ=255)
                        gTGAAAAACAATAAAAACCGCCATTGCAACAACGgc            >  1:4217963/1‑36 (MQ=255)
                        aaaaaaaaAAAAAAAACCCGCCATTGCAACAATGgc            <  1:5309197/33‑1 (MQ=16)
                         tGAAAAACAATAAAAACCGCCATTGCAACAATGgcg           >  1:7850141/1‑36 (MQ=255)
                         aaaaaaaAAAAAAAACCCGCCATTGCAACAATGgcg           <  1:122419/34‑1 (MQ=21)
                          gAAAAACAATAAAAACCGCCATTGCAACAATGgcgc          >  1:7584869/1‑36 (MQ=255)
                          aaaaaaAAAAAAAAACCGCCATTGCAACAATGgcgc          <  1:3108806/35‑1 (MQ=25)
                           aaaaaCAATAAAAACCGCCATTGCAACAATGGCGCa         >  1:2059766/1‑36 (MQ=255)
                            aaaaCAATAAAAACCGCCATTGCAACAATGGCGCaa        >  1:1845099/1‑36 (MQ=255)
                            aaaaCAAAAAaaaccgccgcagattggtatggcaaa        >  1:4908588‑M2/1‑10 (MQ=255)
                             aaaCAATAAAAACCGCCCTTGCAACAATGCCGCaac       >  1:6462075/1‑35 (MQ=25)
                             aaaAAAAAAAAACCGCCATTGCAACAATGGCGCAAt       <  1:3095111/36‑1 (MQ=37)
                              aaCAATAAAAACCGCCATTGCAACAATGGCGCAAtt      >  1:1596554/1‑36 (MQ=255)
                              aaCAATAAAAACCGCCATTGCAACAATGGCGCAAtt      >  1:1741131/1‑36 (MQ=255)
                              aaCAATAAAAACCGCCATTGCAACAATGGCGCAAtt      >  1:5047552/1‑36 (MQ=255)
                               aCAATAAAAACCGCCATTGCAACAATGGCGCAATTc     >  1:819376/1‑36 (MQ=255)
                                 aaTAAAAACCGCCATTGCAACAATGGCGCAATTCgg   >  1:6620795/1‑36 (MQ=255)
                                 aaTAAAAACCGCCATTGCAACAATGGCGCAATTCgg   >  1:1433812/1‑36 (MQ=255)
                                 aaTAAAAACCGCCATTGCAACAATGGCGCAATTCgg   >  1:6525060/1‑36 (MQ=255)
                                  aTAAAAACCGCCATTGCAACAATGGCGCAATTTGGa  >  1:4116640/1‑36 (MQ=255)
                                   |                                  
GACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGA  >  REL606/4266617‑4266686

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG