Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,545,147 T→C 15.2% S265P (TCA→CCA)  nikA → nickel transporter subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,545,1470T→C15.2% 4.4 ‑3.9 ‑0.240793nikAnickel transporter subunit
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (19/14):  total (19/21)
Fisher's exact test strand distribution p-value = 2.02e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.74e-01

GCTTTAGCCAGAATCCGGCTTACCACACCCAACTGTCACAGCCGATCGAAACCGTGATGCTGGCGCTCAAT  >  REL606/3545112‑3545182
                                   |                                   
gNNTTAGCCAGAATCCGGNTTNNCACACCCAACTGt                                     <  1:7150701/36‑1 (MQ=37)
gCTTTAGCCAGAATCCGGCTTACCACACCCAACTGt                                     >  1:5193166/1‑36 (MQ=255)
gCTTTAGCCAGAATCCGGCTTACCACACCCAACTGt                                     <  1:3271111/36‑1 (MQ=255)
gCTTTAGCCAGAATCCGGCTTACCACACCCAACTGt                                     >  1:304973/1‑36 (MQ=255)
 cTTTAGCCAGAATCCGGCTTACCACACCCAACTGTc                                    >  1:2411268/1‑36 (MQ=255)
   ttAGCCAGAATCCGGCTTACCACACCCAACTGTcac                                  <  1:3348900/36‑1 (MQ=255)
   ttAGCCAGAATCCGGCTTACCACACCCAACTGTcac                                  >  1:1370651/1‑36 (MQ=255)
     aGCCAGAATCCGGCTTACCACACCCAACTGTCACAg                                <  1:2628708/36‑1 (MQ=255)
     aGCCAGAATCCGGCTTACCACACCCAACTGTCACAg                                >  1:630028/1‑36 (MQ=255)
      gCCAGAATCCGGCTTACCACACCCAACTGTCACAGc                               <  1:2870569/36‑1 (MQ=255)
      gCCAGAATCCGGCTTACCACACCCAAATGTAACAGc                               >  1:1902369/1‑36 (MQ=37)
         aaaaTCCGGCTTACCACACCCAACTGTCACAGCCGa                            <  1:6415545/34‑1 (MQ=255)
           aacccGGCTTACCACACCCAACGGTCACAGCCGATc                          <  1:443853/33‑1 (MQ=37)
           aaTCCGGCTTACCACACCCAACTGTCACAGCCGATc                          <  1:7021999/36‑1 (MQ=255)
           aaTCCGGCTTACCACACCCAACTGTCACAGCCGATc                          >  1:3350737/1‑36 (MQ=255)
           aaTCCGGCTTACCACACCCAACTGTCACAGCCGATc                          >  1:4115826/1‑36 (MQ=255)
            aTCCGGCTTACCACACCCCACTGTCACAGCCGATcc                         >  1:3271564/1‑35 (MQ=255)
            aTCCGGCTTACCACACCCAACTGTCACACCCTATCg                         >  1:817251/1‑36 (MQ=37)
               cGGCTTACCACACCCAACTGTCACAGCCGATCGaaa                      >  1:7506947/1‑36 (MQ=255)
               cGGCTAACCACACCCAACTGTCACAGCCGATCGaaa                      <  1:1527478/36‑1 (MQ=255)
                ggCTTACCACACCCAACTGTCACAGCCGATCGAAAc                     >  1:2870231/1‑36 (MQ=255)
                  ctgaCCACACCAAACGGTCACAGCCGATCGAAACCg                   <  1:4593894/33‑1 (MQ=25)
                   ttACCACACCCAAGTGTCACAGCCGATCGAAACCGt                  >  1:4327240/1‑36 (MQ=255)
                   ttACCACACAAAACGCGCCCAGCCGATCGAAACCGt                  <  1:6533983/36‑1 (MQ=9)
                   taaCCCCACCAAACGGTCCCAGCCGATCGAAACCGt                  <  1:3495989/34‑1 (MQ=12)
                    ccccACACCAAACTGCCACAGCCGATCGAAACCGTg                 <  1:6812603/34‑1 (MQ=25)
                     acaacacCAAACTGCCCCAGCCGATCGAAACCGTGa                <  1:4567355/33‑1 (MQ=19)
                     aCCACACCCAACTGTCACAGCCGATCGAAACCGTGa                >  1:7384680/1‑36 (MQ=255)
                     aCCACACCCAACTGTCACAGCCGATCGAAACCGCGa                >  1:5183872/1‑36 (MQ=255)
                      ccACCCCAAACTGTCACAGCCGATCGAAACCGTGAt               <  1:5245466/36‑1 (MQ=37)
                       cacaCCAAACTGTCACAGCCGATCGAAACCGTGATg              <  1:4580215/36‑1 (MQ=255)
                       aacacCAAACTGTCACAGCCGATCGAAACCGTGATg              <  1:6674165/35‑1 (MQ=255)
                       aacacCAAACTGTAACAGCCGATCGAAACCGTGATg              <  1:199480/35‑1 (MQ=25)
                        acacCAAACTGCCCCAGCCGATCGAAACCGTGATGc             <  1:7020337/36‑1 (MQ=25)
                         cacCCAACTGTCACAGCCGAGCGAAACCGCGATGCt            >  1:7360836/1‑36 (MQ=37)
                         cacCAAACTGCCCCAGCCGATCGAAACCGTGATGCt            <  1:2488884/36‑1 (MQ=25)
                         cacCAAACGGTCACAGCCGATCGAAACCGTGATGCt            <  1:6690016/36‑1 (MQ=37)
                          acCCAACTGTCACAGCCGATCGAAACCGTGATGCTg           >  1:6368903/1‑36 (MQ=255)
                          acCCAACTGTCACAGCCGATCGAAACCGTGATGCTg           >  1:4600926/1‑36 (MQ=255)
                           cccAACTGTCACAGCCGATCGAAACCGTGATGCTgg          >  1:6981352/1‑36 (MQ=255)
                              aaCTGCCACAGCCGATCGAAACCGTGATGCTGgcgc       <  1:6150614/36‑1 (MQ=255)
                               ccTGCCACAGCCGATCGAAACCGTGATGCTGGCGCt      <  1:5090894/35‑1 (MQ=255)
                                 tGTCACAGCCGATCGAAACCGTGATGCTGGCGCTCa    >  1:3595423/1‑36 (MQ=255)
                                  gTCACAGCCGATCGAAACCGTGATGCTGGCGCTCaa   >  1:7327021/1‑36 (MQ=255)
                                   tCACAGCCGATCGAAACCGTGATGCTGGCGCTCAAt  >  1:6511319/1‑36 (MQ=255)
                                   |                                   
GCTTTAGCCAGAATCCGGCTTACCACACCCAACTGTCACAGCCGATCGAAACCGTGATGCTGGCGCTCAAT  >  REL606/3545112‑3545182

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 13 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG