Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,585,8130C→A 7.7% 0.2 ‑5.9 ‑0.057752yfhRpredicted peptidase
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (38/14):  total (38/19)
Fisher's exact test strand distribution p-value = 2.78e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.47e-01

CCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATG  >  REL606/2585778‑2585848
                                   |                                   
cccATGGCAGCACAGCGAAAAGTTGTATAGCCTGGc                                     >  1:4147985/1‑36 (MQ=255)
 ccATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCa                                    <  1:4147770/36‑1 (MQ=255)
  cATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCaa                                   >  1:1916426/1‑36 (MQ=255)
  cATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCaa                                   >  1:1002072/1‑36 (MQ=255)
    tGGCAGCACAGCGAAAAGTTGTATAGCCTGGCaaaa                                 <  1:5345658/36‑1 (MQ=255)
     ggCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAag                                >  1:381347/1‑36 (MQ=255)
     ggCAGCACAGCGAAAAGCTGTATAGCCTGGCAAAag                                >  1:4168452/1‑36 (MQ=255)
      gcagcaCAGCGAAAAGTTGTATAGCCTGGCAAAaga                               >  1:5496500/1‑36 (MQ=255)
       cagcaCAGCGAAAAGTTGTATAGCCTGGCAAAagag                              <  1:4794460/36‑1 (MQ=255)
       cagcaCAGCGAAAAGTTGTATAGCCTGGCAAAagag                              >  1:5579036/1‑36 (MQ=255)
        agcaCAGCGAAAAGTTGTATAGCCTGGCAAAAGAGc                             >  1:963518/1‑36 (MQ=255)
        agcaCAGCGAAAAGTTGTATAGCCTGGCAAAAGAGc                             >  1:1631632/1‑36 (MQ=255)
        agcaCAGCGAAAAGTTGTATAGCCTGGCAAAAGAGc                             >  1:5232711/1‑36 (MQ=255)
         gcaCAGCGAAAAGTTGTATAGCCTGGCAAAAGAGcc                            >  1:4400398/1‑36 (MQ=255)
          cacaGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCa                           <  1:5678951/36‑1 (MQ=255)
           acaGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCaa                          >  1:4926619/1‑36 (MQ=255)
           acaGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCaa                          >  1:5383683/1‑36 (MQ=255)
           acaGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCaa                          >  1:5665726/1‑36 (MQ=255)
             aGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCaaaa                        <  1:2783004/36‑1 (MQ=255)
             aGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCaaaa                        >  1:4308913/1‑36 (MQ=255)
             aGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCaaaa                        >  1:573803/1‑36 (MQ=255)
              ggaaaaaGTTGAATAGCCGGGAAAAAGAGCCAAAAc                       <  1:1619556/33‑1 (MQ=16)
              gaaaaaaTTTGTATAGCCTGGAAAAAGAGCCAAAAc                       <  1:2360972/33‑1 (MQ=25)
              gCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAAc                       >  1:435355/1‑36 (MQ=255)
               cGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACg                      >  1:2875280/1‑36 (MQ=255)
                gAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACgg                     >  1:1072666/1‑36 (MQ=255)
                  aaaGTTTTATAGCCTGGAAAAAGAGCCAAAACGTtt                   <  1:3796968/36‑1 (MQ=25)
                  aaaGTTGTATAGCCTGGCAAAAGAGCCAAAACGGtt                   <  1:5876419/36‑1 (MQ=255)
                   aaGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTg                  <  1:4798199/36‑1 (MQ=255)
                    aGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGa                 >  1:2983583/1‑36 (MQ=255)
                     gTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGAt                >  1:4759317/1‑36 (MQ=255)
                     gTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGAt                >  1:1677778/1‑36 (MQ=255)
                     gTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGAt                >  1:5205346/1‑36 (MQ=255)
                      ttGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATc               >  1:3065649/1‑36 (MQ=255)
                       tGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATcc              >  1:3883175/1‑36 (MQ=255)
                       tGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATcc              >  1:2532401/1‑36 (MQ=255)
                        gTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCt             <  1:118421/36‑1 (MQ=255)
                        gTATACCCTGGAAAAAGAGCCAAAACGGTTGATCCt             <  1:4297992/36‑1 (MQ=37)
                         tataGCCTGGCAAAAGAGCCAAAACGGTTGATCCTa            <  1:3397795/36‑1 (MQ=255)
                         tataGCCTGGCAAAAGAGCCAAAACGGTTGATCCTa            <  1:5937429/36‑1 (MQ=255)
                          ataGCCTGGCAAAAGAGCCAAAACGGTTGATCCTaa           <  1:3963234/36‑1 (MQ=255)
                           taGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAAt          >  1:537524/1‑36 (MQ=255)
                           taGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAAt          >  1:3396892/1‑36 (MQ=255)
                            aGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATc         >  1:4118772/1‑36 (MQ=255)
                            aGCATGGCAAAAGAGCCAAAACGGTTGATCCTAATc         <  1:5698641/36‑1 (MQ=255)
                             gacTGGAAAAAGAGCCATAACGGTTGATCCTAATcc        <  1:513104/34‑1 (MQ=37)
                             gCCTGGCAAAAGAGCCAAAACGGTTGATCCTCATcc        >  1:251650/1‑36 (MQ=255)
                             gCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATcc        >  1:4869809/1‑36 (MQ=255)
                             gCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATcc        >  1:3731765/1‑36 (MQ=255)
                             gCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATcc        >  1:603769/1‑36 (MQ=255)
                             gCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATcc        >  1:1684602/1‑36 (MQ=255)
                             gCCTGGCAAAAGAGCCAAAACGGTTGATCATAATcc        >  1:1499712/1‑36 (MQ=255)
                               cTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCg      <  1:2882061/36‑1 (MQ=255)
                               cTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCg      >  1:828606/1‑36 (MQ=255)
                                 ggCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGa    >  1:5451149/1‑36 (MQ=255)
                                  gCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGAt   >  1:2597708/1‑36 (MQ=255)
                                  gCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGAt   <  1:3989561/36‑1 (MQ=255)
                                  gCAAAAGAGCCAAAACGGTTGATCCTAATCCAGGAt   >  1:1748543/1‑36 (MQ=255)
                                   cAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATg  >  1:1917220/1‑36 (MQ=255)
                                   cAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATg  <  1:599095/36‑1 (MQ=255)
                                   cAAAAGAGCCAAAACGGTTGATCCTAATCCAGGATg  >  1:530656/1‑36 (MQ=255)
                                   |                                   
CCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATG  >  REL606/2585778‑2585848

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