Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606293,9420G→T16.8% 10.7 ‑10.6 ‑0.062intergenic (+212/‑315)ykgD/ykgEpredicted DNA‑binding transcriptional regulator/predicted oxidoreductase
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (11/0):  ref base (16/35):  total (27/35)
Fisher's exact test strand distribution p-value = 2.57e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.93e-01

TTACTTGATACGTGATTATTGTCGGTGATTTTTTTGTACACCATCAGTGTATATCTCAATAACCCCTGAAT  >  REL606/293907‑293977
                                   |                                   
ttACTTGATACGTGATTATTGTCGGTGATTTTTTTg                                     <  1:5471157/36‑1 (MQ=255)
 tACTTGATACGTGATTATTGTCGGTGATTTTTTTGt                                    <  1:4670968/36‑1 (MQ=255)
 tACTTGATACGTGATTATTGTCGGTGATTTTTTTGt                                    <  1:2115440/36‑1 (MQ=255)
  aCTTGATACGTGATTATTGTCGGTGATTTTTTTTTa                                   >  1:1215164/1‑36 (MQ=255)
  aCTTGATACGTGATTATTGTCGGTGATTTTTTTTTa                                   >  1:3928944/1‑36 (MQ=255)
   cTTGATACGTGATTATTGTCGGTGATTTTTTTTTac                                  >  1:1301169/1‑36 (MQ=255)
    ttGATACGTGATTATTGTCTGTGGTTTTTTTGTaca                                 <  1:418959/36‑1 (MQ=37)
     tGATACGTGATTATTGTCGGTGATTTTTTTTTacac                                >  1:4371147/1‑36 (MQ=255)
      gATACGTGATTATTGTCGGTGGTTTTTTTTTTCAcc                               >  1:4155323/1‑36 (MQ=21)
      gATACGTGATTATTGTCGGTGATTTTTTTTTacaca                               >  1:3153673/1‑35 (MQ=255)
        tACGTGATTATTGTCGGTGATTTTTTTGTACACCAt                             <  1:5449344/36‑1 (MQ=255)
        tACGTGATTATTGTCGGTGATTTTTTTGTACACCAt                             <  1:5901460/36‑1 (MQ=255)
         aCGTGATTATTGTCGGTGATTTTTTTTTACACCATc                            >  1:2710469/1‑36 (MQ=255)
         aCGTGATTATTGTCGGTGATTTTTTTGTACGCCATc                            <  1:742935/36‑1 (MQ=255)
         aCGTGATTATTGTCGGTGATTTTTTTGTACACCATc                            <  1:606688/36‑1 (MQ=255)
          cGTGATTATTGTCGGTGATTTTTTTGTACACCATCa                           <  1:3380989/36‑1 (MQ=255)
           gTGATTATTGTCGGTGATTTTTTTTTTCACCATCAg                          >  1:899151/1‑36 (MQ=37)
           gTGATTATTGTCGGTGATTTTTTTTTACACCATCAg                          >  1:4218753/1‑36 (MQ=255)
            tGATTATTGTCGGTGATTTTTTTGTACACCATCAgt                         <  1:1268058/36‑1 (MQ=255)
            tGATTATTGTCGGTGATTTTTTTGTACACCATCAgt                         <  1:3792041/36‑1 (MQ=255)
             gATTATTGTCGGTGCTTTTTTTTTACNNNCTCCgtg                        >  1:1987626/1‑36 (MQ=12)
              attattGTCGGTGATTTTTTTGTACACCATCAgtgt                       <  1:3724367/36‑1 (MQ=255)
                 attGTCGGTTATTTTTTTGTACTCCATCAGTGtata                    <  1:689310/36‑1 (MQ=37)
                 attGTCGGTGATTTTTTTGTACACCATCAGTGtata                    <  1:245961/36‑1 (MQ=255)
                 attGTCGGTGATTTTTTTGTACACCATCAGTGtata                    <  1:1529705/36‑1 (MQ=255)
                  ttGTCGGTGATTTTTTTGTACACCATCAGTGTatat                   <  1:5478128/36‑1 (MQ=255)
                   tGTCGGTGATTTTTTTTTACACCATCAGTGTATAtc                  >  1:3519962/1‑36 (MQ=255)
                   tGTCGGTGATTTTTTTGTACACCATCAGTGTATAtc                  >  1:5064231/1‑36 (MQ=255)
                   tGTCGGTGATTTTTTTGTACACCATCAGGGTATAtc                  >  1:3828694/1‑36 (MQ=255)
                    gTCGGTGATTTTTTTGTACACCATCAGTGTATAtct                 <  1:4279652/36‑1 (MQ=255)
                    gTCGGTGATTTTTTTGTACACCATCAGTGTATAtct                 <  1:3367697/36‑1 (MQ=255)
                    gTCGGTGATTTTTTTGTACACCATCAGTGTATAtct                 >  1:1873017/1‑36 (MQ=255)
                    gTCGGTGATTTTTTTGTACACCATCAGTGTATAtct                 <  1:2316658/36‑1 (MQ=255)
                     tCGGTGATTTTTTTGTACACCATCAGTGTATAtctc                >  1:3130558/1‑36 (MQ=255)
                     tCGGTGATTTTTTTGTACACCATCAGTGTATAtctc                <  1:3022744/36‑1 (MQ=255)
                     tCGGTGATTTTTTTGTACACCATCAGTGTATAtctc                <  1:3669029/36‑1 (MQ=255)
                     tCGGTGATTTTTTTGTACACCATCAGTGTATAtctc                <  1:5263087/36‑1 (MQ=255)
                     tCGGTGATTTTTTTGTACACCATCAGTGTATAtctc                >  1:4254772/1‑36 (MQ=255)
                       ggTGATTTTTTTGTACACCATCAGTGTATATCTCaa              <  1:2148031/36‑1 (MQ=255)
                       ggTGATTTTTTTGTACACCATCAGTGTATATCTCaa              >  1:2555552/1‑36 (MQ=255)
                       ggTGATTTTTTTGTACACCATCAGTGTATATCTCaa              <  1:5517641/36‑1 (MQ=255)
                       ggTGATTTTTTTGTACACCATCAGTGTATATCTCaa              >  1:1436272/1‑36 (MQ=255)
                       ggTGATTTTTTTGTACACCATCAGTGTATATCTCaa              >  1:666888/1‑36 (MQ=255)
                       ggTGATTTTTTTGTACACCATCAGGGTATATCTCaa              >  1:405700/1‑36 (MQ=255)
                        gTGATTTTTTTGTACACCATCAGTGTATATCTCAAt             <  1:2578318/36‑1 (MQ=255)
                         tGATTTTTTTGTACACCATCAGTGTATATCTCAATa            <  1:246454/36‑1 (MQ=255)
                         tGATTTTTTTGTACACCATCAGTGTATATCTCAATa            >  1:1518263/1‑36 (MQ=255)
                          gATTTTTTTGTACACCATCAGTTTATATCTCAATaa           <  1:526646/36‑1 (MQ=255)
                          gATTTTTTTGTACACCATCAGTGTATATCTCAATaa           >  1:3722879/1‑36 (MQ=255)
                           aTTTTTTTGTACACCATCAGTGTATATCTCAATAAc          <  1:1681804/36‑1 (MQ=255)
                            tttttttGTACACCATCAGTGTATATCTCAATAAcc         >  1:1235959/1‑36 (MQ=255)
                            tttttttGTACACCATCAGTGTATATCTCAATAAcc         <  1:2196491/36‑1 (MQ=255)
                             ttttttGTACACCATCAGTGTATATCTCAATAAccc        <  1:4419766/36‑1 (MQ=255)
                             ttttttGTACACCATCAGTGTATATCTCAATAAccc        <  1:4396191/36‑1 (MQ=255)
                             ttttttGTACACCATCAGTGTATATCTCAATAAccc        >  1:1721195/1‑36 (MQ=255)
                             ttttttGTACACCATCAGTGTATATCTCAATAAccc        <  1:2887366/36‑1 (MQ=255)
                             ttttttGTACACCATCAGTGTATATCTCAATAAccc        <  1:1308261/36‑1 (MQ=255)
                             ttttttGTACACCATCAGTGTATATCTCAATAAccc        <  1:2952195/36‑1 (MQ=255)
                              tttttGTACACCATCAGTGTATATCTCAATAAcccc       >  1:2336995/1‑36 (MQ=255)
                              tttttGTACACCATCAGTGTATATCTCAATAAcccc       >  1:4193250/1‑36 (MQ=255)
                              tttttGTACACCATCAGTGTATATCTCAATAAcccc       <  1:2953294/36‑1 (MQ=255)
                               ttttGTACACCATCAGTGTATATCTCAATAACCCCt      <  1:2662852/36‑1 (MQ=255)
                                  tGTACACCATCAGTGTATATCTCAATAACCCCTGaa   >  1:938740/1‑36 (MQ=255)
                                   gTACACCATCAGTGTATATCTCAATAACCCCTGAAt  >  1:1289257/1‑36 (MQ=255)
                                   gTACACCATCAGCGTATATACCAATAACCCCTGAAt  >  1:1189931/1‑36 (MQ=25)
                                   |                                   
TTACTTGATACGTGATTATTGTCGGTGATTTTTTTGTACACCATCAGTGTATATCTCAATAACCCCTGAAT  >  REL606/293907‑293977

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