Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,928,5130T→A 6.0% 0.1 ‑5.6 0.071297rhotranscription termination factor Rho
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (41/23):  total (42/29)
Fisher's exact test strand distribution p-value = 3.62e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATT  >  REL606/3928479‑3928547
                                  |                                  
gTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa                                   >  1:2156757/1‑36 (MQ=255)
gTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa                                   <  1:2029078/36‑1 (MQ=255)
gTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa                                   >  1:897690/1‑36 (MQ=255)
gTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa                                   >  1:1062159/1‑36 (MQ=255)
 tttCAACCTCCGCACTGGTGATACCATCTCTGGTaa                                  <  1:5577118/36‑1 (MQ=255)
  ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAg                                 >  1:3652901/1‑36 (MQ=255)
  ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAg                                 <  1:7087186/36‑1 (MQ=255)
  ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAg                                 >  1:3985326/1‑36 (MQ=255)
  ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAg                                 <  1:6741143/36‑1 (MQ=255)
   tCAACCTCCGCACTGGTGATACCATCTCTGGTAAGa                                >  1:4617118/1‑36 (MQ=255)
   tCAACCTCCGCACTGGTGATACCATCTCTGGTAAGa                                >  1:7369062/1‑36 (MQ=255)
   tCAACCTCCGCACTGGTGATACCATCTCTGGTAAGa                                <  1:1884993/36‑1 (MQ=255)
    cAACCTCCGCACTGGTGATACCATCTCTGGTAAGAt                               <  1:2615332/36‑1 (MQ=255)
    cAACCTCCGCACGGGTGATACCATCTCTGGTAAGAt                               >  1:4259113/1‑36 (MQ=255)
      aCCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                             <  1:607213/36‑1 (MQ=255)
      aCCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                             >  1:3249638/1‑36 (MQ=255)
      aCCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                             <  1:2175158/36‑1 (MQ=255)
      aCCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                             <  1:1759647/36‑1 (MQ=255)
      aCCTCCGCACTGGTGATACCATCTCTGGTAAGATTc                             <  1:1680227/36‑1 (MQ=255)
       ccTCCGCACTGGTGATACCATCTCTGGTAAGATTcg                            >  1:3932192/1‑36 (MQ=255)
       ccTCCGCACTGGTGATACCATCTCTGGTAAGATTcg                            >  1:4514410/1‑36 (MQ=255)
        cTCCGCACTGGTGATACCATCTCTGGTAAGATTcgc                           <  1:606445/36‑1 (MQ=255)
         tacGCGCTGGTTATACCATCTATGGTAAGATTcgcc                          <  1:697339/34‑1 (MQ=25)
         tCCGCACTGGTGATACCATCTCTGGTAAGATTcgcc                          <  1:1465118/36‑1 (MQ=255)
          ccGCACTGGTGATACCATCTCTGGTAAGATTcgccc                         <  1:5562204/36‑1 (MQ=255)
          ccGCACTGGTGATACCATCTCTGGTAAGATTcgccc                         >  1:1891718/1‑36 (MQ=255)
          ccGCACTGGTGATACCATCTCTGGTAAGATTcgccc                         >  1:6596757/1‑36 (MQ=255)
          ccGCACTGGTGATACCATCTCTGGTAAGATTcgccc                         <  1:3187871/36‑1 (MQ=255)
           cGCACTGGTGATACCATCTCTGGTAAGATTcgcccg                        <  1:5181946/36‑1 (MQ=255)
            gCACTGGTGATACCATCTCTGGTAAGATTcgcccgc                       <  1:2709485/36‑1 (MQ=255)
            gCACTGGTGATACCATCTCTGGTAAGATTcgcccgc                       <  1:7439913/36‑1 (MQ=255)
             cACTGGTGATACCATCTCTGGTAAGATTcgcccgcc                      >  1:5758840/1‑36 (MQ=255)
             cACGGGTGATCCCATCTCTGGTAAGATTcgcccgcc                      <  1:3540279/36‑1 (MQ=37)
              aCTGGTGATACCATCTCTGGTAAGATTCGCcctccg                     >  1:3709159/1‑36 (MQ=255)
              aCTGGTGATACCATCTCTGGTAAGATTCGCccgccg                     >  1:6286886/1‑36 (MQ=255)
              aCTGGTGATACCATCTCTGGTAAGATTCGCccgccg                     >  1:2684922/1‑36 (MQ=255)
              aCTGGTGATACCATCTCTGGTAAGATTCGCccgccg                     >  1:4464289/1‑36 (MQ=255)
              aCTGGTGATACCATCTCTGGTAAGATTCGCccgccg                     >  1:4777102/1‑36 (MQ=255)
              aCTGGTGATACCATCTCTGGGAAGATTCGCccgccg                     >  1:897038/1‑36 (MQ=255)
               cTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGa                    >  1:269466/1‑36 (MQ=255)
                tGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaa                   >  1:1651218/1‑36 (MQ=255)
                tGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaa                   >  1:7698397/1‑36 (MQ=255)
                 ggTGATACCATCTCTGGTAAGATTCGCCCGCCGaaa                  >  1:4215777/1‑36 (MQ=255)
                  gTGATCCCATCTCTGGAAAGATTCGCCCGCCGAaag                 <  1:416974/36‑1 (MQ=25)
                  gTGATACCATCTCTGGTAAGATTCGCCCGCCGAaag                 >  1:7287863/1‑36 (MQ=255)
                  gTGATACCATCTCGGGAAAGATTCGCCCGCCGAaag                 <  1:5255981/36‑1 (MQ=37)
                   tGATACCATCTCTGGTAAGATTCGCCCGCCGAaaga                >  1:1967764/1‑36 (MQ=255)
                   tGATACCATCTCTGGTAAGATTCGCCCGCCGAaaga                <  1:2806798/36‑1 (MQ=255)
                    aaTACCATCTTTGGAATAGTGCGCCCGCCGAaagaa               <  1:1874324/35‑1 (MQ=9)
                     aTACCATCTCTGGTAAGATTCGCCCGCCGAaagaag              >  1:7414583/1‑36 (MQ=255)
                     aTACCATCTCTGGTAAGATTCGCCCGCCGAaagaag              >  1:6937112/1‑36 (MQ=255)
                     aTACCATCTCTGGTAAGATTCGCCCGCCGAaagaag              >  1:2691188/1‑36 (MQ=255)
                     aTACCATCTCTGGTAAGATTCGCCCGCCGAaagaag              >  1:5107/1‑36 (MQ=255)
                      tACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAgg             >  1:3586587/1‑36 (MQ=255)
                       aCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGt            <  1:3787343/36‑1 (MQ=255)
                        ccATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTg           <  1:6673424/36‑1 (MQ=255)
                        ccATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTg           >  1:879896/1‑36 (MQ=255)
                        caaTCTCTGAAAAGATCCGCCCGCCGAAAGAAGGTg           <  1:3052541/34‑1 (MQ=21)
                         cATCTCGGGAAAGATTCGCCCGCCGAAAGAAGGTGa          <  1:652174/36‑1 (MQ=37)
                          aTCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGaa         >  1:6115044/1‑36 (MQ=255)
                           tctcTGGTAAGATTCGCCCGCCGAAAGAAGGTGAAc        >  1:187726/1‑36 (MQ=255)
                           tctcGGGTAAGATTCGCCCGCCGAAAGAAGGTGAAc        <  1:1059548/36‑1 (MQ=255)
                             tctGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGc      >  1:1010143/1‑36 (MQ=255)
                             tctGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGc      >  1:5850560/1‑36 (MQ=255)
                             tctGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGc      >  1:3581132/1‑36 (MQ=255)
                               tGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTa    >  1:1628890/1‑36 (MQ=255)
                               tGGAAAGATTCGCCCGCCGAAAGAAGGTGAACGCTa    <  1:1188685/36‑1 (MQ=255)
                                ggTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAt   >  1:5231171/1‑36 (MQ=255)
                                ggTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAt   >  1:6556235/1‑36 (MQ=255)
                                ggTAAGATTCGCCCGCCGAAAGAAGGGGAACGCTAt   >  1:3819160/1‑36 (MQ=255)
                                 gTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAtt  >  1:6056726/1‑36 (MQ=255)
                                  |                                  
GTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATT  >  REL606/3928479‑3928547

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