Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,171,6820G→T 8.7% 1.7 ‑9.5 ‑0.081286ygjTpredicted inner membrane protein, part of terminus
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (19/53):  total (28/53)
Fisher's exact test strand distribution p-value = 7.66e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.89e-01

GATTGAGAAATCGCTGGCGGTCGATAACGTCTTTGTCTGGCTGATGTTGTTCAGCTATTTCTCCGTTCCG  >  REL606/3171648‑3171717
                                  |                                   
gATTGAGAAATCGCTGGCGGTCGATAACGTCTTTGt                                    <  1:2862130/36‑1 (MQ=255)
 aTTGAGAAATCGCTGGCGGTCGATAACGTCTTTTTc                                   >  1:5886275/1‑36 (MQ=255)
 aTTGAGAAATCGCTGGCGGTCGATAACGTCTTTTTc                                   >  1:7098802/1‑36 (MQ=255)
 aTTGAGAAATCGCTGGCGGTCGATAACGTCTTTGTc                                   <  1:6638348/36‑1 (MQ=255)
 aTTGAGAAATCGCTGGCGGTCGATAACGTCTTTGTc                                   >  1:7335253/1‑36 (MQ=255)
  ttaagaAATCGTTGGCGTTCGATAACGTCTTTGTCt                                  <  1:2693031/33‑1 (MQ=25)
  ttGAGAAATCGCTGGCGGTCGATAACGTCTTTTTCt                                  >  1:6443764/1‑36 (MQ=255)
  ttGAGAAATCGCTGGCGGTCGATAACGTCTTTGTCt                                  <  1:607011/36‑1 (MQ=255)
  ttGAGAAATCGCTGGCGGTCGATAACGTCTTTGTCt                                  <  1:7191652/36‑1 (MQ=255)
  ttGAGAAATCGCTGGCGGTCGATAACGTCTTTGTCt                                  <  1:914529/36‑1 (MQ=255)
   tGAGAAATCGCTGGCGGTCGATAACGTCTTTGTCTg                                 <  1:583314/36‑1 (MQ=255)
   tGAGAAATCGCTGGCGGTCGATAACGTCTTTGTCTg                                 <  1:1494042/36‑1 (MQ=255)
   tGAGAAATCGCTGGCGGTCGATAACGTCTTTGTCTg                                 <  1:5961970/36‑1 (MQ=255)
    gagaAATCGCTGGCGGTCGATAACGTCTTTTTCttg                                >  1:5578080/1‑34 (MQ=37)
    gagaAATCGCTGGCGGTCGATAACGTCTTTGTCTgg                                <  1:7396135/36‑1 (MQ=255)
     agaAATCGCTGGCGGTCGATAACGTCTTTTTCTGGc                               >  1:4900338/1‑36 (MQ=255)
      gaAATCGCTGGCGGTCGATAACGTCTTTGTCTGGCt                              <  1:6222508/36‑1 (MQ=255)
       aaaTCGCTGGCGGTCGATAACGTCTTTGTCTGGCTg                             <  1:4895774/36‑1 (MQ=255)
        aaTCGCTGGCGGTCGATAAGGTTTTTCTTTGGCTGa                            >  1:6623139/1‑36 (MQ=19)
        aaTCGCTGGCGGTCGATAACGTCTTTTTCTGGTTgt                            >  1:3016060/1‑35 (MQ=37)
        aaTCGCTGGCGGTCGATAACGTCTTTGTCTGGCTGa                            <  1:408832/36‑1 (MQ=255)
        aaTCGCTGGCGGTCCCTAAAGTCTTTTTCTGGCttg                            >  1:6301310/1‑34 (MQ=12)
         aTCGCTGGCGGTCGATAACGTCTTTGTCTGGCTGAt                           <  1:5399490/36‑1 (MQ=255)
          tCGCTGGCGGTCGATAACGTCTTTGTCTGGCTGAtg                          <  1:499139/36‑1 (MQ=255)
          tCGCTGGCGGGCGATAACGTCTTTGTCTGGCTGAtg                          >  1:726232/1‑36 (MQ=255)
            gCTGGCGGTCGATAACGTCTTTGTCTGGCTGAtgtt                        <  1:3739953/36‑1 (MQ=255)
            gCTGGCGGTCGATAACGTCTTTGTCTGGCTGAtgtt                        <  1:168334/36‑1 (MQ=255)
             cTGGCGGTCGATAACGTCTTTGTCTGGCTGAtgttg                       <  1:5373620/36‑1 (MQ=255)
             cTGGCGGTCGATAACGTCTTTGTCTGGCTGAtgttg                       <  1:6318746/36‑1 (MQ=255)
              tGGCGGTCGATAACGTCTTTTTCTGGCTGAtgtttt                      >  1:1399122/1‑34 (MQ=37)
              tGGCGGTCGATAACGTCTTTGTCTGGCTGAtgttgt                      <  1:4657177/36‑1 (MQ=255)
              tGGCGGTCGATAACGTCTTTGTCTGGCNNAttttgt                      >  1:7832761/1‑36 (MQ=38)
              tGGCGGTCGATAACGTATTTGTATGGCTGAtgttgt                      >  1:132539/1‑36 (MQ=37)
              gggCGGTCGATAACGTCTTTGTCTGGCTGAtgttgt                      <  1:3173321/35‑1 (MQ=255)
               ggCGGTCGATAACGTCTTTGTCTGGCTGATgttgtt                     <  1:5735996/36‑1 (MQ=255)
               ggCGGTCGATAACGTCTTTGTCTGGCTGATgttgtt                     <  1:674014/36‑1 (MQ=255)
               ggCGGTAGATAACGTCTTTGTCTGGCTGATgttgtt                     <  1:6456541/36‑1 (MQ=255)
                gCGGTCGATAACGTCTTTGTCTGGCTGATGTTGTTc                    <  1:3198418/36‑1 (MQ=255)
                gCGGTCGATAACGTCTTTGTCTGGCTGATGTTGTTc                    <  1:5918382/36‑1 (MQ=255)
                gCGGTCGATAACGTCTTTGTCTGGCTGATGTTGTTc                    <  1:24151/36‑1 (MQ=255)
                gCGGTCGATAACGTCTTTGTCTGGCTGATGTTGTTc                    <  1:1532957/36‑1 (MQ=255)
                 cGGTCGATAACGTCTTTGTCTGGCTGATGTTGTTCa                   <  1:3934686/36‑1 (MQ=255)
                 cGGTCGATAACGTCTTTGTCTGGCTGATGTTGTTCa                   <  1:2506245/36‑1 (MQ=255)
                 cGGTCGATAACGTCTTTGTCTGGCTGATGTTGTTCa                   >  1:432658/1‑36 (MQ=255)
                  ggTCGATAACGTCTTTGTCTGGCTGATGTTGTTCAg                  >  1:1002712/1‑36 (MQ=255)
                   gTCGATAACGTCTTTGTCTGGCTGATGTTGTTCAGc                 <  1:6981052/36‑1 (MQ=255)
                   gTCGATAACGTCTTTGTCTGGCTGATGTTGTTCAGc                 >  1:2399865/1‑36 (MQ=255)
                    tCGATAACGTCTTTGTCTGGCTGATGTTGTTCAGCt                <  1:7732077/36‑1 (MQ=255)
                     cGATAACGTCTTTGTCTGGCTGATGTTGTTCAGCTa               <  1:3529280/36‑1 (MQ=255)
                     cGATAACGTCTTTGTCTGGCTGATGTTGTTAAgttt               >  1:5213716/1‑33 (MQ=25)
                      gATAACGTCTTTGTCTGGCTGATGTTGTTCAGCTAt              >  1:326702/1‑36 (MQ=255)
                       aTAACGTCTTTGTCTGGCTGATGTTGTTCAGCTAtt             <  1:5726889/36‑1 (MQ=255)
                       aTAACGTCTTTGTCTGGCTGATGTTGTTCAGCTAtt             <  1:4151960/36‑1 (MQ=255)
                         aaCGTCTTTGTCTGGCTGATGTTGTTCAGTTATTtc           >  1:5655119/1‑36 (MQ=255)
                         aaCGTCTTTGTCTGGCTGATGTTGTTCAGCTATTtc           <  1:6800871/36‑1 (MQ=255)
                         aaCGTCTTTGTCTGGCTGATGTTGTTCAGCTATTtc           <  1:93770/36‑1 (MQ=255)
                         aaCGTCTTTGTCTGGCTGATGTTGTTAAGCTATTtc           >  1:2932442/1‑36 (MQ=255)
                          aCGTCTTTGTCTGGCTGATGTTGTTCAGCTATTtct          >  1:1907519/1‑36 (MQ=255)
                          aCGTCTTTGTCTGGCTGATGTTGTTCAGCTATTtct          <  1:3652211/36‑1 (MQ=255)
                           cGTCTTTGTCTGGCTGATGTTGTTCAGCTATTtctc         <  1:3694286/36‑1 (MQ=255)
                           cGTCTTTGTCTGGCTGATGTTGTTCAGCTATTtctc         <  1:5663715/36‑1 (MQ=255)
                            gTCTTTGTCTGGCTGATGTTGTTCAGCTATTTCTcc        >  1:1075075/1‑36 (MQ=255)
                            gTCTTTGTCTGGCTGATGTTGTTCAGCTATTTCTcc        <  1:6078454/36‑1 (MQ=255)
                            gTCTTTGTCTGGCTGATGTTGTTCAGCTATTTCTcc        <  1:6965694/36‑1 (MQ=255)
                            gTCTTTGTCTGGCTGATGTTGTTCAGCTATTTCTcc        <  1:684710/36‑1 (MQ=255)
                            gTCTTTGTCTGGCTGATATTGTTCAGCTATTTCTcc        >  1:2334321/1‑36 (MQ=255)
                             tCTTTGTCTGGCTGATGTTGTTCAGCTATTTCTcct       >  1:5038460/1‑35 (MQ=255)
                             tCTTTGTCTGGCTGATGTTGTTCAGCTATTTCTCCg       >  1:3065999/1‑36 (MQ=255)
                             tCTTTGTCTGGCTGATGTTGTTCAGCTATTTCTCCg       >  1:4958295/1‑36 (MQ=255)
                             tCTTTGTCTGGCTGATGTTGTTCAGCTATTTCTCCg       <  1:2426909/36‑1 (MQ=255)
                              cTTTGTCTGGCTGATGTTGTTCAGCTATTTCTCCGt      <  1:295178/36‑1 (MQ=255)
                              cTTTGTCTGGCTGATGTTGTTCAGCTATTTCTCCGt      <  1:5891503/36‑1 (MQ=255)
                               tttGTCTGGCTGATGTTGTTCAGCTATTTCTCCGtt     <  1:6628425/36‑1 (MQ=255)
                               tttGTCTGGCTGATGTTGTTCAGCTATTTCTCCGtt     <  1:3177095/36‑1 (MQ=255)
                               tttGTCTGGCTGATGTTGTTCAGCTATTTCTCCGtt     >  1:2102960/1‑36 (MQ=255)
                               tttGTCTGGCTGATGTTGTTCAGCTATTTCTCCGtt     >  1:4706488/1‑36 (MQ=255)
                               tttGTCTGGCTGATGTTGTTCAGCTATTTCTCCGtt     <  1:5249017/36‑1 (MQ=255)
                                ttGTCTGGCTGATGTTGTTCATCTATTTCTCCGTTc    <  1:4184673/36‑1 (MQ=255)
                                ttGTCTGGCTGATGTAGTTCAGCTATTTCTCCGTTc    <  1:2483667/36‑1 (MQ=255)
                                 tGTCTGGCTGATGTTGTTCAGCTATTTCTCCGTTcc   <  1:289156/36‑1 (MQ=255)
                                 tGTCTGGCTGATGTTGTTCAGCTATTTCTCCGTTcc   <  1:1697298/36‑1 (MQ=255)
                                  gTCTGGCTGATGTTGTTCAGCTATTTCTCCGTTCcg  <  1:5298655/36‑1 (MQ=255)
                                  |                                   
GATTGAGAAATCGCTGGCGGTCGATAACGTCTTTGTCTGGCTGATGTTGTTCAGCTATTTCTCCGTTCCG  >  REL606/3171648‑3171717

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG