Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,226,6880C→A13.4% 4.9 ‑7.2 0.073973malKfused maltose transport subunit, ATP‑binding component of ABC superfamily/regulatory protein
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (0/12):  ref base (32/29):  total (32/41)
Fisher's exact test strand distribution p-value = 7.30e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AAACTCAAATCCATATCCAGATCCCTTCCATTCGTCAAAACCTGGTGTACCGCCAGAACGACGTGGTGTTG  >  REL606/4226653‑4226723
                                   |                                   
aaaCTCAAATCCATATCCAGATCCCTTCCATTCGTc                                     >  1:3957939/1‑36 (MQ=255)
aaaCTCAAATCCATATCCAGATCCCTTCCATTCGTc                                     <  1:314726/36‑1 (MQ=255)
aaaCTCAAATCCATATCCAGATCCCTTCCATTCGTc                                     >  1:2678005/1‑36 (MQ=255)
aaaCTCAAATCCATATCCAGATCCCTTCCATTCGTc                                     <  1:6625877/36‑1 (MQ=255)
aaaCTCAAATCCATATCCAGATCCCTTCCATTCGTc                                     >  1:887183/1‑36 (MQ=255)
aaaCTCAAATCCATATCCAGATCCCTTCCATTCGTc                                     <  1:7821857/36‑1 (MQ=255)
 aaCTCAAATCCATAGCCAGATCCCTTCCATTCGTCa                                    <  1:2111883/36‑1 (MQ=255)
  acccAAAGCCATATCCAGATCCCTTCCATTCGTCaa                                   <  1:2581083/33‑1 (MQ=37)
  aCTCAAATCCATATCCAGATCCCTTCCATTCGTCaa                                   <  1:117350/36‑1 (MQ=255)
  aCTCAAATCCATATCCAGATCCCTTCCATTCGTCaa                                   >  1:5204751/1‑36 (MQ=255)
  aCTCAAATCCATATCCAGATCCCTTCCATTCGTCaa                                   <  1:5992252/36‑1 (MQ=255)
  aCTCAAATCCATATCCAGATCCCTTCCATTCGTCaa                                   >  1:1425627/1‑36 (MQ=255)
   cTCAAATCCATATCCAGATCCCTTCCATTCGTCaaa                                  <  1:7692868/36‑1 (MQ=255)
   cTCAAATCCATATCCAGATCCCTTCCATTCGTCaaa                                  <  1:3443886/36‑1 (MQ=255)
    tCAAATCCATATCCAGATCCCTTCCATTCGTCaaaa                                 <  1:6494459/36‑1 (MQ=255)
    tCAAATCCATATCCAGATCCCTTCCATTCGTCaaaa                                 <  1:6991073/36‑1 (MQ=255)
    tCAAATCCATATCCAGATCCCTTCCATTCGTCaaaa                                 <  1:5980728/36‑1 (MQ=255)
    tCAAATCCATATCCAGATCCCTTCCATTCGTCaaaa                                 >  1:4502672/1‑36 (MQ=255)
     cAAATCCATATCCATATCCCTTCCATTCGTCaactc                                <  1:3946369/36‑4 (MQ=37)
     cAAATCCATATCCAGATCCCTTCCATTCGTCAAAAc                                >  1:6432027/1‑36 (MQ=255)
     cAAATCCATATCCAGATCCCTTCCATTCGTCAAAAc                                <  1:6888138/36‑1 (MQ=255)
     cAAATCCATATCCAGATCCCTTCCATTCGTCAAAAc                                <  1:302758/36‑1 (MQ=255)
      aaaTCCATATCCAGATCCCTTCCATTCGTCAAAAcc                               <  1:6707431/36‑1 (MQ=255)
      aaaTCCATATCCAGATCCCTTCCATTCGTCAAAAcc                               <  1:3100160/36‑1 (MQ=255)
        aTCCATATCCAGATCCCTTCCATTCGTCAAAACCTg                             <  1:3919992/36‑1 (MQ=255)
        aTCCATATCCAGATCCCTTCCATTCGTCAAAACCTg                             >  1:4296045/1‑36 (MQ=255)
         tCCATATCCAGATCCCTTCCATTCGTCAAAACCTgg                            <  1:4833356/36‑1 (MQ=255)
         tCCATATCCAGATCCCTTCCATTCGTCAAAACCTgg                            <  1:1247003/36‑1 (MQ=255)
          ccATATCCAGATCCCTTCCATTCGTCAAAACCTGgt                           >  1:2831715/1‑36 (MQ=255)
          ccATATCCAGATCCCTTCCATTCGTCAAAACCTGgt                           >  1:2130671/1‑36 (MQ=255)
           cATATCCAGATCCCTTCCATTCGTCAAAACCTGgtg                          <  1:721294/36‑1 (MQ=255)
           cATATCCAGATCCCTTCCATTCGTCAAAACCTGgtg                          <  1:500182/36‑1 (MQ=255)
           cATATCCAGATCCCTTCCATTCGTCAAAACCTGgtg                          >  1:6268045/1‑36 (MQ=255)
           cATATCCAGATCCCTTCCATTCGTCAAAACCTGgtg                          >  1:1093876/1‑36 (MQ=255)
           cATATCCAGATCCCTTCCATTCGTCAAAACCTGgtg                          >  1:1563298/1‑36 (MQ=255)
            atatCCAGATCCCTTCCATTCGTCAAAACCTGgtgt                         >  1:77505/1‑36 (MQ=255)
             tatCCAGATCCCTTCCATTCGTCAAAACCTGGTGTa                        <  1:6817127/36‑1 (MQ=255)
             tatCCAGATCCCTTCCATTCGTCAAAACCTGGTGTa                        <  1:704566/36‑1 (MQ=255)
             tatCCAGATCCCTTCCATTCGTCAAAACCTGGTGTa                        <  1:117693/36‑1 (MQ=255)
              atCCAGATCCCTTCCATTCGTCAAAACCTGGTGTAc                       >  1:5206843/1‑36 (MQ=255)
              atCCAGATCCCTTCCATTCGTCAAAACCTGGTGTAc                       >  1:7332883/1‑36 (MQ=255)
               tCCAGATCCCTTCCATTCGTCAAAACCTGGTGTAcc                      >  1:1260235/1‑36 (MQ=255)
                 cAGATCCCTTCCATTCGCCAAAACCTGGTGTAcccc                    >  1:1580819/1‑34 (MQ=38)
                  aGATCCCTTCCATTCGTCAAAACCTGGTGTACCGcc                   >  1:7721541/1‑36 (MQ=255)
                  aGATCCCTTCCATTCGTCAAAACCTGGTGTACCGcc                   >  1:1825892/1‑36 (MQ=255)
                  aGATCCCTTCCATTCGTCAAAACCTGGTGTACCGcc                   >  1:3077202/1‑36 (MQ=255)
                  aGATCCCTTCCATTCGTCAAAACCTGGTGTACCGcc                   >  1:7366861/1‑36 (MQ=255)
                  aGACCCCTTCCATTCGTCAAAACCTGGTGTACCGcc                   <  1:5482561/36‑1 (MQ=255)
                   gATCCCTTCCATTCGTCAAAACCTGGTGTACCGCCa                  >  1:2400707/1‑36 (MQ=255)
                    aTCCCTTCCATTCGTCAAAACCTGGTGTACCGCCAg                 <  1:4852615/36‑1 (MQ=255)
                    aTCCCTTCCATTCGTCAAAACCTGGTGTACCGCCAg                 >  1:5278303/1‑36 (MQ=255)
                     tCCCTTCCATTCGTCAAAACCTGGTGTACCGCCAGa                >  1:2805374/1‑36 (MQ=255)
                     ccccTTCCATCCGTAAAACCCTGGTGTACCGCCAGa                <  1:5661077/35‑1 (MQ=25)
                       ccTTCCATTCGTCAAAACCTGGTGTACCGCCAGAac              >  1:7464042/1‑36 (MQ=255)
                       ccTTCCATTCGTAAAAACCTGGTGTACCGCCAGAac              <  1:6227083/36‑1 (MQ=255)
                       ccTTCCATTCGAAAAACCCGGGTGTACCGCCAGAac              <  1:2866692/36‑1 (MQ=16)
                         ttCCATTCGTCAAAACCTGGTGTACCGCCAGAacga            <  1:7759875/36‑1 (MQ=255)
                         ttCCATTCGTCAAAACCTGGTGTACCGCCAGAacga            >  1:4583107/1‑36 (MQ=255)
                         ttCCATTCAAAAAACCGTGTTTTACCGCCAGAacga            <  1:2904540/36‑1 (MQ=2)
                          tCCATTCGTCAAAACCTGGTGTACCGCCAGAacgac           >  1:4562159/1‑36 (MQ=255)
                          tCCATTCGTAAAAACCTGGTGTACCTCCAGAacgac           <  1:7440099/36‑1 (MQ=37)
                          tCCATTCGCAAAAACCGGGTGTACCGCCAGAacgac           <  1:3975679/36‑1 (MQ=25)
                          ccaaTTCGAAAAACCCTGGTGAACCGCCAGAacgac           <  1:6311738/33‑1 (MQ=16)
                           ccATTCGTCAAAACCTGGTGTACCGCCAGAacgacg          >  1:211717/1‑36 (MQ=255)
                           ccATTCGTCAAAACCTGGTGTACCGCCAGAacgacg          >  1:7479254/1‑36 (MQ=255)
                           ccATTCGAAAAAACCTGGTGTACCGCCAGAacgacg          <  1:1655634/36‑1 (MQ=37)
                            cATTCGTCAAAACCTGGTGTACCGCCAGAACGACgt         >  1:727684/1‑36 (MQ=255)
                            cATTCGTAAAACCCTGGTGTACCGCCAGAACGACgt         <  1:5947526/36‑1 (MQ=37)
                            aaTTCGTCAAAACCTGGTGTACCGCCAGAACGACgt         <  1:7099833/35‑1 (MQ=255)
                             aTTCGTCAAAACCTGGTGTAAAGCCAGAACGACgtg        >  1:2691246/1‑36 (MQ=37)
                             aTTCGAAAAACCCGGGTGTACCGCCAGAACGACgtg        <  1:3764828/36‑1 (MQ=16)
                              ttCGTCAAAACCTGGTGTACCGCCAGAACGACgtgg       <  1:5416020/36‑1 (MQ=255)
                               tCGTCAAAACCTGGTGTACCGCCAGAACGACgtggt      <  1:7100636/36‑1 (MQ=255)
                               tCGTCAAAACCTGGTGTACCGCCAGAACGACgtggt      >  1:4140831/1‑36 (MQ=255)
                               tCGTAAAACCCGGGTGTACCGCCAGAACGACgtggt      <  1:6696410/36‑1 (MQ=25)
                               tCGTAAAAACCTGTTGTACCGCCAGAACGACgtggt      <  1:2225681/36‑1 (MQ=37)
                                cGTCAAAACCTGGTGTACCGCCAGAACGACgtggtg     <  1:7319399/36‑1 (MQ=255)
                                 gTCAAAACCTGGTGTACCGCCAGAACGACGTGgtgt    >  1:1380978/1‑36 (MQ=255)
                                  tCAAAACCTGGTGTACCGCCAGAACGACGTGGTGtt   >  1:4020736/1‑36 (MQ=255)
                                   cAAATCCTGGTGTGCCGCCAGAACGACGTGGTGTTg  <  1:7459304/36‑1 (MQ=37)
                                   cAAAACCTGGTGTACCGCCAGAACGACGTGGTGTTg  >  1:2909943/1‑36 (MQ=255)
                                   cAAAACCTGGTGTACCGCCAGAACGACGTGGTGTTg  >  1:6520378/1‑36 (MQ=255)
                                   |                                   
AAACTCAAATCCATATCCAGATCCCTTCCATTCGTCAAAACCTGGTGTACCGCCAGAACGACGTGGTGTTG  >  REL606/4226653‑4226723

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG