Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,253,643 +G 6.1% intergenic (‑157/‑251) napF ← / → eco ferredoxin‑type protein, predicted role in electron transfer to periplasmic nitrate reductase (NapA)/ecotin precursor

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,253,6431.→G 6.1% 2.6 ‑2.7 ‑0.033intergenic (‑157/‑251)napF/ecoferredoxin‑type protein, predicted role in electron transfer to periplasmic nitrate reductase (NapA)/ecotin precursor
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (7/24):  total (9/24)
Fisher's exact test strand distribution p-value = 6.82e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.65e-01

GTAAAACCTGTAGCAGGATCAAAAGAGTGGC.GGGAAGCGTGGCAAAAACGGGCTTTTGCTCACA  >  REL606/2253613‑2253676
                               |                                 
gTAAAACCTGTAGCAGGATCAAAAGAGTGGC.ggggg                              >  1:5390766/1‑34 (MQ=255)
  aaaaCCTGTAGCAGGATCAAAAGAGTGGC.GGGAAGc                            <  1:6147739/36‑1 (MQ=255)
  aaaaCCTGTAGCAGGATCAAAAGAGTGGC.GGGAAGc                            <  1:5732647/36‑1 (MQ=255)
   aaaCCTGTAGCAGGATCAAAAGAGTGGC.GGGGAggg                           >  1:5038741/1‑34 (MQ=38)
     aCCTGTAGCAGGATCAAAGGAGTCGC.GGGAAGCGTg                         <  1:4751570/36‑1 (MQ=37)
     aCCTGTAGCAGGATCAAAAGAGTGGC.GGGAAGCGTg                         <  1:3226121/36‑1 (MQ=255)
      ccTGTAGCAGGATCAAAAGAGTGGC.GGGAAGCGTgg                        <  1:7939384/36‑1 (MQ=255)
       ctctAGCAGGGTCAAAAGAGTGGC.GGGAAGCGTGGc                       <  1:1234065/33‑1 (MQ=38)
        tGTAGCAGGATCAAAAGAGTGGC.GGGGAGAGTgggg                      >  1:3198162/1‑34 (MQ=37)
         gTAGCAGGATCAAAAGAGTGGC.GGGAAGCGTGGCaa                     <  1:1499316/36‑1 (MQ=255)
          tAGCAGGATCAAAAGAGTGGC.GGGAAGCGTGGCaaa                    <  1:2363672/36‑1 (MQ=255)
           aGCAGGATCAAAAGAGTGGC.GGGAAGCGTGGCaaaa                   <  1:7380922/36‑1 (MQ=255)
              aGGATCAAAAGAGTGGC.GGGGAGGGGGGGaaaaaag                >  1:525350/1‑34 (MQ=16)
              aGGATCAAAAGAGTGGC.GGGAAGCGTGGCAAAAACg                <  1:576311/36‑1 (MQ=255)
               ggATCAAAAGAGTGGCGGGGAAGGGGGGAaaaaaag                >  1:5556268/1‑34 (MQ=12)
               ggATCAAAAGAGAGGCGGGGAAGGGAGGGaaaaagg                >  1:5137773/1‑34 (MQ=2)
                gATCAAAAGAGTGGC.GGGAAGCGTGGCAAAAACggg              <  1:2831963/36‑1 (MQ=255)
                gATCAAAAGAGTGGC.GGGAAGCGTGGCAAAAACggg              <  1:3938473/36‑1 (MQ=255)
                gATCAAAAGAGTAAC.GGGAAGCGTGGCAAAAACggg              <  1:6759277/36‑1 (MQ=37)
                 aTCAAAAGAGTGGC.GTGAATCTTGGCAAAAACGGGc             <  1:2863248/36‑1 (MQ=21)
                 aTCAAAAGAGTGGC.GGGAAGCGTGGCAAAAACGGGc             <  1:6307025/36‑1 (MQ=255)
                  tCAAAAGAGTGGC.GGGAAGCGGGGCAAAAAAggggt            >  1:7690411/1‑34 (MQ=25)
                      aaGAGTGGC.GGGGAGCGTGGGAAAAACGGGGttttt        >  1:669359/1‑35 (MQ=25)
                      aaGAGTGGC.GGGAAGCGTGGCAAAAACGGGCTTTTg        <  1:64759/36‑1 (MQ=255)
                      aaGAGTGGC.GGGAAGCGTGGCAAAAACGGGCTTTTg        <  1:7395272/36‑1 (MQ=255)
                       agagTGGC.GGGAAGCGTGGCAAAAACGGGCTTTTGc       <  1:2488083/36‑1 (MQ=255)
                       agagTGGC.GGGAAGCGTGGCAAAAACGGGCTTTTGc       <  1:329048/36‑1 (MQ=255)
                          gTGGC.GGGAAGCGTGGCAAAAACGGGCTTTTGCTca    <  1:1554190/36‑1 (MQ=255)
                          gTGGC.GGGAAGCGTGGCAAAAACGGGCTTTTGCTca    <  1:5737606/36‑1 (MQ=255)
                          gTGGC.GGGAAGCGTGGCAAAAACGGGCTTTTGCTca    <  1:5743550/36‑1 (MQ=255)
                          gTGGC.GGGAAGCGTGGCAAAAACGGGCTTTTGCTca    <  1:2887881/36‑1 (MQ=255)
                            ggC.GGGAAGCGTGGCAAAAACGGGCTTTTGCTcaca  >  1:1933389/1‑36 (MQ=255)
                            ggC.GGGAAGCGTGGCAAAAACGGGCTTTTGCTcaca  <  1:5249617/36‑1 (MQ=255)
                            ggC.GGGAAGCGTGGCAAAAACGGGCTTTTGCTcaca  >  1:6766942/1‑36 (MQ=255)
                               |                                 
GTAAAACCTGTAGCAGGATCAAAAGAGTGGC.GGGAAGCGTGGCAAAAACGGGCTTTTGCTCACA  >  REL606/2253613‑2253676

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 18 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG