Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,603,6880G→C 2.9% 0.1 0.0 ‑0.568835yhjDconserved inner membrane protein
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (1/1):  ref base (30/36):  total (31/37)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.15e-01

TTCTGGCTGCTATTGGTTTTGAAGTGATTAAAATCGTGATGACCTACACCCTGCCATCGTTGATGAAATCC  >  REL606/3603653‑3603723
                                   |                                   
ttCTGGCTGCTATTGGTTTTGAAGTGATTAAAATCg                                     <  1:1019494/36‑1 (MQ=255)
ttCTGGATGCTATTGGTTTTGAAGTGATTAAAATCg                                     <  1:2353592/36‑1 (MQ=255)
 tCTGGCTGCTATTGGTTTTGAAGTGATTAAAATCGt                                    <  1:4697696/36‑1 (MQ=255)
 tCTGGCTGCTATGGGTTTTGAAGTGATTAAAAGCGt                                    <  1:1596466/36‑1 (MQ=37)
  cTGGCTGCTATTGGTTTTGAAGTGATTAAAATCGtg                                   <  1:5097637/36‑1 (MQ=255)
   tGGCTGCTATTGGTTTTGAAGTGATTAAAATCGtga                                  <  1:5254352/36‑1 (MQ=255)
    ggCTGCTATTGGTTTTGAAGTGATTAAAATCGtgat                                 >  1:7218678/1‑36 (MQ=255)
    ggCTGCTATTGGTTTTGAAGTGATTAAAATCGtgat                                 <  1:6788251/36‑1 (MQ=255)
    ggCTGCTATTGGTTTTGAAGTGATTAAAATCGtgat                                 >  1:244429/1‑36 (MQ=255)
    ggCTGCTATTGGTTTTGAAGTGATTAAAATCGtgat                                 >  1:6063529/1‑36 (MQ=255)
     gctgctATTGGTTTTGAAGTGATTAAAATCGtgatg                                >  1:6386684/1‑36 (MQ=255)
      ctgctATTGGTTTTGAAGTGATTAAAATCGtgatga                               >  1:6166911/1‑36 (MQ=255)
       tgctATTGGTTTTGAAGTGATTAAAATCGTGATGAc                              <  1:4152338/36‑1 (MQ=255)
         ctATTGGTTTTGAAGTGATTAAAATCGTGGTGACCt                            <  1:4367423/36‑1 (MQ=255)
         ctATTGGTTTTGAAGTGATTAAAATCGTGATGACCt                            <  1:1506262/36‑1 (MQ=255)
          tATTGGTTTTGAAGTGATTAAAATCGTGATGACCTa                           >  1:7862072/1‑36 (MQ=255)
           aTTGGTTTTGAAGTGATTAAAATCGTGATGACCTac                          <  1:5110408/36‑1 (MQ=255)
           aTTGGTTTTGAAGTGATTAAAATCGTGATGACCTac                          >  1:3312050/1‑36 (MQ=255)
           aTTGGTTTTGAAGTGATTAAAATCGTGATGACCTac                          <  1:7594339/36‑1 (MQ=255)
           aTTGGTTTTGAAGTGATTAAAATCGTGATGACCTac                          >  1:7342311/1‑36 (MQ=255)
           aTTGGTTTTGAAGTGATTAAAATCCTGATGACCTac                          <  1:3733804/36‑1 (MQ=255)
            ttGGTTTTGAAGTGATTAAAATCGTGATGACCTaca                         >  1:1165289/1‑36 (MQ=255)
            ttGGTTTTGAAGTGATTAAAATCGTGATGACCTaca                         <  1:5497914/36‑1 (MQ=255)
             tGGTTTTGAAGTGATTAAAATCGTGATGACCTacac                        >  1:3437805/1‑36 (MQ=255)
             tGGTTTTGAAGTGATTAAAATCGTGATGACCTacac                        >  1:6079678/1‑36 (MQ=255)
              ggTTTTGAAGTGATTAAAATCGTGATGACCTACAcc                       <  1:460164/36‑1 (MQ=255)
              ggTTTTGAAGTGATTAAAATCGTGATGACCTACAcc                       >  1:7354744/1‑36 (MQ=255)
              ggTTTTGAAGTGATTAAAATCGTGATGACCTACAcc                       <  1:3896514/36‑1 (MQ=255)
               gTTTTGAAGTGATTAAAATCGTGATGACCTACAccc                      <  1:4551521/36‑1 (MQ=255)
                ttttGAAGTGATTAAAATCGTGATGACCTACACCCt                     >  1:4494694/1‑36 (MQ=255)
                ttttGAAGTGATTAAAATCGTGATGACCTACACCCt                     <  1:7809611/36‑1 (MQ=255)
                 tttGAAGTGATTAAAATCGTGATGACCTACACCCTg                    >  1:257438/1‑36 (MQ=255)
                  ttGAAGTGATTAAAATCGTGATGACCTACACCCTGc                   >  1:6244580/1‑36 (MQ=255)
                  ttGAAGTGATTAAAATCGTGATGACCTACACCCTGc                   >  1:502945/1‑36 (MQ=255)
                  ttGAAGTGATTAAAATCGTGATGACCTACACCCTGc                   <  1:4887141/36‑1 (MQ=255)
                  ttGAAGTGATTAAAATCGTGATGACCTACACCCTGc                   <  1:4723373/36‑1 (MQ=255)
                  ttGAAGTGATTAAAATCGTGATGACCTACACCCTGc                   <  1:230783/36‑1 (MQ=255)
                   tGAAGTGATTAAAATCGTGATGACCTACACCCTGcc                  >  1:770744/1‑36 (MQ=255)
                    gAAGTGATTAAAATCGTGATGACCTACACCCTGCCa                 <  1:4910822/36‑1 (MQ=255)
                    gAAGTGATTAAAATCGTGATGACCTACACCCTGCCa                 >  1:2403058/1‑36 (MQ=255)
                     aaGTGATTAAAATCGTGATGACCTACACCCTGCCAt                >  1:4945201/1‑36 (MQ=255)
                     aaGTGATTAAAATCGTGATGACCTACACCCTGCCAt                <  1:4711438/36‑1 (MQ=255)
                      aGTGATTAAAATCGTGATTACCTACACCCTGCCATc               <  1:6151588/36‑1 (MQ=255)
                      aGTGATTAAAATCGTGATGACCTACACCCTGCCATc               <  1:5805714/36‑1 (MQ=255)
                      aGTGATTAAAATCGTGATGAACTACACCCTGCAATc               >  1:1916023/1‑36 (MQ=37)
                      aGTGATTAAAATCCTGATGACCTACACCCTGCCATc               >  1:7175574/1‑36 (MQ=255)
                       gTGATTAAAATCGTGATGACCTACACCCTGCCATCg              >  1:6713595/1‑36 (MQ=255)
                        tGATTAAAATCGTGATGGCCTACACCCTGCCATCGt             <  1:1428364/36‑1 (MQ=255)
                        tGATTAAAATCGTGATGACCTACACCCTGCCATCGt             <  1:6366913/36‑1 (MQ=255)
                         gATTAAAATCGTGATGACCTACACCCTGCCATCGtt            <  1:5784217/36‑1 (MQ=255)
                         gATTAAAATCGTGATGACCTACACCCTGCCATCGtt            <  1:5362268/36‑1 (MQ=255)
                          aTTAAAATCGTGATTACCTACACCCTGCCATCGTtg           <  1:6760389/36‑1 (MQ=255)
                          aTTAAAATCGTGATGACCTACACCCTGCCATCGTtg           >  1:4620846/1‑36 (MQ=255)
                           ttAAAATCGTGATGACCTACACCCTGCCATCGTtga          <  1:4912069/36‑1 (MQ=255)
                            tAAAATCGTGATGACCTACACCCTGCCATCGTtgat         <  1:6950536/36‑1 (MQ=255)
                            tAAAATCGTGATGACCTACACCCTGCCATCGTtgat         <  1:5794303/36‑1 (MQ=255)
                            tAAAATCGTGATGACCTACACCCTGCCATCGTtgat         <  1:2305026/36‑1 (MQ=255)
                            tAAAATCGTGATGACCTACACCCTGCCATCGTtgat         <  1:2296018/36‑1 (MQ=255)
                             aaaaTCGTGATGACCTACACCCTGCCATCGTtgatg        >  1:243968/1‑36 (MQ=255)
                             aaaaTCGTGATGACCTACACCCTGCCATCGTtgatg        >  1:516144/1‑36 (MQ=255)
                              aaaTCGTGATGACCTATACCCTGCCATCGTtgatga       >  1:7600835/1‑36 (MQ=255)
                               aaacGTGATGTCCTACGCCCTGCCATCGTTGATGaa      <  1:4588605/33‑1 (MQ=25)
                               aaTCGTGATGACCTAGACCCTGACATCGTTGATGaa      >  1:5962480/1‑36 (MQ=38)
                               aaTCGTGATGACCTACACCCTGCCATCGTTGATGaa      <  1:4170541/36‑1 (MQ=255)
                               aaTCGTGATGACCTACACCCTGCCATCGTTGATGaa      >  1:3959124/1‑36 (MQ=255)
                                aTCGTGATGACCTACACCCTGCCATCGTTGATGaaa     >  1:812636/1‑36 (MQ=255)
                                aTCGTGATGACCTACACCCTGCCATCGTTGATGaaa     >  1:3637137/1‑36 (MQ=255)
                                aTCGTGATGACCTACACCCTCCCATCGTTGATGaaa     >  1:2691244/1‑36 (MQ=255)
                                  cGTGATGACCTACACCCTGCCATCGTTGATGAAATc   <  1:229060/36‑1 (MQ=255)
                                  cGTGATGACCTACACCCTGCCATCGTTGATGAAATc   <  1:1211364/36‑1 (MQ=255)
                                   gTGATGACCTACACCCTGCCATCGTTGATGAAATcc  >  1:4057266/1‑36 (MQ=255)
                                   |                                   
TTCTGGCTGCTATTGGTTTTGAAGTGATTAAAATCGTGATGACCTACACCCTGCCATCGTTGATGAAATCC  >  REL606/3603653‑3603723

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 37 ≤  ATCG