Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,248,9570A→T94.1% 5.3 ‑0.5 ‑0.1682292infBtranslation initiation factor IF‑2
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (35/29):  ref base (3/1):  total (38/30)
Fisher's exact test strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.67e-01

AATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTG  >  REL606/3248924‑3248988
                                 |                               
aaTCAGGTTATAGATGACGGAGTAGTAACGCAGTTc                               >  1:5038181/1‑36 (MQ=255)
aaTCAGGTTATAGATGACGGAGTAGTAACGCAGTTc                               <  1:3977609/36‑1 (MQ=255)
aaTCAGGTTATAGATGACGGAGTAGTAACGCAGTTc                               >  1:1830446/1‑36 (MQ=255)
 aTCAGGTTATAGATGACGGAGTAGTAACGCAGTTcc                              <  1:7905356/36‑1 (MQ=255)
 aTCAGGTTATAGATGACGGAGTAGTAACGCAGTTcc                              <  1:3814745/36‑1 (MQ=255)
  tCAGGTTATAGATGACGGAGTAGTAACGCAGTTCCa                             >  1:6000604/1‑36 (MQ=255)
  tCAGGTTATAGATGACGGAGTAGTAACGCAGTTCCa                             <  1:4744245/36‑1 (MQ=255)
  tCAGGTTATAGATGACGGAGTAGTAACGCAGTTCCa                             <  1:2796548/36‑1 (MQ=255)
   cAGGTTATAGATGACGGAGTAGTAACGCAGTTCCAg                            <  1:948895/36‑1 (MQ=255)
    aGGTTATAGATGACGGAGTAGTAACGCAGTTCCAgg                           >  1:7262606/1‑36 (MQ=255)
    aGGTTATAGATGACGGAGTAGTAACGCAGTTCCAgg                           >  1:1687959/1‑36 (MQ=255)
     ggTTATAGATGACGGAGTAGTAACGCAGTTCCAGGc                          <  1:1386575/36‑1 (MQ=255)
     ggTTATAGATGACCGAGTAGTAACGCAGTTCCAGGc                          >  1:5301696/1‑36 (MQ=37)
      gTTATAGATGACGGAGTAGTAACGCAGTTCCAGGCt                         <  1:3310878/36‑1 (MQ=255)
      gTTATAGATGACGGAGTAGTAACGCAGTTCCAGGCt                         <  1:1119912/36‑1 (MQ=255)
       ttATAGATGACGGTGTAGTAACGCAGTTCCAGGCtt                        >  1:1947231/1‑36 (MQ=37)
       ttATAGATGACGGAGTAGTAACGCAGTTCCAGGCtt                        <  1:5868738/36‑1 (MQ=255)
       ttATAGATGACGGAGTAGTAACGCAGTTCCAGGCtt                        >  1:2887406/1‑36 (MQ=255)
        tataGATGACGGAGTAGTAACGCAGTTCCAGGCttt                       >  1:4164725/1‑36 (MQ=255)
         ataGATGACGGAGTAGTAACGCAGTTCCAGGCTTTc                      >  1:1194672/1‑36 (MQ=255)
         ataGATGACGGAGTAGTAACGCAGTTCCAGGCTTTc                      >  1:2471370/1‑36 (MQ=255)
          taGATGACGGAGTAGTAACGCAGTTCCAGGCTTTcc                     >  1:2208465/1‑36 (MQ=255)
          taGATGACGGAGTAGTAACGCAGTTCCAGGCTTTcc                     <  1:4655723/36‑1 (MQ=255)
          taGATGACGGAGTAGTAACGCAGTTCCAGGCTTTcc                     >  1:1709943/1‑36 (MQ=255)
          taGATGACGGAGTAGTAACGCAGTTACAGGCtttac                     >  1:3203192/1‑34 (MQ=37)
           aGATGACGGAGTAGTAACGCAGTTCCAGGCTTTCCg                    >  1:5526357/1‑36 (MQ=255)
           aGATGACGGAGTAGTAACGCAGTTCCAGGCTTTCCg                    <  1:1420748/36‑1 (MQ=255)
           aGATGACGGAGTAGTAACGCAGTTCCAGGCTTTCCg                    >  1:2676534/1‑36 (MQ=255)
           aGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCg                    <  1:6018193/36‑1 (MQ=255)
           aGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCg                    >  1:1519609/1‑36 (MQ=255)
            gATGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGc                   >  1:6561418/1‑36 (MQ=255)
             aTGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCt                  <  1:1743063/36‑1 (MQ=255)
              tGACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCtt                 <  1:2142475/36‑1 (MQ=255)
               gACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTc                <  1:833531/36‑1 (MQ=255)
               gACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTc                <  1:4563393/36‑1 (MQ=255)
               gACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTc                <  1:5634435/36‑1 (MQ=255)
               gACGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTc                >  1:7245502/1‑36 (MQ=255)
                aCGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCa               <  1:3687842/36‑1 (MQ=255)
                 cGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCaa              <  1:6576746/36‑1 (MQ=255)
                 cGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCaa              >  1:3936414/1‑36 (MQ=255)
                 cGGAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCaa              >  1:5603684/1‑36 (MQ=255)
                  ggAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCAAt             <  1:3409598/36‑1 (MQ=255)
                  ggAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCAAt             <  1:7799286/36‑1 (MQ=255)
                  ggAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAAt             >  1:7596404/1‑36 (MQ=255)
                  ggAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAAt             >  1:2860048/1‑36 (MQ=255)
                   gAGTAGTAACGCAGTTCCAGGCTTTCCGCTTCAATc            <  1:393839/36‑1 (MQ=255)
                    agtagtAACGCAGTTCCAGGCTTTCCGCTTCAATCa           >  1:6281333/1‑36 (MQ=255)
                    agtagtAACGCAGTTCCAGGCTTTCCGCTTCAATCa           >  1:2000420/1‑36 (MQ=255)
                    agtagtAACGCAGTTCCAGGCTTTCCGCTTCAATCa           >  1:6375920/1‑36 (MQ=255)
                     gtagtaACGCAGTTCCAGGCTTTCCGCTTCAATCAc          <  1:5641799/36‑1 (MQ=255)
                     gtagtaACGCAGTTCCAGGCTTTCCGCTTCAATCAc          >  1:6310182/1‑36 (MQ=255)
                     gtagtaACGCAGTTCCAGGCTTTCCGCTTCAATCAc          <  1:6212873/36‑1 (MQ=255)
                      tagtaACGCAGTTCCAGGCTTTCCGCTTCAATCACt         >  1:7390900/1‑36 (MQ=255)
                      tagtaACGCAGTTCCAGGCTTTCCGCTTCAATCACt         <  1:963781/36‑1 (MQ=255)
                      tagtaACGCAGTTCCAGGCTTTCCGCTTCAATCACt         >  1:1336099/1‑36 (MQ=255)
                        gtaACGCAGTTCCAGGCTTTCCGCTTCAATCACttt       >  1:2871128/1‑36 (MQ=255)
                        gtaACGCAGTTCCAGGCTTTCCGCTTCAATCACttt       >  1:157202/1‑36 (MQ=255)
                        gtaACGCAGTTCCAGGCTTTCCGCTTCAATCACttt       <  1:3264520/36‑1 (MQ=255)
                         taccGCAGTTCCAGGCTTTCCGCTTCAATCCCTTTa      <  1:1476201/33‑1 (MQ=25)
                         taACGCAGTTCCAGGCTTTCCGCTTCAATCACTTTa      <  1:1323356/36‑1 (MQ=255)
                         taACGCAGTTCCAGGCTTTCCGCTTCAATCACTTTa      >  1:1709501/1‑36 (MQ=255)
                          aaCGCAGTTCCAGGCTTTCCGCTTCAATCACTTTAc     >  1:6363275/1‑36 (MQ=255)
                          aaCGCAGTTCCAGGCTTTCCGCTTCAATCACTTTAc     >  1:7174050/1‑36 (MQ=255)
                          aaCGCAGTTCCAGGCTTTCCGCTTCAACCACTTTAc     >  1:5197225/1‑36 (MQ=37)
                            cGCAGTTCCAGGCTTTCCGCTTCAATCACTTTACGt   <  1:5865343/36‑1 (MQ=255)
                            cGCAGTTCCAGGCTTTCCGCTTCAATCACTTTACGt   >  1:6602764/1‑36 (MQ=255)
                             gCAGTTCCAGGCTTTCCGCTTCAATCACTTTACGTg  >  1:2784484/1‑36 (MQ=255)
                             gCAGTTCCAGGCTTTCCGCTTCAATCACTTTACGTg  >  1:3740267/1‑36 (MQ=255)
                                 |                               
AATCAGGTTATAGATGACGGAGTAGTAACGCAGATCCAGGCTTTCCGCTTCAATCACTTTACGTG  >  REL606/3248924‑3248988

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 37 ≤  ATCG