Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,758,1371.→C 2.6% 0.2 ‑1.7 0.076coding (943/1263 nt)ygbKhypothetical protein
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (0/2):  ref base (44/30):  total (44/32)
Fisher's exact test strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GCCAGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.CCGCCAGCACTGACGCATTGGCAGCAATTCAACAG  >  REL606/2758102‑2758172
                                     |                                   
gCCAGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCa                                      <  1:7369003/36‑1 (MQ=255)
gCCAGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCa                                      <  1:7059614/36‑1 (MQ=255)
gCCAGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCa                                      >  1:1957972/1‑36 (MQ=255)
gCCAGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCa                                      >  1:785828/1‑36 (MQ=255)
 ccAGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.c                                    >  1:3102893/1‑36 (MQ=255)
 ccAGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.c                                    <  1:1092964/36‑1 (MQ=255)
  cAGGAAAGTGTACTTGCTCCACTGGTTTTTT.CCA.cc                                   >  1:2738010/1‑36 (MQ=255)
  cAGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.cc                                   >  1:7879575/1‑36 (MQ=255)
  cAGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.cc                                   >  1:4990204/1‑36 (MQ=255)
  cAGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.cc                                   >  1:6022696/1‑36 (MQ=255)
  cAGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.cc                                   <  1:5476418/36‑1 (MQ=255)
  cAGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.cc                                   <  1:7269472/36‑1 (MQ=255)
  cAGGAAAGTGTACTTGCTCCAATGGTTTTTG.CCaaa                                    >  1:5111974/1‑34 (MQ=38)
   aGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.CCg                                  >  1:3718586/1‑36 (MQ=255)
   aGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.CCg                                  <  1:7796119/36‑1 (MQ=255)
   aGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.CCg                                  <  1:5848084/36‑1 (MQ=255)
    ggAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.CCGc                                 >  1:604024/1‑36 (MQ=255)
    ggAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.CCGc                                 <  1:6752621/36‑1 (MQ=255)
    ggAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.CCGc                                 <  1:7100333/36‑1 (MQ=255)
     gAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.CCGcc                                >  1:6053122/1‑36 (MQ=255)
     gAAAGTGTACTTGCTCAACTGGTTCTTG.CCA.CCGcc                                <  1:5396010/36‑1 (MQ=38)
      aaaGTGTCCTTCCCCCACTGGTTTTTG.CCA.CCGCCa                               <  1:4026809/36‑1 (MQ=25)
      aaaGTGTACTTGCTCCACTGGTTTTTG.CCA.CCGCCa                               >  1:2777339/1‑36 (MQ=255)
      aaaGTGTACTTGCTCCACTGGTTTTTG.CCA.CCGCCa                               >  1:5019014/1‑36 (MQ=255)
       aaGTGTACTTGCTCCACTGGTTTTTT.CCA.CCGCCAg                              <  1:5152303/36‑1 (MQ=255)
       aaGTGTACTTGCTCCACTGGTTTTTG.CCA.CCGCCAg                              <  1:6691398/36‑1 (MQ=255)
       aaGTGTACTTGCTCCACTGGTTTTTG.CCA.CCGCCAg                              <  1:3690591/36‑1 (MQ=255)
        aGTGTACTTGCTCCACTGGTTTTTT.TCA.ACTCCAGc                             >  1:4761678/1‑36 (MQ=16)
        aGTGTACTTGCTCCACTGGTTTTTG.CCA.CCGCCAGc                             >  1:844401/1‑36 (MQ=255)
        aGTGTACTTGCTCCACTGGTTTTTG.CCA.CCGCCAGc                             >  1:7440733/1‑36 (MQ=255)
        aGTGTACTTGCTCCACTGGTTTTTG.CCA.CCGCCAGc                             >  1:7329049/1‑36 (MQ=255)
        aGTGTACTTGCTCCACTGGTTTTTG.CCA.CCGCCAGc                             >  1:4673006/1‑36 (MQ=255)
        aGTGTACTTGCTCCACTGGGTTTTG.CCA.CCGCAAGc                             >  1:1762984/1‑36 (MQ=37)
         gtgtACTTGCTCCACTGTTTTTTC.CCA.CCGCCAGCa                            <  1:7791472/36‑1 (MQ=37)
          tgtACTTGCTCCACTGGTTTTTG.CCA.CCGCCAGCAc                           >  1:2730053/1‑36 (MQ=255)
           gtACTTGCTCCACTGGTTTTTG.CCA.CCGCCAGCACt                          >  1:7442381/1‑36 (MQ=255)
            tccTTGCTCCACTGTTTTTTC.CCC.CCGCCAGCACTg                         <  1:3427310/34‑1 (MQ=21)
            tACTTGCTCCACTGGTTTTTG.CCA.CCGCCAGCACTg                         >  1:8218231/1‑36 (MQ=255)
             tcTTGCTCCCCTGTTTTTTG.CCA.CCGCCAGCACTGa                        <  1:1740011/35‑1 (MQ=37)
             aCTTGCTCCACTGGTTTTTG.CCA.CCGCCAGCACTGa                        >  1:3538705/1‑36 (MQ=255)
             aCTTGCTCCACTGGTTTTTG.CCA.CCGCCAGCACTGa                        >  1:2804205/1‑36 (MQ=255)
             aCTTGCTCCACTGGTTTTTG.CCA.CCGCCAGCACTGa                        >  1:2698054/1‑36 (MQ=255)
             aCTTGCTCCACTGGTTTTTG.CCA.ACGCCAGCGCTGa                        >  1:694887/1‑36 (MQ=37)
             aCTTGCTCCACTGGCTTTTG.CCA.CCGCCACCACTGa                        >  1:7490813/1‑36 (MQ=37)
               ttGCTCCACTGGTTTTTG.CCA.CCGCCAGCACTGACg                      >  1:1974329/1‑36 (MQ=255)
               ttGCTCCACTGGTTTTTG.CCA.CCGCCAGCACTGACg                      >  1:5402118/1‑36 (MQ=255)
               ttGCTCCACTGGTTTTTG.CCA.CCGCCAGCACTGACg                      >  1:1716668/1‑36 (MQ=255)
                tGCTCCACTGGTTTTTT.CCA.CCGCCAGCACTGAAGc                     >  1:7806353/1‑36 (MQ=37)
                  cTCAACTGGTTTTTG.CCA.CCGCCAGCACTGACGCAt                   >  1:5879853/1‑36 (MQ=255)
                   cccACTGGTTTTTC.CCA.CCGCCAGCACTGACGCAtt                  <  1:7171687/35‑1 (MQ=255)
                    ccACTGGTTTTTG.CCA.CCGCCAGCACTGACGCATTg                 >  1:3445610/1‑36 (MQ=255)
                    ccACTGGTTTTTG.CCA.CCGCCAGCACTGACGCATTg                 >  1:7061200/1‑36 (MQ=255)
                    ccACTGGTTTTTG.CCA.CCGCCAGCACTGACGCATTg                 >  1:4282943/1‑36 (MQ=255)
                    ccACTGGTTTTTG.CCA.CCGCCAGCACTGACGCATTg                 >  1:389637/1‑36 (MQ=255)
                    ccACTGGTTTTTG.CCA.CCGCCAGCACTGACGCATTg                 >  1:5683339/1‑36 (MQ=255)
                      aCTGTTTTTTC.CCC.CCGCCAGCACTGACGCATTGgc               <  1:5125436/36‑1 (MQ=25)
                      aCTGTTTTTGC.CCA.CCGCCAGCACTGACGCATTGgc               <  1:8023002/36‑1 (MQ=21)
                      aCTGGTTTTTG.CCA.CCGCCAGCACTGACGCATTGgc               >  1:4836144/1‑36 (MQ=255)
                      aCTGGTTTTTG.CCA.CCGCCAGCACTGACGCATTGgc               >  1:2787258/1‑36 (MQ=255)
                      aCTGGTTTTTG.CCA.CCGCCAGCACTGACGCATTGgc               >  1:3377321/1‑36 (MQ=255)
                       cTGTTTTTTC.CCA.CCGCCAGCACTGACGCATTGgca              <  1:3329908/36‑1 (MQ=37)
                       cTGGTTTTTG.CCA.CCGCCAGCACTGACGCATTGgca              >  1:3695396/1‑36 (MQ=255)
                       cTGGTTTTTG.CCA.CCGCCAGCACTGACGCATTGgca              >  1:3022067/1‑36 (MQ=255)
                       cTGGTTTTTG.CCA.CCGCCAGCACTGACGCATTGgca              >  1:1722114/1‑36 (MQ=255)
                        tGGTTTTTG.CCA.CCGCCAGCACTGACGCATTGgcag             >  1:3531549/1‑36 (MQ=255)
                        cggTTTTTGCCCA.CCGCCAGCACTGACGCATTGgca              <  1:4282758/35‑1 (MQ=25)
                        aggTTTTTC.CCCCCCGCCAGCACTGACGCATTGgca              <  1:3436260/35‑1 (MQ=16)
                         gttttttG.CCA.CCGCCAGCACTGACGCATTGgcagc            <  1:8023639/34‑1 (MQ=255)
                         gttttttG.CAA.CCGCCAGCACTGACGCATTGgcagc            <  1:199102/34‑1 (MQ=37)
                         ggTTTTTG.CGA.CCGCCAGCACTGACGCATTGgaagc            >  1:458462/1‑36 (MQ=37)
                         ggTTTTTC.CCCCCCGCCAGCACTGACGCATTGgcag             <  1:1379711/36‑1 (MQ=19)
                          ttttttC.CCA.CCGCCAGCACTGACGCATTGgcagca           <  1:2500505/35‑1 (MQ=255)
                          ttttttC.CCA.CCGCCAGCACTGACGCATTGgcagca           <  1:1600736/35‑1 (MQ=255)
                          tttttgG.CCA.CCGCCAGCACTGACGCATTGgcagca           <  1:973383/35‑1 (MQ=255)
                          gTTTTTG.CCA.CCGCCAGCACTGACGCATTGgcagca           >  1:3135944/1‑36 (MQ=255)
                           tttttG.CCA.CCGCCAGCACTGACGCATTGGCAGCaa          >  1:2140398/1‑36 (MQ=255)
                           tttttG.CCA.CCGCCAGCACTGACGCATTGGCAGCaa          >  1:7724962/1‑36 (MQ=255)
                            ttttC.CCA.CCGCCAGCACTGACGCATTGGCAGCAAt         <  1:7563098/36‑1 (MQ=255)
                             tttC.CCA.CCGCCAGCACTGACGCATTGGCAGCAAtt        <  1:8327417/36‑1 (MQ=255)
                              ttG.CCA.CCGCCAGCACTGACGCATTGGCAGCAATTc       >  1:1123633/1‑36 (MQ=255)
                               tccccA.CCGCCAGCACTGACGCATTGGCAGCAATTc       <  1:4315405/33‑1 (MQ=255)
                               tccccA.CCGCCAGCACTGACGCATTGGCAGCAATTc       <  1:5628387/33‑1 (MQ=255)
                                g.CCA.CCGCCAGCACTGACGCATTGGCAGCAATTCaa     >  1:8273832/1‑36 (MQ=255)
                                g.CCA.CCGCCAGCACTGACGCATTGGCAGCAATTCaa     >  1:4701631/1‑36 (MQ=255)
                                 cccA.CCGCCAGCACTGACGCATTGGCAGCAATTCaa     <  1:296269/35‑1 (MQ=255)
                                  ccA.CCGCCAGCACTGACGCATTGGCAGCAATTCAAc    <  1:7264802/36‑1 (MQ=255)
                                  ccA.CCGCCAGCACTGACGCATTGGCAGCAATTCAAc    >  1:3509921/1‑36 (MQ=255)
                                   cA.CCGCCAGCACTGACGCATTGGCAGCAATTCAACa   <  1:7954192/36‑1 (MQ=255)
                                   cA.CCGCCAGCACTGACGCATTGGCAGCAATTCAACa   <  1:67406/36‑1 (MQ=255)
                                    a.CCGCCAGCACTGACGCATTGGCAGCAATTCAACAg  >  1:2774849/1‑36 (MQ=255)
                                    a.CCGCCAGCACTGACGCATTGGCAGCAATTCAACAg  >  1:3803887/1‑36 (MQ=255)
                                    a.CCGCCAGCACTGACGCATTGGCAGCAATTCAACAg  >  1:1148192/1‑36 (MQ=255)
                                     |                                   
GCCAGGAAAGTGTACTTGCTCCACTGGTTTTTG.CCA.CCGCCAGCACTGACGCATTGGCAGCAATTCAACAG  >  REL606/2758102‑2758172

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG