Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,513,999 +T 5.4% coding (2743/3741 nt) narZ ← nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,513,9991.→T 5.4% 2.7 ‑1.8 ‑0.037coding (2743/3741 nt)narZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  new base (0/2):  ref base (22/13):  total (22/15)
Fisher's exact test strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.97e-01

GAAATCGATTTCATCCTGGGTATTCCACGAAATCCC.TTT.ACCGCCATTGCCAAGTTTGTCCATCAATGGCC  >  REL606/1513964‑1514034
                                    |                                    
gAAATCGATTTCATCCTGGGTATTCCACGAAATccc                                       <  1:410939/36‑1 (MQ=255)
gAAATCGATTTCATCCTGGGTATTCCACGAAATccc                                       <  1:3537561/36‑1 (MQ=255)
gAAATCGATTTCATCCTGGGGATTCCACGAAATccc                                       >  1:244641/1‑36 (MQ=255)
gAAATCGATTTCATCATGGGTATTCCACGAAATccc                                       >  1:2653122/1‑36 (MQ=255)
 aaaTCGATTTCATCCTGGGTATTCCACGAAATCCC.t                                     <  1:6561766/36‑1 (MQ=255)
 aaaTCGATTTCATCCTGGGTATTCCACGAAATCCC.t                                     <  1:2778448/36‑1 (MQ=255)
  aaTCGATTTCATCCTGGGTATTCCACGAAATCCC.tt                                    >  1:249437/1‑36 (MQ=255)
   aTCGATTTCATCCTGGGTATTCCACGAAATCCC.ttt                                   <  1:6303002/36‑1 (MQ=255)
   aTCGATTTCATCCTGGGTATTCCACGAAATCCC.ttt                                   <  1:7496869/36‑1 (MQ=255)
   aTCGATTTCATCCTGGGTATTCCACGAAATCCC.ttt                                   <  1:4689320/36‑1 (MQ=255)
   aTCGATTTCATCCTGGGTATTCCACGAAATCCC.ttt                                   >  1:7902204/1‑36 (MQ=255)
    tCGATTTCATCCTGGGTATTCCACGAAATCCC.TTT.a                                 >  1:4782529/1‑36 (MQ=255)
    tCGATTTCATCCTGGGTATTCCACGAAATCCC.TTT.a                                 <  1:2927691/36‑1 (MQ=255)
      gATTTCATCCTGGGTATTCCACGAAATCCC.TTT.Acc                               <  1:5889843/36‑1 (MQ=255)
      gATTTCATCCTGGGTATTCCACGAAATCCC.TTT.Acc                               >  1:3967981/1‑36 (MQ=255)
      gATTTCATCCTGGGTATTCCACGAAATACC.TTT.Acc                               >  1:5742603/1‑36 (MQ=255)
       aTTTCATCCTGGGTATTCCACGAAATCCC.TTT.ACCg                              <  1:2336942/36‑1 (MQ=255)
       aTTTCATCCTGGGTATTCCACGAAATCCC.TTT.ACCg                              <  1:5927106/36‑1 (MQ=255)
         ttCATCCTGGGTATTCCACGAAATCCC.TTT.ACCGcc                            >  1:3736079/1‑36 (MQ=255)
         ttCATCCTGGGTATTCCACGAAATCCC.TTT.ACCGcc                            >  1:8034956/1‑36 (MQ=255)
           cATCCTGGGTATTCCACGAAATCCC.TTT.ACCGCCAt                          >  1:3431670/1‑36 (MQ=255)
           aaTCCTGGGCATTCCACGAAATCCC.TTT.ACCGCCAt                          >  1:4469303/2‑36 (MQ=255)
            aTCCTGGGTATTCCACGAAATCCC.TTT.ACCGCCAtt                         >  1:7999934/1‑36 (MQ=255)
            aTCCTGGGTATTCCACGAAATCCC.TTT.ACCGCCAtt                         >  1:8252022/1‑36 (MQ=255)
              cgtGGTTATCCAACGAGATCCC.TTT.ACCGCCATTGc                       <  1:1769769/34‑1 (MQ=16)
              cgggggTATCCCACAAAATCCT.TTT.CCCGCCATTGc                       <  1:111469/33‑1 (MQ=16)
              ccTGGGTATTCCACGAAATCCC.TTT.ACCGCCATTGc                       >  1:7456556/1‑36 (MQ=255)
               cTGGGTATTCCACGAAATCCC.TTT.ACCGCCATTGcc                      >  1:2334340/1‑36 (MQ=255)
               cTGGGTATTCCACGAAATCCC.TTT.ACCGCCATTGcc                      >  1:6282445/1‑36 (MQ=255)
                 ggttATTACACAAAACCTT.TTG.ACCGCCATTGCCaa                    <  1:2979680/33‑1 (MQ=0)
                 gggTATTCCACGAAATCCC.TTT.ATCGCCTTTGCCaa                    >  1:5647734/1‑36 (MQ=37)
                 gggTATTCCAAAAACCCCT.TTT.CCCGCCATTGCCaa                    <  1:3594701/36‑1 (MQ=9)
                 gggTATCCCACAAAACCCT.TTT.ACCGCCATTGCCaa                    <  1:4214684/36‑1 (MQ=19)
                  ggTTTTCCAAAAAATCCCTTTT.ACCGCCATTGCCaa                    <  1:5760674/36‑1 (MQ=9)
                    tATTCCACGAAATCCC.TTT.ACCGCCATTGCCAAGtt                 >  1:907274/1‑36 (MQ=255)
                    tATTCCACGAAATCCC.TTT.ACCGCCATTGCCAAGtt                 >  1:7530732/1‑36 (MQ=255)
                    tATTCCACGAAATCCC.TTT.ACCGCCATTGCCAAGtt                 >  1:3519850/1‑36 (MQ=255)
                        ccACGAAATCCC.TTT.ACCGCCATTGCCAAGTTTGTc             >  1:4450835/1‑36 (MQ=255)
                        ccACGAAATCCC.TTT.ACCGCCATTGCCAAGTTTGTc             >  1:7719911/1‑36 (MQ=255)
                           caaaaTCCT.TTT.ACCGCCATTGCCAAGTTTGTCCAt          <  1:7412645/34‑1 (MQ=37)
                           cGAAATCCC.TTT.ACCGCCATTGCCAAGTTTGTCCAt          >  1:3845411/1‑36 (MQ=255)
                           cGAAATCCC.TTT.ACCGCCATTGCCAAGTTTGTCCAt          >  1:7971697/1‑36 (MQ=255)
                           aaaaaTCCCTTTT.CCCGCCATTGCCAAGTTTGTCCa           <  1:5989076/34‑1 (MQ=21)
                            gAAATCCC.TTT.ACCGCCATTGCCAAGTTTGTCCATc         >  1:7845702/1‑36 (MQ=255)
                            caaaACCC.TTTAACCGCCATTGCCAAGTTTGTCCAt          <  1:760909/35‑1 (MQ=25)
                             aaaCCCC.TTT.ACCGCCATTGCCAAGTTTGTCCATCa        <  1:7033001/36‑1 (MQ=255)
                             aaaCCCC.TTT.ACCGCCATTGCCAAGTTTGTCCATCa        <  1:6853634/36‑1 (MQ=255)
                              accccc.TTT.CCCGCCATTGCCAAGTTTGTCCATCaa       <  1:4953478/33‑1 (MQ=37)
                               aTCCC.TTT.ACCGCCATTGCCAAGTTTGTCCATCAAt      >  1:2270246/1‑36 (MQ=255)
                                 ccc.TTTAACCGCCATTGCCAAGTTTGTCCATCAATg     <  1:7959743/36‑1 (MQ=37)
                                 ccc.TTC.CCCGCCATTGCCAAGTTTGTCCATCAATgg    <  1:2062710/36‑1 (MQ=37)
                                  cc.TTT.ACCGCCATTGCCAAGTTTGTCCATCAATGGc   >  1:7310014/1‑36 (MQ=255)
                                   c.TTT.ACCGCCATTGCCAAGTTTGTCCATCAATGGcc  >  1:4977858/1‑36 (MQ=255)
                                    |                                    
GAAATCGATTTCATCCTGGGTATTCCACGAAATCCC.TTT.ACCGCCATTGCCAAGTTTGTCCATCAATGGCC  >  REL606/1513964‑1514034

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