Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,136,6290G→C11.8% 3.1 ‑6.8 0.074319uppPundecaprenyl pyrophosphate phosphatase
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (0/10):  ref base (35/29):  total (35/39)
Fisher's exact test strand distribution p-value = 1.14e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GTTAAACAATGACTTAATCGTGTCGTGGAACAACAGCCCCAGTACCACCGCCGGAATCATCCCCAGCAAA  >  REL606/3136594‑3136663
                                   |                                  
gTTAAACAATGACTTAATCGTGTCGTGGAACAACAg                                    >  1:4098714/1‑36 (MQ=255)
gTTAAACAATGACTTAATCGTGTCGTGGAACAACAg                                    >  1:5989578/1‑36 (MQ=255)
gTTAAACAATGACTTAATCGTGTCGTGGAACAACAg                                    >  1:6638980/1‑36 (MQ=255)
 ttAAACAATGACTTAATCGTGTCGTGGAACAACAGc                                   >  1:5852923/1‑36 (MQ=255)
 ttAAACAATGACTTAATCGTGTCGTGGAACAACAGc                                   <  1:3923779/36‑1 (MQ=255)
 ttAAACAATGACTTAATCGTGTCGTGGAACAACAGc                                   <  1:7809208/36‑1 (MQ=255)
  tAAACAATGACTTAATCGTGTCGTGGAACAACAGcc                                  >  1:4870969/1‑36 (MQ=255)
  tAAACAATGACTTAATCGTGTCGTGGAACAACAGcc                                  <  1:3760578/36‑1 (MQ=255)
   aaaCAATGACTTAATCGTGTCGTGGAACAACAGccc                                 <  1:2956495/36‑1 (MQ=255)
   aaaCAATGACTTAATCGTGTCGTGGAACAACAGccc                                 <  1:7738110/36‑1 (MQ=255)
    aaCAATGACTTAATCGTGTCGTGGAACAACAGcccc                                >  1:5230280/1‑36 (MQ=255)
    aaCAATGACTTAATCGTGTCGTGGAACAACAGcccc                                <  1:7050217/36‑1 (MQ=255)
     aCAATGACTTAATCGTGTCGTGGAACAACAGCCCCa                               <  1:4133417/36‑1 (MQ=255)
     aCAATGACTTAATCGTGTCGTGGAACAACAGCCCCa                               >  1:3252734/1‑36 (MQ=255)
     aCAATGACGTAATCGTGTCGTGGAACAACAGCCCCa                               >  1:7145747/1‑36 (MQ=255)
      cAATGACTTAATCGTGTCGTGGAAGAACAGCCCCAg                              <  1:1500824/36‑1 (MQ=255)
      cAATGACTTAATCGTGTCGTGGAACAACAGCCCCAg                              >  1:5374720/1‑36 (MQ=255)
      cAATGACTTAATCGTGTCGTGGAACAACAGCCCCAg                              >  1:868522/1‑36 (MQ=255)
       aaTGACTTAATCGTGTCGTGGAACAACAGCCCCAGt                             >  1:1288461/1‑36 (MQ=255)
       aaTGACTTAATCGTGTCGTGGAACAACAGCCCCAGt                             <  1:5803891/36‑1 (MQ=255)
       aaTGACTTAATCGTGTCGTGGAACAACAGCCCCAGt                             <  1:1458474/36‑1 (MQ=255)
        aTGACTTAATCGTGTCGTGGAACAACAGCCCCAGTa                            <  1:1085255/36‑1 (MQ=255)
        aTGACTTAATCGTGTCGTGGAACAACAGCCCCAGTa                            >  1:4323293/1‑36 (MQ=255)
         tggcTTAATCGTGTCGTGGAACAACAGCCCCAGTac                           <  1:5498252/33‑1 (MQ=255)
         tGACTTAATCGTGTCGTGGAACAACAGCCCCAGTac                           <  1:3672270/36‑1 (MQ=255)
         tGACTTAATCGTGTCGTGGAACAACAGCCCCAGTac                           <  1:5265487/36‑1 (MQ=255)
         tGACTTAATCGTGTCGTGGAACAACAGCCCCAGTac                           <  1:1526084/36‑1 (MQ=255)
         tGACTTAATCGTGTCGTGGAACAACAGCCCCAGTac                           <  1:5347728/36‑1 (MQ=255)
         tGACTTAATCGTGTCGTGGAACAACAGCCCCAGTac                           <  1:2166016/36‑1 (MQ=255)
         tGACTTAATCGTGTCGTGGAACAACAGCCCCAGTac                           <  1:6892232/36‑1 (MQ=255)
           aCTTAATCGTGTCGTGGAACAACAGCCCCAGTacca                         >  1:7406899/1‑36 (MQ=255)
           aCTTAATCGTGTCGTGGAACAACAGCCCCAGTacca                         >  1:1509211/1‑36 (MQ=255)
             ttAATCGTGTCGTGGAACAACAGCCCCAGTaccacc                       <  1:7805924/36‑1 (MQ=255)
             ttAATCGTGTCGTGGAACAACAGCCCCAGTaccacc                       <  1:4679966/36‑1 (MQ=255)
             taaaTCGTGTCGTGGAACAACAGCCCCAGTaccacc                       <  1:4602271/34‑1 (MQ=255)
              tAATCGTGTCGTGTAACAACAGCCCCAGTCCCAccg                      <  1:5215277/36‑1 (MQ=37)
              tAATCGTGTCGTGGAACAACAGCCCCAGTACCAccg                      <  1:5286297/36‑1 (MQ=255)
               aaTCGTGTCGTGTAACAACAGCCCCAGTACCAccgc                     <  1:5919715/36‑1 (MQ=255)
               aaTCGTGTCGTGGAACAACAGCCCCAGTACCAccgc                     >  1:6303706/1‑36 (MQ=255)
               aaTCGTGTCGTGGAACAACAGCCCCAGTACCAccgc                     >  1:3984553/1‑36 (MQ=255)
                aTCGTGTCGTGGAACAACAGCCCCAGTACCAccgcc                    >  1:4503467/1‑36 (MQ=255)
                aTCGTGTCGTGGAACAACAGCCCCAGTACCAccgcc                    >  1:1115682/1‑36 (MQ=255)
                aTCGTGTCGTGGAAAAACAGCCCCAGTACCAccgcc                    <  1:1688634/36‑1 (MQ=255)
                aTCGTGTCGTGAAAAAACAGCCCCAGTACCAccgcc                    <  1:991633/36‑1 (MQ=37)
                 tcgtgtcgtgGAACAACAGCCCCAGTACCAccgccg                   >  1:5116978/1‑36 (MQ=255)
                 tcgtgtcgtgGAACAACAGCCCCAGTACCAccgccg                   >  1:1797960/1‑36 (MQ=255)
                  cgtgtcgtgGAACAACAGCCCCAGTACCACCGCCgg                  >  1:2828808/1‑36 (MQ=255)
                   gtgtcgtgGAACAACAGCCCCAGTACCACCGCCGGa                 >  1:7220419/1‑36 (MQ=255)
                   gtgtcgtgGAACAACAGCCCCAGTACCACCGCCGGa                 >  1:6809480/1‑36 (MQ=255)
                    tgtcgtgGAACAACAGCCCCAGTACCACCGCCGGaa                >  1:6988092/1‑36 (MQ=255)
                    tgtcgtgGAACAACAGCCCCAGTACCACCGCCGGaa                >  1:3453596/1‑36 (MQ=255)
                     gtcgtgGAACAACAGCCCCAGTACCACCGCCGGAat               >  1:5800058/1‑36 (MQ=255)
                     gtcgggGAACAACAGCCCCAGTACCACCGCCGGAat               <  1:678882/36‑1 (MQ=255)
                       cgtgGAACAACAGCCCCAGTACCACCGCCGGAatca             >  1:2992174/1‑36 (MQ=255)
                       cgtgGAACAACAGCCCCAGTACCACCGCCGGAatca             >  1:6194926/1‑36 (MQ=255)
                       cgtgGAACAACAGCCCCAGTACCACCGCCGGAatca             >  1:4131103/1‑36 (MQ=255)
                       cgtgAAACAACACCCCCAGTACCACCGCCGGAatca             <  1:2915895/36‑1 (MQ=38)
                        gtgGAACAACAGCCCCAGTACCACCGCCGGAatcat            >  1:2723366/1‑36 (MQ=255)
                          ggAACAACAGCCCCAGTACCACCGCCGGAATCATcc          >  1:4098557/1‑36 (MQ=255)
                          ggAACAACAGCCCCAGTACCACCGCCGGAATCATcc          >  1:6494079/1‑36 (MQ=255)
                           aaacaacAGCCCCAGTACCACCGCCGGAATCATccc         <  1:4699676/35‑1 (MQ=255)
                            aacaccACCCCAACTCCCACCGCCGGAATCATcccc        <  1:4768733/36‑1 (MQ=16)
                            aacaacAGCCCCAGTACCACCGCCGGAATCATcccc        >  1:2232078/1‑36 (MQ=255)
                            aacaacACCCCCAGTACCACCGCCGGAATCATcccc        <  1:2314648/36‑1 (MQ=255)
                            aacaacACCCCCAGTACCACCGCCGGAATCATcccc        <  1:4757652/36‑1 (MQ=255)
                            aacaacACCCCCAAGGCCACCGCCGGAATCATcccc        <  1:541416/36‑1 (MQ=19)
                            aaaaacaCCCCCAGCCCCCCCGCCGGAATCATcccc        <  1:4087094/33‑1 (MQ=12)
                             ccaacaCCCCCATTACCACCGCCGGAATCATCCCCa       <  1:1876333/35‑1 (MQ=37)
                             caaacaGCCCCAGTACCACCGCCGGAATCATCCCCa       <  1:274305/34‑1 (MQ=255)
                             acaacaGCCCCAGTACCACCGCCGGAATCATCCCCa       >  1:7038/1‑36 (MQ=255)
                             acaacaCCCCAAGTACCACCGCCGGAATCATCCCCa       <  1:1400367/36‑1 (MQ=38)
                             aaaacaCCCCCATTACCACCGCCGGAATCATCCCCa       <  1:6988538/34‑1 (MQ=37)
                             aaaacaCCCCCACTGCCCCCGCCGGAATCATCCCCa       <  1:1803113/34‑1 (MQ=16)
                              caacaGCCCCAGTACCACAGCCGGAATCATCCCCAg      >  1:725395/1‑36 (MQ=255)
                                 caGCCCCAGTACCACCGCCGGAATCATCCCCAGCaa   >  1:730022/1‑36 (MQ=255)
                                 caGCCCCAGGACCACCGCCGGAATCATCCCCAGCaa   >  1:5453561/1‑36 (MQ=255)
                                  aGCCCCAGTACCACCGCCGGAATCATCCCCAGCaaa  >  1:2402792/1‑36 (MQ=255)
                                   |                                  
GTTAAACAATGACTTAATCGTGTCGTGGAACAACAGCCCCAGTACCACCGCCGGAATCATCCCCAGCAAA  >  REL606/3136594‑3136663

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