Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,545,8110A→T 6.7% 0.2 ‑8.4 0.0661457nikAnickel transporter subunit
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (15/44):  total (22/44)
Fisher's exact test strand distribution p-value = 2.19e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GACATTCTGACCCGTCTGCATGACGAGGCGGTTTATCTGCCTATCAGTTACATCTCAATGATGGTGGTAT  >  REL606/3545777‑3545846
                                  |                                   
gACATTCTGACCCGTCTGCATGACGAGGCGGTTTAt                                    <  1:4930259/36‑1 (MQ=255)
 aCATTCTGACCCGTCTGCATGACGAGGCGGTTTATc                                   <  1:1355630/36‑1 (MQ=255)
 aCATTCTGACCCGTCTGCATGACGAGGCGGTTTATc                                   <  1:351427/36‑1 (MQ=255)
  cATTCTGACCCGTCTGCATGACGAGGCGGTTTATCt                                  <  1:410394/36‑1 (MQ=255)
  cATTCTGACCCGTCTGCATGACGAGGCGGTTTATCt                                  <  1:3205793/36‑1 (MQ=255)
   aTTCTGACCCGTCTGCATGACGAGGCGGTTTTTCtt                                 >  1:3804282/1‑35 (MQ=255)
   aTTCTGACCCGTCTGCATGACGAGGCGGTTTTTCtt                                 >  1:2646985/1‑35 (MQ=255)
   aTTCTGACCCGTCTGCATGACGAGGCGGTTTATCTg                                 <  1:2462470/36‑1 (MQ=255)
     tCTGACCCGTCTGCATGACGAGGCGGTTTATCTGcc                               <  1:784205/36‑1 (MQ=255)
     tCTGACCCGTCTGCATGACGAGGCGGTTTATCTGcc                               <  1:5001965/36‑1 (MQ=255)
     gcTGACCCGTCTGCATGACGAGGCGGTTTATCTGcc                               <  1:370130/35‑1 (MQ=255)
       tGACCCGTCTGCATTGCGAGGCGGTTTATCTGCCTa                             <  1:6491533/36‑1 (MQ=37)
       tGACCCGTCTGCATGACGAGGCGGTTTATCTTCCtt                             >  1:7189171/1‑35 (MQ=255)
       tGACCCGTCTGCATGACGAGGCGGTTTATCTGCCTa                             <  1:6422476/36‑1 (MQ=255)
       tGACCCGTCTGCATGACGAGGCGGTTTATCTGCCTa                             <  1:2891516/36‑1 (MQ=255)
        gACCCGTCTGCATGACGAGGCGGTTTATCTGCCTAt                            <  1:4401132/36‑1 (MQ=255)
         aCCCGTCTGCATGACGAGGCGGTTTATCTGCCTATc                           <  1:762181/36‑1 (MQ=255)
         aCCCGTCTGCATGACGAGGCGGTTTATCTGCCTATc                           <  1:5350391/36‑1 (MQ=255)
          cccGTCTGCATGACGAGGCGGTTTTTCTGCCTATCa                          >  1:7135315/1‑36 (MQ=255)
          cccGTCTGCATGACGAGGCGGTTTATCTGCCTATCa                          <  1:3932917/36‑1 (MQ=255)
           ccGTCTGCATGACGAGGCGGTTTATCTTCCTATCat                         >  1:1248849/1‑35 (MQ=255)
           ccGTCTGCATGACGAGGCGGTTTATCTGCCTATCAg                         <  1:3472065/36‑1 (MQ=255)
           ccGTCTGCATGACGAGGCGGTTTATCTGCCTATCAg                         <  1:3852164/36‑1 (MQ=255)
             gTCTGCATGACGAGGCGGTTTTTCTTCCTATCCGtt                       >  1:2209859/1‑36 (MQ=25)
             gTCTGCATGACGAGGCGGTTTATTTTCCTATCAGtt                       >  1:2661881/1‑36 (MQ=37)
             gTCTGCATGACGAGGCGGTTTATCTGCCTATCAGtt                       <  1:412421/36‑1 (MQ=255)
              tCTGCATGACGAGGCGGTTTATCTGCCTATCAGTTa                      <  1:7181602/36‑1 (MQ=255)
              tCTGCATGACGAGGCGGTTTATCTGCCTATCAGTTa                      <  1:6411922/36‑1 (MQ=255)
                tGCATGACGAGGCGGTTTTTCTTCCTTTTAGTTGCa                    >  1:5551070/1‑36 (MQ=16)
                tGCATGACGAGGCGGTTTATCTGCCTATCAGTTACa                    <  1:1009999/36‑1 (MQ=255)
                tGCATGACGAGGCGGTTTATCTGCCTATCAGTTACa                    <  1:4263807/36‑1 (MQ=255)
                 gCATGACGAGGCGGTTTATCTGCCTATCAGTTACAt                   <  1:4204091/36‑1 (MQ=255)
                  cATGACGAGGCGGTTTATCTTCCTATCAGTTACAtc                  >  1:4985042/1‑36 (MQ=255)
                  cATGACGAGGCGGTTTATCTGCCTATCAGTTACAtc                  <  1:6162348/36‑1 (MQ=255)
                  cATGACGAGGCGGTTTATCTGCCTATCAGTTACAtc                  <  1:445797/36‑1 (MQ=255)
                  cATGACGAGGCGGTTTATCTGCCTATCAGTTACAtc                  <  1:5743558/36‑1 (MQ=255)
                  cATGACGAGGCGGTTTATCTGCCTATCAGTTACAtc                  <  1:1358966/36‑1 (MQ=255)
                  cATGACGAGGCGGCTTTTCTTCCTATTAGTTACAtc                  >  1:2751623/1‑36 (MQ=19)
                   aTGACGAGGCGGTTTATCTTCCTATCAGTTACAtct                 >  1:1325102/1‑36 (MQ=255)
                    tGACGAGTCGGTTTATCTGCCTCTCAGTTACAtctc                <  1:5591060/36‑1 (MQ=38)
                    tGACGAGGCGGTTTATCTGCCTATCAGTTACAtctc                <  1:5433448/36‑1 (MQ=255)
                     gACGAGGCGGTTTATCTGCCTATCAGTTACATCTCa               <  1:2102378/36‑1 (MQ=255)
                     gACGAGGCGGGTTTTTTTCCTATCAGTTACAtctct               >  1:5523549/1‑35 (MQ=16)
                      aCGAGGCGGTTTATCTGCCTATCAGTTACATCTCaa              <  1:6342019/36‑1 (MQ=255)
                      aCGAGGCGGTTTAGCTGCCTATCAGTTACATCTCaa              >  1:3049454/1‑36 (MQ=255)
                       cGAGGCGGTTTATCTGCCTATCAGTTACATCTCAat             <  1:4299994/36‑1 (MQ=255)
                       cGAGGCGGTTTATCTGCCTATCAGTTACATCTCAat             <  1:3671934/36‑1 (MQ=255)
                       cGAGGCGGTTTATCTGCCTATCAGTTACATCTCAat             <  1:5413633/36‑1 (MQ=255)
                         aGGCGGTTTATCTTCCTATCAGTTAAATCTTAatga           >  1:550112/1‑36 (MQ=37)
                         aGGCGGTTTATCTGCCTATCAGTTACATCTCAatga           <  1:7411051/36‑1 (MQ=255)
                         aGGCGGTTTATCTGCCTATCAGTTACATCTCAatga           >  1:6766812/1‑36 (MQ=255)
                          ggCGGTTTATCTGCCTATCAGTTACATCTCAatgat          <  1:1041964/36‑1 (MQ=255)
                          ggCGGTTTATCTGCCTATCAGTGACATCTCAatgat          <  1:2568294/36‑1 (MQ=255)
                           gCGGTTTATCTGCCTATCAGTTACATCTCAatgatg         >  1:3989344/1‑36 (MQ=255)
                            cGGTTTATCTGCCTATCAGTTACATCTCAATGAtgg        <  1:2473122/36‑1 (MQ=255)
                            cGGTTTATCTGCCTATCAGTTACATCTCAATGAtgg        >  1:676098/1‑36 (MQ=255)
                            cGGTTTATCTGCCTATCAGTTACATCTCAATGAtgg        >  1:5862830/1‑36 (MQ=255)
                              gTTTATCTGCCTATCATTTACATCTCAATGAtggtg      >  1:7246806/1‑36 (MQ=255)
                              gTTTATCTGCCTATCAGTTACATCTCAATGAtggtg      >  1:1995831/1‑36 (MQ=255)
                              gTTTATCTGCCTATCAGTTACATCTCAATGAtggtg      <  1:5152823/36‑1 (MQ=255)
                               tttATCTGCCTATCAGTTACATCTCAATGAtggtgg     <  1:1086588/36‑1 (MQ=255)
                               tttATCTGCCTATCAGTTACATCTCAATGAtggtgg     >  1:7486436/1‑36 (MQ=255)
                                ttATCTGCCTATCAGTTACATCTCAATGATggtggt    <  1:5655393/36‑1 (MQ=255)
                                 tATCTGCCTATCAGTTACATCTCAATGATGGTGGTa   >  1:1679226/1‑36 (MQ=255)
                                 tATCTGCCTATCAGTTACATCTCAATGATGGTGGTa   <  1:4567317/36‑1 (MQ=255)
                                 tATCTGCCTATCAGTTACATCTCAATGATGGTGGTa   <  1:526880/36‑1 (MQ=255)
                                  aTCTGCCTATCAGTTACATCTCAATGATGGTGGTAt  <  1:5485192/36‑1 (MQ=255)
                                  |                                   
GACATTCTGACCCGTCTGCATGACGAGGCGGTTTATCTGCCTATCAGTTACATCTCAATGATGGTGGTAT  >  REL606/3545777‑3545846

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