Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,266,649 C→A 9.7% intergenic (‑345/‑48) 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,6490C→A 9.7% 5.1 ‑2.8 0.053intergenic (‑345/‑48)acs/nrfAacetyl‑coenzyme A synthetase/nitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  new base (2/4):  ref base (35/12):  total (37/16)
Fisher's exact test strand distribution p-value = 6.02e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTC  >  REL606/4266615‑4266683
                                  |                                  
aTGACAGTGTATGTGCGGTTACTGTCGTGAAAAACa                                   <  1:782258/36‑1 (MQ=255)
aTGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACa                                   >  1:5795566/1‑36 (MQ=255)
aTGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACa                                   >  1:2389429/1‑36 (MQ=255)
aTGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACa                                   >  1:7266178/1‑36 (MQ=255)
aTGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACa                                   >  1:7194884/1‑36 (MQ=255)
 tGACAGTGTAGGTGCGGTTACTTTCGTGAAAAACaa                                  >  1:2129985/1‑36 (MQ=255)
 tGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACaa                                  <  1:400565/36‑1 (MQ=255)
 tGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACaa                                  >  1:6163886/1‑36 (MQ=255)
 tGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACaa                                  <  1:6998740/36‑1 (MQ=255)
  gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                 >  1:5463755/1‑36 (MQ=255)
  gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                 >  1:7596362/1‑36 (MQ=255)
  gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                 >  1:3732881/1‑36 (MQ=255)
  gACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAt                                 <  1:3959198/36‑1 (MQ=255)
    cAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa                               >  1:2285525/1‑36 (MQ=255)
    cAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa                               <  1:4828138/36‑1 (MQ=255)
    cAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATaa                               >  1:7925611/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:2262089‑M2/1‑34 (MQ=255)
      gtgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaa                             <  1:7278335‑M2/36‑3 (MQ=255)
       tgtAGGTGCGGTTACTGTCGTGAAAATCAATataaa                            >  1:3913464/1‑36 (MQ=37)
       tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaaa                            <  1:2036716‑M2/36‑4 (MQ=255)
       tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaaa                            <  1:2729121‑M2/36‑4 (MQ=255)
       tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAaaa                            >  1:3354270‑M2/1‑33 (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:1215817/1‑36 (MQ=255)
                ggTTACTGTCGTGAAAAACAATAAAAACCGCCATTg                   >  1:6793296/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:1741131/1‑36 (MQ=255)
                                aaCAATAAAAACCGCCATTGCAACAATGGCGCAAtt   >  1:1596554/1‑36 (MQ=255)
                                aaCAATAAAAACCGCCATTGCAACAATGGCGCAAtt   >  1:5047552/1‑36 (MQ=255)
                                 aCAATAAAAACCGCCATTGCAACAATGGCGCAATTc  >  1:819376/1‑36 (MQ=255)
                                  |                                  
ATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTC  >  REL606/4266615‑4266683

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 22 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG