Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,069,853 G→T 9.3% G202G (GGC→GGA pitB ← phosphate transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,069,8530G→T 9.3% 2.0 ‑4.5 0.049606pitBphosphate transporter
Reads supporting (aligned to +/- strand):  new base (5/1):  ref base (11/32):  total (16/33)
Fisher's exact test strand distribution p-value = 1.09e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CGCAATACGCGTCCAGAATGGCGG.TTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTG  >  REL606/3069820‑3069888
                                  |                                   
cgcAATACGCGTCCAGAATGGCGG.TTTACGTTTGcc                                   <  1:5404221/36‑1 (MQ=255)
 gcAATACGCGTCCAGAATGGCGG.TTTACGTTTGCCt                                  <  1:7504612/36‑1 (MQ=255)
  cAATACGCGTCCAGAATGGCGG.TTTACGTTTGCCtt                                 <  1:5237793/36‑1 (MQ=255)
   aaTACGCGTCCAGAATGGCGG.TTTACGTTTGCCttt                                <  1:3392821/36‑1 (MQ=255)
    aTACGCGTCCAGAATGGCGG.TTTCCGTTTGCCtttt                               <  1:5627722/36‑1 (MQ=255)
    aTACGCGTCCAGAATGGCGG.TTTACGTTTGCCtttt                               <  1:2538853/36‑1 (MQ=255)
    aTACGCGTCCAGAATGGCGG.TTTACGTTTGCCtttt                               <  1:390158/36‑1 (MQ=255)
     tACGCGTCCAGAATGGCGG.TTTACGTTTGCCttttt                              <  1:6687599/36‑1 (MQ=255)
      aCGCGTCCAGAATGGCGG.TTTACGTTTGCCtttttt                             <  1:3707255/36‑1 (MQ=255)
      aCGCGTCCAGAATGGCGG.GTTTTGTTTTTTtttttt                             >  1:4049199/1‑36 (MQ=2)
         cgTCCAGAATGGCGG.GTTTTGTTTTTTTTTTTTTtt                          >  1:4758178/1‑36 (MQ=2)
          gTCCAGAATGGCGG.CTTTTGTTTTGTTTTTTTTttt                         >  1:4650377/1‑36 (MQ=9)
           tCCAGCATGGCGG.TTTACTTTTGCCTTTTTTCtttt                        <  1:7407445/36‑1 (MQ=37)
           tCCAGAATGGCGG.TTTACGTTTGCCTTTTTTCtttt                        <  1:6027/36‑1 (MQ=255)
            ccAGAATGGCGG.TTTTGTTTTTTCTTTTTTTttttt                       >  1:270877/1‑36 (MQ=2)
            ccAGAATGGCGG.TTTGCGTTTGCCTTTTTTCttttt                       <  1:6783654/36‑1 (MQ=255)
            ccAGAATGGCGG.TTTACGTTTGCCTTTTTTCttttt                       <  1:1365717/36‑1 (MQ=255)
             cAGAATGGCGG.TTTACGTTTGCCTTTTTTCTTTTTg                      <  1:4013338/36‑1 (MQ=255)
              aGAGTGGCGG.TTTCCGTTTGCCTTTTTTCTTTTTGc                     <  1:5476319/36‑1 (MQ=37)
              aGAATGGCGG.TTTACGTTTGCCTTTTTTCTTTTTGc                     <  1:6380643/36‑1 (MQ=255)
               gAATGGCGGCTTTAATTTTGCCTTTTTTTttttttt                     >  1:3151108/1‑34 (MQ=9)
                 aTGGCGG.TTTACGTTTGCCTTTTTTCTTTTTGCGAt                  <  1:4543398/36‑1 (MQ=255)
                 aTGGCGG.TTTACGTTTGCCTTTTTTCTTTTTGCGAt                  <  1:5454737/36‑1 (MQ=255)
                 aTGGCGG.TTTACGTTTGCCTTTTTTCTTTTTGCGAt                  <  1:1581483/36‑1 (MQ=255)
                   ggCGG.TTTTCGTTTTCCTTTTTTTTTTTTGCGATCt                >  1:7749724/1‑36 (MQ=25)
                      gg.TTTACGTTTGCCTTTTTTCTTTTTGCGATCTTcc             <  1:6554384/36‑1 (MQ=255)
                      gg.TTTACGTTTGCCTTTTTTCTTTTTGCGATCTTcc             <  1:7681526/36‑1 (MQ=255)
                      gg.TTTACGTTTGCCTTTTTTCTTTTTGCGATCTTcc             <  1:4811064/36‑1 (MQ=255)
                       g.TTTACGTTTGCCTTTTTTCTTTTTTCGATCTTCCg            >  1:1056003/1‑36 (MQ=255)
                       g.TTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCg            >  1:29681/1‑36 (MQ=255)
                          ttACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGa          >  1:7308634/1‑36 (MQ=255)
                          ttACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGa          <  1:7063326/36‑1 (MQ=255)
                           tACGTTTGCCTTTTTTCTTTTTTCGATCTTCCGGaa         >  1:3501312/1‑36 (MQ=255)
                            aCGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAAt        >  1:4870107/1‑36 (MQ=255)
                             cGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATg       <  1:4510119/36‑1 (MQ=255)
                             cGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATg       <  1:2424765/36‑1 (MQ=255)
                              gTTTTCCTTTTTTCTTTTTGCGATCTTCCGGAATGc      <  1:2709304/36‑1 (MQ=255)
                              gTTTGCCTTTTTTCTTTTTGCTATCTTCCGGAATGc      >  1:3179518/1‑36 (MQ=255)
                              gTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGc      >  1:7598504/1‑36 (MQ=255)
                              gTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGc      >  1:7591180/1‑36 (MQ=255)
                              gTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGc      <  1:7286687/36‑1 (MQ=255)
                              gTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGc      <  1:2781570/36‑1 (MQ=255)
                              gTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGc      <  1:4818813/36‑1 (MQ=255)
                              gTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGc      >  1:5220906/1‑36 (MQ=255)
                              gTTTGCCTTTTTTCTTTTTGCGATCGTCCGGAATGc      >  1:1531780/1‑36 (MQ=255)
                               tttGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCg     <  1:5522492/36‑1 (MQ=255)
                                 tGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGt   <  1:3235307/36‑1 (MQ=255)
                                 tGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGt   <  1:4126581/36‑1 (MQ=255)
                                 tGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGt   <  1:7810867/36‑1 (MQ=255)
                                  gCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTg  <  1:4602011/36‑1 (MQ=255)
                                  |                                   
CGCAATACGCGTCCAGAATGGCGG.TTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTG  >  REL606/3069820‑3069888

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