Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL60666,6110G→T14.3% 3.1 ‑20.6 0.0801974polBDNA polymerase II
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (17/0):  ref base (13/50):  total (30/50)
Fisher's exact test strand distribution p-value = 1.18e-09
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATAGCTCCTGCTGAAACT  >  REL606/66576‑66645
                                   |                                  
aTGGTTTCGCGTACATATTCCTGATATGGCTCGTTg                                    <  1:5337129/36‑1 (MQ=255)
aTGGTTTCGCGTACATATTCCTGATATGGCTCGTTg                                    >  1:2396183/1‑36 (MQ=255)
 tGGTTTCGCGTACATATTCCTGATATGGCTCGTTGc                                   <  1:3353013/36‑1 (MQ=255)
  ggTTTCGCGTACATATTCCTGATATGGCTCGTTTCg                                  >  1:7119053/1‑36 (MQ=255)
   gTTTCGCGTACATATTCCTGATATGGCTCGTTTCgg                                 >  1:6295490/1‑36 (MQ=255)
    tttCGCGTACATATTCCTGATATGGCTCGTTGCGGa                                <  1:5682333/36‑1 (MQ=255)
    tttCGCGTACATATTCCTGATATGGCTCCTTTCggg                                >  1:2894382/1‑35 (MQ=37)
     ttCGCGTACATATTCCTGATATGGCTCGTTTCggta                               >  1:2627160/1‑34 (MQ=38)
     ttCGCGTACATATTCCTGATATGGCTCGTTGCGGaa                               <  1:6826930/36‑1 (MQ=255)
     ttCGCGTACATATTCCTGATATGGCTCGTTGCGGaa                               <  1:2605224/36‑1 (MQ=255)
        gcgTACATATTCCTTATATGGCTCGTTGCGGAAGAt                            <  1:6937716/36‑1 (MQ=255)
        gcgTACATATTCCTGATATGGCTCGTTTCGGAAGAt                            >  1:986880/1‑36 (MQ=255)
        gcgTACATATTCCTGATATGGCTCGTTGCGGAAGAt                            <  1:4861465/36‑1 (MQ=255)
        gcgTACATATTCCTGATATGGCTCGTTGCGGAAGAt                            <  1:6452768/36‑1 (MQ=255)
        gcgTACATATTCCTGATATGGCTCGTTGCGGAAGAt                            <  1:7104856/36‑1 (MQ=255)
        gcgTACATATTCCTGATATGGCTCGTTGCGGAAGAt                            >  1:5481938/1‑36 (MQ=255)
        gcatACATATTCCTGATATGGCTAGTTGCGGAAGAt                            <  1:7262210/33‑1 (MQ=38)
         cgTACATATTCCTGATATGGCTCGTTGCGGAAGATg                           <  1:6581383/36‑1 (MQ=255)
         cgTACATATTCCTGATATGGCTCGTTGCGGAAGATg                           <  1:5343191/36‑1 (MQ=255)
          gTACATATTCCTGATATGGCTCGTTGCGGAAGATgc                          <  1:4071355/36‑1 (MQ=255)
          gTACATATTCCTGATATGGCTACTTTAGGGAGATgc                          >  1:5403081/1‑36 (MQ=16)
           tACATATTCCTGATATGGCTCGTTTCGGAAGAttgg                         >  1:4742404/1‑33 (MQ=37)
           tACATATTCCTGATATGGCTCGTTTCGGAAGATgcg                         >  1:2815411/1‑36 (MQ=255)
           tACATATTCCTGATATGGCTCGTTGCGGAAGATgcg                         <  1:4772442/36‑1 (MQ=255)
             cATATTCCTGATATGGCTCGTTTCGGTAGATTCGCa                       >  1:872977/1‑36 (MQ=25)
             cATATTCCTGATATGGCTCGTTTCGGGAGATTCGCa                       >  1:56316/1‑36 (MQ=37)
             cATATTCCTGATATGGCTCGTTGCGGAAGATGCGCa                       <  1:5581238/36‑1 (MQ=255)
             cATATTCCTGATATGGCTCGTTGCGGAAGATGCGCa                       <  1:4923695/36‑1 (MQ=255)
             cATATTCCTGATATGGCTCGTTGCGGAAGATGCGCa                       <  1:2739038/36‑1 (MQ=255)
              atatTCCTGATATGGCTCGTTGCGGAAGATGCGCAg                      <  1:2372697/36‑1 (MQ=255)
               tatTCCTGATATGGCTCGTTGCGGTAGATGTGCAgg                     >  1:1107693/1‑36 (MQ=37)
               tatTCCTGATATGGCTCGTTGCGGAAGATGCGCAgg                     <  1:6576141/36‑1 (MQ=255)
               tatTCCTGATATGGCTCGTTGCGGAAGATGCGCAgg                     <  1:3736828/36‑1 (MQ=255)
               tatTCCTGATATGGCTCCCTTCGGGCGATGTGCAgg                     >  1:4006021/1‑36 (MQ=9)
                ataccTGATATGGCTAGTTGCGGAAGATGCGCAGGt                    <  1:4948107/33‑1 (MQ=37)
                atTCCTGATATGGCTCTTTTCGGCATATGCGCAGGt                    >  1:3560091/1‑36 (MQ=19)
                atTCCTGATATGGCTCGTTTCGGAAGATGCGCAGGt                    >  1:4871868/1‑36 (MQ=255)
                atTCCTGATATGGCTCGTTGCGGAAGATGCGCCGGt                    >  1:6135434/1‑36 (MQ=255)
                atTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGt                    <  1:7179074/36‑1 (MQ=255)
                atTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGt                    <  1:2298345/36‑1 (MQ=255)
                 ttCCTGATATGGCTCGTTTCGGAAGATGCGCAggga                   >  1:1185761/1‑34 (MQ=37)
                 ttCCTGATATGGCTCGTTGCGGGAGATGCGCAGGta                   >  1:556517/1‑36 (MQ=255)
                 ttCCTGATATGGCTCGTTGCGGAAGATGCGCAGGta                   <  1:1961547/36‑1 (MQ=255)
                 ttCCTGATATGGCTCCTTGCGGGAGATGCGCCggga                   >  1:3754827/1‑34 (MQ=21)
                  tCCTGATATGGCTCGTTTCGGAAGATGCGCAGGtat                  >  1:5083431/1‑36 (MQ=255)
                  tCCTGATATGGCTCGTTGCGGAAGATGCGCAGGtat                  <  1:5871088/36‑1 (MQ=255)
                  tCCTGATATGGCTCGTTGCGGAAGATGCGCAGGtat                  >  1:7460088/1‑36 (MQ=255)
                  tCCTGATATGGCTCGTTGCGGAAGATGCGCAGGtat                  <  1:3076064/36‑1 (MQ=255)
                  tCCTGATATGGCTCGCTTCGGGAGGTGCGCTGGtat                  >  1:5601512/1‑36 (MQ=16)
                  tCCTGATATGGCTCCCTGCGGAAGATGCGCAGGtat                  >  1:6253392/1‑36 (MQ=25)
                   ccTGATATGGCTCGTTGCGGAAGATGCGCAGGtata                 <  1:3514122/36‑1 (MQ=255)
                   ccTGATATGGCTCGTTGCGGAAGATGCGCAGGtata                 <  1:5491689/36‑1 (MQ=255)
                   ccTGATATGGCTCGTTGCGGAAGATGCGCAGGtata                 <  1:781524/36‑1 (MQ=255)
                   ccTGATATGGCTCGTTGCGGAAGATGCGCAGGtata                 <  1:4666252/36‑1 (MQ=255)
                   ccTGATATGGCTCGTTGCGGAAGATGCGCAGGtata                 <  1:6744392/36‑1 (MQ=255)
                   ccTGATATGGCTCGTTGCGGAAGATGCGCAGGtata                 <  1:7480701/36‑1 (MQ=255)
                   ccTGATATGGCTCCTTTCGGGAGATTGGCAGGgatt                 >  1:4435985/1‑35 (MQ=9)
                    cTGATATGGCTCGTTGCGGAAGATGCGCAGGTATAg                <  1:5424781/36‑1 (MQ=255)
                     tGATATGGCTCGTTGCGGAAGATGCGCAGGTATAGc               <  1:7562612/36‑1 (MQ=255)
                       ctatGGCTCGTTGCGGAAGATGCGCAGGTATAGCTc             <  1:5984617/35‑1 (MQ=255)
                       atatGGCTCGTTGCGGAAGATGCGCAGGTATAGCTc             <  1:448322/36‑1 (MQ=255)
                       atatGGCTCGTTGCGGAAGATGCGCAGGTATAGCTc             >  1:7746216/1‑36 (MQ=255)
                        tatGGCTCGTTGCGGAAGATGCGCAGGTATAGCTcc            <  1:6295813/36‑1 (MQ=255)
                         atGGCTCGTTGCGGAAGATGCGCAGGTATAGCTCct           <  1:5293799/36‑1 (MQ=255)
                         atGGCTCGTTGCGGAAGATGCGCAGGTATAGCTCct           <  1:4396817/36‑1 (MQ=255)
                          tGGCTCGTTGCGGCAGATGCGCAGGTATAGCTCctg          >  1:2673056/1‑36 (MQ=255)
                          tGGCTCGTTGCGGAAGATGCGCAGGTATAGCTCctg          <  1:5853064/36‑1 (MQ=255)
                           ggCTCGTTGCGGAAGATGCGCAGGTATAGCTCctgc         <  1:4194134/36‑1 (MQ=255)
                            gCTCGTTGCTGAAGATGCGCAGTTATAGCTCctgct        <  1:7508249/36‑1 (MQ=38)
                            gCTCGTTGCGGAAGATGCGCAGGTATAGCTCctgct        >  1:2334352/1‑36 (MQ=255)
                             cTCGTTGCGGAAGATGCGCAGGTATAGCTCctgctg       >  1:4021643/1‑36 (MQ=255)
                              tCGTTGCGGAAGATGCGCAGGTATAGCTCCTGCTGa      <  1:7143000/36‑1 (MQ=255)
                                gTTGCGGAAGATGCGCAGGTATAGCTCCTGCTGaaa    >  1:4223063/1‑36 (MQ=255)
                                gTTGCGGAAGATGCGCAGGTATAGCTCCTGCTGaaa    <  1:5162120/36‑1 (MQ=255)
                                gTTGCGGAAGATGCGCAGGTATAGCTCCTGCTGaaa    <  1:3357040/36‑1 (MQ=255)
                                gTTGCGGAAGATGCGCAGGTATAGCTCCTGCTGaaa    >  1:28188/1‑36 (MQ=255)
                                gTTGAGGAAGATGCGCAGGTATAGCTCCTGCTGaaa    <  1:7689632/36‑1 (MQ=255)
                                 ttGCGGAAGATGCGCAGGTATAGCTCCTGCTGAAAc   <  1:4642057/36‑1 (MQ=255)
                                 ttGCGGAAGATGCGCAGGTATAGCTCCTGCTGAAAc   <  1:906258/36‑1 (MQ=255)
                                 ttGCGGAAGATGCGCAGGTATAGCTCCTGCTGAAAc   <  1:3384532/36‑1 (MQ=255)
                                  tGCGGAAGATGCGCAGGTATAGCTCCTGCTGAAAct  <  1:2071704/36‑1 (MQ=255)
                                  tGCGGAAGATGCGCAGGTATAGCTCCTGCTGAAAct  <  1:11904/36‑1 (MQ=255)
                                   |                                  
ATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATAGCTCCTGCTGAAACT  >  REL606/66576‑66645

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 18 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG