Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,768,2030G→T 9.5% 2.5 ‑10.4 0.085intergenic (‑5/+45)ygbE/cysCconserved inner membrane protein/adenylylsulfate kinase
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (10/0):  ref base (23/52):  total (33/52)
Fisher's exact test strand distribution p-value = 2.96e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TTGTTAGTGTAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCC  >  REL606/2768168‑2768238
                                   |                                   
ttGTTAGTGTAATGTTATGGCTATTACGCATATTTg                                     <  1:42803/36‑1 (MQ=255)
ttGTTAGTGTAATGTTATGGCTATTACGCATATTTg                                     <  1:3601768/36‑1 (MQ=255)
  gTTAGTGTAATGTTATGGCTATTACGCCAATTTTAt                                   >  1:2556261/1‑36 (MQ=25)
  gTTAGTGTAATGTTATGGCTATTACGCATATTTGAt                                   <  1:3583497/36‑1 (MQ=255)
  gTTAGTGTAATGTTATGGCTATTACGCATATTTGAt                                   <  1:4630183/36‑1 (MQ=255)
  gTTAGTGTAATGTTATGGCTATTACGCATATTTGAt                                   <  1:5118570/36‑1 (MQ=255)
   ttAGTGTAATGTTATGGCTATTACGCATATTTTATc                                  >  1:3651341/1‑36 (MQ=255)
   ttAGTGTAATGTTATGGCTAATACCCATATTTTATc                                  >  1:875373/1‑36 (MQ=25)
    tAGTGTAATGTTATGGCTATTTCGCATATTTGATcc                                 <  1:3617659/36‑1 (MQ=255)
    tAGTGTAATGTTATGGCTATTACGCATATTTTAttc                                 >  1:804270/1‑34 (MQ=37)
    tAGTGTAATGTTATGGCTATTACGCATATTTGATcc                                 <  1:385517/36‑1 (MQ=255)
    tAGTGTAATGTTATGGCTATTACGCATATTTGATcc                                 <  1:5801844/36‑1 (MQ=255)
    tAGTGTAATGTTATGGCTATTACGCATATTTGATcc                                 <  1:2956918/36‑1 (MQ=255)
    tAGTGTAATGTTATGGCTATTAAGCATATTTTATca                                 >  1:5356116/1‑35 (MQ=37)
     aGTGTAATGTTATGGCTATTACGCATATTTGATCCt                                <  1:5437385/36‑1 (MQ=255)
     aGTGTAATGTTATGGCTATTACGCATATTTGATCCt                                <  1:1018564/36‑1 (MQ=255)
     aGTGTAATGTTATGGCTATTAAGCATATTTTATCCt                                >  1:530074/1‑36 (MQ=37)
     aGTGTAATGTGATGGCTATTACGCATATTTGATCCt                                <  1:5507213/36‑1 (MQ=255)
      gtgtAATGTTATGGCTATTAAGCATATTTTATTCtt                               >  1:2570012/1‑35 (MQ=21)
       tgtAATGTTATGGCTATTACGCATATTTGATCCTGt                              <  1:5487230/36‑1 (MQ=255)
       tgtAATGTTATGGCTATTACGCATATTTGATCCTGt                              <  1:4545347/36‑1 (MQ=255)
       tgtAATGTTATGGCTATTACGCATATTTGATCCTGt                              <  1:1085631/36‑1 (MQ=255)
       tgtAATGTTATGGCTATTACGCATATTTGATCCTGt                              <  1:4758978/36‑1 (MQ=255)
       tgtAATGTTATGGCTAATACTCATATTTTATTCttt                              >  1:2863667/1‑34 (MQ=16)
       tgtAATGTTATGGCTAATAAGCATATTTGATCCttt                              >  1:5465812/1‑34 (MQ=37)
        gtAATGTTATGGCTATTACGCATATTTGATCCTGtt                             <  1:1694459/36‑1 (MQ=255)
        gtAATGTTATGGCTATTACGCATATTTGATCCTGtt                             <  1:4704975/36‑1 (MQ=255)
        gtAATGTTATGGCTAATTCGCATATATGATCCTTtt                             >  1:1434244/1‑36 (MQ=19)
         tAATGTTATGGCTATTACGCATATTTGATCCTGTTa                            <  1:2474385/36‑1 (MQ=255)
         tAATGTTATGGCTATTACGCATATTTGATCCTGTTa                            <  1:2087476/36‑1 (MQ=255)
          aaTGTTATGGCTATTACGCATATTTTATCCTTtttc                           >  1:224363/1‑34 (MQ=25)
          aaTGTTATGGCTATTACGCATATTTGATCCTGTTAc                           <  1:60961/36‑1 (MQ=255)
           aTGTTATGGCTATTACGCATATTTGATCCTGTTACt                          <  1:4298412/36‑1 (MQ=255)
           aTGTTATGGCTATTACGCATATTTGATCCTGTTACt                          <  1:1890381/36‑1 (MQ=255)
           aTGTTATGGCTATTACGCATATTTGATCCTGTTACt                          <  1:2262206/36‑1 (MQ=255)
            tGTTATGGCTATTACGCATATTTGCTCCTGTTTCtt                         <  1:5767720/36‑1 (MQ=37)
             gTTATGGCTATTACGCATATTTGATCCTGTTACttt                        <  1:383783/36‑1 (MQ=255)
             gTTATGGCTATTACGCATATTTGATCCTGTTACttt                        <  1:4472580/36‑1 (MQ=255)
             gTTATGGCTATTACGCATATTTGATCCTGTTACttt                        <  1:4711959/36‑1 (MQ=255)
             gTTATGGCTATTACGCATATTTGATCCTGTTACttt                        <  1:2217971/36‑1 (MQ=255)
              ttATGGCTATTACGCATATTTGATCCTGTTACTTTg                       <  1:645285/36‑1 (MQ=255)
               tATGGCTATTACGCATATTTTATCCTGTTACtttta                      >  1:753336/1‑34 (MQ=38)
               tATGGCTATTACGCATATTTGATCCTGTTACTTTGa                      <  1:421654/36‑1 (MQ=255)
                aTGGCTATTACGCATATTTGATCCTGTTACTTTGAc                     <  1:3249875/36‑1 (MQ=255)
                aTGGCTATTACGCATATTTGATCCTGTTACTTTGAc                     <  1:5643338/36‑1 (MQ=255)
                aTGGCTATTACGCATATTTGATCCTGTTACTTTGAc                     >  1:4164678/1‑36 (MQ=255)
                aTGGCTAATGCGCATATTTGATCCTGTTACTTTGAc                     <  1:546172/36‑1 (MQ=37)
                 tGGCTATTACGCATATTTGATCCTGTTACTTTGACg                    >  1:4019160/1‑36 (MQ=255)
                   gCTATTACGCATATTTGATCCTGTTACTTTGACGaa                  >  1:3977870/1‑36 (MQ=255)
                   gCTATTACGCATATTTGATCCTGTTACTTTGACGaa                  <  1:1333918/36‑1 (MQ=255)
                   gCTATTACGCATATTTGATCCTGTTACTTTGACGaa                  <  1:3378148/36‑1 (MQ=255)
                   gCTATTACGCATATTTGATCCTGTTACTTTGACGaa                  >  1:2759003/1‑36 (MQ=255)
                   gCTATTACGCATATTTGATCCTGTTAATTTGACGaa                  >  1:4023568/1‑36 (MQ=255)
                    cTATTACGCATATTTGATCCTGTTACTTTGACGAAc                 >  1:824952/1‑36 (MQ=255)
                    cTATTACGCATATTTGATCCTGTTACTTTGACGAAc                 <  1:1218515/36‑1 (MQ=255)
                    cTATTACGCATATTTGATCCTGTTACTTTGACGAAc                 <  1:214090/36‑1 (MQ=255)
                     tATTACGCATATTTGATCCTGTTACTTTGACGAAcc                <  1:964447/36‑1 (MQ=255)
                     tATTACGCATATTTGATCCTGTTACTTTGACGAAcc                <  1:1032148/36‑1 (MQ=255)
                     tATTACGCATATTTGATCCTGTTACTTTGACGAAcc                <  1:66345/36‑1 (MQ=255)
                      aTTACGCATATTTGCTCCTGTTACTTTGACGAACCg               <  1:5256241/36‑1 (MQ=255)
                      aTTACGCATATTTGATCCTGTTACTTTGACGAACCg               >  1:4177002/1‑36 (MQ=255)
                      aTTACGCATATTTGATCCTGTTACTTTGACGAACCg               >  1:1885279/1‑36 (MQ=255)
                       ttACGCATATTTGATCCTGTTACTTTGACGAACCgg              >  1:2803240/1‑36 (MQ=255)
                       ttACGCATATTTGATCCTGTTACTTTGACGAACCgg              >  1:2269533/1‑36 (MQ=255)
                        tACGCATATTTGATCCTGTTACTTTGACGAACCggg             >  1:870501/1‑36 (MQ=255)
                        tACGCATATTTGATCCTGTTACTTTGACGAAACggg             >  1:4853575/1‑36 (MQ=255)
                         aCGCATATTTGATCCTGTTACTTTGACGAACCGGGc            >  1:2509976/1‑36 (MQ=255)
                           gCATATTTGATCCTGTTACTTTTACGAACCGGGCAt          <  1:2951062/36‑1 (MQ=255)
                           gCATATTTGATCCTGTTACTTTGACGAACCGGGCAt          <  1:3619391/36‑1 (MQ=255)
                           gCATATTTGATCCTGTTACTTTGACGAACCGGGCAt          >  1:5559265/1‑36 (MQ=255)
                           gCATATTTGATCCTGTTACTTTGACGAACCGGGCAt          <  1:4602067/36‑1 (MQ=255)
                           gCATATTTGATCCTGTTACTGTGACGAAACGGGGAt          >  1:2377032/1‑36 (MQ=37)
                            cATATTTGATCCTGTTACTTTGACGAACCGGGCATg         <  1:885511/36‑1 (MQ=255)
                            cATATTTGATCCTGTTACTTTGACGAACCGGGCATg         <  1:3002551/36‑1 (MQ=255)
                            cATATTTGATCCTGTTACTTTGACGAACCGGGCATg         >  1:5689302/1‑36 (MQ=255)
                            cATATTTGATCCTGTTACTTTGACGAACCGGGCATg         <  1:4332479/36‑1 (MQ=255)
                             atatTTGATCCTGTTCCTTTGACGAACCGGGCATgg        <  1:4571546/36‑1 (MQ=255)
                             atatTTGATCCTGTTACTTTGACGAACCGGGCATgg        >  1:3689177/1‑36 (MQ=255)
                             atatTTGATCCTGTTACTTTGACGAACCGGGCATgg        >  1:284445/1‑36 (MQ=255)
                              tatTTGATCCTGTTACTTTGACTAACCGGGCATGGa       <  1:1968238/36‑1 (MQ=255)
                               atTTGATCCTGTTACTTTGACGAACCGGGCATGGaa      <  1:5818526/36‑1 (MQ=255)
                               atTTGATCCTGTTACTTTGACGAACCGGGCATGGaa      >  1:2340161/1‑36 (MQ=255)
                               atTTGATCCTGTTACTTTGACGAACCGGGCATGGaa      >  1:275846/1‑36 (MQ=255)
                                tttGATCCTGTTACTTTGACGAACCGGGCATGGaaa     <  1:1471334/36‑1 (MQ=255)
                                tttGATCCTGTTACTTTGACGAACCGGGCATGGaaa     >  1:5191111/1‑36 (MQ=255)
                                 ttGATTCTGTTACTTTGACTACCCGGGCATGGAAAc    <  1:1498496/36‑1 (MQ=25)
                                 ttGATCCTGTTACTTTGACGAACCGGGCATGGAAAc    <  1:3420353/36‑1 (MQ=255)
                                   gATCCTGTTACTTTGACGAACCGGGCATGGAAAccc  <  1:1564766/36‑1 (MQ=255)
                                   gATCCTGTTACTTTGACGAACCGGGCATGGAAAccc  >  1:604849/1‑36 (MQ=255)
                                   gATCCTGTTACTTTGACGAACCGGGCATGGAAAccc  >  1:4348973/1‑36 (MQ=255)
                                   gATCCTGTTACTTTGACGAACCGGGCATGGAAAccc  <  1:2841623/36‑1 (MQ=255)
                                   gATCCTGTTACTTTGACGAACCGGGCATGGAAAccc  >  1:5055288/1‑36 (MQ=255)
                                   gATCCTGTTACTTTGACGAACCGGGCATGGAAAccc  <  1:3858514/36‑1 (MQ=255)
                                   gATCCTGTTACTTTGACGAACCGGGCATGGAAAccc  >  1:896148/1‑36 (MQ=255)
                                   |                                   
TTGTTAGTGTAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCC  >  REL606/2768168‑2768238

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